Skip to content

Commit

Permalink
Modify the getNodalFieldSolution function in FEI.
Browse files Browse the repository at this point in the history
Previous to this change, the nodal loop would continue if the node was
a shared node and was not attached to the local element. This continue
was causing issues in parallel when the nodal data was being synched across
processor after a FETI solve when the run involved element death.
This created many issues and support tickets we have been working through for a
long time.

This was a great team effort with Alan, Kendall, Mike, and Patrick.
Alan and Kendall will work on pushing this change to Trilinos.
  • Loading branch information
alanw0 authored and jwillenbring committed Nov 16, 2017
1 parent 377094b commit 014ae92
Showing 1 changed file with 0 additions and 6 deletions.
6 changes: 0 additions & 6 deletions packages/fei/base/fei_FEDataFilter.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -2156,12 +2156,6 @@ int FEDataFilter::getNodalFieldSolution(int fieldID,
//next loop iteration.
if (!hasField) continue;

std::vector<GlobalID>& shNodeIDs = nodeCommMgr.getSharedNodeIDs();
int shIndex = fei::binarySearch(nodeID, &shNodeIDs[0], shNodeIDs.size());
if (shIndex > -1) {
if (!(problemStructure_->isInLocalElement(nodeNumber))) continue;
}

int dof_id = problemStructure_->getFieldDofMap().get_dof_id(fieldID, 0);

int offset = fieldSize*i;
Expand Down

0 comments on commit 014ae92

Please sign in to comment.