tension_function.cpp

Go to the documentation of this file.

/*****************************************************************************
** Includes
*****************************************************************************/

#include <iostream>
#include <string>
#include <gtest/gtest.h>
#include <ecl/formatters/floats.hpp>
#include <ecl/formatters/strings.hpp>
#include "../../include/ecl/geometry/tension_function.hpp"

/*****************************************************************************
** Using
*****************************************************************************/

using std::string;
using ecl::RightAlign;
using ecl::Format;
using ecl::TensionFunction;

/*****************************************************************************
** Tests
*****************************************************************************/

TEST(TensionFunction,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);

    std::cout << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << "                    Constructors" << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << std::endl;

    TensionFunction function;
    function = TensionFunction::Interpolation(2.0,1.0,2.0,3.0,2.0,3.0);

    std::cout << "f(tau,x) =" << std::endl;
    std::cout << function << std::endl;

    std::cout << string_format("t=0.01");
    std::cout << string_format("t=1");
    std::cout << string_format("t=2");
    std::cout << string_format("t=3");
    std::cout << string_format("t=4");
    std::cout << string_format("t=10") << std::endl;
    for ( int i = 0; i <= 50; ++i ) {
        std::cout << format(function(0.01,2.0 + 0.02*i));
        for ( int j = 1; j < 5; j+=1 ) {
            std::cout << format(function(1.0*j,2.0 + 0.02*i));
        }
        std::cout << format(function(10.0,2.0 + 0.02*i));
        std::cout << std::endl;
    }

    std::cout << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << "                    Derivatives" << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << std::endl;

    std::cout << string_format("t=0.01");
    std::cout << string_format("t=1");
    std::cout << string_format("t=2");
    std::cout << string_format("t=3");
    std::cout << string_format("t=4");
    std::cout << string_format("t=10") << std::endl;
    for ( int i = 0; i <= 50; ++i ) {
        std::cout << format(function.derivative(0.01,2.0 + 0.02*i));
        for ( int j = 1; j < 5; j+=1 ) {
            std::cout << format(function.derivative(1.0*j,2.0 + 0.02*i));
        }
        std::cout << format(function.derivative(10.0,2.0 + 0.02*i));
        std::cout << std::endl;
    }


    std::cout << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << "                  2nd Derivatives" << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << std::endl;

    std::cout << string_format("t=0.01");
    std::cout << string_format("t=1");
    std::cout << string_format("t=2");
    std::cout << string_format("t=3");
    std::cout << string_format("t=4");
    std::cout << string_format("t=10") << std::endl;
    for ( int i = 0; i <= 50; ++i ) {
        std::cout << format(function.dderivative(0.01,2.0 + 0.02*i));
        for ( int j = 1; j < 5; j+=1 ) {
            std::cout << format(function.dderivative(1.0*j,2.0 + 0.02*i));
        }
        std::cout << format(function.dderivative(10.0,2.0 + 0.02*i));
        std::cout << std::endl;
    }


    std::cout << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << "                      Passed" << std::endl;
    std::cout << "***********************************************************" << std::endl;
    std::cout << std::endl;

    testing::InitGoogleTest(&argc,argv);
    return RUN_ALL_TESTS();
}