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))