def rotateStress2D(stressmat,phi,degoption=False):
    A = Mmath.getRotMatA(phi,degoption)
    return np.einsum('kj,ij',A,stressmat)
def rotateStress2D(stressmat, phi, degoption=False):
    A = Mmath.getRotMatA(phi, degoption)
    return np.einsum('kj,ij', A, stressmat)
def rotateStrain2D(strainmat,phi,degoption=False):
    A = Mmath.getRotMatA(phi,degoption)
    R = Mmath.getR()
    RAR = np.einsum('ij,jk,kl',R,A,spla.inv(R))
    return np.einsum('kj,ij',RAR,strainmat)
def rotateStrain2D(strainmat, phi, degoption=False):
    A = Mmath.getRotMatA(phi, degoption)
    R = Mmath.getR()
    RAR = np.einsum('ij,jk,kl', R, A, spla.inv(R))
    return np.einsum('kj,ij', RAR, strainmat)