anchor.diff:        ENDDO   
         
         ! D2_53, D2_63, D2_64 
  -      p%D2_63 = MATMUL( p%PhiM, p%MMB )
  +      !p%D2_63 = p%PhiRb_TI - MATMUL( p%PhiM, p%MMB ) 
  +      CALL LAPACK_GEMM( 'N', 'N', 1.0_ReKi, p%PhiM, p%MMB, 0.0_ReKi, p%D2_63, ErrStat2, ErrMsg2 ); if(Failed()) return;
         p%D2_63 = p%PhiRb_TI - p%D2_63
   
         !p%D2_64 = MATMUL( p%PhiM, p%PhiM_T )
  
positive.code:       p%D1_13 = p%MBB - p%D1_13
  
        ! TODO cant use LAPACK due to type conversions FEKi->ReKi
        !p%D1_14 = MATMUL( p%MBM, p%PhiM_T ) - MATMUL( TI_transpose, TRANSPOSE(PHiRb))  
        p%D1_14 = MATMUL( TI_transpose, TRANSPOSE(PHiRb)) 
        !CALL LAPACK_GEMM( 'T', 'T', 1.0_ReKi, p%TI,   PHiRb,  0.0_ReKi, p%D1_14, ErrStat2, ErrMsg2 ); if(Failed()) return   ! p%D1_14 = MATMUL( TRANSPOSE(TI), TRANSPOSE(PHiRb))  
        CALL LAPACK_GEMM( 'N', 'T', 1.0_ReKi, p%MBM, p%PhiM, -1.0_ReKi, p%D1_14, ErrStat2, ErrMsg2 ); if(Failed()) return   ! p%D1_14 = MATMUL( p%MBM, TRANSPOSE(p%PhiM) ) - p%D1_14 
  
     
  
negative.code:    !p%D1_15=-TI_transpose  !this is 6x6NIN
     IF ( p%nDOFM > 0 ) THEN ! These values don't exist for nDOFM=0; i.e., p%nDOFM == 0
        ! TODO cant use LAPACK due to type conversions FEKi->ReKi
        p%MBM = MATMUL( TI_transpose, MBmb )  ! NOTE: type conversion
        !CALL LAPACK_gemm( 'T', 'N', 1.0_ReKi, p%TI, MBmb, 0.0_ReKi, p%MBM, ErrStat2, ErrMsg2); if(Failed()) return
        !p%CBM = MATMUL( TRANSPOSE(p%TI), CBMb )    != CBMt
        !CALL LAPACK_gemm( 'T', 'N', 1.0_ReKi, p%TI, CBMb, 0.0_ReKi, p%CBM, ErrStat2, ErrMsg2); if (Failed()) return
  