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

Go to the documentation of this file.

subroutine symbolic(valsizeMax)

use global

use mf

use DFPORT

implicit none



integer indSize
integer ptrSize
integer valSize,valSizeMax

real(8), allocatable :: JacVal(:),luval(:)
integer, allocatable :: JacColInd(:),JacRowPtr(:),iw(:),uptr(:)


integer numberOfColors




!###########################################################################################################################
!detect numberOfColors
valsize=1
indSize=valSize
ptrSize=globalcurrentsetlength+1

valSizeMax=valsize

allocate(JacVal(valSizeMax),JacColInd(valSizeMax),JacRowPtr(ptrSize),uptr(globalcurrentsetlength))

call colorJacobianCRS(JacVal,JacColInd,JacRowPtr,globalcurrentsetlength,1,valSize,indSize,uptr,1,numberOfColors)


deallocate(JacVal,JacColInd,JacRowPtr,uptr)
!###########################################################################################################################
! calculate the colorbased jacobian
valsize=globalcurrentsetlength*numberOfColors ! limit of jacVal allocation
indSize=valSize
ptrSize=globalcurrentsetlength+1

valSizeMax=valsize


allocate(JacVal(valSizeMax),JacColInd(valSizeMax),JacRowPtr(ptrSize),uptr(globalcurrentsetlength))

jacVal=0.d0
JacColInd=0
JacRowPtr=0
uptr=0

tempy=y



allocate(yBackup(size(y)))
yBackup=0.d0
yBackup=y



call random_seed()
CALL RANDOM_NUMBER(y)

tempy=y


myModel%physics(currentPhysics)%evalOrder=>jacDiscOrder

call colorJacobianCRS(JacVal,JacColInd,JacRowPtr,globalcurrentsetlength,1,valSize,indSize,uptr,0,numberOfColors)
!###########################################################################################################################


myModel%physics(currentPhysics)%evalOrder=>myModel%physics(currentPhysics)%discOrder


y=yBackup
tempy=y
deallocate(yBackup)



deallocate(JacVal,JacColInd,JacRowPtr,uptr)
!###########################################################################################################################
! calculate the colorbased jacobian
!valsize=globalcurrentsetlength*numberOfColors ! limit of jacVal allocation
indSize=valSize
ptrSize=globalcurrentsetlength+1

valSizeMax=valsize
valSizeMax=floor(valsize*jacSF)

!all of these can be perfomed in symbolic




end subroutine symbolic

---