tspline.cpp

/*************************************************************************************
 * MechSys - A C++ library to simulate (Continuum) Mechanical Systems                *
 * Copyright (C) 2005 Dorival de Moraes Pedroso <dorival.pedroso at gmail.com>       *
 * Copyright (C) 2005 Raul Dario Durand Farfan  <raul.durand at gmail.com>           *
 *                                                                                   *
 * This file is part of MechSys.                                                     *
 *                                                                                   *
 * MechSys is free software; you can redistribute it and/or modify it under the      *
 * terms of the GNU General Public License as published by the Free Software         *
 * Foundation; either version 2 of the License, or (at your option) any later        *
 * version.                                                                          *
 *                                                                                   *
 * MechSys is distributed in the hope that it will be useful, but WITHOUT ANY        *
 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A   *
 * PARTICULAR PURPOSE. See the GNU General Public License for more details.          *
 *                                                                                   *
 * You should have received a copy of the GNU General Public License along with      *
 * MechSys; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, *
 * Fifth Floor, Boston, MA 02110-1301, USA                                           *
 *************************************************************************************/

#include <iostream>
#include <iomanip>
#include <cmath>

#include "util/numstreams.h"
#include "numerical/spline.h"

using Util::_6_3;
using std::cout;
using std::endl;

int main(int argc, char **argv) try
{
    int  n = 21;
    REAL X[] = { 0.9, 1.3, 1.9 , 2.1, 2.6, 3.0, 3.9, 4.4 , 4.7 , 5.0, 6.0 , 7.0, 8.0 , 9.2 , 10.5, 11.3, 11.6, 12.0, 12.6, 13.0, 13.3  };
    REAL Y[] = { 1.3, 1.5, 1.85, 2.1, 2.6, 2.7, 2.4, 2.15, 2.05, 2.1, 2.25, 2.3, 2.25, 1.95,  1.4,  0.9,  0.7,  0.6,  0.5,  0.4,  0.25 };
    Spline S(X,Y,n);
    cout << S << endl;

    cout << _6_3()<< "x" << _6_3()<< "y" << endl;
    for (double x=0.9; x<13.3; x+=0.1)
        cout << _6_3()<< x << _6_3()<< S.Eval(x) << endl;

    return 0;
}
catch (Exception * e)
{
    e->Cout();
    if (e->IsFatal()) {delete e; exit(1);}
    delete e;
}
catch (...)
{
    std::cout << "Some exception (...) ocurred\n";
}

---