anchor.diff:        if( nspace==2 ) then
           XJ(1:2,1:2)=MATMUL( elecoord(1:2,1:nn), deriv(1:nn,1:2) )
           DET=XJ(1,1)*XJ(2,2)-XJ(2,1)*XJ(1,2)
  -        if( det==0.d0 ) stop "Math error in GetGlobalDeriv! Determinaut==0.0"
  +        if( det==0.d0 ) stop "Math error in GetGlobalDeriv! Determinant==0.0"
           DUM=1.d0/DET
           XJI(1,1)= XJ(2,2)*DUM
           XJI(1,2)=-XJ(1,2)*DUM
  
positive.code:            -jacobian(3,1)*jacobian(2,2)*jacobian(1,3)                                         &
             -jacobian(2,1)*jacobian(1,2)*jacobian(3,3)                                         &
             -jacobian(1,1)*jacobian(3,2)*jacobian(2,3)
          if( det==0.d0 ) stop "Math error in getJacobain! Determinaut==0.0"
  ! INVERSION OF JACOBIAN
          dum=1.d0/det
          inverse(1,1)=DUM*( jacobian(2,2)*jacobian(3,3)-jacobian(3,2)*jacobian(2,3) )
  
negative.code:            -XJ(3,1)*XJ(2,2)*XJ(1,3)                                             &
             -XJ(2,1)*XJ(1,2)*XJ(3,3)                                             &
             -XJ(1,1)*XJ(3,2)*XJ(2,3)
          if( det==0.d0 ) stop "Math error in GetGlobalDeriv! Determinant==0.0"
  ! INVERSION OF JACOBIAN
          DUM=1.d0/DET
          XJI(1,1)=DUM*( XJ(2,2)*XJ(3,3)-XJ(3,2)*XJ(2,3) )
  