Python GaussQuadratures.quadrature方法代碼示例

本文整理匯總了Python中pymor.tools.quadratures.GaussQuadratures.quadrature方法的典型用法代碼示例。如果您正苦於以下問題：Python GaussQuadratures.quadrature方法的具體用法？Python GaussQuadratures.quadrature怎麽用？Python GaussQuadratures.quadrature使用的例子？那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在類pymor.tools.quadratures.GaussQuadratures的用法示例。

在下文中一共展示了GaussQuadratures.quadrature方法的7個代碼示例，這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚，您的評價將有助於係統推薦出更棒的Python代碼示例。

示例1: init

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
    def __init__(self):
        def tensor_points(P):
            PP0, PP1 = np.array(np.meshgrid(P, P))
            return np.array((PP0.ravel(), PP1.ravel())).T

        def tensor_weights(W):
            return np.dot(W[:, np.newaxis], W[np.newaxis, :]).ravel()
        self._quadrature_points = [tensor_points(GaussQuadratures.quadrature(npoints=p + 1)[0])
                                   for p in xrange(GaussQuadratures.maxpoints())]
        self._quadrature_weights = [tensor_weights(GaussQuadratures.quadrature(npoints=p + 1)[1])
                                    for p in xrange(GaussQuadratures.maxpoints())]
        self._quadrature_npoints = np.arange(1, GaussQuadratures.maxpoints() + 1) ** 2
        self._quadrature_orders = GaussQuadratures.orders
        self._quadrature_order_map = GaussQuadratures.order_map

開發者ID:maltedschumacher，項目名稱:pyMor，代碼行數:16，代碼來源:referenceelements.py

示例2: init

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
 def __init__(self, flux, flux_derivative, gausspoints=5, intervals=1):
     self.flux = flux
     self.flux_derivative = flux_derivative
     self.gausspoints = gausspoints
     self.intervals = intervals
     self.build_parameter_type(inherits=(flux, flux_derivative))
     points, weights = GaussQuadratures.quadrature(npoints=self.gausspoints)
     points = points / intervals
     points = ((np.arange(self.intervals, dtype=np.float)[:, np.newaxis] * (1 / intervals))
               + points[np.newaxis, :]).ravel()
     weights = np.tile(weights, intervals) * (1 / intervals)
     self.points = points
     self.weights = weights

開發者ID:JuliaBru，項目名稱:pymor，代碼行數:15，代碼來源:fv.py

示例3: quadrature

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
 def quadrature(self, order=None, npoints=None, quadrature_type='default'):
     if quadrature_type == 'default' or quadrature_type == 'gauss':
         P, W = GaussQuadratures.quadrature(order, npoints)
         return P[:, np.newaxis], W
     else:
         raise NotImplementedError('quadrature_type must be "default" or "gauss"')

開發者ID:tobiasleibner，項目名稱:pymor，代碼行數:8，代碼來源:referenceelements.py

示例4: test_points

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
 def test_points(self):
     for order in GaussQuadratures.orders:
         P, _ = GaussQuadratures.quadrature(order)
         assert float_cmp_all(P, np.sort(P))
         assert 0.0 < P[0]
         assert P[-1] < 1.0

開發者ID:lucas-ca，項目名稱:pymor，代碼行數:8，代碼來源:tools.py

示例5: test_weights

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
 def test_weights(self):
     for order in GaussQuadratures.orders:
         _, W = GaussQuadratures.quadrature(order)
         assert float_cmp(sum(W), 1)

開發者ID:lucas-ca，項目名稱:pymor，代碼行數:6，代碼來源:tools.py

示例6: test_points

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
 def test_points(self):
     for order in GaussQuadratures.orders:
         P, _ = GaussQuadratures.quadrature(order)
         np.testing.assert_array_equal(P, np.sort(P))
         self.assertLess(0.0, P[0])
         self.assertLess(P[-1], 1.0)

開發者ID:BarbaraV，項目名稱:pymor，代碼行數:8，代碼來源:tools.py

示例7: test_weights

# 需要導入模塊: from pymor.tools.quadratures import GaussQuadratures [as 別名]
# 或者: from pymor.tools.quadratures.GaussQuadratures import quadrature [as 別名]
 def test_weights(self):
     for order in GaussQuadratures.orders:
         _, W = GaussQuadratures.quadrature(order)
         self.assertAlmostEqual(sum(W), 1)

開發者ID:BarbaraV，項目名稱:pymor，代碼行數:6，代碼來源:tools.py

注：本文中的pymor.tools.quadratures.GaussQuadratures.quadrature方法示例由純淨天空整理自Github/MSDocs等開源代碼及文檔管理平台，相關代碼片段篩選自各路編程大神貢獻的開源項目，源碼版權歸原作者所有，傳播和使用請參考對應項目的License；未經允許，請勿轉載。

示例1: __init__