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!
Best,
--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
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.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/CAO2HPGW8067gaFaYx2QjVHc719xdUVSh%3Dyg2UdhrFa4U38cn5Q%40mail.gmail.com.
