예제 #1
0
 def test_hermfromroots(self):
     res = herm.hermfromroots([])
     assert_almost_equal(trim(res), [1])
     for i in range(1, 5):
         roots = np.cos(np.linspace(-np.pi, 0, 2 * i + 1)[1::2])
         pol = herm.hermfromroots(roots)
         res = herm.hermval(roots, pol)
         tgt = 0
         assert_(len(pol) == i + 1)
         assert_almost_equal(herm.herm2poly(pol)[-1], 1)
         assert_almost_equal(res, tgt)
예제 #2
0
 def test_hermfromroots(self) :
     res = herm.hermfromroots([])
     assert_almost_equal(trim(res), [1])
     for i in range(1,5) :
         roots = np.cos(np.linspace(-np.pi, 0, 2*i + 1)[1::2])
         pol = herm.hermfromroots(roots)
         res = herm.hermval(roots, pol)
         tgt = 0
         assert_(len(pol) == i + 1)
         assert_almost_equal(herm.herm2poly(pol)[-1], 1)
         assert_almost_equal(res, tgt)
예제 #3
0
 def test_hermroots(self):
     assert_almost_equal(herm.hermroots([1]), [])
     assert_almost_equal(herm.hermroots([1, 1]), [-.5])
     for i in range(2, 5):
         tgt = np.linspace(-1, 1, i)
         res = herm.hermroots(herm.hermfromroots(tgt))
         assert_almost_equal(trim(res), trim(tgt))
예제 #4
0
 def test_hermroots(self) :
     assert_almost_equal(herm.hermroots([1]), [])
     assert_almost_equal(herm.hermroots([1, 1]), [-.5])
     for i in range(2,5) :
         tgt = np.linspace(-1, 1, i)
         res = herm.hermroots(herm.hermfromroots(tgt))
         assert_almost_equal(trim(res), trim(tgt))