1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
"""Test the multilinear interpolation Cython implementation."""
import unittest
import numpy as np
from geotiepoints.multilinear import MultilinearInterpolator
ARR1 = np.array([[0.99819991, 0.46872418, 0.39356122, 0.33170364, 0.79171217,
0.31300346, 0.69075206, 0.47492703, 0.2012821, 0.57347047,
0.12934437, 0.7384143, 0.84650654, 0.63113452, 0.35514309,
0.36773267, 0.42887193, 0.83856559, 0.1081537, 0.33634562,
0.25445117, 0.30121727, 0.19697695, 0.92937056, 0.40359487,
0.10237384, 0.37803665, 0.94699248, 0.08045698, 0.01366914],
[0.72748552, 0.70872219, 0.36883461, 0.52914895, 0.22308535,
0.4444687, 0.52394334, 0.52870835, 0.31756298, 0.13776131,
0.11812231, 0.46974149, 0.12318789, 0.76525517, 0.97328814,
0.6580273, 0.93059119, 0.32785305, 0.57033161, 0.80133526,
0.4311177, 0.44957946, 0.81073879, 0.79356296, 0.77565555,
0.90520185, 0.76064422, 0.78609587, 0.43915797, 0.50745485]])
RES1 = np.array([[1.25149566, 0.86801702, 0.59233148, 0.66310018, 0.86093089,
0.59448266, 0.88145882, 0.72057176, 0.44395844, 0.62692691,
0.22956686, 0.89308545, 0.8823242, 1.01691582, 1.06050602,
0.78964244, 1.04244001, 0.93481853, 0.61242017, 0.90454307,
0.5570492, 0.59214105, 0.87052131, 1.24339837, 0.9031915,
0.93282679, 0.88212567, 1.25090173, 0.47821939, 0.51308049],
[1.13891573, 0.79536539, 0.54753759, 0.61188808, 0.82629032,
0.55155224, 0.8025709, 0.64724795, 0.42423377, 0.60563713,
0.22485215, 0.82022895, 0.86317281, 0.92860053, 1.00514142,
0.73270743, 0.9756435, 0.88385772, 0.59570712, 0.85230456,
0.52319791, 0.55035212, 0.83992216, 1.12698577, 0.84054343,
0.91689093, 0.82346408, 1.13391799, 0.46731605, 0.5109711]])
def assertNumpyArraysEqual(self, other):
if self.shape != other.shape:
raise AssertionError("Shapes don't match")
if not np.allclose(self, other):
raise AssertionError("Elements don't match!")
class TestMultilinearInterpolator(unittest.TestCase):
"""Class for unit testing the multilinear interpolation method
"""
def setUp(self):
pass
def test_multilinear_interp(self):
"""Test the multilinear interpolation"""
smin = [-1, -1]
smax = [1, 1]
orders = [5, 5]
f = lambda x: np.vstack([
np.sqrt(x[0, :]**2 + x[1, :]**2),
np.power(x[0, :]**3 + x[1, :]**3, 1.0 / 3.0)
])
interp = MultilinearInterpolator(smin, smax, orders)
interp.set_values(f(interp.grid))
result = interp(ARR1)
# exact_values = f(ARR1)
assertNumpyArraysEqual(result, RES1)
def tearDown(self):
"""Clean up"""
return
def suite():
"""The suite for Multilinear Interpolator"""
loader = unittest.TestLoader()
mysuite = unittest.TestSuite()
mysuite.addTest(loader.loadTestsFromTestCase(TestMultilinearInterpolator))
return mysuite
if __name__ == '__main__':
unittest.main()
|
