ContinuousQuadratureDataTransfer error

[Junxiang Wang]

Dear Deal.II community,

I am using the CellDataStorage class to record my history state variable and try to interpolate it onto my refined mesh. Therefore, the history data can be accessible for my analysis of the refined mesh.

However, I checked the data on the quadrature point on the refined mesh, It is not consistent with the data obtained by fe_values[introspection.extractors.heat_diffu].get_function_values(

 old_relevant_solution, old_timestep_heat_values); 

Can you give me some hints to debug it?

""""

 if (mesh_changed)

{
heat_pre = old_timestep_heat_values[q];

 cout << "data history" << std::endl;

cout << lqph[q]->get_initial("Temperature", 0.0) << std::endl;

cout << "solution vector" << std::endl;

}

"""

Thanks a lot in advance.

Best regards,

Junxiang.

[Wolfgang Bangerth]

Junxiang:

> I am using the CellDataStorage class to record my history state variable and
> try to interpolate it onto my refined mesh. Therefore, the history data can be
> accessible for my analysis of the refined mesh.

You might want to compare to compare with the implementation in this new code
gallery program:

https://dealii.org/developer/doxygen/deal.II/code_gallery_Phase_field_fracture_model_in_3D.html

> However, I checked the data on the quadrature point on the refined mesh, It is
> not consistent with the data obtained by
> fe_values[introspection.extractors.heat_diffu].get_function_values(
>  old_relevant_solution, old_timestep_heat_values);
>
> Can you give me some hints to debug it?

What does "not consistent" mean? What exactly is it that you expect, and what
exactly is it that you get?

Best
Wolfgang