On 07/19/2017 04:09 PM, Jaekwang Kim wrote:
because my dof_handler has doc for pressure, size does not match...
how can I avoid this problem ?
You just need to pass the *entire* solution vector, not just the first
block. Otherwise, DataOut cannot find out which elements of the vector
correspond to the degrees of freedom enumerated in your DoFHandler.
You will also have to change this:
Magnitude<dim>::compute_derived_quantities_vector ([...])
{
[...]
Assert(uh[i].size() == 2, ExcDimensionMismatch (uh[i].size(), 2));
[...]
The size will now be 3 (or probably dim+1), even though you are only
using the zeroth and first component of uh in computing the magnitude.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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.