From a1095c6a4fccbbacc9036bdd606bc6080a59858c Mon Sep 17 00:00:00 2001
From: Andreas Naumann <andreas.naumann@tu-dresden.de>
Date: Mon, 7 Jun 2010 13:47:48 +0000
Subject: [PATCH] removed debug code
---
AMDiS/src/ProblemVec.cc | 6 ------
AMDiS/src/SubAssembler.cc | 10 ----------
AMDiS/src/ZeroOrderAssembler.cc | 9 ---------
3 files changed, 25 deletions(-)
diff --git a/AMDiS/src/ProblemVec.cc b/AMDiS/src/ProblemVec.cc
index f76dafca..7824e8c2 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 ad6a82da..d1a70e05 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 351679d6..e4b85bd7 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");
}
--
GitLab