Dear All,
I have a question. I want to define a tensor with components like the below
expression, where eta is a scalar function in cylindrical coordinate
system. I defined this component in dealii like this:
FE_Q<dim> fe
const QGauss<dim> quadrature_formula(3)
FEValues<dim> fe_values (fe, quadrature_formula,
update_gradients|
update_quadrature_points);
const unsigned int n_q_points =
quadrature_formula.size();
std::vector< Tensor<1, dim> >
eta_gradients(n_q_points);
std::vector< Tensor<2, dim> >
eta_hessians(n_q_points);
fe_values.get_function_gradients(eta_gradients);
fe_values.get_function_hessians(eta_hessians);
for(unsigned int q_point = 0; q_point <
n_q_points; ++q_point)
const double comp11 = std::pow(eta, 2)*
std::pow((1-eta), 2)+(eta_hessians[q_point][0] +
eta_hessians[q_point][1]) - eta_gradients[q_point][0] *
eta_gradients[q_point][0];
I do really appreciate it if you let me know whether this part of code for
the following component is true or not.
.
<https://lh3.googleusercontent.com/-ki-xrwAfxYc/V8oOoAsTsmI/AAAAAAAAACo/uWArz5VCC1UsAWC4s9EcvnKVzMzp2P_5ACLcB/s1600/2.PNG>
--
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.
