On 12/15/2016 06:11 PM, Jiaqi ZHANG wrote:
std::cout<<"testing quadrature points:\n";
const Point<dim> v0 = cell->vertex(0);
Quadrature<dim> q_point(v0);

Quadrature objects are defined on the reference cell [0,1]^d. So you will need to pass a point in reference coordinates to the quadrature object, but cell->vertex(0) only gives you a point in real space coordinates. You will have to transform it back into the reference coordinate system.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           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.

Reply via email to