曼哈顿距离 Fortran

Question

我想根据公式计算Fortran中2个数组之间的曼哈顿距离：

d = Sum(|P(i)-R(i)|)

所以我在 Fortran 中做了一个代码：

function DistM(v, u, dim)
    integer dim
    real(8) v(dim), u(dim), DistM
    DistM=sum(abs(v-u))
end function DistM

我用这个来调用这个函数：

Coeff=DistM(tempvector1,tempvector2, dim)

但这似乎不起作用（我没有得到任何回报）。我还尝试取 2 的幂，然后对其进行 sqrt，但它会卡在 sqrt 上（如果我分步运行它）。 我尝试了另一个功能并且那个有效（见下文），但是这个不起作用:(:

function Roznica(v, u, dim)
        integer dim
        real(8) v(dim), u(dim), Rozn
        Rozn=sum((v-u)**2)/dim
end function Roznica

谁有想法？

Answer 1

我会像这样简单地编写一个曼哈顿距离函数（使用假定形状的数组，您不再需要像您那样使用自动数组，前提是该函数具有显式接口）：

pure function L1(v)
  real, intent(in) :: v(:)  ! <- shape assumed from actual argument
  real :: L1
  L1 = sum(abs(v))
end function

然后，如果您有两个向量，您只需将函数称为 foo = L1(p - q)。