def test_A_bar_k(self):
for eps in [5., 10.]:
pts = np.linspace(0, 1, 21)[1:-1]
natural_kernel = laplacian_natural_kernel.LaplacianNaturalKernel(eps, True)
for x in pts:
for y in pts:
res = natural_kernel.A_bar_k(x, y)
res_quad, err = A_bar_k_quadrature(x, y, eps)
np.testing.assert_almost_equal(res, res_quad, err_msg='Failed at ({}, {}), epsilon={}'.format(x, y, eps))
def test_I_4(self):
# exclude points directly on the boundary
for eps in [5., 10.]:
pts = np.linspace(0, 1, 21)[1:-1]
natural_kernel = laplacian_natural_kernel.LaplacianNaturalKernel(eps, True)
for x in pts:
for y in pts[(pts >= x)]:
res = natural_kernel.I_4(x, y)
res_quad, err = I_4_quadrature(x, y, eps)
np.testing.assert_almost_equal(res, res_quad, err_msg='Failed at ({}, {}), epsilon={}'.format(x, y, eps))