GradHSolver.h

Go to the documentation of this file.

#pragma once
 
 
#include "sph/solvers/AsymmetricSolver.h"
 
NAMESPACE_SPH_BEGIN
 
class IAsymmetricDerivative : public IDerivative {
public:
    virtual void evalAsymmetric(const Size idx,
        ArrayView<const Size> neighs,
        ArrayView<const Vector> gradi,
        ArrayView<const Vector> gradj) = 0;
};
 
class IAsymmetricTerm : public Polymorphic {
public:
    virtual void setAsymmetricDerivatives(Array<AutoPtr<IAsymmetricDerivative>>& asymmetricDerivatives) = 0;
};
 
// class EquationHolder : public EquationHolder3
 
 
class GradHSolver : public AsymmetricSolver {
private:
    Array<AutoPtr<IAsymmetricDerivative>> asymmetricDerivatives;
 
    struct SecondThreadData {
        Array<Vector> grads;
    };
 
    ThreadLocal<SecondThreadData> secondData;
 
public:
    GradHSolver(IScheduler& scheduler,
        const RunSettings& settings,
        const EquationHolder& basicTerms,
        Array<AutoPtr<IAsymmetricTerm>>&& asymmetricTerms);
 
    virtual void create(Storage& storage, IMaterial& material) const override;
 
protected:
    virtual void loop(Storage& storage, Statistics& stats) override;
};
 
 
NAMESPACE_SPH_END