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

Fixed: Clear reaction force component after it being added when a group

Browse Source 
of nodes is specified, such that nodes that are member in more than one
set are not added twice
This commit is contained in:
Knut Morten Okstad 2022-03-29 16:18:05 +02:00
parent 3fe7c9842f
commit 3ce4b7f369
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/LinAlg/SAM.C
View File

@ -976,9 +976,13 @@ Real SAM::getReaction (int dir, const RealArray& rf, const IntVec* nodes) const
if (!nodes || std::find(nodes->begin(),nodes->end(),i+1) != nodes->end())
{
int idof = madof[i]+dir-2;
if (idof < madof[i+1]-1)
if (msc[idof] < 0 && -msc[idof] <= (int)rf.size())
retVal += rf[-msc[idof]-1];
int ipr = -msc[idof]-1;
if (idof < madof[i+1]-1 && ipr >= 0 && ipr < (int)rf.size())
{
retVal += rf[ipr];
if (nodes) // clear this force component to avoid it being added twice
const_cast<RealArray&>(rf)[ipr] = 0.0;
}
}
return retVal;