mirror of
https://codeberg.org/vcbferreira/NuFI_deal.ii
synced 2026-08-12 14:33:18 +02:00
59 lines
1.7 KiB
C++
59 lines
1.7 KiB
C++
#ifndef NUFI_SOLVER_H
|
|
#define NUFI_SOLVER_H
|
|
|
|
#include <boost/qvm/mat_access.hpp>
|
|
#include <cmath>
|
|
#include <deal.II/base/point.h>
|
|
#include <deal.II/base/tensor.h>
|
|
#include <deal.II/numerics/vector_tools.h>
|
|
#include <vector>
|
|
|
|
#include "nufi/fields.h" //dont remove
|
|
#include "nufi/grids.h"
|
|
#include "nufi/parameters.h"
|
|
#include "nufi/poisson_problem.h"
|
|
|
|
using namespace dealii;
|
|
|
|
class NuFISolver {
|
|
public:
|
|
NuFISolver();
|
|
|
|
void run();
|
|
std::vector<double>
|
|
eval_rho(unsigned int n, std::vector<double> &x,
|
|
const std::vector<GridStructure<1>> &grid_struct,
|
|
const std::vector<SolutionSnapshot<1>> &phi_history,
|
|
const unsigned int Nv = Parameters::NV) const;
|
|
std::vector<double>
|
|
eval_ftilda(unsigned int, std::vector<double> &x, double u,
|
|
const std::vector<GridStructure<1>> &grid_struct,
|
|
const std::vector<SolutionSnapshot<1>> &phi_history) const;
|
|
std::vector<double>
|
|
eval_f(unsigned int n, std::vector<double> &x, double u,
|
|
const std::vector<GridStructure<1>> &grid_struct,
|
|
const std::vector<SolutionSnapshot<1>> &phi_history) const;
|
|
|
|
std::vector<double>
|
|
eval_rho_points(unsigned int n, const std::vector<Point<1>> &points,
|
|
const std::vector<GridStructure<1>> &grid_struct,
|
|
const std::vector<SolutionSnapshot<1>> &phi_history,
|
|
const unsigned int Nv) const;
|
|
|
|
private:
|
|
unsigned int Nt = std::floor(Parameters::TMAX / Parameters::DT);
|
|
unsigned int Nx = Parameters::CALC_NX;
|
|
|
|
double Lx = Parameters::LX;
|
|
|
|
double x_min = Parameters::X_DOMAIN_LEFT;
|
|
double x_max = Parameters::X_DOMAIN_RIGHT;
|
|
|
|
std::vector<double> rho;
|
|
|
|
unsigned int order;
|
|
|
|
PoissonProblem<1> poisson;
|
|
};
|
|
#endif
|