Hi all,

I'm solving a phase field fracture problem and need to store the history of 
elastic energy at each Gauss points. I wonder if there is anyway to 
transfer( or interpolate) the energy stored at Gauss points of the old mesh 
to the new mesh during mesh refinement operation. The snippet code below 
shows how I store history of energy (called 'HistoryField' variable).

for (auto &cell : dof_handler.active_cell_iterators())
{
if (cell->is_locally_owned())
{
PointHistory<dim> *local_quadrature_points_history =
reinterpret_cast<PointHistory<dim> *>(cell->user_pointer());
for(int q=0; q< quadrature_formula.size(); q++)
{
local_quadrature_points_history[q].HistoryField = 0.0;
local_quadrature_points_history[q].max_Ybar = 0.0;

}
}
}

-- 
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.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/a1feed80-74a9-458d-bd05-82a3ca30be24n%40googlegroups.com.

Reply via email to