demona: D:/demona/2.0/manual/doc/2-31set/problems/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 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=100
00047 
00048                 call RunNewton
00049 
00050 
00051 
00052         case(2)
00053 
00054                 inexactNewton=off
00055 
00056                 maxNewtonIter=5
00057 
00058                 maxPicardIter=1000
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 print*,' '
00096 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00097 print*,'MODEL SOLVER ENDS'
00098 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00099 print*,'[][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][][]'
00100 print*,' '
00101 
00102 
00103 return
00104 
00105 
00106 
00107 
00108 do myP=1,size(pVec)
00109 
00110         currentPhysics=pVec(myP)
00111 
00112         allocate(myModel%physics(myP)%dVec(1))
00113 
00114         dVec=>myModel%physics(myP)%dVec
00115 
00116         dVec(1)=1
00117 
00118 
00119 
00120 end do
00121 
00122 
00123 
00124 do myD=1,size(dVec)
00125 
00126         currentDomain=dVec(myD)
00127 
00128 end do
00129 
00130 
00131 
00132 print*,'modelsolver set'
00133 
00134 
00135 
00136 
00137 
00138 
00139 
00140 end subroutine modelSolver