anchor.diff:        CALL AllocAry( Temp , nL    , nL , 'Temp' , ErrStat2 , ErrMsg2); if(Failed()) return
         CALL AllocAry( MU   , nM_out, nL , 'Mu'   , ErrStat2 , ErrMsg2); if(Failed()) return
         MU   = TRANSPOSE(PhiL)
  -      Temp = MATMUL( MU, MLL )
  -      MU(1:nM_Out,1:nM_Out) = MATMUL( Temp, PhiL )
  +      !Temp = MATMUL( MU, MLL )
  +      CALL LAPACK_gemm( 'N', 'N', 1.0_FeKi, MU, MLL, 0.0_FeKi, Temp  , ErrStat2, ErrMsg2); if(Failed()) return
  +      !MU(1:nM_Out,1:nM_Out) = MATMUL( Temp, PhiL )
  +      CALL LAPACK_gemm( 'N', 'N', 1.0_FeKi, Temp, PhiL, 0.0_FeKi, MU  , ErrStat2, ErrMsg2); if(Failed()) return
         DEALLOCATE(Temp)
         ! PhiL = MATMUL( PhiL, MU2 )  !this is the nondimensionalization (MU2 is diagonal)   
         DO I = 1, nM_out
  
positive.code:          MBM = MATMUL( MRL, PhiL(:,1:nM) ) + MBM    !This had PhiM      
        ENDIF
        
        KBB = MATMUL(KRL, PhiR)   
        KBB = KBB + KRR
  
        if (present(CRR)) then
  
negative.code: SUBROUTINE CraigBamptonReduction_FromPartition( MRR, MLL, MRL, KRR, KLL, KRL, nR, nL, nM, nM_Out,&
                       MBB, MBM, KBB, PhiL, PhiR, OmegaL, ErrStat, ErrMsg,&
                       CRR, CLL, CRL, CBB, CBM, CMM)
     USE NWTC_LAPACK, only: LAPACK_getrs, LAPACK_getrf, LAPACK_gemm
     INTEGER(IntKi),         INTENT(  in)  :: nR
     INTEGER(IntKi),         INTENT(  in)  :: nL
     INTEGER(IntKi),         INTENT(  in)  :: nM_Out
  