GAS/GAS_SPHForces.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 Lait
 *      Side Effects Software Inc
 *      123 Front Street West, Suite 1401
 *      Toronto, Ontario
 *      Canada   M5J 2M2
 *      416-504-9876
 *
 * NAME:        GAS_SPHForces.h ( GAS Library, C++)
 *
 * COMMENTS:
 */

#ifndef __GAS_SPHForces__
#define __GAS_SPHForces__

#include "GAS_API.h"

#include <GU/GU_Detail.h>

#include <UT/UT_ThreadedAlgorithm.h>

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

class GU_Detail;
class GAS_SPH;

class GAS_API GAS_SPHForces : public GAS_SubSolver
{
public:
    GET_DATA_FUNC_S(GAS_NAME_GEOMETRY, GeometryName);
    GET_DATA_FUNC_F("gasconstant", GasConstant);
    GET_DATA_FUNC_B("dopressure", DoPressure);
    GET_DATA_FUNC_B("doviscosity", DoViscosity);
    GET_DATA_FUNC_B("dotension", DoTension);
    GET_DATA_FUNC_B("computev", ComputeVelocity);

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

    /// Moves the fields given by FieldName by the given velocity
    /// field.  Moves the points given by GeometryName likewise.
    virtual bool         solveGasSubclass(SIM_Engine &engine,
                                SIM_Object *obj,
                                SIM_Time time,
                                SIM_Time timestep);

    bool                 shouldMultiThread() { return true; }

protected:
    void                 calculatePressure(GU_Detail *gdp, fpreal gasconstant,
                                int i = 0, int n = -1);
    void                 accumulatePressureForce(GU_Detail *gdp, GAS_SPH *sph,
                                int i = 0, int n = -1);
    void                 accumulateViscousForce(GU_Detail *gdp, GAS_SPH *sph,
                                int i = 0, int n = -1);
    void                 accumulateSurfaceForce(GU_Detail *gdp, GAS_SPH *sph,
                                int i = 0, int n = -1);

    THREADED_METHOD2(GAS_SPHForces, shouldMultiThread(), calcPressure,
                                GU_Detail *, gdp, fpreal, gasconstant)
    void                 calcPressurePartial(GU_Detail *gdp, fpreal gasconstant,
                                const UT_JobInfo &info);

    THREADED_METHOD4(GAS_SPHForces, shouldMultiThread(), calcForces,
                                GU_Detail *, gdp, bool, dopressure,
                                bool, doviscosity, bool, dotension)
    void                 calcForcesPartial(GU_Detail *gdp,
                                bool dopressure, bool diviscosity,
                                bool dotension, const UT_JobInfo &info);

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

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

private:
    static const SIM_DopDescription     *getDopDescription();

    DECLARE_STANDARD_GETCASTTOTYPE();
    DECLARE_DATAFACTORY(GAS_SPHForces,
                        GAS_SubSolver,
                        "Gas SPH Forces",
                        getDopDescription());
};

#endif

---