GAS/GAS_ParticlePressure.h

Go to the documentation of this file.

/*
 * PROPRIETARY INFORMATION.  This software is proprietary to
 * Side Effects Software Inc., and is not to be reproduced,
 * transmitted, or disclosed in any way without written permission.
 *
 * Produced by:
 *      Jeff Chadwick
 *      Side Effects Software Inc
 *      123 Front Street West, Suite 1401
 *      Toronto, Ontario
 *      Canada   M5J 2M2
 *      416-504-9876
 *
 * NAME:        GAS_ParticlePressure.h ( GAS Library, C++)
 *
 * COMMENTS:    This node is intended to produce a pressure-like
 *              force acting on the particles of a particle fluid
 *              simulation.  This force is fairly similar to the
 *              SPH style pressure calculation in GAS_SPHForces,
 *              but uses an idea referred to as "double density
 *              relaxation" (see paper repository).
 *
 */

#ifndef __GAS_ParticlePressure__
#define __GAS_ParticlePressure__

#include "GAS_API.h"

#include "GAS_SubSolver.h"
#include "GAS_Utils.h"

#include <UT/UT_PtrArray.h>
#include <UT/UT_RefArray.h>
#include <UT/UT_ThreadedAlgorithm.h>

class GEO_Point;
class GU_Detail;
class GAS_SPH;
class UT_FloatArray;

class GAS_API GAS_ParticlePressure : public GAS_SubSolver
{
public:
    GET_DATA_FUNC_S(GAS_NAME_GEOMETRY, GeometryName);
    GET_DATA_FUNC_F("gasconstant",  GasConstant);
    GET_DATA_FUNC_F("repulsionstrength", RepulsionStrength);
    GET_DATA_FUNC_B("dopressure",   DoPressure);
    GET_DATA_FUNC_B("computev",     ComputeVelocity);

    bool                 shouldMultiThread() { return true; }

protected:
    explicit             GAS_ParticlePressure(const SIM_DataFactory *factory);
    virtual             ~GAS_ParticlePressure();

    /// Computes a pressure force for the given particle field
    /// and either adjusts the velocity directly, or just
    /// adds to the force variable.
    virtual bool         solveGasSubclass(SIM_Engine &engine,
                                    SIM_Object *obj,
                                    SIM_Time time,
                                    SIM_Time timestep);

protected:
    void                 calculatePressureForce(GU_Detail *gdp,
                                    GAS_SPH *sph, fpreal gasconstant,
                                    fpreal repulsionstrength,
                                    SIM_Time timestep);

    /// A function that calculates the special density values required
    /// by this pressure computation.
    THREADED_METHOD5(GAS_ParticlePressure, shouldMultiThread(),
                                calcPressures,
                                GU_Detail *, gdp,
                                UT_FloatArray &, pressures,
                                UT_FloatArray &, nearpressures,
                                fpreal, gasconstant,
                                fpreal, repulsionstrength)
    void                 calcPressuresPartial(GU_Detail *gdp,
                                UT_FloatArray &pressures,
                                UT_FloatArray &nearpressures,
                                fpreal gasconstant,
                                fpreal repulsionstrength,
                                const UT_JobInfo &info);

    THREADED_METHOD6(GAS_ParticlePressure, shouldMultiThread(),
                                calcPressureForce,
                                GU_Detail *, gdp,
                                const UT_FloatArray &, pressures,
                                const UT_FloatArray &, nearpressures,
                                fpreal, gasconstant,
                                fpreal, repulsionstrength,
                                SIM_Time, timestep)
    void                 calcPressureForcePartial(GU_Detail *gdp,
                                const UT_FloatArray &pressures,
                                const UT_FloatArray &nearpressures,
                                fpreal gasconstant,
                                fpreal repulsionstrength,
                                SIM_Time timestep,
                                const UT_JobInfo &info);

    THREADED_METHOD2(GAS_ParticlePressure, shouldMultiThread(), applyForces,
                                GU_Detail *, gdp, SIM_Time, timestep)
    void                 applyForcesPartial(GU_Detail *gdp, SIM_Time timestep,
                                const UT_JobInfo &info);

    THREADED_METHOD1(GAS_ParticlePressure, shouldMultiThread(), clearForces,
                                GU_Detail *, gdp)
    void                 clearForcesPartial(GU_Detail *gdp,
                                const UT_JobInfo &info);

private:
    static const SIM_DopDescription     *getDopDescription();

    // Store a private list of point lists.
    static UT_RefArray<UT_PtrArray<GEO_Point *> > ourPointLists;

    DECLARE_STANDARD_GETCASTTOTYPE();
    DECLARE_DATAFACTORY(GAS_ParticlePressure,
                        GAS_SubSolver,
                        "Gas Particle Pressure",
                        getDopDescription());
};

#endif

---