Почему этот код F # так медленно?

Реализация Левенштейна в С# и F #. Версия С# в 10 раз быстрее для двух строк около 1500 символов. С#: 69 мс, F # 867 мс. Зачем? Насколько я могу судить, они делают то же самое? Не имеет значения, является ли это сборкой Release или Debug.

EDIT: Если кто-то приходит сюда специально для реализации Edit Distance, он сломан. Рабочий код здесь.

С#:

private static int min3(int a, int b, int c)
{
   return Math.Min(Math.Min(a, b), c);
}

public static int EditDistance(string m, string n)
{
   var d1 = new int[n.Length];
   for (int x = 0; x < d1.Length; x++) d1[x] = x;
   var d0 = new int[n.Length];
   for(int i = 1; i < m.Length; i++)
   {
      d0[0] = i;
      var ui = m[i];
      for (int j = 1; j < n.Length; j++ )
      {
         d0[j] = 1 + min3(d1[j], d0[j - 1], d1[j - 1] + (ui == n[j] ? -1 : 0));
      }
      Array.Copy(d0, d1, d1.Length);
   }
   return d0[n.Length - 1];
}

F #

let min3(a, b, c) = min a (min b c)

let levenshtein (m:string) (n:string) =
   let d1 = Array.init n.Length id
   let d0 = Array.create n.Length 0
   for i=1 to m.Length-1 do
      d0.[0] <- i
      let ui = m.[i]
      for j=1 to n.Length-1 do
         d0.[j] <- 1 + min3(d1.[j], d0.[j-1], d1.[j-1] + if ui = n.[j] then -1 else 0)
      Array.blit d0 0 d1 0 n.Length
   d0.[n.Length-1]

Ответ 1

Проблема заключается в том, что функция min3 скомпилирована как общая функция, которая использует общее сравнение (я думал, что это использует только IComparable, но на самом деле это более сложно - он будет использовать структурное сравнение для типов F #, и это справедливо сложная логика).

> let min3(a, b, c) = min a (min b c);;
val min3 : 'a * 'a * 'a -> 'a when 'a : comparison

В версии С# эта функция не является общей (она принимает только int). Вы можете улучшить версию F #, добавив аннотации типов (чтобы получить то же, что и на С#):

let min3(a:int, b, c) = min a (min b c)

... или делая min3 как inline (в этом случае он будет специализироваться на int при использовании):

let inline min3(a, b, c) = min a (min b c);;

Для случайной строки str длиной 300 я получаю следующие числа:

> levenshtein str ("foo" + str);;
Real: 00:00:03.938, CPU: 00:00:03.900, GC gen0: 275, gen1: 1, gen2: 0
val it : int = 3

> levenshtein_inlined str ("foo" + str);;
Real: 00:00:00.068, CPU: 00:00:00.078, GC gen0: 0, gen1: 0, gen2: 0
val it : int = 3