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

Go to the documentation of this file.

subroutine generateIndexMatrix

use global

implicit none



integer i,j,k,t

integer myI,myJ,myK,myT

integer sumI,sumJ,sumK

real(8) myVal

integer,pointer :: myNx,myNy



integer sumNodes,sumDof,sum

integer, pointer :: myNodeMatrix(:,:,:),myNodeVectorArray(:,:)
integer, pointer :: myDofMatrix(:,:,:,:),myDofVectorArray(:,:)

real(8), allocatable :: resVec(:)

!real(8), allocatable :: setTol(:) ! to be general,we can define a tolerance for each unk, if F is scaled, this might not be necessary






integer,pointer :: minI,minJ,maxI,maxJ,maxUnk,myLocalDof









!allocate(myNodeMatrix(fullNx,fullNy,2))
!allocate(myDofMatrix(fullNx,fullNy,4,myLocalDof))


myNodeMatrix=>myModel%physics(currentPhysics)%domain(currentDomain)%NodeMatrix
myDofMatrix=>myModel%physics(currentPhysics)%domain(currentDomain)%DofMatrix

myNodeMatrix=0
myDofMatrix=0

!-----------------------------------------------------------------------------------------
! declarations on the top will change for myModel struct
!-----------------------------------------------------------------------------------------

sumNodes=1
sumDof=1

! if set is on, perhaps first unk's should be called. but no, numbering is unk1,unk2,unk1+1,unk2+1 based
! later set can be generated based on the unk


! this is valid only for fd domains
! finite volume domains should be formed using minI maxI formation


minI=>myModel%physics(currentPhysics)%domain(currentDomain)%minI
maxI=>myModel%physics(currentPhysics)%domain(currentDomain)%maxI
minJ=>myModel%physics(currentPhysics)%domain(currentDomain)%minJ
maxJ=>myModel%physics(currentPhysics)%domain(currentDomain)%maxJ


myNx=>myModel%physics(currentPhysics)%domain(currentDomain)%nx
myNy=>myModel%physics(currentPhysics)%domain(currentDomain)%ny



myLocalDof=>myModel%physics(currentPhysics)%dof



if (myModel%physics(currentPhysics)%dFX.eq.dFV) print*,'skip function is neeeded for no-compatible FV mode'



                        do j=minJ,maxJ

                                do i=minI,maxI

                                        ! skip function here  ! should be user input in the model field
                                        ! if skip_is_yes cycle

                                        
                                

                                        do t=1,myLocalDof
                                                
                                                myDofMatrix(i,j,t,relative)=sumDof

                                                sumDof=sumDof+1


                                        end do

                                end do

                        end do


                        if ((poissonCompatibility.eqv.on).and.(myModel%physics(currentPhysics)%dFX.eq.dFV)) then

                                sumDof=1
                                myDofMatrix=0

                                        ! fill the dofMatrix
                                        do j=1,myNy-1

                                                do i=1,myNx-1

                                                        ! skip function here  ! should be user input in the model field
                                                        ! if skip_is_yes cycle
                                                

                                                        do t=1,myLocalDof
                                                                
                                                                myDofMatrix(i,j,t,relative)=sumDof

                                                                sumDof=sumDof+1


                                                        end do

                                                end do

                                        end do

                                        if (myModel%physics(currentPhysics)%domain(currentDomain)%numberOfGhostCells.gt.0) then

                                                ! first ghost cell domain
                                                
                                                ! bottom and top ghost cells
                                                do j=0,myNy,myNy

                                                        do i=1,myNx-1

                                                                do t=1,myLocalDof

                                                                        myDofMatrix(i,j,t,relative)=sumDof
                                                                        sumDof=sumDof+1

                                                                end do


                                                        end do

                                                end do



                                                ! left and right ghost cells
                                                do i=0,myNx,myNx

                                                        do j=1,myNy-1

                                                                do t=1,myLocalDof

                                                                        myDofMatrix(i,j,t,relative)=sumDof
                                                                        sumDof=sumDof+1

                                                                end do

                                                        end do

                                                end do


                                        end if


                                        if (myModel%physics(currentPhysics)%domain(currentDomain)%numberOfGhostCells.gt.1) then

                                                ! second ghost cell domain

                                                ! bottom and top ghost cells
                                                do j=-1,myNy+1,myNy+2

                                                        do i=1,myNx-1

                                                                do t=1,myLocalDof

                                                                        myDofMatrix(i,j,t,relative)=sumDof
                                                                        sumDof=sumDof+1

                                                                end do


                                                        end do

                                                end do



                                                ! left and right ghost cells
                                                do i=-1,myNx+1,myNx+2

                                                        do j=1,myNy-1

                                                                do t=1,myLocalDof
                                                                        myDofMatrix(i,j,t,relative)=sumDof
                                                                        sumDof=sumDof+1
                                                                end do


                                                        end do
                                                end do

                                        end if






                        end if


                        do j=1,myNy

                                do i=1,myNx

                                        ! skip function here  ! should be user input in the model field
                                        ! if skip_is_yes cycle


                                        !print*,sumNodes
                                
                                        myNodeMatrix(i,j,relative)=sumNodes


                                        sumNodes=sumNodes+1

                                end do

                        end do


