solver/problem_param.cpp

Go to the documentation of this file.

/****************************************
 * INCLUDES 
 ****************************************/

#include "problem_param.h"

#include <cmath> //cos,sin

/****************************************
 * FUNCTIONS 
 ****************************************/

problem_parameters::problem_parameters (
    const int N_,
    const double Alpha,
    const double Beta,
    const double Gamma,
    const double regularization)
{
    N = N_;

    i2Q[0] = 1/(2*(Beta/2));
    i2Q[1] = 1/(2*(Alpha/2));
    i2Q[2] = 1/(2*regularization);

    i2P = 1/(2 * (Gamma/2));

    spar = new state_parameters[N];
}



problem_parameters::~problem_parameters()
{
    if (spar != NULL)
        delete [] spar;
}



void problem_parameters::set_state_parameters (
    const double* T_,
    const double* h_,
    const double h_initial,
    const double* angle)
{
    for (int i = 0; i < N; i++)
    {
        spar[i].cos = cos(angle[i]);
        spar[i].sin = sin(angle[i]);

        if (i == 0)
        {
            spar[i].A6 = T_[i]*T_[i]/2 - (h_[0] - h_initial);
        }
        else
        {
            spar[i].A6 = T_[i]*T_[i]/2 - (h_[i] - h_[i-1]);
        }

        spar[i].T = T_[i];
        spar[i].h = h_[i];

        spar[i].B[2] = T_[i];
        spar[i].B[1] = T_[i]*T_[i]/2;
        spar[i].B[0] = spar[i].B[1]*T_[i]/3 - h_[i]*T_[i];

        spar[i].A3 = T_[i];
    }
}