def test_multiexp(self):
A = multisym(rnd.randn(self.k, self.m, self.m))
e = np.zeros((self.k, self.m, self.m))
for i in range(self.k):
e[i] = expm(A[i])
np_testing.assert_allclose(multiexp(A, sym=True), e)
def test_multiexp(self):
A = multisym(rnd.randn(self.k, self.m, self.m))
e = np.zeros((self.k, self.m, self.m))
for i in range(self.k):
e[i] = expm(A[i])
np_testing.assert_allclose(multiexp(A, sym=True), e)
def test_multiexp_singlemat(self):
# A is a positive definite matrix
A = rnd.randn(self.m, self.m)
A = A + A.T
np_testing.assert_allclose(multiexp(A, sym=True), expm(A))
def test_multiexp_singlemat(self):
# A is a positive definite matrix
A = rnd.randn(self.m, self.m)
A = A + A.T
np_testing.assert_allclose(multiexp(A, sym=True), expm(A))