SummationSolver.h

Go to the documentation of this file.

#pragma once
 
 
#include "sph/solvers/SymmetricSolver.h"
 
NAMESPACE_SPH_BEGIN
 
template <Size Dim>
class SummationSolver : public SymmetricSolver<Dim> {
private:
    Size maxIteration;
    Float targetDensityDifference;
    bool adaptiveH;
    Array<Float> rho, h;
 
    LutKernel<Dim> densityKernel;
 
    using ThreadData = typename SymmetricSolver<Dim>::ThreadData;
 
public:
    SummationSolver(IScheduler& scheduler,
        const RunSettings& settings,
        const EquationHolder& additionalEquations = {});
 
    SummationSolver(IScheduler& scheduler,
        const RunSettings& settings,
        const EquationHolder& additionalEquations,
        AutoPtr<IBoundaryCondition>&& bc);
 
    virtual void create(Storage& storage, IMaterial& material) const override;
 
private:
    virtual void beforeLoop(Storage& storage, Statistics& stats) override;
 
    virtual void sanityCheck(const Storage& storage) const override;
};
 
NAMESPACE_SPH_END