Removed memory routines (ms.f'
This commit is contained in:
parent
71e45d3763
commit
a6958a3b2f
@ -125,12 +125,6 @@ c
|
||||
ifg(ii0)=jjb
|
||||
90 continue
|
||||
100 continue
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('ib',i1hi+1,0)
|
||||
call memacct('b',kb,0)
|
||||
call memacct('jb',kb,0)
|
||||
c
|
||||
return
|
||||
c
|
||||
@ -213,11 +207,6 @@ c===> i0 is an f-variable
|
||||
60 continue
|
||||
ib(i0hi+1)=kb
|
||||
if(kb.gt.ndimb) go to 9902
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('b',kb,0)
|
||||
call memacct('jb',kb,0)
|
||||
return
|
||||
c
|
||||
c===> error messages
|
||||
@ -288,11 +277,6 @@ c call outa(k,5)
|
||||
20 continue
|
||||
30 continue
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('b',ib(i0hi+1),0)
|
||||
call memacct('jb',ib(i0hi+1),0)
|
||||
c
|
||||
c compress the matrix.
|
||||
c
|
||||
kb=ib(i0lo)
|
||||
|
||||
@ -99,12 +99,6 @@ c
|
||||
jcnbhi=0
|
||||
jvalxx=jvalmx+1
|
||||
nspts=0
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('icg',jvalxx,0)
|
||||
call memacct('ifg',jvalxx,0)
|
||||
call memacct('ib',ihi+ishift,0)
|
||||
c
|
||||
if(ncolor.eq.-1) then
|
||||
c
|
||||
|
||||
@ -158,14 +158,6 @@ c
|
||||
endif
|
||||
70 continue
|
||||
imax(k)=ic
|
||||
c write(6,9000) k
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('a',ka,0)
|
||||
call memacct('ja',ka,0)
|
||||
call memacct('ifg',ichi,0)
|
||||
|
||||
return
|
||||
c
|
||||
c===> error messages
|
||||
|
||||
@ -26,9 +26,7 @@ globals.h
|
||||
iinvert.f
|
||||
matrix.c
|
||||
matrix.h
|
||||
memoryx.f
|
||||
misc.f
|
||||
ms.f
|
||||
opdfn.f
|
||||
outprt.f
|
||||
pcg.c
|
||||
|
||||
@ -129,10 +129,5 @@ c
|
||||
do 100 i=ilo,ihi
|
||||
if(icg(i).gt.0) jb(ib(i))=i
|
||||
100 continue
|
||||
c
|
||||
c set permanent memory usage
|
||||
c
|
||||
call memacct('b',ib(ihi+1),1)
|
||||
call memacct('jb',ib(ihi+1),1)
|
||||
return
|
||||
end
|
||||
|
||||
@ -66,32 +66,6 @@ c
|
||||
ib(1)=1
|
||||
k=1
|
||||
c
|
||||
c set array limits
|
||||
c
|
||||
call memlimit('a',ndima)
|
||||
call memlimit('ja',ndima)
|
||||
call memlimit('ia',ndimu)
|
||||
call memlimit('iu',ndimu)
|
||||
call memlimit('ip',ndimu)
|
||||
call memlimit('icg',ndimu)
|
||||
call memlimit('ifg',ndimu)
|
||||
call memlimit('ib',ndimu)
|
||||
call memlimit('b',ndimb)
|
||||
call memlimit('jb',ndimb)
|
||||
c
|
||||
c set initial data usage (temporary)
|
||||
c
|
||||
call memacct('a',ia(imax(1)+1),1)
|
||||
call memacct('ja',ia(imax(1)+1),1)
|
||||
call memacct('ia',imax(1)+1,1)
|
||||
call memacct('iu',imax(1),1)
|
||||
call memacct('ip',imax(1),1)
|
||||
call memacct('icg',0,1)
|
||||
call memacct('ifg',0,1)
|
||||
call memacct('ib',0,1)
|
||||
call memacct('b',0,1)
|
||||
call memacct('jb',0,1)
|
||||
c
|
||||
c load ifg into ifc
|
||||
c
|
||||
c do 10 i=1,imax(1)
|
||||
@ -102,17 +76,6 @@ c
|
||||
call trunc(1,imin,imax,a,ia,ja)
|
||||
c
|
||||
call symm(1,1,imin,imax,a,ia,ja,icg,ifg)
|
||||
c
|
||||
c set fine level data usage (permanent)
|
||||
c
|
||||
call memacct('a',ia(imax(1)+1),1)
|
||||
call memacct('ja',ia(imax(1)+1),1)
|
||||
call memacct('ia',imax(1)+1,1)
|
||||
call memacct('iu',imax(1),1)
|
||||
call memacct('ip',imax(1),1)
|
||||
call memacct('ib',imax(1)+1,1)
|
||||
c
|
||||
call msample
|
||||
c
|
||||
if(levels.le.1) then
|
||||
close(6)
|
||||
@ -140,7 +103,6 @@ cveh call stats(k-1,levels,idump,
|
||||
cveh * nun,imin,imax,a,ia,ja,iu,ip,icg,ifg,
|
||||
cveh * b,ib,jb,ipmn,ipmx,iv,xp,yp)
|
||||
|
||||
c call memoryx(k-1,imin,imax,ia,ib,ipmn,ipmx)
|
||||
c
|
||||
c=> test for coarsest grid
|
||||
c
|
||||
@ -156,20 +118,6 @@ c
|
||||
call trunc(k,imin,imax,a,ia,ja)
|
||||
|
||||
call symm(k,isymm,imin,imax,a,ia,ja,icg,ifg)
|
||||
c
|
||||
c set final data usage (permanent)
|
||||
c
|
||||
call memacct('a',ia(imax(k)+1),1)
|
||||
call memacct('ja',ia(imax(k)+1),1)
|
||||
call memacct('ia',imax(k)+1,1)
|
||||
call memacct('iu',imax(k),1)
|
||||
call memacct('ip',imax(k),1)
|
||||
call memacct('icg',imax(k-1),1)
|
||||
call memacct('ib',imax(k)+1,1)
|
||||
call memacct('b',ib(imax(k-1)+1),1)
|
||||
call memacct('jb',ib(imax(k-1)+1),1)
|
||||
c
|
||||
call msample
|
||||
c
|
||||
go to 20
|
||||
c
|
||||
@ -182,7 +130,6 @@ c compute & print statistics after coarsening
|
||||
* b,ib,jb,ipmn,ipmx,iv,xp,yp)
|
||||
|
||||
|
||||
call msample
|
||||
c
|
||||
c test operators
|
||||
c
|
||||
@ -195,7 +142,6 @@ c=> compress unneeded entries from interpolation/restriction
|
||||
c
|
||||
c call crushb(levels,imin,imax,icg,b,ib,jb)
|
||||
c
|
||||
call mreport(1)
|
||||
|
||||
close(6)
|
||||
|
||||
|
||||
@ -375,12 +375,6 @@ c
|
||||
do 30 i=i0lo,i0hi
|
||||
ja(ia(i))=i
|
||||
30 continue
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('a',ka,0)
|
||||
call memacct('ja',ka,0)
|
||||
call memacct('ia',i0hi+ishift+1,0)
|
||||
c
|
||||
return
|
||||
c
|
||||
|
||||
@ -125,12 +125,6 @@ c
|
||||
ifg(ii0)=jjb
|
||||
90 continue
|
||||
100 continue
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('ib',i1hi+1,0)
|
||||
call memacct('b',kb,0)
|
||||
call memacct('jb',kb,0)
|
||||
c
|
||||
return
|
||||
c
|
||||
@ -213,11 +207,6 @@ c===> i0 is an f-variable
|
||||
60 continue
|
||||
ib(i0hi+1)=kb
|
||||
if(kb.gt.ndimb) go to 9902
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('b',kb,0)
|
||||
call memacct('jb',kb,0)
|
||||
return
|
||||
c
|
||||
c===> error messages
|
||||
@ -288,11 +277,6 @@ c call outa(k,5)
|
||||
20 continue
|
||||
30 continue
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('b',ib(i0hi+1),0)
|
||||
call memacct('jb',ib(i0hi+1),0)
|
||||
c
|
||||
c compress the matrix.
|
||||
c
|
||||
kb=ib(i0lo)
|
||||
|
||||
@ -99,12 +99,6 @@ c
|
||||
jcnbhi=0
|
||||
jvalxx=jvalmx+1
|
||||
nspts=0
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('icg',jvalxx,0)
|
||||
call memacct('ifg',jvalxx,0)
|
||||
call memacct('ib',ihi+ishift,0)
|
||||
c
|
||||
if(ncolor.eq.-1) then
|
||||
c
|
||||
|
||||
@ -158,14 +158,6 @@ c
|
||||
endif
|
||||
70 continue
|
||||
imax(k)=ic
|
||||
c write(6,9000) k
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('a',ka,0)
|
||||
call memacct('ja',ka,0)
|
||||
call memacct('ifg',ichi,0)
|
||||
|
||||
return
|
||||
c
|
||||
c===> error messages
|
||||
|
||||
@ -26,9 +26,7 @@ globals.h
|
||||
iinvert.f
|
||||
matrix.c
|
||||
matrix.h
|
||||
memoryx.f
|
||||
misc.f
|
||||
ms.f
|
||||
opdfn.f
|
||||
outprt.f
|
||||
pcg.c
|
||||
|
||||
@ -129,10 +129,5 @@ c
|
||||
do 100 i=ilo,ihi
|
||||
if(icg(i).gt.0) jb(ib(i))=i
|
||||
100 continue
|
||||
c
|
||||
c set permanent memory usage
|
||||
c
|
||||
call memacct('b',ib(ihi+1),1)
|
||||
call memacct('jb',ib(ihi+1),1)
|
||||
return
|
||||
end
|
||||
|
||||
@ -66,32 +66,6 @@ c
|
||||
ib(1)=1
|
||||
k=1
|
||||
c
|
||||
c set array limits
|
||||
c
|
||||
call memlimit('a',ndima)
|
||||
call memlimit('ja',ndima)
|
||||
call memlimit('ia',ndimu)
|
||||
call memlimit('iu',ndimu)
|
||||
call memlimit('ip',ndimu)
|
||||
call memlimit('icg',ndimu)
|
||||
call memlimit('ifg',ndimu)
|
||||
call memlimit('ib',ndimu)
|
||||
call memlimit('b',ndimb)
|
||||
call memlimit('jb',ndimb)
|
||||
c
|
||||
c set initial data usage (temporary)
|
||||
c
|
||||
call memacct('a',ia(imax(1)+1),1)
|
||||
call memacct('ja',ia(imax(1)+1),1)
|
||||
call memacct('ia',imax(1)+1,1)
|
||||
call memacct('iu',imax(1),1)
|
||||
call memacct('ip',imax(1),1)
|
||||
call memacct('icg',0,1)
|
||||
call memacct('ifg',0,1)
|
||||
call memacct('ib',0,1)
|
||||
call memacct('b',0,1)
|
||||
call memacct('jb',0,1)
|
||||
c
|
||||
c load ifg into ifc
|
||||
c
|
||||
c do 10 i=1,imax(1)
|
||||
@ -102,17 +76,6 @@ c
|
||||
call trunc(1,imin,imax,a,ia,ja)
|
||||
c
|
||||
call symm(1,1,imin,imax,a,ia,ja,icg,ifg)
|
||||
c
|
||||
c set fine level data usage (permanent)
|
||||
c
|
||||
call memacct('a',ia(imax(1)+1),1)
|
||||
call memacct('ja',ia(imax(1)+1),1)
|
||||
call memacct('ia',imax(1)+1,1)
|
||||
call memacct('iu',imax(1),1)
|
||||
call memacct('ip',imax(1),1)
|
||||
call memacct('ib',imax(1)+1,1)
|
||||
c
|
||||
call msample
|
||||
c
|
||||
if(levels.le.1) then
|
||||
close(6)
|
||||
@ -140,7 +103,6 @@ cveh call stats(k-1,levels,idump,
|
||||
cveh * nun,imin,imax,a,ia,ja,iu,ip,icg,ifg,
|
||||
cveh * b,ib,jb,ipmn,ipmx,iv,xp,yp)
|
||||
|
||||
c call memoryx(k-1,imin,imax,ia,ib,ipmn,ipmx)
|
||||
c
|
||||
c=> test for coarsest grid
|
||||
c
|
||||
@ -156,20 +118,6 @@ c
|
||||
call trunc(k,imin,imax,a,ia,ja)
|
||||
|
||||
call symm(k,isymm,imin,imax,a,ia,ja,icg,ifg)
|
||||
c
|
||||
c set final data usage (permanent)
|
||||
c
|
||||
call memacct('a',ia(imax(k)+1),1)
|
||||
call memacct('ja',ia(imax(k)+1),1)
|
||||
call memacct('ia',imax(k)+1,1)
|
||||
call memacct('iu',imax(k),1)
|
||||
call memacct('ip',imax(k),1)
|
||||
call memacct('icg',imax(k-1),1)
|
||||
call memacct('ib',imax(k)+1,1)
|
||||
call memacct('b',ib(imax(k-1)+1),1)
|
||||
call memacct('jb',ib(imax(k-1)+1),1)
|
||||
c
|
||||
call msample
|
||||
c
|
||||
go to 20
|
||||
c
|
||||
@ -182,7 +130,6 @@ c compute & print statistics after coarsening
|
||||
* b,ib,jb,ipmn,ipmx,iv,xp,yp)
|
||||
|
||||
|
||||
call msample
|
||||
c
|
||||
c test operators
|
||||
c
|
||||
@ -195,7 +142,6 @@ c=> compress unneeded entries from interpolation/restriction
|
||||
c
|
||||
c call crushb(levels,imin,imax,icg,b,ib,jb)
|
||||
c
|
||||
call mreport(1)
|
||||
|
||||
close(6)
|
||||
|
||||
|
||||
@ -375,12 +375,6 @@ c
|
||||
do 30 i=i0lo,i0hi
|
||||
ja(ia(i))=i
|
||||
30 continue
|
||||
c
|
||||
c set temp memory usage
|
||||
c
|
||||
call memacct('a',ka,0)
|
||||
call memacct('ja',ka,0)
|
||||
call memacct('ia',i0hi+ishift+1,0)
|
||||
c
|
||||
return
|
||||
c
|
||||
|
||||
Loading…
Reference in New Issue
Block a user