1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
#! /usr/bin/env python
from __future__ import print_function
from openturns import *
ref = NumericalMathFunction("x", "sin(x)")
size = 12
locations = NumericalPoint(size)
values = NumericalPoint(size)
derivatives = NumericalPoint(size)
# Build locations/values/derivatives with non-increasing locations
for i in range(size):
locations[i] = 10.0 * i * i / (size - 1.0) / (size - 1.0)
values[i] = ref([locations[i]])[0]
derivatives[i] = ref.gradient([locations[i]])[0, 0]
evaluation = PiecewiseHermiteEvaluationImplementation(
locations, values, derivatives)
print("evaluation=", evaluation)
# Check the values
for i in range(2 * size):
x = [-1.0 + 12.0 * i / (2.0 * size - 1.0)]
print("f( %.12g )=" % x[0], evaluation(x), ", ref=", ref(x))
|
