[deal.II] get_function_2nd_derivatives vs get_function_hessians

[misabel]

Dear DealII's users:

I am working with two versions of DealII: version 6.0.0 and version 6.1.0.

For both cases I use:

std::vector< std::vector< Tensor< 2, DIMENSION  > > > second_derivatives;
second_derivatives.resize(n_face_q_points);
for (unsigned int i=0;i<n_face_q_points;i++)
                second_derivatives[i].resize(DIMENSION+1);

In order to get the second derivatives of a vector with the first version
(6.0.0), I must use
“fe_face_values.get_function_2nd_derivatives(Bvector_solution,
second_derivatives)”. Using the second version (6.1.0), I have to use
“fe_face_values.get_function_hessians(Bvector_solution,second_derivatives)”

I thought the results I would obtain would be the same ones, but using
“get_function_2nd_derivatives()”, my program breaks up.

I use the variable “second_derivatives” in the following way for both cases:

cell_rhs(i) += viscosity*fe_face_values.JxW(q_point)*
                (fe_face_values.shape_value(i, q_point) *
                 ((fe_face_values.normal_vector(q_point)(0)*
                (derivadas_segundas[q_point][0][0][0]+
                 derivadas_segundas[q_point][1][1][0]))
                +
                (fe_face_values.normal_vector(q_point)(1)*
                (derivadas_segundas[q_point][0][0][1]+
derivadas_segundas[q_point[1[1[1])))
                );

I would be grateful if you could give me some advice.
Thanks in advance and best regards!
Isa


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii