demona: D:/demona/2.0/manual/doc/2-31set/problems/bratu/1/solver/modelSolver.f90 Source File

00001 subroutine modelSolver
00002 
00003 use global
00004 
00005 use model
00006 
00007 !use mf
00008 
00009 implicit none
00010 
00011 
00012 integer i
00013 
00014 integer myP,myD
00015 
00016 real(8) myEps
00017 
00018 real(8),allocatable :: testVec(:)
00019 
00020 print*,' '
00021 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00022 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00023 print*,'MODEL SOLVER STARTS'
00024 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00025 print*,' '
00026 
00027 
00028 !recursive call of modelSolver ???
00029 
00030 
00031 call solverParameters !(mySolver)
00032 
00033 select case(mySolver)
00034 
00035 
00036         case(0)
00037 
00038                 call RunHybrid
00039 
00040 
00041         case(1)
00042 
00043                 inexactNewton=on
00044 
00045                 maxNewtonIter=300
00046 
00047                 call RunNewton
00048 
00049 
00050 
00051         case(2)
00052 
00053                 inexactNewton=off
00054 
00055                 maxNewtonIter=5
00056 
00057                 maxPicardIter=100
00058 
00059                 call RunPicard
00060 
00061 
00062 
00063 
00064         case(3)
00065 
00066                 inexactNewton=off
00067 
00068                 maxNewtonIter=5
00069 
00070                 call gs
00071 
00072         case(4)
00073 
00074                 inexactNewton=on
00075 
00076                 maxNewtonIter=200
00077 
00078                 call RunMS
00079 
00080 
00081 
00082 
00083         case default
00084         print*,'invalid mySolver'
00085         print*,'demona stop in modelSolver!!'
00086         stop
00087         
00088 
00089 
00090 
00091 end select 
00092 
00093 
00094 print*,' '
00095 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00096 print*,'MODEL SOLVER ENDS'
00097 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00098 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00099 print*,' '
00100 
00101 
00102 return
00103 
00104 
00105 
00106 
00107 do myP=1,size(pVec)
00108 
00109         currentPhysics=pVec(myP)
00110 
00111         allocate(myModel%physics(myP)%dVec(1))
00112 
00113         dVec=>myModel%physics(myP)%dVec
00114 
00115         dVec(1)=1
00116 
00117 
00118 
00119 end do
00120 
00121 
00122 
00123 do myD=1,size(dVec)
00124 
00125         currentDomain=dVec(myD)
00126 
00127 end do
00128 
00129 
00130 
00131 print*,'modelsolver set'
00132 
00133 
00134 
00135 
00136 
00137 
00138 
00139 end subroutine modelSolver