2 ///////////////////////////////////////////////////////////////////////////
3 // This program is free software: you can redistribute it and/or modify //
4 // it under the terms of the version 3 of the GNU General Public License //
5 // as published by the Free Software Foundation. //
7 // This program is distributed in the hope that it will be useful, but //
8 // WITHOUT ANY WARRANTY; without even the implied warranty of //
9 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU //
10 // General Public License for more details. //
12 // You should have received a copy of the GNU General Public License //
13 // along with this program. If not, see <http://www.gnu.org/licenses/>. //
15 // Written and (C) by Francois Fleuret //
16 // Contact <francois.fleuret@idiap.ch> for comments & bug reports //
17 ///////////////////////////////////////////////////////////////////////////
29 int main(int argc, char **argv) {
31 // X will be the "identity mapping"
32 Mapping X = Mapping::X;
34 // A undefined mapping
37 // Now contains some complex stuff
38 f = sin(X * X - (X * 3.0) + 0.5);
41 cout << f * f << endl;
43 // fp will be the first derivative
44 Mapping fp = f.derivative();
46 // fs the second derivative
47 Mapping fs = fp.derivative();
49 // Compute the second derivative at 5
50 cout << fs(5.0) << endl;