From c59cb365b00ed633e4afb4f14d343b2e0bde99ef Mon Sep 17 00:00:00 2001 From: Hrvoje Jasak Date: Thu, 26 May 2011 11:29:40 +0100 Subject: [PATCH] Fix template issue --- .../interpolations/RBFInterpolation/RBFInterpolation.C | 4 +++- .../RBFInterpolation/RBFInterpolationTemplates.C | 2 +- 2 files changed, 4 insertions(+), 2 deletions(-) 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);