D:/demona/2.0/manual/doc/2-31set/problems/test/1/solver/modelSolver.f90

Go to the documentation of this file.

subroutine modelSolver

use global

use model

!use mf

implicit none


integer i

integer myP,myD

real(8) myEps

real(8),allocatable :: testVec(:)






print*,pVec




! here we will form our first template solver.
! dvec's will be created considering the physics of the problem
! if the problem is partitioned in domains
! then 






if (demonaSolver.eq.dsNewton) then



        print*,'newton is called'
        
        
        call newton
        
        
        
        !



end if





stop



return











stop




do myP=1,size(pVec)

        currentPhysics=pVec(myP)

        allocate(myModel%physics(myP)%dVec(1))

        dVec=>myModel%physics(myP)%dVec

        dVec(1)=1



end do



do myD=1,size(dVec)

        currentDomain=dVec(myD)

end do



print*,'modelsolver set'


!< should be called for all pvec and dvec in the set   <---------------------
call set

!print*,currentSetlength



allocate(testVec(currentSetlength))

testVec=0.d0


do i=1,currentSetlength

        testVec(i)=i


end do


myEps=dsqrt(epsilon(0.d0))



!< use setVec instead. setVec is combination of local setVecs   <---------------------

tempy(myModel%physics(currentPhysics)%domain(currentDomain)%setVec)=tempy(myModel%physics(currentPhysics)%domain(currentDomain)%setVec)+myeps*testVec
!tempy(setVec)=tempy(setVec)+myeps*testVec

!print*,tempy



!print*,discrete(0)
!print*,discrete(1)

print*,(discrete(1)-discrete(0))/myeps

!print*,matVec(testVec)



end subroutine modelSolver

---