Commit 70c8d9cc authored by Thomas Witkowski's avatar Thomas Witkowski
Browse files

* NOTHING

parent 6577015f
...@@ -88,7 +88,7 @@ namespace AMDiS { ...@@ -88,7 +88,7 @@ namespace AMDiS {
neighbourCoord_[i][0] = macroNeighbour->coord[1]; neighbourCoord_[i][0] = macroNeighbour->coord[1];
neighbourCoord_[i][1] = macroNeighbour->coord[2]; neighbourCoord_[i][1] = macroNeighbour->coord[2];
} else { } else {
ERROR_EXIT("should not happen!\n"); ERROR_EXIT("Should not happen!\n");
} }
neighbourCoord_[i][2] = oppCoord_[i]; neighbourCoord_[i][2] = oppCoord_[i];
...@@ -102,15 +102,28 @@ namespace AMDiS { ...@@ -102,15 +102,28 @@ namespace AMDiS {
neighbourCoord_[i][0] = macroNeighbour->coord[2]; neighbourCoord_[i][0] = macroNeighbour->coord[2];
neighbourCoord_[i][1] = macroNeighbour->coord[0]; neighbourCoord_[i][1] = macroNeighbour->coord[0];
} else { } else {
ERROR_EXIT("should not happen!\n"); ERROR_EXIT("Should not happen!\n");
} }
neighbourCoord_[i][2] = oppCoord_[i]; neighbourCoord_[i][2] = oppCoord_[i];
break; break;
case 2:
if (*(macroNeighbour->getElement()->getDOF(2)) == *(element_->getDOF(0))) {
neighbourCoord_[i][0] = macroNeighbour->coord[2];
neighbourCoord_[i][1] = macroNeighbour->coord[1];
} else if (*(macroNeighbour->getElement()->getDOF(2)) == *(element_->getDOF(1))) {
neighbourCoord_[i][0] = macroNeighbour->coord[0];
neighbourCoord_[i][1] = macroNeighbour->coord[2];
} else {
ERROR_EXIT("Should not happen!\n");
}
break;
default: default:
std::cout << "------------- Error --------------" << std::endl; std::cout << "------------- Error --------------" << std::endl;
std::cout << " Element index = " << element_->getIndex() << "\n\n"; std::cout << " Neighbour counter = " << i << "\n";
std::cout << " Element index = " << element_->getIndex() << "\n\n";
for (int j = 0; j < neighbours; j++) { for (int j = 0; j < neighbours; j++) {
if (mel->getNeighbour(j)) { if (mel->getNeighbour(j)) {
std::cout << " Neighbour " << j << ": " std::cout << " Neighbour " << j << ": "
......
...@@ -1359,6 +1359,8 @@ namespace AMDiS { ...@@ -1359,6 +1359,8 @@ namespace AMDiS {
// and now partition the mesh // and now partition the mesh
partitionMesh(adaptInfo); partitionMesh(adaptInfo);
#if 0
// 2009-03-31 THOMAS, DO NOT REMOVE THIS BLOCK
TraverseStack stack; TraverseStack stack;
ElInfo *elInfo = stack.traverseFirst(mesh, -1, Mesh::CALL_LEAF_EL); ElInfo *elInfo = stack.traverseFirst(mesh, -1, Mesh::CALL_LEAF_EL);
int nLeaves = 0; int nLeaves = 0;
...@@ -1377,8 +1379,9 @@ namespace AMDiS { ...@@ -1377,8 +1379,9 @@ namespace AMDiS {
} }
mesh->setNumberOfLeaves(nLeaves); mesh->setNumberOfLeaves(nLeaves);
#endif
#if 0 #if 1
globalRefineOutOfPartition(adaptInfo); globalRefineOutOfPartition(adaptInfo);
refineOverlap(adaptInfo); refineOverlap(adaptInfo);
......
...@@ -25,6 +25,7 @@ ...@@ -25,6 +25,7 @@
#include "Lagrange.h" #include "Lagrange.h"
#include "Flag.h" #include "Flag.h"
#include "TraverseParallel.h" #include "TraverseParallel.h"
#include "VtkWriter.h"
namespace AMDiS { namespace AMDiS {
...@@ -573,6 +574,8 @@ namespace AMDiS { ...@@ -573,6 +574,8 @@ namespace AMDiS {
{ {
FUNCNAME("ProblemVec::estimate()"); FUNCNAME("ProblemVec::estimate()");
VtkWriter::writeFile(solution->getDOFVector(0), "test.vtu");
clock_t first = clock(); clock_t first = clock();
#ifdef _OPENMP #ifdef _OPENMP
......
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