Commit 1b0e1b16 authored by Thomas Witkowski's avatar Thomas Witkowski
Browse files

Do not assemble zero operators.

parent 51092637
......@@ -274,7 +274,8 @@ namespace AMDiS {
std::vector<Operator*>::iterator it = operators.begin();
std::vector<double*>::iterator factorIt = operatorFactor.begin();
for (; it != operators.end(); ++it, ++factorIt)
if ((*it)->getNeedDualTraverse() == false)
if ((*it)->getNeedDualTraverse() == false &&
(*factorIt == NULL || **factorIt != 0.0))
(*it)->getElementMatrix(elInfo, elementMatrix, *factorIt ? **factorIt : 1.0);
if (factor != 1.0)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment