diff --git a/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolation.C b/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolation.C index d7b24e1a7..d82cfa85b 100644 --- a/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolation.C +++ b/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolation.C @@ -49,7 +49,7 @@ void Foam::RBFInterpolation::calcB() const // Determine inverse of boundary connectivity matrix label polySize(4); - if(!polynomials_) + if (!polynomials_) { polySize = 0; } @@ -152,6 +152,8 @@ void Foam::RBFInterpolation::calcB() const // Collect ALL control points from ALL CPUs // Create an identical inverse for all CPUs + Info<< "Inverting RBF motion matrix" << endl; + BPtr_ = new scalarSquareMatrix(A.LUinvert()); } diff --git a/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolationTemplates.C b/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolationTemplates.C index e23b7ca59..f37264fb2 100644 --- a/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolationTemplates.C +++ b/src/OpenFOAM/interpolations/RBFInterpolation/RBFInterpolationTemplates.C @@ -72,7 +72,7 @@ Foam::tmp > Foam::RBFInterpolation::interpolate // 3) Return displacements using tresult() const label nControlPoints = controlPoints_.size(); - const scalarSquareMatrix& mat = B(); + const scalarSquareMatrix& mat = this->B(); // Determine interpolation coefficients Field alpha(nControlPoints, pTraits::zero);