SymmetricSolver.h

Go to the documentation of this file.

#pragma once
 
 
#include "sph/equations/Derivative.h"
#include "sph/equations/EquationTerm.h"
#include "sph/kernel/Kernel.h"
#include "thread/ThreadLocal.h"
#include "timestepping/ISolver.h"
 
NAMESPACE_SPH_BEGIN
 
class IBoundaryCondition;
class ISymmetricFinder;
 
template <Size Dim>
class SymmetricSolver : public ISolver {
protected:
    struct ThreadData {
        DerivativeHolder derivatives;
 
        Array<NeighbourRecord> neighs;
 
        Array<Size> idxs;
 
        Array<Vector> grads;
    };
 
    IScheduler& scheduler;
 
    ThreadLocal<ThreadData> threadData;
 
    EquationHolder equations;
 
    AutoPtr<IBoundaryCondition> bc;
 
    AutoPtr<ISymmetricFinder> finder;
 
    LutKernel<Dim> kernel;
 
public:
    SymmetricSolver(IScheduler& scheduler,
        const RunSettings& settings,
        const EquationHolder& eqs,
        AutoPtr<IBoundaryCondition>&& bc);
 
    SymmetricSolver(IScheduler& scheduler, const RunSettings& settings, const EquationHolder& eqs);
 
    ~SymmetricSolver();
 
    virtual void integrate(Storage& storage, Statistics& stats) override;
 
    virtual void create(Storage& storage, IMaterial& material) const override;
 
protected:
    virtual void loop(Storage& storage, Statistics& stats);
 
    virtual void beforeLoop(Storage& storage, Statistics& stats);
 
    virtual void afterLoop(Storage& storage, Statistics& stats);
 
    virtual RawPtr<const IBasicFinder> getFinder(ArrayView<const Vector> r);
 
    virtual void sanityCheck(const Storage& storage) const;
};
 
NAMESPACE_SPH_END