D:/demona/2.0/manual/doc/2-31set/processor/processor.f90

Go to the documentation of this file.

subroutine processor

use global
use model
use myFunctions
!use processorParameters

implicit none


integer i

integer myP,myD

integer dummy,dummySkip

integer startIndex,endIndex


integer realTimeIter
integer physicsIter
integer physicsGroup
integer pseudoTimeIter
integer domainIter
integer domainGroup


integer,pointer :: IM(:,:,:,:)


real(8) ptNorm,ptNormPrev,ptNormIni
real(8) serNormCurrent,serNormPrevious

real(8) dTauTest,dtauBackup

!switching between the solvers
! solver rule!!! fas-fmg-pseudotransiet2newton



!################################################
!< create physics on the set
allocate(mymodel%pVec(numberOfPhysics))  !<-------------------------!!!!!!!!!!!
!allocate(mymodel%pVec(1))
!allocate(mymodel%pVec(2))
!allocate(mymodel%pVec(3))


pVec=>mymodel%pVec



!do i=1,numberOfPhysics
do i=1,size(pVec)
        pVec(i)=i
end do

myModel%endIndex=0



!################################################
!################################################
!< create domains on the set
do myP=1,size(pVec)

        currentPhysics=pVec(myP)
        
        myModel%endIndex(currentPhysics)=1
        myModel%physics(currentPhysics)%endIndex=0


        numberOfDomains=>mymodel%physics(currentPhysics)%numberOfdomains

        
        allocate(mymodel%physics(currentPhysics)%dVec(numberOfDomains))        !<-------------------------!!!!!!!!!!!
        !allocate(mymodel%physics(currentPhysics)%dVec(1))



        dVec=>mymodel%physics(currentPhysics)%dVec


        do i=1,size(dVec)

                dVec(i)=i

                myModel%physics(currentPhysics)%endIndex(i)=1



        end do



end do
!################################################
!< create physics on the set


!print*,mymodel%physics(1)%dVec
910 FORMAT((E18.10))

911 FORMAT((I6,1X,F20.12))


if (readInitialGuess.eqv.yes) then

        Open (4400,File='x.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
        !Open (4400,File='x1000.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
        !Open (4400,File='x1.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
        !Open (4400,File='x100.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
                
        read(4400,910)yPTO

        close(4400)

end if




call currentSet


myModel%physics(1)%physicalParameters(1)=myModel%physics(1)%physicalParametersFinal(1)

call continuation

call postprocessor
!stop

!test=1

!call modelSolver


!second model solver call!!!


!print*,'press any key to start newton'
!pause
!mySolver=1



y=yPTO


Open (4402,File='normPT.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
Open (4404,File='dtau.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering

dt=0.0001d0

dtau=0.01d0

!dtau=0.01d0

gssolver=off



tempy=yPTO




yPTO=y ! so zero norm for the time derivatives

serNormCurrent=dsqrt(dot_product(discrete(0),discrete(0)))

yPTO=tempy

dummyskip=0

serNormPrevious=serNormCurrent


call solverParameters


ptNormPrev=1.d0
ptNormIni=1.d0


do i=1,maxPseudoTimeIter
        
        print*,'==================i',i
        newtonsuccess=0
        call modelSolver
        !stop

        yPTO=y-yPTO
        ptNorm=dsqrt(dot_product(yPTO,yPTO)) !/globalCurrentSetLength

        write(4402,911)i,ptNorm
        write(4404,911)i,dtau

        print*,'pseudo time norm=',ptNorm

        if (i.eq.1) ptNormIni=ptNorm

        yPTO=y ! both initialization and zero time derivative residual for serNormCurrent

        serNormCurrent=dsqrt(dot_product(discrete(0),discrete(0)))
        
        !print*,'==============================='
        !print*,serNormPrevious
        !print*,serNormCurrent
        !print*,serNormPrevious/serNormCurrent
        !print*,'==============================='       

        print*,'old dtau=',dtau
        
        !dtau=ser(serNormPrevious,serNormCurrent,dtau,1.d0)
        !dtau=ser(serNormPrevious,serNormCurrent,dtau,(serNormCurrent/serNormPrevious)**20.d0)

        !dtau=ser(ptNormPrev,ptNorm,dtau,1.d0)

        ptNormPrev=ptNorm

        print*,'new dtau=',dtau


        !if (ptNorm.le.1e-5) then
        if (ptNorm/ptNormIni.le.1e-5) then
        

                print*,'pt converged'
                exit

        end if

        serNormPrevious=serNormCurrent


        if (mod(i,10).eq.0) call postprocessor
        !pause



end do

close(4402)
close(4404)






Open (4400,File='x.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
        
        
! later modify it with parameter success or physics or model success.
if (newtonSuccess.eq.1) WRITE(4400,910)y 

close(4400)


return



!call NewtonSet

!IM(:,:,:,flag)=1

call currentSet



call DBGcheckFlag(0)


call DBGsqIni(0)


call DBGcheckBC(0)


call DBGtestDiscrete(0,1)


!stop






return






do realTimeIter=1,maxRealTimeIter

        

        do physicsIter=1,maxPhysicsIter

                
                do physicsGroup=1,maxPhysicsGroup
                        
                        pVec=>myModel%pg(physicsGroup)%pVec

                        ! meaning physics on pVec will be solved together

                        if (allocated(myModel%pVec)) deallocate(myModel%pVec)

                        allocate(myModel%pVec(size(pVec)))

                        pVec=>myModel%pVec
                        pVec=myModel%pg(physicsGroup)%pVec




                        do pseudoTimeIter=1,maxPseudoTimeIter

                                
                                !< maxDomainIter might be modified with physicsrule

                                
                                call physicsRule



                                do domainIter=1,maxDomainIter


                                        do domainGroup=1,maxDomainGroup





                                                call modelSolver




                                        end do






                                end do





                                !!call modelSolver



                        end do



                end do




        end do





end do






!call modelSolver






!modelsolver
!physicssolver
!domainsolver.




end subroutine processor

---