EquilibriumSolver.h

Go to the documentation of this file.

#pragma once
 
 
#include "math/SparseMatrix.h"
#include "objects/wrappers/Outcome.h"
#include "sph/equations/EquationTerm.h"
#include "sph/solvers/SymmetricSolver.h"
 
NAMESPACE_SPH_BEGIN
 
class IGravity;
 
class EquilibriumEnergySolver {
    IScheduler& scheduler;
 
    AutoPtr<IGravity> gravity;
 
    AutoPtr<IBasicFinder> finder;
 
    SymmetrizeSmoothingLengths<LutKernel<3>> kernel;
 
    Size boundaryThreshold;
 
public:
    EquilibriumEnergySolver(IScheduler& scheduler,
        const RunSettings& settings,
        AutoPtr<IGravity>&& gravity,
        const Size boundaryThreshold);
 
    ~EquilibriumEnergySolver();
 
    Outcome solve(Storage& storage, Statistics& stats);
};
 
class EquilibriumStressSolver {
private:
    IScheduler& scheduler;
 
    AutoPtr<ISymmetricFinder> finder;
 
    SymmetrizeSmoothingLengths<LutKernel<3>> kernel;
 
    SymmetricSolver<DIMENSIONS> equationSolver;
 
    Size boundaryThreshold;
 
    SparseMatrix matrix;
 
public:
    EquilibriumStressSolver(IScheduler& scheduler,
        const RunSettings& settings,
        const EquationHolder& equations);
 
    ~EquilibriumStressSolver();
 
    Outcome solve(Storage& storage, Statistics& stats);
 
    void create(Storage& storage, IMaterial& material);
};
 
NAMESPACE_SPH_END