Page MenuHomec4science

main.cpp
No OneTemporary

File Metadata

Created
Sun, Oct 20, 07:47

main.cpp

#include <iostream>
#include "VectorialODE.h"
#include "Input.h"
#include "InputConstruction.h"
#include "Solution.h"
#include "Bashforth.h"
#include "BashforthFirstStep.h"
#include "BashforthSecondStep.h"
#include "BashforthThirdStep.h"
#include "BashforthFourthStep.h"
int main() {
const double time_step = 1;
const int dimension= 2;
const int order= 4;
const int number_steps= 10;
const int write_output_timestep= 1;
vector<vector<double>> initial_condition_matrix = {{0, 2.9, 15.1, 61.2},{0, 1.7, 5.3, 14.4}};
vector<vector<double>> changing_matrix= {{1.1, 0},{0.1, 0}};
vector<vector<double>> past_step_matrix= {{1,0.44367,-0.60632,-0.98167,-0.26475,0.74675,0.92737,0.07613,-0.85982,-0.83907},
{1,0.44367,-0.60632,-0.98167,-0.26475,0.74675,0.92737,0.07613,-0.85982,-0.83907}};
Input input(time_step, dimension, order, number_steps, write_output_timestep, initial_condition_matrix, changing_matrix, past_step_matrix); //Declaration of an input object.
int columns =number_steps/write_output_timestep;
vector<vector<double>> solution_matrix (dimension, vector<double>(columns,0));
Solution solution(dimension, columns); //Declaration of an output object
BashforthFourthStep vectorialODE(input, solution);
vectorialODE.SolveVectorialODE();
}

Event Timeline