anchor.diff:      real(kind=kreal) :: tempv(3)
       real(kind=kreal) :: temps
     
  -		call eigen3(tensor,eigval,princnormal)
  -	  
  +    call eigen3(tensor,eigval,princnormal)
  +    
       if (eigval(1)<eigval(2)) then
         temps=eigval(1)
         eigval(1)=eigval(2)
  
positive.code:     complex(kind=kreal):: x1,x2,x3
      real(kind=kreal):: rtemp
      real(kind=kreal) :: mat(3,4)   
  		integer :: i
      s11 = tensor(1)
      s22 = tensor(2)
      s33 = tensor(3)
  
negative.code:   end subroutine transformation
  
    SUBROUTINE get_principal (tensor, eigval, princmatrix)
    
      implicit none
      integer i,j
      real(kind=kreal) :: tensor(1:6)
  