!print*,'make the vectors absolute'

myModel%physics(currentPhysics)%domain(currentDomain)%totalNode=sumNodes-1

myModel%physics(currentPhysics)%domain(currentDomain)%totalDof=sumDof-1

myModel%physics(currentPhysics)%totalDof=myModel%physics(currentPhysics)%domain(currentDomain)%totalDof+myModel%physics(currentPhysics)%totalDof

myModel%physics(currentPhysics)%totalNode=myModel%physics(currentPhysics)%domain(currentDomain)%totalNode+myModel%physics(currentPhysics)%totalNode

! zero will be replaced with sum of all previouns unknowns.





myNodeMatrix(:,:,absolute)=myNodeMatrix(:,:,relative)+modelNode
myDofMatrix(:,:,:,absolute)=myDofMatrix(:,:,:,relative)+modelDof
!myNodeMatrix(:,:,absolute)=myNodeMatrix(:,:,relative)
!myDofMatrix(:,:,:,absolute)=myDofMatrix(:,:,:,relative)

myModel%physics(currentPhysics)%domain(currentDomain)%relativeIndex=modelDof
myModel%physics(currentPhysics)%domain(currentDomain)%relativeIndexNode=modelNode


modelDof=myModel%physics(currentPhysics)%domain(currentDomain)%totalDof+modelDof
modelNode=myModel%physics(currentPhysics)%domain(currentDomain)%totalNode+modelNode




! vectorization should be indexmatrix independent
! i.e. numbering independent.
! the size of the vector is already known
! so relative index will be read for the index number
! later, that index will be filled only. so entries
! will not be added one-by-one in order but as random entries



!do j=fullNy,1,-1
!       print*,myDofMatrix(:,j,1,relative)
!end do



return



print*,'now create a residual vector'

allocate(resVec(size(myDofVectorArray,dim=1)))

resVec=0.d0


resVec=0.4d0

resVec(1)=0.6d0
resVec(3)=0.6d0
resVec(11)=0.6d0
resVec(13)=0.6d0
resVec(15)=0.6d0
resVec(17)=0.6d0
resVec(25)=0.6d0
resVec(27)=0.6d0
resVec(43)=0.6d0
resVec(45)=0.6d0
resVec(47)=0.6d0

resVec(2)=2.1d0
resVec(10)=2.1d0
resVec(12)=2.1d0
resVec(16)=2.1d0
resVec(22)=2.1d0
resVec(24)=2.1d0
resVec(26)=2.1d0
resVec(38)=2.1d0


!print*,resVec

print*,'flag the values'


do i=1,size(resVec)


        myVal=resVec(i)
        myT=myDofVectorArray(i,4)


        if (myVal.ge.setTol(myT)) then

                myI=myDofVectorArray(i,2)
                myJ=myDofVectorArray(i,3)

                myDofMatrix(myI,myJ,myT,flag)=1

        end if


end do



do j=fullNy,1,-1

        print*,myDofMatrix(:,j,1,flag)


end do



print*,'--------------------'


do j=-1,-2

        print*,'woink'

end do



print*,'nodes are flagged'






return





! might use allocate reallocate instead of a double run of the same loop
! need a jvec and ivalues as construct





print*,'each time a set is solved, new set should be found using all values'
print*,'so we need a default set function'




end subroutine generateIndexMatrix

---