main.cc

   1
   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.                         //
   6 //                                                                       //
   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.                              //
  11 //                                                                       //
  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/>.  //
  14 //                                                                       //
  15 // Written and (C) by Francois Fleuret                                   //
  16 // Contact <francois.fleuret@idiap.ch> for comments & bug reports        //
  17 ///////////////////////////////////////////////////////////////////////////
  18
  19 #include <iostream>
  20 #include <fstream>
  21 #include <cmath>
  22 #include <stdio.h>
  23 #include <stdlib.h>
  24
  25 #include "mappings.h"
  26
  27 using namespace std;
  28
  29 int main(int argc, char **argv) {
  30
  31   // X will be the "identity mapping"
  32   Mapping X = Mapping::X;
  33
  34   // A undefined mapping
  35   Mapping f;
  36
  37   // Now contains some complex stuff
  38   f = sin(X * X - (X * 3.0) + 0.5);
  39
  40   // Display it
  41   cout << f * f << endl;
  42
  43   // fp will be the first derivative
  44   Mapping fp = f.derivative();
  45
  46   // fs the second derivative
  47   Mapping fs = fp.derivative();
  48
  49   // Compute the second derivative at 5
  50   cout << fs(5.0) << endl;
  51
  52 }