anchor.diff:                 temp6(:) = 0.0D0
                  temp6(1:3) = y%BldMotion%TranslationAcc(1:3,temp_id2)
                  temp6(4:6) = y%BldMotion%RotationAcc(1:3,temp_id2)
  -               temp6(:) = MATMUL(MoTens,temp6)
  +               temp6(1:3) = MATMUL(TRANSPOSE(p%GlbRot),temp6(1:3))
  +               temp6(4:6) = MATMUL(TRANSPOSE(p%GlbRot),temp6(4:6))
              ENDIF
              OtherState%Acc(temp_id+1:temp_id+6) = temp6(1:6)
          ENDDO
  
positive.code:    SUBROUTINE BD_GenerateDynamicElementGA2(uuN0,uuNf,vvNf,aaNf,                 &
                                             Stif0,Mass0,gravity,u,damp_flag,beta,&
                                             elem_total,node_elem,dof_node,ngp,   &
                                             StifK,RHS,MassM,DampG)
  
     REAL(ReKi),        INTENT(IN   ):: uuN0(:,:)
     REAL(ReKi),        INTENT(IN   ):: uuNf(:)
  
negative.code:    Eref = 0.0D0
  
     DO i=1,niter
         CALL BD_GenerateDynamicElementGA2(uuN0,uuNf,vvNf,aaNf,                 &
                                           Stif0,Mass0,gravity,u,damp_flag,beta,&
                                           elem_total,node_elem,dof_node,ngp,   &
                                           StifK,RHS,MassM,DampG,&
                                           ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         StifK = MassM + coef(7) * DampG + coef(8) * StifK
         DO j=1,node_total
             temp_id = (j-1)*dof_node
  