主题:cgl_lgs请指教!
你好,冒昧打扰,
10.4日本人发帖关于内存不足的问题,当时建议改为动态分配,、
但是更改后还是同样的‘memory allocation failed'.
所以请问,是不是我得矩阵就是太大了,我这台电脑不足以运行?
我把更改后的部分贴上来,请专家们给看看:
parameter(np=9000)
c integer n,np,indx(9000)
c real*8 a(np,np),a2(np,np),y(np,np),prod(np,np),d
c character result*8
c---------------Dynamic allocation---------------------
integer n,np
integer,allocatable:: indx(:)
real*8,allocatable:: a(:,:),a2(:,:),y(:,:),prod(:,:)
real*8 d
character result*8
allocate(indx(np))
allocate(a(np,np),a2(np,np),y(np,np),prod(np,np))
c_________________________________________________________
open(unit=8,file='matrix.dat',status='unknown')
open(unit=10,file='prepared.dat',status='old')
c
所出现的问题提示在附件里显示。
谢谢大家
10.4日本人发帖关于内存不足的问题,当时建议改为动态分配,、
但是更改后还是同样的‘memory allocation failed'.
所以请问,是不是我得矩阵就是太大了,我这台电脑不足以运行?
我把更改后的部分贴上来,请专家们给看看:
parameter(np=9000)
c integer n,np,indx(9000)
c real*8 a(np,np),a2(np,np),y(np,np),prod(np,np),d
c character result*8
c---------------Dynamic allocation---------------------
integer n,np
integer,allocatable:: indx(:)
real*8,allocatable:: a(:,:),a2(:,:),y(:,:),prod(:,:)
real*8 d
character result*8
allocate(indx(np))
allocate(a(np,np),a2(np,np),y(np,np),prod(np,np))
c_________________________________________________________
open(unit=8,file='matrix.dat',status='unknown')
open(unit=10,file='prepared.dat',status='old')
c
所出现的问题提示在附件里显示。
谢谢大家