def test_computeRJF_part0():
"""Test computeRJF() with random 3x3 ... 10x10 matrices.
"""
np.random.seed(1234)
for d in range(3,10):
A = np.random.random((d,d))
J,T,_ = LinearDS.computeRJF(A)
Ac = T.dot(A).dot(np.linalg.inv(T))
np.testing.assert_almost_equal(np.linalg.norm(J-Ac, 'fro'), 0)
def test_computeRJF_part0():
"""Test computeRJF() with random 3x3 ... 10x10 matrices.
"""
np.random.seed(1234)
for d in range(3,10):
A = np.random.random((d,d))
J,T,_ = LinearDS.computeRJF(A)
Ac = T.dot(A).dot(np.linalg.inv(T))
np.testing.assert_almost_equal(np.linalg.norm(J-Ac, 'fro'), 0, 2)
def test_computeRJF_part1():
"""Test computeRJF() with MATLAB example.
"""
# MATLAB test example, see 'doc jordan' in MATLAB
A = np.asarray([[1,-3,-2],
[-1,1,-1],
[2,4,5]])
J,T,_ = LinearDS.computeRJF(A)
ref = np.asarray([[3,0,0],
[0,2,0],
[0,0,2]])
np.testing.assert_almost_equal(np.linalg.norm(J-ref,'fro'), 0)
def test_computeRJF_part1():
"""Test computeRJF() with MATLAB example.
"""
# MATLAB test example, see 'doc jordan' in MATLAB
A = np.asarray([[1,-3,-2],
[-1,1,-1],
[2,4,5]])
J,T,_ = LinearDS.computeRJF(A)
ref = np.asarray([[3,0,0],
[0,2,0],
[0,0,2]])
np.testing.assert_almost_equal(np.linalg.norm(J-ref,'fro'), 0, 2)
