/
hermite.cpp
95 lines (86 loc) · 1.95 KB
/
hermite.cpp
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
85
86
87
88
89
90
91
92
93
94
95
#include "hermite.h"
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
using namespace std;
/*!
Constructs a new Hermite polynomial of a given degree.
@param degree Degree of the polynomial.
*/
Hermite::Hermite()
{
}
/*!
Evaluate the polynomial at the given value based on the polynomial's degree.
@param x Value at which to evaluate the polynomial.
@returns Value of the polynomial at the given x
*/
double Hermite::evaluate(int degree, double x) {
switch(degree) {
case 0:
return 1;
break;
case 1:
return 2*x;
break;
case 2:
return 4*x*x - 2;
break;
case 3:
return 8*x*x*x - 12;
break;
case 4:
return 16*x*x*x*x - 48*x*x + 12;
break;
default:
cerr << __PRETTY_FUNCTION__ << ": Hermite polynomial of unknown degree called. Degree was: " << degree << endl;
throw 999;
return 0;
}
}
double Hermite::derivative(int degree, double x) {
switch(degree) {
case 0:
return 0;
break;
case 1:
return 2;
break;
case 2:
return 4 * 2 * x;
break;
case 3:
return 8 * 3 * x * x;
break;
case 4:
return 16 * 4 * x * x * x - 48 * 2 * x;
break;
default:
cerr << __PRETTY_FUNCTION__ << ": Hermite polynomial of unknown degree called. Degree was: " << degree << endl;
throw 999;
return 0;
}
}
double Hermite::doubleDerivative(int degree, double x) {
switch(degree) {
case 0:
return 0;
break;
case 1:
return 0;
break;
case 2:
return 4 * 2;
break;
case 3:
return 8 * 3 * 2 * x;
break;
case 4:
return 16 * 4 * 3 * x * x - 48 * 2;
break;
default:
cerr << __PRETTY_FUNCTION__ << ": Hermite polynomial of unknown degree called. Degree was: " << degree << endl;
throw 999;
return 0;
}
}