Thanks a lot! I get the same answer by using quadrature.point(q), and I get
the similar answer after I transfer the real space point into reference
space( [0, 1]x[0, 1] ), and then use fe.shape_value(reference_point).
Best
Yidong
在 2018年5月7日星期一 UTC+8下午2:37:22，Wolfgang Bangerth写道：
>
> 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: [email protected]
> <javascript:>
> 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 [email protected].
For more options, visit https://groups.google.com/d/optout.
