-
Notifications
You must be signed in to change notification settings - Fork 35
/
speed_comparison.py
54 lines (36 loc) · 1.13 KB
/
speed_comparison.py
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
from interpolation.splines.eval_cubic_numba import vec_eval_cubic_spline_3
from interpolation.splines.filter_cubic import filter_coeffs
from interpolation.splines.multilinear_numba import multilinear_interpolation
from interpolation.splines.misc import mlinspace
import numpy
K = 50
d = 3
N = 10**6
# N = 100
a = numpy.array([0.0,0.0,0.0])
b = numpy.array([1.0,1.0,1.0])
orders = numpy.array([K,K,K],dtype=int)
V = numpy.random.random(orders)
C = filter_coeffs(a,b,orders,V)
X = numpy.random.random((N,3))
res = numpy.zeros(N)
res2 = res.copy()
vec_eval_cubic_spline_3(a,b,orders,C,X,res)
multilinear_interpolation(a,b,orders,V,X,res)
import time
t1 = time.time()
vec_eval_cubic_spline_3(a,b,orders,C,X,res)
t2 = time.time()
multilinear_interpolation(a,b,orders,V,X,res2)
t3 = time.time()
print("Cubic: {}".format(t2-t1))
print("Linear: {}".format(t3-t2))
# assert(abs(res-res2).max()<1e-10)
# scipy
from scipy.interpolate import RegularGridInterpolator
pp = [numpy.linspace(a[i],b[i],orders[i]) for i in range(3)]
rgi = RegularGridInterpolator(pp, V)
t1 = time.time()
rgi(X)
t2 = time.time()
print("Scipy (linear): {}".format(t2-t1))