def test_eval_U_grid1(self):
log = logging.getLogger('TestCPsf.test_eval_U_grid1')
enz = CPsf(8)
r_sp = float(abs(normal(size=1, ))[0])
ph_sp = float(normal(size=1, )[0])
f_sp = float(normal(size=1, )[0])
enz.make_pol_grid([r_sp], [ph_sp], [f_sp])
beta = normal(size=enz.czern.nk) + 1j * normal(size=enz.czern.nk)
a = enz.eval_grid(beta, 0, 0, 0)
b = enz.U(beta, r_sp, ph_sp, f_sp)
log.debug('r = {:e}, ph = {:e}, f = {:e}'.format(r_sp, ph_sp, f_sp))
log.debug('a {:s} {}'.format(str(type(a)), a))
log.debug('b {:s} {}'.format(str(type(b)), b))
log.debug('{}'.format(abs(a - b)))
self.assertTrue(abs(a - b) < self.max_enorm)
def test_eval_U_grid1(self):
log = logging.getLogger("TestCPsf.test_eval_U_grid1")
enz = CPsf(8)
r_sp = float(abs(normal(size=1))[0])
ph_sp = float(normal(size=1)[0])
f_sp = float(normal(size=1)[0])
enz.make_pol_grid([r_sp], [ph_sp], [f_sp])
beta = normal(size=enz.czern.nk) + 1j * normal(size=enz.czern.nk)
a = enz.eval_grid(beta, 0, 0, 0)
b = enz.U(beta, r_sp, ph_sp, f_sp)
log.debug("r = {:e}, ph = {:e}, f = {:e}".format(r_sp, ph_sp, f_sp))
log.debug("a {:s} {}".format(str(type(a)), a))
log.debug("b {:s} {}".format(str(type(b)), b))
log.debug("{}".format(abs(a - b)))
self.assertTrue(abs(a - b) < self.max_enorm)
