> Could someone please help me understand where I go wrong. The error can be 
> easily reproduced by setting /fePolyDegree >1 in the attached MWE.

Paras -- it's a bug in deal.II: The assertion you run into is a bit 
overzealous and flags something that is perfectly legitimate.

I'll fix the actual problem later this weekend, but here is a workaround for 
you: Where you call the build_patches() function, replace
    dataOutPrimary.build_patches(qMapping, nSubdivisions);
by
    dataOutPrimary.build_patches(qMapping, nSubdivisions,
                                 DataOut<dim>::curved_inner_cells);
and your code should work (doing essentially the same as before -- maybe even 
in a more correct way).

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.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/09d9297d-fc3d-14a6-5978-135b2e1ba3b6%40colostate.edu.

Reply via email to