On 05/07/2018 02:20 PM, Zhao Yidong wrote:
2.Iterate on every quadrature point in every cell, store them. (These points
are in real space right?)
|
for(unsignedq=0;q<quadrature.size();++q){
Point<dim>point =fe_values.quadrature_point(q);
|
3. Use fe.shape_value() to get shape value at every degree of a cell,
store_point is the point I store in the 2nd step.
|
doubleshapeValue =fe.shape_value(i,store_point);// shape value
|
This is your problem. 'store_point' lives in real space, but 'fe.shape_value'
wants a point on the reference cell. You need to call
fe.shape_value (quadrature.point(q));
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.