Как повысить точность регулярного гринтерполятора - PullRequest
Новая
       21

Как повысить точность регулярного гринтерполятора

0 голосов
/ 16 мая 2019

Я создал файл h5 аналитической функции на обычном скелете. Чтобы получить производные этой функции, я использовал numpy.gradient здесь. Но проблема в том, что у меня очень плохой результат для регулярного сетчатого интерполятора, особенно для интерполяции градиентов. Я сравнил свой результат, полученный с помощью Regulargridinterpolator, с аналитическим результатом, и ошибка достаточно велика. Как я могу получить лучшую точность от Regulargridinterpolator или какой интерполятор обеспечивает лучшую точность для данных 3D регулярной сетки?

Я создал файл h5 следующим образом: 

import numpy as np
from numpy import gradient
import h5py
from scipy.interpolate import RegularGridInterpolator


def f(x,y,z):
    return ( -1 / np.sqrt(x**2 + y**2 + z**2))

#grid
x = np.linspace(0, 100, 32)  # since the boxsize is 320 Mpc/h
y = np.linspace(0, 100, 32)
z = np.linspace(0, 100, 32)


mesh_data = phi_an(*np.meshgrid(x, y, z, indexing='ij', sparse=True))


#create h5 file
h5file = h5py.File('analytic.h5', 'w')
h5file.create_dataset('/x', data=x)
h5file.create_dataset('/y', data=y)
h5file.create_dataset('/z', data=z)
h5file.create_dataset('/mesh_data', data=mesh_data)

h5file.close()

и код интерполяции приведен ниже: 

interpolator_f = RegularGridInterpolator((x, y, z), data[...])
interpolator_df_dx = RegularGridInterpolator((x, y, z), np.gradient(data[...], 100.0/32.0, axis=0))
...