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

Go to the documentation of this file.

subroutine picard

use global

use mf

use model

implicit none


integer picardIter
integer i,j


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

integer,allocatable :: picardSuccessVec(:)

real(8),allocatable :: myV(:)

real(8),allocatable :: compareVec(:)

!real(8) :: damping=0.005
!real(8),allocatable :: damping(:)


real(8),allocatable :: oldVec(:)

real(8) picardTol

real(8) picardNorm


real(8) myP,myD



do i=1,5*5*2


        !y(i)=i

end do





picardTol=1e-5




!Open (4400,File='x.dat',Status='unknown')  !<--- physics and domain and subdomain based numbering
        
        
        
!910 FORMAT((E18.10))


!read(4400,910)y
!close(4400)




!IM=>myModel%physics(currentPhysics)%domain(currentDomain)%dofMatrix

dof=>myModel%physics(currentPhysics)%dof  ! later will be replaced with picardDof to extract max dof of picard with pVec,dVEc


allocate(picardSuccessVec(dof))


damping=>myModel%physics(currentPhysics)%damping

!setDof=1 ! picard solves only for one unk


pVec=>myModel%pVec

do picardIter=1,maxPicardIter

        print*,'#########################'
        print*,'picard Iter:',picardIter

        picardSuccessVec=0


        !if (mod(picardIter,100).eq.0 ) then

                !damping(2)=damping(2)*2.d0

        !end if

                
        !do i=1,picardDof
        do i=1,dof

                print*,'-----------'
                !print*,'picardDof:',i
                print*,'dof:',i
                print*,'        '

                currentDof=i
                
                
                do myP=1,size(pVec)

                        currentPhysics=pVEc(myP)

                        dVec=>myModel%physics(currentPhysics)%dVec

                        do myD=1,size(dVec)

                                currentDomain=dVec(myD)

                                IM=>myModel%physics(currentPhysics)%domain(currentDomain)%dofMatrix
                                
                                IM(:,:,:,flag)=0
                                IM(:,:,i,flag)=1

                        end do

                end do


                call currentSet


                !print*,setVec
                !stop





                !if (damping(i).ne.1.d0) then

                allocate(oldVec(globalCurrentSetLength))
                
                oldVec=y(setVec)

                !end if
                
                
                
                !print*,globalCurrentSetLength
                !stop
                
                !print*,setvec

                !allocate(myV(globalCurrentSetLength))
                !myV=1.d0


                !myV=matvec(myV)

                !do j=1,5*5*2


                !       print*,j,myV(j)

                !end do


                !stop


                !print*,myModel%physics(currentPhysics)%domain(currentDomain)%setDof
                !stop


                
                call Newton

                !print*,size(setVec),size(oldVec)

                !pause

                picardNorm=dsqrt(dot_product(y(setVec)-oldVec,y(setVec)-oldVec)/globalcurrentSetLength)
                
                print*,'picardNorm=',picardNorm

                if (picardNorm.le.picardTol) then

                        !picardSuccessVec(i)=newtonSuccess ! if newton return with no iterations
                        picardSuccessVec(i)=1

                end if



                !if (damping(i).ne.1.d0) then 
                
                y(setVec)=damping(i)*y(setVec)+oldVec*(1.d0-damping(i))

                deallocate(oldVec)

                !end if


                !print*,setVec
                !stop


                !call tecFD

                !print*,y
                !stop



                

        end do

        call postProcessor


        !pause
        
        
        !stop


        if (minval(picardSuccessVec).eq.1) then

                print*,'picard converged at iter:',picardIter
                exit

        end if

        


end do




end subroutine picard

---