00001 subroutine RunHybrid
00002
00003 use global
00004
00005 implicit none
00006
00007 integer myP,myD
00008
00009 integer,pointer :: IM(:,:,:,:)
00010
00011 print*,'HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH'
00012 print*,'HYBRID STARTS FOR',myHybrid
00013 print*,' '
00014
00015 select case(myHybrid)
00016
00017
00018 case(1)
00019
00020 inexactNewton=off
00021
00022 maxNewtonIter=5
00023
00024 maxPicardIter=100
00025
00026 call RunPicard
00027
00028
00029 call globalSet
00030
00031 inexactNewton=on
00032
00033 maxNewtonIter=20
00034
00035 call RunNewton
00036
00037 case(2)
00038
00039 call RunGS
00040
00041 call globalSet
00042
00043 call RunNewton
00044
00045
00046
00047
00048 case default
00049
00050 print*,'invalid myHybrid'
00051 print*,'demona stop in RunHybrid!!'
00052 stop
00053
00054
00055
00056 end select
00057
00058
00059 print*,' '
00060 print*,'HYBRID ENDS FOR',myHybrid
00061 print*,'HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH'
00062
00063
00064
00065
00066 end subroutine RunHybrid
