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