一个matmul的问题

little_kar

请教大家一个问题，我在用matmul这个函数的时候，ZKM=ZKM+matmul(PST,dk,PS)  !总体刚度矩阵!
ZFM=ZFM+matmul(PST,FV)  !总体载荷向量!
其中，PS是3*90000的二维数组，dk是3*3的，PST是PS的转置，FV是3*1的二维数组，
然后调试的时候，下面有一个警告和一个错误，分别是这样的：
警告：The number of arguments is incompatible with intrinsic procedure, assume 'external'.   [MATMUL]
错误：This name does not have a type, and must have an explicit type.   [MATMUL]
请问怎么回事，怎么修改呢？小弟有点急，谢谢大家

fcode

数组是 real 类型的，但是数组的大小（编译器是用整数表示数组的大小）
是数组的大小超出了integer的范围

little_kar

fcode 发表于 2017-2-27 18:23
integer4 都表达不了数组的大小，那说明数组的确是太大太大了。用 integer8 也不行。（因为你内存就那么大~ ...

我用的都是real，也就是说，只能吧数组变小了，比如变成100*100的就可以了？

fcode

integer4 都表达不了数组的大小，那说明数组的确是太大太大了。用 integer8 也不行。（因为你内存就那么大~~）

pasuka

little_kar 发表于 2017-2-27 16:36
按照你说的  以解决，多谢多谢，额，又遇到一个问题啊：forrtl: severe (179): Cannot allocate array -  ...

使用huge函数看看是否超过整型范围
https://gcc.gnu.org/onlinedocs/gcc-6.3.0/gfortran/HUGE.html#HUGE

li913

little_kar 发表于 2017-2-27 16:36
按照你说的  以解决，多谢多谢，额，又遇到一个问题啊：forrtl: severe (179): Cannot allocate array -  ...

表示数组大小的变量，其值太大，溢出了。可以试试用integer(8).

little_kar

li913 发表于 2017-2-27 09:46
matmul 只有两个参数。matmul(PST,dk,PS)  改为  matmul(matmul(PST,dk),PS) 。理论上可以，但如果矩阵太大 ...

按照你说的  以解决，多谢多谢，额，又遇到一个问题啊：forrtl: severe (179): Cannot allocate array - overflow on array size calculation，这个不懂啊，是个断点

little_kar

matmul（matmul（PST，dk），PS）

little_kar

好的 多谢，我试试

li913

matmul 只有两个参数。matmul(PST,dk,PS)  改为  matmul(matmul(PST,dk),PS) 。理论上可以，但如果矩阵太大，实际情况就不好说了。

fcode

请给出更多的代码。尽量是当前程序单元的全部。