OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

bugfix for mixed spline fields: ip must be cleared before calling scatterInd for the second basis

Browse Source 
git-svn-id: http://svn.sintef.no/trondheim/IFEM/trunk@1349 e10b68d5-8a6e-419e-a041-bce267b0401d
This commit is contained in:
kmo 2011-12-17 14:14:20 +00:00 committed by Knut Morten Okstad
parent d5036ab66d
commit 554594ed2b
4 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
src/ASM/SplineField2D.C
View File

@ -118,6 +118,7 @@ bool SplineField2D::gradFE (const FiniteElement& fe, Vector& grad) const
}
dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
ip.clear();
ASMs2D::scatterInd(basis->numCoefs_u(),basis->numCoefs_v(),
basis->order_u(),basis->order_v(),
spline.left_idx,ip);

1
src/ASM/SplineField3D.C
View File

@ -123,6 +123,7 @@ bool SplineField3D::gradFE (const FiniteElement& fe, Vector& grad) const
}
dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
ip.clear();
ASMs3D::scatterInd(basis->numCoefs(0),basis->numCoefs(1),basis->numCoefs(2),
basis->order(0),basis->order(1),basis->order(2),
spline.left_idx,ip);

1
src/ASM/SplineFields2D.C
View File

@ -125,6 +125,7 @@ bool SplineFields2D::gradFE (const FiniteElement& fe, Matrix& grad) const
}
dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
ip.clear();
ASMs2D::scatterInd(basis->numCoefs_u(),basis->numCoefs_v(),
basis->order_u(),basis->order_v(),
spline.left_idx,ip);

1
src/ASM/SplineFields3D.C
View File

@ -130,6 +130,7 @@ bool SplineFields3D::gradFE(const FiniteElement& fe, Matrix& grad) const
}
dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
ip.clear();
ASMs3D::scatterInd(basis->numCoefs(0),basis->numCoefs(1),basis->numCoefs(2),
basis->order(0),basis->order(1),basis->order(2),
spline.left_idx,ip);