當前位置: 首頁>>代碼示例>>Python>>正文


Python legendre.legroots方法代碼示例

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


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

示例1: _getNodes

# 需要導入模塊: from numpy.polynomial import legendre [as 別名]
# 或者: from numpy.polynomial.legendre import legroots [as 別名]
def _getNodes(self):
        """
        Computes Gauss-Lobatto integration nodes.

        Calculates the Gauss-Lobatto integration nodes via a root calculation of derivatives of the legendre
        polynomials. Note that the precision of float 64 is not guarantied.

        Copyright by Dieter Moser, 2014

        Returns:
            np.ndarray: array of Gauss-Lobatto nodes
        """

        M = self.num_nodes
        a = self.tleft
        b = self.tright

        roots = leg.legroots(leg.legder(np.array([0] * (M - 1) + [1], dtype=np.float64)))
        nodes = np.array(np.append([-1.0], np.append(roots, [1.0])), dtype=np.float64)

        nodes = (a * (1 - nodes) + b * (1 + nodes)) / 2

        return nodes 
開發者ID:Parallel-in-Time,項目名稱:pySDC,代碼行數:25,代碼來源:gauss_lobatto.py

示例2: test_legroots

# 需要導入模塊: from numpy.polynomial import legendre [as 別名]
# 或者: from numpy.polynomial.legendre import legroots [as 別名]
def test_legroots(self):
        assert_almost_equal(leg.legroots([1]), [])
        assert_almost_equal(leg.legroots([1, 2]), [-.5])
        for i in range(2, 5):
            tgt = np.linspace(-1, 1, i)
            res = leg.legroots(leg.legfromroots(tgt))
            assert_almost_equal(trim(res), trim(tgt)) 
開發者ID:Frank-qlu,項目名稱:recruit,代碼行數:9,代碼來源:test_legendre.py

示例3: test_legroots

# 需要導入模塊: from numpy.polynomial import legendre [as 別名]
# 或者: from numpy.polynomial.legendre import legroots [as 別名]
def test_legroots(self) :
        assert_almost_equal(leg.legroots([1]), [])
        assert_almost_equal(leg.legroots([1, 2]), [-.5])
        for i in range(2, 5) :
            tgt = np.linspace(-1, 1, i)
            res = leg.legroots(leg.legfromroots(tgt))
            assert_almost_equal(trim(res), trim(tgt)) 
開發者ID:ktraunmueller,項目名稱:Computable,代碼行數:9,代碼來源:test_legendre.py


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