def test_stacking(self):
a1 = ch.Ch(np.arange(10).reshape(2, 5))
b1 = ch.Ch(np.arange(20).reshape(4, 5))
c1 = ch.vstack((a1, b1))
c1_check = np.vstack((a1.r, b1.r))
residuals1 = (c1_check - c1.r).ravel()
a2 = ch.Ch(np.arange(10).reshape(5, 2))
b2 = ch.Ch(np.arange(20).reshape(5, 4))
c2 = ch.hstack((a2, b2))
c2_check = np.hstack((a2.r, b2.r))
residuals2 = (c2_check - c2.r).ravel()
self.assertFalse(np.any(residuals1))
self.assertFalse(np.any(residuals2))
d0 = ch.array(np.arange(60).reshape((10, 6)))
d1 = ch.vstack((d0[:4], d0[4:]))
d2 = ch.hstack((d1[:, :3], d1[:, 3:]))
tmp = d2.dr_wrt(d0).todense()
diff = tmp - np.eye(tmp.shape[0])
self.assertFalse(np.any(diff.ravel()))
def test_stacking(self):
a1 = ch.Ch(np.arange(10).reshape(2,5))
b1 = ch.Ch(np.arange(20).reshape(4,5))
c1 = ch.vstack((a1,b1))
c1_check = np.vstack((a1.r, b1.r))
residuals1 = (c1_check - c1.r).ravel()
a2 = ch.Ch(np.arange(10).reshape(5,2))
b2 = ch.Ch(np.arange(20).reshape(5,4))
c2 = ch.hstack((a2,b2))
c2_check = np.hstack((a2.r, b2.r))
residuals2 = (c2_check - c2.r).ravel()
self.assertFalse(np.any(residuals1))
self.assertFalse(np.any(residuals2))
d0 = ch.array(np.arange(60).reshape((10,6)))
d1 = ch.vstack((d0[:4], d0[4:]))
d2 = ch.hstack((d1[:,:3], d1[:,3:]))
tmp = d2.dr_wrt(d0).todense()
diff = tmp - np.eye(tmp.shape[0])
self.assertFalse(np.any(diff.ravel()))
def test_stacking(self):
a1 = ch.Ch(np.arange(10).reshape(2,5))
b1 = ch.Ch(np.arange(20).reshape(4,5))
c1 = ch.vstack((a1,b1))
c1_check = np.vstack((a1.r, b1.r))
residuals1 = (c1_check - c1.r).ravel()
a2 = ch.Ch(np.arange(10).reshape(5,2))
b2 = ch.Ch(np.arange(20).reshape(5,4))
c2 = ch.hstack((a2,b2))
c2_check = np.hstack((a2.r, b2.r))
residuals2 = (c2_check - c2.r).ravel()
self.assertFalse(np.any(residuals1))
self.assertFalse(np.any(residuals2))
