implicit.hpp

Go to the documentation of this file.

#pragma once
 
#include "../particles.hpp"
#include "../refvalues.hpp"
#include "dirichlet_boundary_condition_generator/interface.hpp"
#include "interface.hpp"
#include "pressure_poisson_equation.hpp"
 
#include <Eigen/Sparse>
#include <memory>
#include <vector>
 
namespace PressureCalculator {
 
class Implicit : public Interface {
public:
    std::vector<double> calc(Particles& particles) override;
    ~Implicit() override;
 
    Implicit(
        int dimension,
        double particleDistance,
        double re_forGradient,
        double re_forLaplacian,
        double dt,
        double compressibility,
        double relaxationCoefficient,
        std::unique_ptr<DirichletBoundaryConditionGenerator::Interface>&& dirichletBoundaryConditionGenerator
    );
 
private:
    Particles particles;
    std::vector<double> pressure; 
    std::unique_ptr<DirichletBoundaryConditionGenerator::Interface> dirichletBoundaryConditionGenerator;
    PressurePoissonEquation pressurePoissonEquation;
 
    void removeNegativePressure();
};
 
} // namespace PressureCalculator