3

我正在尝试为数学算法（共轭梯度法）运行代码。在这样做时，我输入了一个双精度矩阵，在我的序言中如此定义。编译时，我收到以下错误：

A=RESHAPE((/ 0,8,0,4,26,8,0,17.5,0,0,0,17.5,0,2.5,-8,4,0,2.5,0,-5,26,0,-8,-5,0 
                        1
Error: Element in INTEGER(4) array constructor at (1) is REAL(4)
make: FTranProjectBuilder: Error: Execution exited with code 2 
*** [cg_main.o] Error 1

我在程序中定义了矩阵的定义是这样给出的（数组定义是我程序的第一个操作）：

PROGRAM cg_main 
IMPLICIT NONE 

INTEGER,PARAMETER                     ::d=5 !use a parameter for the dimensions (simple)
DOUBLE PRECISION,DIMENSION(d,d)       ::A !matrix
INTEGER,DIMENSION(2)                  ::order2 = (/ 2, 1 /) !matrix reshape order

[MORE DECLARATIONS HERE]

A=RESHAPE((/ 0,8,0,4,26,8,0,17.5,0,0,0,17.5,0,2.5,-8,4,0,2.5,0,-5,26,0,-8,-5,0 /),(/d,d/), order2) !specify dxd matrix

[MORE CODE HERE]

END PROGRAM

该代码在我的矩阵输入中没有十进制数字的情况下工作，但似乎与我的小数无关，我不知道为什么。

4

1 回答 1

2

前面的回答是不正确的，标准对数组构造函数的处理方式很清楚。请参阅 Fortran 2008 C4103。

问题形式的数组构造函数中的所有元素必须具有相同的声明类型和相同类型的类型参数。此处的编译器会抱怨，因为将实数表达式与那些整数表达式混合在一起违反了此约束。

仅仅使第一个元素成为实数/双精度元素是不够的。也许编译器正在使用有关“数组构造函数”的错误消息integer(4)作为“第一个元素是integer(4)但不是所有其他元素是”的简写。

为了解决这个问题，使用这种形式的数组构造函数，您应该使所有元素具有相同的类型/种类（可能是实数/双精度）。

正如其他答案中所述（正确），数组构造函数还有另一种形式：

[ real ::0,8,0,4,26,8,0,17.5,0,0,0,17.5 ... ]

它将所有元素转换为指定的类型。你甚至可以使用这种形式来创建一个大小为零的数组：

 [ real :: ]
于 2017-10-05T08:26:50.877 回答