Esempio n. 1
0
 def test_clenshaw__empty(self):
     self.assertEquals(clenshaw(np.array([]), np.array([])).size, 0)
     self.assertEquals(clenshaw(np.array([]), np.array([1.])).size, 0)
     self.assertEquals(clenshaw(np.array([1.]), np.array([])).size, 0)
     self.assertEquals(clenshaw(self.pts, np.array([])).size, 0)
     self.assertEquals(clenshaw(np.array([]), self.pts).size, 0)
     self.assertNotEquals(clenshaw(np.array([.1]), np.array([.1])).size, 0)
Esempio n. 2
0
 def test_clenshaw__empty(self):
     self.assertEquals(clenshaw(np.array([]), np.array([])).size, 0)
     self.assertEquals(clenshaw(np.array([]), np.array([1.])).size, 0)
     self.assertEquals(clenshaw(np.array([1.]), np.array([])).size, 0)
     self.assertEquals(clenshaw(self.pts, np.array([])).size, 0)
     self.assertEquals(clenshaw(np.array([]), self.pts).size, 0)
     self.assertNotEquals(clenshaw(np.array([.1]), np.array([.1])).size, 0)
Esempio n. 3
0
def evalTester(method, fun, evalpts, chebpts):

    x = evalpts
    xk = chebpts
    fvals = fun(xk)

    if method is bary:
        vk = Chebtech2._barywts(fvals.size)
        a = bary(x, fvals, xk, vk)
        tol_multiplier = 1e0

    elif method is clenshaw:
        ak = Chebtech2._vals2coeffs(fvals)
        a = clenshaw(x, ak)
        tol_multiplier = 2e1

    b = fun(evalpts)
    n = evalpts.size
    tol = tol_multiplier * scaled_tol(n)

    return infNormLessThanTol(a, b, tol)
Esempio n. 4
0
def evalTester(method, fun, evalpts, chebpts):

    x = evalpts
    xk = chebpts
    fvals = fun(xk)

    if method is bary:
        vk = Chebtech2._barywts(fvals.size)
        a = bary(x, fvals, xk, vk)
        tol_multiplier = 1e0

    elif method is clenshaw:
        ak = Chebtech2._vals2coeffs(fvals)
        a = clenshaw(x, ak)
        tol_multiplier = 2e1

    b = fun(evalpts)
    n = evalpts.size
    tol = tol_multiplier * scaled_tol(n)

    return infNormLessThanTol(a, b, tol)
Esempio n. 5
0
 def __call__clenshaw(self, x):
     return clenshaw(x, self.coeffs)
Esempio n. 6
0
 def test_clenshaw__scalar_input(self):
     for x in self.xx:
         self.assertTrue(np.isscalar(clenshaw(x, self.ak)))
     self.assertFalse(np.isscalar(clenshaw(xx, self.ak)))
Esempio n. 7
0
 def test_clenshaw__scalar_input(self):
     for x in self.xx:
         self.assertTrue(np.isscalar(clenshaw(x,self.ak)))
     self.assertFalse(np.isscalar(clenshaw(xx,self.ak)))
Esempio n. 8
0
 def __call__clenshaw(self, x):
     return clenshaw(x, self.coeffs)