36namespace cellCellStencils
54 label nD = donorCcs.
size();
76 d[j] = donorCcs[j] -
sample;
80 if (
mag(d[j]) < ROOTVSMALL)
86 LSQw[j] = 1.0/
magSqr(d[j]);
113 for (label i = 0; i < 3; i++)
120 W -= Row[j]*LSQw[j]*dw[i];
124 if (
mag(W) < SMALL*
mag(saveDiag))
136 - ATAinvAT[0][j]*LSQw[j]*dw[0]
137 - ATAinvAT[1][j]*LSQw[j]*dw[1]
138 - ATAinvAT[2][j]*LSQw[j]*dw[2]
static const Foam::dimensionedScalar A("", Foam::dimPressure, 611.21)
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
Minimal example by using system/controlDict.functions:
void setSize(const label n)
Alias for resize()
A templated (M x N) rectangular matrix of objects of <Type>, containing M*N elements,...
Singular value decomposition of a rectangular matrix.
scalarRectangularMatrix VSinvUt() const
Return the matrix product V S^(-1) U^T (the pseudo inverse)
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
void size(const label n)
Older name for setAddressableSize.
const Cmpt & z() const
Access to the vector z component.
const Cmpt & y() const
Access to the vector y component.
const Cmpt & x() const
Access to the vector x component.
Calculation of interpolation stencils.
Inverse-distance-weighted interpolation stencil.
virtual bool update()
Update stencils. Return false if nothing changed.
Least-squares-weighted interpolation stencil.
virtual ~leastSquares()
Destructor.
virtual void stencilWeights(const point &sample, const pointList &donorCcs, scalarList &weights) const
Calculate lsq weights for single acceptor.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
dimensioned< Type > T() const
Return transpose.
Mesh data needed to do the Finite Volume discretisation.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
static constexpr const zero Zero
Global zero (0)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
#define forAll(list, i)
Loop across all elements in list.