Dear all, I was trying to write a roots finding method along cell faces in 2D, so I need to know the solution values at any point of cell faces. Is this the right way to use Quadrature <https://www.dealii.org/8.4.1/doxygen/deal.II/classQuadrature.html>< dim >:: Quadrature <https://www.dealii.org/8.4.1/doxygen/deal.II/classQuadrature.html>(const Point <https://www.dealii.org/8.4.1/doxygen/deal.II/classPoint.html>< dim > & point)?
std::cout<<"testing quadrature points:\n"; const Point<dim> v0 = cell->vertex(0); Quadrature<dim> q_point(v0); FEValues<dim> fe_p_values (mapping(),fe, q_point, update_values | update_quadrature_points ); fe_p_values.reinit(cell); const std::vector<Point<dim>> &q_v0 = fe_p_values.get_quadrature_points(); std::cout<<" v0 = "<<v0 <<" qv0 = "<< q_v0[0]<<"\n\n"; But it gives me the following results: testing quadrature points: v0 = 0 0 qv0 = 0 0 testing quadrature points: v0 = 0.015625 0 qv0 = 0.0158691 0 testing quadrature points: v0 = 0 0.015625 qv0 = 0 0.0158691 testing quadrature points: v0 = 0.015625 0.015625 qv0 = 0.0158691 0.0158691 testing quadrature points: v0 = 0.03125 0 qv0 = 0.0317383 0 testing quadrature points: v0 = 0.046875 0 qv0 = 0.0476074 0 testing quadrature points: v0 = 0.03125 0.015625 qv0 = 0.0317383 0.0158691 Why can I not get the same points? Any help would be appreciated! Thanks in advance. Jiaqi -- 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.