Commit f8684306 authored by Stefan's avatar Stefan
Browse files

add pressure:solver to par

parent 633d1d55
......@@ -4,7 +4,7 @@ c Note: Keys have to be in captial letters
c
integer PARDICT_NKEYS
parameter(PARDICT_NKEYS = 100)
parameter(PARDICT_NKEYS = 101)
character*132 pardictkey(PARDICT_NKEYS)
data
......@@ -105,3 +105,4 @@ c
& pardictkey(98) / 'GENERAL:MAXNUMELEMENTS' /
& pardictkey(99) / 'GENERAL:MAXNUMHISTORYPOINTS' /
& pardictkey(100)/ 'GENERAL:FILTERMODES' /
& pardictkey(101)/ 'PRESSURE:SOLVER' /
......@@ -436,6 +436,20 @@ c set parameters
call finiparser_getDbl(d_out,'mesh:firstBCFieldIndex',ifnd)
if(ifnd .eq. 1) param(33) = int(d_out)
call finiparser_getString(c_out,'pressure:solver',ifnd)
if (ifnd .eq. 1) then
call capit(c_out,132)
if (index(c_out,'GMRES') .eq. 1) then
param(42) = 0
else if (index(c_out,'CGFLEX') .eq. 1) then
param(42) = 2
else
write(6,*) 'value: ',trim(c_out)
write(6,*) 'not supported for pressure:solver!'
goto 999
endif
endif
call finiparser_getString(c_out,'pressure:preconditioner',ifnd)
if (ifnd .eq. 1) then
call capit(c_out,132)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment