I am getting an exception that I don't understand. I am trying to write a function that, given a point in your domain, finds the active cell which contains the point.
There is a function GridTools::find_active_cell_around_point that already does this, by the way.
The error is fairly deep in the stacktrace, and is an exception of type ExcAssignmentOfUnusedObject. The code actually runs for a while successfully before reaching this error, so it is clearly not every time that the function is called this problem arises.
This sounds like a memory corruption issue, as if the triangulation is changing while your function is running. Are you using multiple threads?
With these sort of problems it is really hard to find what's going on without seeing the complete code. Can you reduce it to something small?
Step31::BoussinesqFlowProblem<2>::run()
It's good to see that something is using this program :-) Are you using a Lagrangian advection method with it?
Best W. ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
