levenshtein

Levenshtein edit distance

  1. levenshtein
  2. Documentation
    1. 8-bit Values Only
      1. levenshtein-distance/byte
        1. Usage
      2. levenshtein-distance/transpose-byte
        1. Usage
    2. Parameterized Edit Distance
      1. Usage
      2. levenshtein-distance/sequence
    3. Parameterized Edit Distance with Trailing
      1. Usage
      2. levenshtein-distance/vector*
    4. Interface Implementation Modules
      1. Cost Implementation Modules
      2. Sequence Implementation Modules
    5. Edit Operators
      1. Usage
      2. levenshtein-operator
      3. make-levenshtein-operator
      4. levenshtein-operator?
      5. clone-levenshtein-operator
      6. levenshtein-operator-ref
      7. levenshtein-operator-set!
      8. levenshtein-operator-delete!
      9. levenshtein-operator-reset
      10. levenshtein-operator=?
      11. levenshtein-insert-operator?
      12. levenshtein-delete-operator?
      13. levenshtein-substitute-operator?
      14. levenshtein-required-operator?
      15. levenshtein-required-operators
      16. levenshtein-base-operators
      17. levenshtein-extended-operators
    6. Path Iterator
      1. Usage
      2. levenshtein-path-iterator
    7. Path Matrix Print
      1. Usage
      2. print-levenshtein-matrix
  3. Bugs & Limitations
  4. Examples
  5. Requirements
  6. Author
  7. Version history
  8. License

Documentation

Levenshtein is a collection of procedures providing various forms of the Levenshtein edit distance calculation.

The Levenshtein edit distance has been used for areas as diverse as soil sample and language dialect analysis. Not just for text strings.

8-bit Values Only

Performs edit distance calculation for byte strings & blobs. All return the total edit cost.

levenshtein-distance/byte

Usage
(import levenshtein-byte)
[procedure] (levenshtein-distance/byte SOURCE TARGET)

Calculates the edit distance from the SOURCE to the TARGET. All costs are unitary.

levenshtein-distance/transpose-byte

Usage
(import levenshtein-transpose-byte)
[procedure] (levenshtein-distance/transpose-byte SOURCE TARGET)

Calculates the edit distance from the SOURCE to the TARGET, taking into account the Transpose operation. All costs are unitary.

By using the Transpose operation the total edit cost is not at least the difference of the sizes of the two strings.

Parameterized Edit Distance

A functor implementing an edit distance algorithm parameterized by a cost and sequence operation modules.

Performs edit distance calculation for sequences.

Usage

(import trunk/levenshtein-sequence-functor)

levenshtein-distance/sequence

