Hi James, perhaps I should clarify one more thing: For hanging nodes the number and position of neighbors at a given face of a cell is the relevant thing: If there is none, then we are at the boundary and do not have hanging nodes, if there is one neigbor then everything is implemented (and only important in the hp-case). If there are four (in 3D), then locally the refinement is isotropic and this is also implemented. We only face a new situation if the number of neighbors is two or three, these cases are not implemented. If you use your aspect ratio algorithm on a hyper-rectangle, all cells are split in the same way and you end up with a mesh where each cell has one neighbor at each face (except for the boundary , of course). Thus this case works. On a general mesh like the one you created externally, however, neighbors will be refined differently and thus you will end up with the cases that are not implemented.
Best, Tobias _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
