Просто возиться с F #, и я пытался создать основную функцию Лагранжа Интерполяция на основе этой C# версии (скопированного из записи вики C++):Переписывая C# код в F #
double Lagrange(double[] pos, double[] val, double desiredPos)
{
double retVal = 0;
for (int i = 0; i < val.Length; ++i)
{
double weight = 1;
for (int j = 0; j < val.Length; ++j)
{
// The i-th term has to be skipped
if (j != i)
{
weight *= (desiredPos - pos[j])/(pos[i] - pos[j]);
}
}
retVal += weight * val[i];
}
return retVal;
}
Лучшее, что я мог бы прийти с помощью моего ограниченного знания F # и функционального программирования было:
let rec GetWeight desiredPos i j (pos : float[]) weight =
match i with
| i when j = pos.Length -> weight
| i when i = j -> GetWeight desiredPos i (j+1) pos weight
| i -> GetWeight desiredPos i (j+1) pos (weight * (desiredPos - pos.[j])/(pos.[i] - pos.[j]))
let rec Lagrange (pos : float[]) (vals : float[]) desiredPos result counter =
match counter with
| counter when counter = pos.Length -> result
| counter -> Lagrange pos vals desiredPos (result + (GetWeight desiredPos counter 0 pos 1.0)* vals.[counter]) (counter+1)
Может кто-то обеспечить лучшее/опрятнее F # версии на основе же C# код?
Я думаю, что это хороший пример того, когда императивный код легче читать и поддерживать, чем функциональный. – gradbot