I currently have the solution to a PDE stored in a Vector<double> solution. I need to compute the solution at some Point<2> loc. I accomplish this using the function
VectorTools::point_value(dof_handler, solution, loc, val); which fills val with the value I care about. Similarly, I get the gradient with point_gradient. This accomplishes what I need but if I have a lot of points like loc and have to do this every timestep then it quickly becomes a computational bottleneck. If I am not refining the mesh, is there a more efficient why of accomplishing this? I assume by storing the specific cell that loc is in? -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/fff649ee-8b5b-428b-95b9-bc7603a8afd6n%40googlegroups.com.
