/usr/include/libwildmagic/Wm5Fluid2Db.h is in libwildmagic-dev 5.13-1+b2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 | // Geometric Tools, LLC
// Copyright (c) 1998-2014
// Distributed under the Boost Software License, Version 1.0.
// http://www.boost.org/LICENSE_1_0.txt
// http://www.geometrictools.com/License/Boost/LICENSE_1_0.txt
//
// File Version: 5.0.1 (2010/10/01)
#ifndef WM5FLUID2DB_H
#define WM5FLUID2DB_H
#include "Wm5PhysicsLIB.h"
#include "Wm5Vector2.h"
// This class is an implementation based on the ideas in "Real-Time Fluid
// Dynamics for Games", by Jos Stam (GDC 2003 Proceedings), but the code
// is faster in that it processes the source, diffusion, and advection
// terms all at the same time. The only Gauss-Seidel iteration is in the
// Poisson solver.
namespace Wm5
{
template <typename Real>
class WM5_PHYSICS_ITEM Fluid2Db
{
public:
// Construction and destruction.
Fluid2Db (Real x0, Real y0, Real x1, Real y1, Real dt, Real denViscosity,
Real velViscosity, int imax, int jmax, int numGaussSeidelIterations,
bool densityDirichlet);
virtual ~Fluid2Db ();
// Member access.
inline Real GetX0 () const;
inline Real GetY0 () const;
inline Real GetX1 () const;
inline Real GetY1 () const;
inline Real GetDt () const;
inline Real GetDx () const;
inline Real GetDy () const;
inline Real GetTime () const;
inline int GetIMax () const;
inline int GetJMax () const;
inline const Real* GetX () const;
inline const Real* GetY () const;
inline Real** GetDensity () const;
inline Vector2<Real>** GetVelocity () const;
// Derived classes must provide initialization at time 0.
virtual Real InitialDensity (Real x, Real y,
int i, int j) = 0;
virtual Vector2<Real> InitialVelocity (Real x, Real y,
int i, int j) = 0;
// Derived classes must provide source terms. The density has a
// source-sink function. The velocity has a force (impulse) function.
virtual Real SourceDensity (Real t, Real x, Real y,
int i, int j) = 0;
virtual Vector2<Real> SourceVelocity (Real t, Real x, Real y,
int i, int j) = 0;
void Initialize ();
void DoSimulationStep ();
protected:
// Update the boundary values after the interior values have been
// assigned.
void DirichletBoundaryZero (Real** data);
void NeumannBoundaryZero (Real** data);
void UpdateDensityBoundary ();
void UpdateVelocityBoundary ();
// Swap the buffer pointers.
void SwapDensityBuffers ();
void SwapVelocityBuffers ();
// Constructor inputs.
Real mX0, mY0, mX1, mY1;
Real mDt;
Real mDenViscosity, mVelViscosity;
int mIMax, mJMax;
int mNumGaussSeidelIterations;
bool mDensityDirichlet;
// Derived quantities;
int mIMaxM1, mJMaxM1, mIMaxP1, mJMaxP1, mKMaxP1, mNumPixels;
Real mDx, mDy;
Real mDxDx, mDyDy;
Real mHalfDivDx, mHalfDivDy;
Real mDtDivDx, mDtDivDy, mDtDivDxDx, mDtDivDyDy;
Real mEpsilon0, mEpsilonX, mEpsilonY;
Real mDenLambdaX, mDenLambdaY;
Real mVelLambdaX, mVelLambdaY;
Real mDenGamma0, mDenGammaX, mDenGammaY;
Real mVelGamma0, mVelGammaX, mVelGammaY;
// Current simulation time.
Real mTime;
// Lattice cell centers.
Real* mX;
Real* mY;
// State arrays.
Real** mDensity0;
Real** mDensity1;
Vector2<Real>** mVelocity0;
Vector2<Real>** mVelocity1;
Real** mDivergence;
Real** mPoisson;
};
#include "Wm5Fluid2Db.inl"
}
#endif
|