Jane, [...] > > Now, > VectorTools::interpolate(dof_handler, exact_solution, solution); > > gives me a dimension mismatch error - I'm assuming that this is due to the > fact that RT elements can't distinguish between the x and y components. > I would expect that RT elements are not interpolatory and expect a corresponding error message. What exactly is the error message you get in debug mode?
> > I then tried adding a component mask in various ways to put the velocity > and pressure solutions in separately, but I still get the same error. > Doing a component mask with solution.block(0) separately also gives me a > number of dofs mismatch error. > > Is there any way around this? The above interpolation works if i have Fe_Q > elements (like in step-22) so I've narrowed it down to the fact that it is > because of the function space, but I need the function spaces as they are. > > I can't seem to find another way around this using the interpolation > function, and hope that there is actually a simple way around this? > Try to use VectorTools::project <https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceVectorTools.html#af282602c17b8eb6afe47ef56fc417ecb> instead. Best, Daniel -- 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. For more options, visit https://groups.google.com/d/optout.
