#!/usr/bin/env python3

import numpy as np
import ndpolator

a1 = np.array([0.0])

ndp = ndpolator.Ndpolator(basic_axes=(a1,))

grid = np.empty((len(a1), 1))

for i, x in enumerate(a1):
    grid[i, 0] = 1.0

ndp.register(table='main', associated_axes=None, grid=grid)

query_pts = [[0.0]]

interps = ndp.ndpolate(table='main', query_pts=query_pts, extrapolation_method='linear')

print("a1 = ", a1)
print("query_pts = ", query_pts)
print("interps['interps'] = ", interps['interps'])