anchor.diff:           CALL AllocAry( Temp , nM, nL , 'Temp' , ErrStat2 , ErrMsg2); if(Failed()) return
            CALL AllocAry( MU   , nM, nL , 'Mu'   , ErrStat2 , ErrMsg2); if(Failed()) return
            MU  = TRANSPOSE(PhiL(1:nL, 1:nM))
  -         Temp = MATMUL( MU, CLL )
  -         CMM = MATMUL( Temp, PhiL(1:nL, 1:nM) )
  +         !Temp = MATMUL( MU, CLL )
  +         CALL LAPACK_gemm( 'N', 'N', 1.0_FeKi, MU, CLL, 0.0_FeKi, Temp  , ErrStat2, ErrMsg2); if(Failed()) return
  +         !   CMM = MATMUL( Temp, PhiL(1:nL, 1:nM) )
  +         CALL LAPACK_gemm( 'N', 'N', 1.0_FeKi, Temp, PhiL, 0.0_FeKi, CMM  , ErrStat2, ErrMsg2); if(Failed()) return
            DEALLOCATE(MU)
            DEALLOCATE(Temp)
         endif
  
positive.code:          
        PhiR_T_MLL = TRANSPOSE(PhiR)
        PhiR_T_MLL = MATMUL(PhiR_T_MLL, MLL)
        MBB = MATMUL(MRL, PhiR)
        MBB = MRR + MBB + TRANSPOSE( MBB ) + MATMUL( PhiR_T_MLL, PhiR )
           
        IF ( nM == 0) THEN
  
negative.code:          MBM = MATMUL( MRL, PhiL(:,1:nM) ) + MBM    !This had PhiM      
        ENDIF
        
        !KBB = MATMUL(KRL, PhiR)   
        CALL LAPACK_gemm( 'N', 'N', 1.0_FeKi, KRL, PhiR, 0.0_FeKi, KBB  , ErrStat2, ErrMsg2); if(Failed()) return
        KBB = KBB + KRR
  
        if (present(CRR)) then
  