CFD/hypre
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

placed tolerance into code

Browse Source
This commit is contained in:
vhenson 1997-02-12 22:37:21 +00:00
parent f420af2741
commit c18a360965
2 changed files with 18 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
seq_linear_solvers/amg/amg_fdrive/amg_fdrive.f
View File

@ -10,7 +10,6 @@ C@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C
IMPLICIT REAL*8 (A-H,O-Z)
C
parameter(ndima=600000,ndimu=250000)
DIMENSION A(ndima),JA(ndima),IA(ndimu),U(ndimu),F(ndimu)
CHARACTER*15 FYSMP,FRHS,FINITU
@ -21,14 +20,14 @@ C READ INPUT FILES WITH MATRIX, RIGHT-HAND SIDE, PROBLEM SPECS,
C INITIAL GUESS.
C
FYSMP = 'AMG2.in.ysmp'
FRHS = 'AMG2.in.rhs'
FINITU='AMG2.in.initu'
FYSMP = 'AMG.in.ysmp'
FRHS = 'AMG.in.rhs'
FINITU='AMG.in.initu'
CVEH
C READ IN MATRIX FROM AMG2.in.ysmp
C READ IN MATRIX FROM AMG.in.ysmp
CVEH
open (8,FILE=FYSMP,STATUS='OLD')
read (8,*) junk
@ -41,7 +40,7 @@ CVEH
CVEH
C READ IN RHS FROM AMG2.in.rhs
C READ IN RHS FROM AMG.in.rhs
CVEH
open (8,FILE=FRHS,STATUS='OLD')
read (8,*) junk
@ -56,7 +55,7 @@ CVEH
CVEH
C READ IN INITIAL GUESS FROM AMG2.in.initu
C READ IN INITIAL GUESS FROM AMG.in.initu
CVEH
open (8,FILE=FINITU,STATUS='OLD')
read (8,*) junk
@ -69,10 +68,12 @@ CVEH
read (8,*) (u(j), j=1,nv)
close (8)
tol = 1.0e-7
call amg_Initialize(data,0)
call amg_SetLogging(3, "AMG2.runlog", data)
call amg_SetLogging(3, "AMG.runlog", data)
call amg_Setup(isterr,a,ia,ja,nv,data)
print *, 'Setup error flag = ', isterr

17
seq_ls/amg/amg_fdrive/amg_fdrive.f
View File

@ -10,7 +10,6 @@ C@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C
IMPLICIT REAL*8 (A-H,O-Z)
C
parameter(ndima=600000,ndimu=250000)
DIMENSION A(ndima),JA(ndima),IA(ndimu),U(ndimu),F(ndimu)
CHARACTER*15 FYSMP,FRHS,FINITU
@ -21,14 +20,14 @@ C READ INPUT FILES WITH MATRIX, RIGHT-HAND SIDE, PROBLEM SPECS,
C INITIAL GUESS.
C
FYSMP = 'AMG2.in.ysmp'
FRHS = 'AMG2.in.rhs'
FINITU='AMG2.in.initu'
FYSMP = 'AMG.in.ysmp'
FRHS = 'AMG.in.rhs'
FINITU='AMG.in.initu'
CVEH
C READ IN MATRIX FROM AMG2.in.ysmp
C READ IN MATRIX FROM AMG.in.ysmp
CVEH
open (8,FILE=FYSMP,STATUS='OLD')
read (8,*) junk
@ -41,7 +40,7 @@ CVEH
CVEH
C READ IN RHS FROM AMG2.in.rhs
C READ IN RHS FROM AMG.in.rhs
CVEH
open (8,FILE=FRHS,STATUS='OLD')
read (8,*) junk
@ -56,7 +55,7 @@ CVEH
CVEH
C READ IN INITIAL GUESS FROM AMG2.in.initu
C READ IN INITIAL GUESS FROM AMG.in.initu
CVEH
open (8,FILE=FINITU,STATUS='OLD')
read (8,*) junk
@ -69,10 +68,12 @@ CVEH
read (8,*) (u(j), j=1,nv)
close (8)
tol = 1.0e-7
call amg_Initialize(data,0)
call amg_SetLogging(3, "AMG2.runlog", data)
call amg_SetLogging(3, "AMG.runlog", data)
call amg_Setup(isterr,a,ia,ja,nv,data)
print *, 'Setup error flag = ', isterr