Commit 00fc1298 authored by Thomas Witkowski's avatar Thomas Witkowski

da

parent 43c4b5c2
......@@ -216,16 +216,8 @@ namespace AMDiS {
MeshStructure elementStructure;
elementStructure.init(structureCode, codeSize);
if (macroInMesh.count(elIndex)) {
MSG("READ MACRO EL %d\n", elIndex);
if (MPI::COMM_WORLD.Get_rank() == 52 && elIndex == 19777) {
debug::writeMacroElementIndexMesh(mesh, "test");
}
if (macroInMesh.count(elIndex))
elementStructure.fitMeshToStructure(mesh, refManager, false, elIndex);
}
if (nValueVectors > 0) {
uint32_t nValuesPerVector = 0;
......
......@@ -717,12 +717,16 @@ namespace AMDiS {
{
FUNCNAME("PetscSolverFeti::testWirebasketEdge()");
return true;
if (meshDistributor->getMesh()->getDim() == 2)
return true;
if (meshDistributor->getIntBoundary().getDegreeOwn(edge) != 3)
return false;
return true;
Element *el = edge.el;
int i0 = el->getVertexOfEdge(edge.ithObj, 0);
int i1 = el->getVertexOfEdge(edge.ithObj, 1);
......@@ -786,7 +790,7 @@ namespace AMDiS {
if (it->rankObj.subObj == FACE)
allFaces[it.getRank()].insert(it->rankObj);
#if 0
std::set<BoundaryObject> allMyEdges;
TraverseStack stack;
ElInfo *elInfo = stack.traverseFirst(meshDistributor->getMesh(), 0, Mesh::CALL_EL_LEVEL | Mesh::FILL_BOUND);
......@@ -818,6 +822,7 @@ namespace AMDiS {
}
}
}
#endif
int nEmptyFaces = 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