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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.