diff --git a/AMDiS/src/ProblemVec.cc b/AMDiS/src/ProblemVec.cc index f76dafcaa79cf6dca4eb31b08bdadb73c0fe9c0f..7824e8c208b5035596f5d581745a5876cd7ebce8 100644 --- a/AMDiS/src/ProblemVec.cc +++ b/AMDiS/src/ProblemVec.cc @@ -1118,16 +1118,10 @@ namespace AMDiS { it != op->getAuxFeSpaces().end(); ++it) { if ((*it)->getMesh() != componentSpaces[i]->getMesh() || (*it)->getMesh() != componentSpaces[j]->getMesh()) { - // std::cout<<"auxFE "<<k<<" in position ("<<i<<","<<j<<") has different mesh\n"<<std::endl; op->setNeedDualTraverse(true); break; } } - /* std::cout<<"after adding MatOp to ("<<i<<","<<j<<"):\n"; - traverseInfo.updateStatus(); - for(int i2=0;i2<nComponents;i2++) - for(int j2=0;j2<nComponents;j2++) - std::cout<<"traverseInfo.getStatus("<<i2<<","<<j2<<"):"<<traverseInfo.getStatus(i2,j2)<<std::endl;*/ } diff --git a/AMDiS/src/SubAssembler.cc b/AMDiS/src/SubAssembler.cc index ad6a82dad73122a9f034d5dfb2cad04a3a66b76b..d1a70e05af215588cad2daa8a0248ca4b56bec78 100644 --- a/AMDiS/src/SubAssembler.cc +++ b/AMDiS/src/SubAssembler.cc @@ -229,17 +229,7 @@ namespace AMDiS { WorldVector<double> *values = &(gradientsAtQPs[vec]->values[0]); gradientsAtQPs[vec]->valid = true; - bool hasNan=false; - for(int i=0;i<localQuad->getNumPoints() && !hasNan;i++) - for(int j=0;j<Global::getGeo(WORLD) && !hasNan ;j++) - hasNan=isnan((values[i])[j]); - TEST_EXIT_DBG(!hasNan)("gradient contains a nan\n"); vec->getGrdAtQPs(smallElInfo, largeElInfo, localQuad, NULL, values); - hasNan=false; - for(int i=0;i<localQuad->getNumPoints() && !hasNan;i++) - for(int j=0;j<Global::getGeo(WORLD) && !hasNan ;j++) - hasNan=isnan((values[i])[j]); - TEST_EXIT_DBG(!hasNan)("gradient contains a nan\n"); return values; } diff --git a/AMDiS/src/ZeroOrderAssembler.cc b/AMDiS/src/ZeroOrderAssembler.cc index 351679d6fb51a12d22d72e54ba662b4ac1bfef67..e4b85bd719b547945ac9ea799f23424db07a064e 100644 --- a/AMDiS/src/ZeroOrderAssembler.cc +++ b/AMDiS/src/ZeroOrderAssembler.cc @@ -99,10 +99,6 @@ namespace AMDiS { for (termIt = terms[myRank].begin(); termIt != terms[myRank].end(); ++termIt) (static_cast<ZeroOrderTerm*>((*termIt)))->getC(elInfo, nPoints, c); - bool hasNan=false; - for(int i=0;i<nPoints && !hasNan ;i++) - hasNan=isnan(c[i]); - TEST_EXIT(!hasNan)("there is a nan in at least one position\n"); if (symmetric) { TEST_EXIT_DBG(nCol == nRow)("nCol != nRow, but symmetric assembling!\n"); @@ -138,11 +134,6 @@ namespace AMDiS { } } } - hasNan=false; - for(unsigned i=0;i<num_rows(mat) && !hasNan;i++) - for(unsigned j=0;j<num_cols(mat) && !hasNan; j++) - hasNan=isnan(mat[i][j]); - TEST_EXIT(!hasNan)("calculated an element matrix containing a nan\n"); }