MINLOC是否适用于从索引0开始的数组？ （Fortran 90/95）

Question

使用C一段时间后，我回到Fortran并在我的代码中将数组从索引0分配到N：

real(kind=dp), dimension(:), allocatable :: a 
allocate(a(0:50))

我需要找到数组最小绝对值的索引，所以我使用MINLOC，并检查这个我将它与MINVAL进行比较：

minloc(abs(a(:)))
minval(abs(a))

MINLOC的结果是索引42，但MINVAL的结果对应41。以下是输出中的相关部分：

Index i    a(i) 

39         0.04667    
40         0.02222    
41         0.00222           !This was clearly the minimum value
42         0.02667

MINLOC = 42
MINVAL = 0.00222

我认为这与Fortran内部没有正确处理索引为0的数组有关，因为以这种方式声明数组不是标准的Fortran样式（但它仍然是允许的！）。

任何人都可以确认或提供解决方法吗？

Answer 1

您的数组a确实从索引0开始，但您没有使用它。您搜索了至少一个数组abs(a(:))。此匿名数组表达式从1开始，因为所有数组都默认执行。

但即使您使用a，结果也是一样的，并且与数组参数传递在Fortran中的工作方式一致。

Fortran标准明确指出：

  

返回的i下标位于1到ei的范围内，其中ei是   ARRAY的范围。如果ARRAY的大小为零，则全部   结果的元素为零。

如果使用假定的形状参数，下限不会自动与数组一起传递。例如，如果您有自己的功能

  function f(arg)
    real :: arg(:)

无论调用代码中的实际参数在何处开始，

arg始终为1。

您可以将其更改为从其他值开始

  function f(arg)
    real :: arg(-42:)

它将从该值开始编入索引。

Answer 2

有两种简单的方法来处理调整从minloc（）获得的索引的复杂性：一种是简单地为所有索引添加lbound（） - 1，另一种是使用具有从1开始的索引的数组指针。示例代码可能如下所示：

(keep-indexed (fn [idx v] (if (pos? v) idx)) input-vector)