def test_gauss_hermite_3d(self): if self.writeData: p,w = gq.gauss_hermite([3,3,4]) np.savez('data/test_gauss_hermite_3d.npz',p=p,w=w) data = helper.load_test_npz('test_gauss_hermite_3d.npz') p,w = gq.gauss_hermite([3,3,4]) np.testing.assert_almost_equal(p, data['p']) np.testing.assert_almost_equal(w, data['w'])
def test_gauss_hermite_1d_int_arg(self): if self.writeData: p,w = gq.gauss_hermite(7) np.savez('data/test_gh1d_7pts.npz',p=p,w=w) data = helper.load_test_npz('test_gh1d_7pts.npz') p,w = gq.gauss_hermite(7) np.testing.assert_almost_equal(p, data['p']) np.testing.assert_almost_equal(w, data['w'])
def av_quadrature_rule(avmap, N): """Get a quadrature rule on the space of active variables. Parameters ---------- avmap : ActiveVariableMap a domains.ActiveVariableMap N : int the number of quadrature nodes in the active variables Returns ------- Yp : ndarray quadrature nodes on the active variables Yw : ndarray quadrature weights on the active variables See Also -------- integrals.quadrature_rule """ m, n = avmap.domain.subspaces.W1.shape if isinstance(avmap.domain, UnboundedActiveVariableDomain): NN = [int(np.floor(np.power(N, 1.0/n))) for i in range(n)] Yp, Yw = gq.gauss_hermite(NN) elif isinstance(avmap.domain, BoundedActiveVariableDomain): if n == 1: Yp, Yw = interval_quadrature_rule(avmap, N) else: Yp, Yw = zonotope_quadrature_rule(avmap, N) else: raise Exception('There is a problem with the avmap.domain.') return Yp, Yw
def gauss_hermite_design(N): """Tensor product Gauss-Hermite quadrature points. Parameters ---------- N : int[] contains the number of points per dimension in the tensor product design Returns ------- design : ndarray N-by-m matrix that contains the design points """ design = gauss_hermite(N)[0] return design
def test_gauss_hermite_3d(self): p,w = gq.gauss_hermite([3,3,4])
def test_gauss_hermite_2d(self): p,w = gq.gauss_hermite([3,3])
def test_gauss_hermite_1d_int_arg(self): p,w = gq.gauss_hermite(7)
def test_gauss_hermite_1d_array_arg(self): p,w = gq.gauss_hermite([7])
def test_gauss_hermite_3d(self): data = helper.load_test_npz('test_gauss_hermite_3d.npz') p,w = gq.gauss_hermite([3,3,4]) np.testing.assert_equal(p, data['p']) np.testing.assert_equal(w, data['w'])
def test_gauss_hermite_1d_int_arg(self): data = helper.load_test_npz('test_gh1d_7pts.npz') p,w = gq.gauss_hermite(7) np.testing.assert_equal(p, data['p']) np.testing.assert_equal(w, data['w'])
def test_gauss_hermite_3d(self): p, w = gq.gauss_hermite([3, 3, 4])
def test_gauss_hermite_2d(self): p, w = gq.gauss_hermite([3, 3])
def test_gauss_hermite_1d_int_arg(self): p, w = gq.gauss_hermite(7)
def test_gauss_hermite_1d_array_arg(self): p, w = gq.gauss_hermite([7])