On Thu, Jul 23, 2020 at 5:14 PM Daniel Arndt <d.arndt.m...@gmail.com> wrote:

> You can do similarly,
> Quadrature q(fe.get_unit_support_points());
> FEValues fe_values (..., q, update_q_points);
> for (const auto& cell)
>   ...
>   points = fe_values.get_quadrature_points();
>   fe_values.get_function_values(values);
>   fe_values.get_function_gradients(gradients);
>   for (unsigned int i=0; i<n_dofs_per_cell;  ++i)
>     v(local_dof_indices[i]) = f(points[i],  values(i), gradients(i));
> Yes! Declaring a quadrature on support points is a great solution. Thanks
a lot!

Stay put, practice social distancing, and be safe!


--Xuefeng Li, (504)865-3340(phone)
   Like floating clouds, the heart rests easy
   Like flowing water, the spirit stays free
      Loyola University New Orleans
   New Orleans, Louisiana (504)865-2051(fax)

The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
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 

Reply via email to