Hello everyone!
How to obtain the value of each quadrature point locally with reference to
q_index, (I made this: std::cout << fe_values.quadrature_point <<, in
step-3 and step-8, since q_index is looping there), when loops each cell,
if my 1D PDE's right-hand-side is a function of x, f(x)?
In step-3 and step-8, I set to Q1 element by switching a parameter in Class
Constructor. When I std::cout << q_index <<, I saw 2 indices {0, 1}
returned, does that mean 2 quadrature points are used? I also std::cout <<
fe_values.JxW(q_index) <<, that gave me 0.005 corresponding to both
indices, so I guess 0.5 for quadrature weights, since I meshed 10 cells in
1D domain of (0, 0.1). However, when I std::cout <<
fe_values.quadrature_point <<, I did not see expected values {0, 1} locally
for weights {0.5, 0.5}, when loops each cell.
I must be wrong in guessing fe_values.JxW(q_index) and
fe_values.quadrature_point. Does step-3 or step-8 use 2 quadrature points
for Q1 element for 1D PDE, and how to obtain the value of each quadrature
point locally with reference to q_index, if I want to build f(x), in terms
of quadrature point locally?
Thank you all!
Judy
--
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/42ec8193-5353-4b81-b6b6-da5cae320c08n%40googlegroups.com.
