Fix template issue
This commit is contained in:
parent
7e994026ee
commit
c59cb365b0
2 changed files with 4 additions and 2 deletions
|
@ -49,7 +49,7 @@ void Foam::RBFInterpolation::calcB() const
|
||||||
// Determine inverse of boundary connectivity matrix
|
// Determine inverse of boundary connectivity matrix
|
||||||
label polySize(4);
|
label polySize(4);
|
||||||
|
|
||||||
if(!polynomials_)
|
if (!polynomials_)
|
||||||
{
|
{
|
||||||
polySize = 0;
|
polySize = 0;
|
||||||
}
|
}
|
||||||
|
@ -152,6 +152,8 @@ void Foam::RBFInterpolation::calcB() const
|
||||||
// Collect ALL control points from ALL CPUs
|
// Collect ALL control points from ALL CPUs
|
||||||
// Create an identical inverse for all CPUs
|
// Create an identical inverse for all CPUs
|
||||||
|
|
||||||
|
Info<< "Inverting RBF motion matrix" << endl;
|
||||||
|
|
||||||
BPtr_ = new scalarSquareMatrix(A.LUinvert());
|
BPtr_ = new scalarSquareMatrix(A.LUinvert());
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
|
@ -72,7 +72,7 @@ Foam::tmp<Foam::Field<Type> > Foam::RBFInterpolation::interpolate
|
||||||
// 3) Return displacements using tresult()
|
// 3) Return displacements using tresult()
|
||||||
|
|
||||||
const label nControlPoints = controlPoints_.size();
|
const label nControlPoints = controlPoints_.size();
|
||||||
const scalarSquareMatrix& mat = B();
|
const scalarSquareMatrix& mat = this->B();
|
||||||
|
|
||||||
// Determine interpolation coefficients
|
// Determine interpolation coefficients
|
||||||
Field<Type> alpha(nControlPoints, pTraits<Type>::zero);
|
Field<Type> alpha(nControlPoints, pTraits<Type>::zero);
|
||||||
|
|
Reference in a new issue