[procedure] (levenshtein-distance/sequence SOURCE TARGET [:insert-cost INSERT-COST] [#:delete-cost DELETE-COST] [#:substitute-cost SUBSTITUTE-COST] (#:get-work-vector GET-WORK-VECTOR) [#:elm-eql ELM-EQL] [#:limit-cost LIMIT-COST]) -> (or false cost)
SOURCE
string.
TARGET
string.
INSERT-COST
number, default 1.
DELETE-COST
number, default 1.
SUBSTITUTE-COST
number, default 1.
ELM-EQL
procedure; (-> object object boolean), default eqv?. The equality predicate.
GET-WORK-VECTOR
procedure, default make-vector.
LIMIT-COST
number or #f, default #f. Quit when cost over limit & return #f.

The SOURCE & TARGET must be the same type, which the instantiating sequence module supports.

Note that the element comparison procedure is passed via the argument list, and not via the sequence implementation module. Annoying when using strings but useful when using vectors.

Parameterized Edit Distance with Trailing

A functor implementing an edit distance algorithm parameterized by a cost operation module.

Performs edit distance calculation for vectors. Allows definition of new edit operations. Will keep track of edit operations performed; the path-matrix. Primarily a toy.

Usage

(import levenshtein-vector-functor)

levenshtein-distance/vector*

[procedure] (levenshtein-distance/vector* SOURCE TARGET [EDIT-OPER ...] [#:elm-eql ELM-EQL] [#:operations? OPERATIONS] [#:limit-cost LIMIT-COST]) -> (or false cost) path-matrix

Calculates the edit distance from the source vector SOURCE to the target vector TARGET. Returns the total edit cost or (values <total edit cost> <performed operations matrix>).

SOURCE
vector.
TARGET
vector.
EDIT-OPER
levenshtein-operator. Edit operation definitions to apply. Defaults are the basic Insert, Delete, and Substitute.
ELM-EQL
procedure; (-> object object boolean), default char=?. The equality predicate.
OPERATIONS
boolean. Include the matrix of edit operations performed? Default #f.
LIMIT-COST
(or number false), default #f. Quit when cost over limit & return #f.

Interface Implementation Modules

Cost Implementation Modules

For use with levenshtein-distance/sequence & levenshtein-distance/vector*.

Sequence Implementation Modules

For use with levenshtein-distance/sequence.

Edit Operators

Edit operation specification. A set of base operations is predefined, but may be overridden. The base set is identified by the keys Insert, Delete, Substitute, and Transpose. A printer and reader are provided for edit operations.

Usage

(import levenshtein-operators)

levenshtein-operator

[record] levenshtein-operator
[procedure] (levenshtein-operator-key OPER) -> {{symbol}}
[procedure] (levenshtein-operator-name OPER) -> {{string}}
[procedure] (levenshtein-operator-cost OPER) -> {{number}}
[procedure] (levenshtein-operator-above OPER) -> {{fixnum}}
[procedure] (levenshtein-operator-left OPER) -> {{fixnum}}

make-levenshtein-operator

[procedure] (make-levenshtein-operator KEY NAME COST ABOVE LEFT)

Returns a new edit operator.

KEY
symbol. Key for the operation.
NAME
string. Describes the operation.
COST
number. The cost of the operation.
ABOVE
non-negative-fixnum. How far back in the source.
LEFT
non-negative-fixnum. How far back in the target.

levenshtein-operator?

[procedure] (levenshtein-operator? OBJECT)

Is the OBJECT a levenshtein operator?

clone-levenshtein-operator

[procedure] (clone-levenshtein-operator EDIT-OPERATION [#:key KEY] [#:name NAME] [#:cost COST] [#:above ABOVE] [#:left LEFT])

Returns a duplicate of the EDIT-OPERATION, with field values provided by the optional keyword arguments. EDIT-OPERATION may be the key of the already defined edit operation.

levenshtein-operator-ref

[procedure] (levenshtein-operator-ref KEY)

Get the definition of an edit operation.

levenshtein-operator-set!

[procedure] (levenshtein-operator-set! EDIT-OPERATION)

Define an edit operation.

levenshtein-operator-delete!

[procedure] (levenshtein-operator-delete! EDIT-OPERATION)

Removes the EDIT-OPERATION definition. EDIT-OPERATION may be the KEY of the already defined edit operation.

levenshtein-operator-reset

[procedure] (levenshtein-operator-reset)

Restore defined edit operations to the base set.

levenshtein-operator=?

[procedure] (levenshtein-operator=? A B)

Are the levenshtein-operator A & levenshtein-operator B equal for all fields?

levenshtein-insert-operator?

levenshtein-delete-operator?

levenshtein-substitute-operator?

levenshtein-required-operator?

[procedure] (levenshtein-insert-operator? OBJECT)
[procedure] (levenshtein-delete-operator? OBJECT)
[procedure] (levenshtein-substitute-operator? OBJECT)
[procedure] (levenshtein-required-operator? OBJECT)

levenshtein-required-operators

[procedure] (levenshtein-required-operators) -> (vector-of levenshtein-operator)

levenshtein-base-operators

[procedure] (levenshtein-base-operators) -> (vector-of levenshtein-operator)

levenshtein-extended-operators

[procedure] (levenshtein-extended-operators OPERLIST) -> (vector-of levenshtein-operator)
OPERLIST
(list-of levenshtein-operator) ; operations to add to the required

Path Iterator

Usage

(import levenshtein-path-iterator)

levenshtein-path-iterator

[procedure] (levenshtein-path-iterator PATH-MATRIX) -> (-> (or false list))

Creates an optimal edit distance operation path iterator over the performed operations matrix PATH-MATRIX. The matrix is usually the result of an invocation of (levenshtein-distance/vector* ... operations: #t).

Each invocation of the iterator will generate a list of the form: ((cost source-index target-index levenshtein-operator) ...). The last invocation will return #f.

Path Matrix Print

Usage

(import levenshtein-print)
[procedure] (print-levenshtein-matrix PATH-MATRIX)

Displays a readable representation of the PATH-MATRIX on the current-output-port.

Bugs & Limitations

C5_CACHE="$HOME/.cache/chicken-install/levenshtein"
C5_HOME=$(csi -n -R chicken.platform -p '(chicken-home)')
cp "$C5_CACHE/levenshtein-cost-*.scm" \
   "$C5_CACHE/levenshtein-sequence-*.scm" "$C5_HOME"

Examples

(import levenshtein-byte)
(levenshtein-distance/byte "ctas" "cats")
;=> 2
(import levenshtein-transpose-byte)
(levenshtein-distance/transpose-byte "ctas" "cats")
;=> 1
;until R⁷RS
(define (string->vector s) (list->vector (string->list s)))
(import levenshtein-sequence-functor levenshtein-cost-fixnum levenshtein-sequence-vector)
(module levenshtein-sequence-fixnum-vector =
  (levenshtein-sequence-functor levenshtein-cost-fixnum levenshtein-sequence-vector))
(import (prefix levenshtein-sequence-fixnum-vector vsfx:))
(vsfx:levenshtein-distance/sequence
  (string->vector "ctas") (string->vector "cats")
  #:elm-eql char=?)
;=> 2
(import levenshtein-operators)
(import levenshtein-vector-functor levenshtein-cost-fixnum)
(module levenshtein-vector-fixnum =
  (levenshtein-vector-functor levenshtein-cost-fixnum))
(import (prefix levenshtein-vector-fixnum vcfx:))
(vcfx:levenshtein-distance/vector*
  (string->vector "ctas") (string->vector "cats")
  #:elm-eql char=?
  ;NOTE if supplying any additional operations
  ;then must supply all operations beyond those
  ;required (insert & delete). While substitute
  ;is a base operation, it is not necessary.
  (levenshtein-operator-ref 'Substitute)
  (levenshtein-operator-ref 'Transpose))
;=> 1
;=> #f
;=> ; 2 values
(import levenshtein-print)
;only useful for the path-matrix
(import levenshtein-vector-functor levenshtein-cost-fixnum)
(module levenshtein-vector-fixnum = (levenshtein-vector-functor levenshtein-cost-fixnum))
(import (prefix levenshtein-vector-fixnum fixn:))
;NOTE default #:elm-eql keyword parameter is char=?
(let-values (((cost pm)
               (fixn:levenshtein-distance/vector*
                 (string->vector "ctas") (string->vector "cats")
                 operations: #t)) )
  (print cost)
  (print-levenshtein-matrix pm) )
;=> 2
;=> ( 0 Substitute) ( 1     Delete) ( 2     Delete) ( 3     Delete)
;=> ( 1     Insert) ( 1 Substitute) ( 1 Substitute) ( 2     Delete)
;=> ( 2     Insert) ( 1 Substitute) ( 2 Substitute) ( 2 Substitute)
;=> ( 3     Insert) ( 2     Insert) ( 2 Substitute) ( 2 Substitute)
; Instantiate the distance measure algorithm
(import levenshtein-path-iterator)
(import levenshtein-vector-functor levenshtein-cost-fixnum)
(module levenshtein-vector-fixnum = (levenshtein-vector-functor levenshtein-cost-fixnum))
(import (prefix levenshtein-vector-fixnum vcfx:))

(define iter
  (levenshtein-path-iterator
    (vcfx:levenshtein-distance/vector*
      (string->vector "YWCQPGK") (string->vector "LAWYQQKPGKA")
      operations: #t)))
; ignoring interpreter feedback & we know the distance is 6
(define r0 (iter))
(define t0 r0)
(define r1 (iter))
(define r2 (iter))
(define r3 (iter))
(define r4 (iter))
(define r5 (iter))
(iter)
; r0 now has #f, since the iterator finishes by returning to the initial caller,
; which is the body of '(define r0 (iter))', thus re-binding r0. However, t0 has
; the original returned value.

Requirements

[check-errors]] [vector-lib]] [srfi-1]] [srfi-13]] [srfi-63]] [srfi-69]] [utf8]] [miscmacros]] record-variants

[test]] [test-utils]]

Author

Kon Lovett

Version history

2.4.2
Fix levenshtein-operator printing.
2.4.1
Correctly type cost interface constants. sequence-string is sequence-utf8.
2.4.0
Operations extend required insert & delete. Add levenshtein-extended-operators & others.
2.3.0
Imported implementation modules. levenshtein-sequence-string is levenshtein-sequence-utf8.
2.2.6
Deactivate functors tests (see Bugs & Limitations).
2.2.7
Activate functors tests.
2.2.6
Fix total-cost result.
2.2.5
.
2.2.4
.
2.2.3
.
2.2.2
Fix levenshtein-distance/vector* type. Add tests,
2.2.1
Fix levenshtein-sequence-vector,
2.2.0
Use record-variants,
2.1.3
Fix levenshtein-path-iterator,
2.1.2
.
2.0.0
Chicken 5 release.
1.0.3
Added types. Re-flow.
1.0.2
Added an "egg tag".
1.0.1
Drop "format-compiler".
1.0.0
Chicken 4 release.

License

Copyright (c) 2012-2024, Kon Lovett. All rights reserved.

Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the Software), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED ASIS, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.