Dear Katie, internally FEFieldFunction uses FindCellAroundPoint. The exception that is thrown indicates that that function failed to locate the point.
Can you reproduce this behavior in a 10 lines snippet of code? Thanks, Luca. -- Luca Heltai <[email protected]> http://people.sissa.it/~heltai/ Scuola Internazionale Superiore di Studi Avanzati Phone: +39 040 3787 449, Office: 732 -- There are no answers, only cross references Begin forwarded message: > From: Katie Leonard <[email protected]> > Date: June 15, 2011 3:36:02 PM GMT+02:00 > To: "[email protected]" <[email protected]> > Subject: Adaptive meshes and Interpolation using FEFieldFunction deal.II > > Dear Luca, > > I am using the FEFieldFunction to interpolate a solution from one grid to > another. When I don't have any refinement on the grid (i.e. no hanging > nodes) then the interpolation works fine. However, if I introduce refinement > then I get the error: > > The violated condition was: > best_cell.first.state() == IteratorState::valid > The name and call sequence of the exception was: > ExcPointNotFound<dim>(p) > Additional Information: > The point <-0.0468750 -0.0375000> could not be found inside any of the > subcells of a coarse grid cell. > but this point is within the triangulation. > > > Currently I am using a solutions with switch-like functions (functions of > tanh). When I reduce the steepness of the solution then I can re-introduce > the adaptivity up to a certain point. > > I was wondering why this is? I assumed that FEFieldfunction could be used > with adaptive meshes. > > Many thanks, > > K Leonard
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
