demona: D:/demona/2.0/manual/doc/2-31set/problems/nc_sf_vor/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 
00032 call solverParameters()
00033 
00034 
00035 
00036 
00037 select case(mySolver)
00038 
00039 
00040         case(0)
00041 
00042                 call RunHybrid
00043 
00044 
00045         case(1)
00046 
00047                 inexactNewton=on
00048 
00049                 maxNewtonIter=20
00050 
00051                 call RunNewton
00052 
00053 
00054 
00055         case(2)
00056 
00057                 inexactNewton=off
00058 
00059                 maxNewtonIter=5
00060 
00061                 maxPicardIter=1000
00062 
00063                 call RunPicard
00064 
00065 
00066 
00067 
00068         case(3)
00069 
00070                 inexactNewton=off
00071 
00072                 maxNewtonIter=5
00073 
00074                 call gs
00075 
00076         case(4)
00077 
00078                 inexactNewton=on
00079 
00080                 maxNewtonIter=20
00081 
00082                 call RunMS
00083 
00084 
00085 
00086 
00087         case default
00088         print*,'invalid mySolver'
00089         print*,'demona stop in modelSolver!!'
00090         stop
00091         
00092 
00093 
00094 
00095 end select 
00096 
00097 
00098 
00099 print*,'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
00100 print*,'lambda=0.5'
00101 print*,'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
00102 
00103 
00104 
00105 print*,' '
00106 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00107 print*,'MODEL SOLVER ENDS'
00108 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00109 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00110 print*,' '
00111 
00112 
00113 return
00114 
00115 
00116 
00117 
00118 
00119 
00120 
00121 
00122 
00123 do myP=1,size(pVec)
00124 
00125         currentPhysics=pVec(myP)
00126 
00127         allocate(myModel%physics(myP)%dVec(1))
00128 
00129         dVec=>myModel%physics(myP)%dVec
00130 
00131         dVec(1)=1
00132 
00133 
00134 
00135 end do
00136 
00137 
00138 
00139 do myD=1,size(dVec)
00140 
00141         currentDomain=dVec(myD)
00142 
00143 end do
00144 
00145 
00146 
00147 print*,'modelsolver set'
00148 
00149 
00150 
00151 
00152 
00153 
00154 
00155 end subroutine modelSolver