def METHOD_NAME():
    class Rel(sk.LagrangianScleronomousR):
        pass
    r = Rel()
    j = np.array([[1, 2, 3], [4, 5, 6]])
    r.setJachqPtr(j)
    # C is transposed()
    r.C()
    assert np.max(r.C() - np.array([[1, 2, 3], [4, 5, 6]])) == 0.0
    assert np.max(r.C() - np.array([[0, 2, 3], [4, 5, 6]])) == 1.0
    r.setJachqPtr(r.C())
    r.C()
    assert np.max(r.C() - np.array([[1, 2, 3], [4, 5, 6]])) == 0.0
    assert np.max(r.C() - np.array([[0, 2, 3], [4, 5, 6]])) == 1.0