Till, without seeing the code in question, it's hard to say what's happening. However, does it make a difference if you use reordering *and* Triangulation::create_triangulation (not the _compatibility version)? Or do you mean here...
> boundary_vertices,BEM_cells, subcelldata); // (empty subcelldata object) > (I've tried without reordering and .create_triangulation(...) first.) ...that you've tried all four permutations and they all don't work? My best guess is that the triangulation is already invalid. One way to test this is to create a surface (dim,dim-1) triangulation from your once refined cube, refine the surface triangulation a couple of times, and output it to a ucd or vtk file. If you visualize this, does it look alright? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
