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

Go to the documentation of this file.

!subroutine setRule(resVec,set,aveNorm)
subroutine setRule(resVec,norm,normRMS)


use global


implicit none


real(8) resVec(globalcurrentsetlength)
integer set(globalcurrentsetlength)
integer i

real(8) myval,norm,normRMS

integer myT,myI,myJ,myII

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

integer dummy

dof=>myModel%physics(currentPhysics)%dof



if (allocated(setTol)) deallocate(setTol)

allocate(setTol(dof))






!setTol=1e-2
!setTol=min(1e-2,1e-2*normF)
setTol=norm*1e-2
!setTol(2)=1.d0


setTol=1e-8
!setTol=0.d0

!print*,(norm/normRMS)


!setTOL=(normRMS/(norm/normRMS))**1.d0

!print*,setTol
!stop

!print*,resvec(1)


!if (resvec(1).eq.settol(1)) stop




myDofVectorArray=>myModel%physics(currentPhysics)%domain(currentDomain)%dofVector

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

dummy=0

do i=1,globalcurrentsetlength

        !print*,'********'


        !myII=setVec(i)
        myII=i

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

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


        if (dabs(myVal).gt.setTol(myT)) then

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

                dummy=dummy+1

        else

                !print*,'++'

                !print*,myI,myJ,myT,setVEc(i)

                IM(myI,myJ,myT,flag)=0

        end if


end do


!print*,'values kept=',dummy


!pause








end subroutine setRule

---