> Does some has a hint? Not really. A typical problem in this context is that queries like cell->is_inside(p) are badly posed in floating point arithmetic if the point p may be close to the boundary of the cell. In this case, depending on roundoff, the answer may be true or false.
Can you send me the simplest form of the program that shows this problem (in private mail)? I may take a look at it then. Best W. PS: On an unrelated note, you may want to add a 'break' to the if-clause in your for-loop -- presumably once you've found a cell, you are not interested in looking through the rest of cell_list any more. -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
