anchor.diff:           MBM = MATMUL( MRL, PhiL(:,1:nM) ) + MBM    !This had PhiM      
         ENDIF
         
  -      KBB = MATMUL(KRL, PhiR)   
  +      !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
  
positive.code:       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 )
        DEALLOCATE(Temp)
        ! PhiL = MATMUL( PhiL, MU2 )  !this is the nondimensionalization (MU2 is diagonal)   
        DO I = 1, nM_out
  
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
  