smooth_linear_splines.cpp

Go to the documentation of this file.

 
/*****************************************************************************
** Includes
*****************************************************************************/
 
#include <iostream>
#include <string>
#include <gtest/gtest.h>
#include <ecl/containers/array.hpp>
#include <ecl/formatters/floats.hpp>
#include <ecl/formatters/strings.hpp>
#include "../../include/ecl/geometry/smooth_linear_spline.hpp"
 
/*****************************************************************************
** Using
*****************************************************************************/
 
using std::cout; using std::endl;
using std::string;
using ecl::Array;
using ecl::RightAlign;
using ecl::Format;
using ecl::SmoothLinearSpline;
 
/*****************************************************************************
** Tests
*****************************************************************************/
 
TEST(SmoothLinearSplinesTests,allEggsInOneBasket) {
        // Haven't got around to running this properly through gtests yet.
        SUCCEED();
}
/*****************************************************************************
** Main program
*****************************************************************************/
 
int main(int argc, char **argv) {
 
    Format<string> string_format; string_format.width(8); string_format.align(RightAlign);
    Format<double> format; format.width(8); format.precision(2); format.align(RightAlign);
 
    /*********************
    ** Initialisation
    **********************/
    Array<double> x_set(6);
    Array<double> y_set(6);
    x_set << 0.0, 1.0, 2.0, 3.0, 4.0, 5.0;
    y_set << 1.0, 2.0, 1.0, 3.0, 4.0, 4.0;
    const double a_max = 10.0;
 
    cout << endl;
    cout << "***********************************************************" << endl;
    cout << "           Smooth Linear Spline Constructor" << endl;
    cout << "***********************************************************" << endl;
    cout << endl;
 
    SmoothLinearSpline spline = SmoothLinearSpline::Interpolation(x_set, y_set, a_max);
//    cout << spline_1 << endl;
 
    cout << endl;
    cout << "***********************************************************" << endl;
    cout << "                      Output" << endl;
    cout << "***********************************************************" << endl;
    cout << endl;
 
    int n = 200;
    cout << string_format("x  ");
    cout << string_format("y  ");
    cout << string_format("y' ");
    cout << string_format("y''") << endl;
    for ( int i = 0; i <= n; ++i ) {
        double x = x_set[0] + i*(x_set.back()-x_set.front())/n;
        cout << format(x);
        cout << format(spline(x));
        cout << format(spline.derivative(x));
        cout << format(spline.dderivative(x)) << endl;
    }
 
    cout << endl;
    cout << "***********************************************************" << endl;
    cout << "                      Passed" << endl;
    cout << "***********************************************************" << endl;
    cout << endl;
 
    testing::InitGoogleTest(&argc,argv);
    return RUN_ALL_TESTS();
}