ComputeMultiplier<dim>::get_names() const
{
std::vector<std::string> solution_names;
- solution_names.push_back ("Gradient norm");
- solution_names.push_back ("Lagrange multiplier");
+ solution_names.push_back ("Gradient_norm");
+ solution_names.push_back ("Lagrange_multiplier");
return solution_names;
}
triangulation.execute_coarsening_and_refinement();
dof_handler.distribute_dofs(fe);
# if DEAL_II_VERSION_GTE(9, 7, 0)
+ present_solution.reinit(dof_handler.n_dofs());
solution_transfer.interpolate(present_solution);
# else
Vector<double> tmp(dof_handler.n_dofs());
);
};
- auto kappa_2 = [=](double T, double lambda){
+ auto kappa_2 = [=](double T, double /*lambda*/){
return -problem.k * std::exp(-problem.theta / T) * Heaviside_func(T - problem.T_ign);
};
setup_system(/*initial_step=*/ false);
# if DEAL_II_VERSION_GTE(9, 7, 0)
+ current_solution.reinit(dof_handler.n_dofs());
solution_transfer.interpolate(current_solution);
# else
Vector<double> tmp(dof_handler.n_dofs());
# if DEAL_II_VERSION_GTE(9, 7, 0)
// stress
+ for (unsigned int i=0; i<dim; ++i)
+ for (unsigned int j=0; j<dim; ++j)
+ {
+ history_stress_field[i][j].reinit(history_dof_handler.n_dofs());
+ }
+
history_stress_field_transfer0.interpolate(history_stress_field[0]);
if ( dim > 1)
{
}
// strain
+ for (unsigned int i=0; i<dim; ++i)
+ for (unsigned int j=0; j<dim; ++j)
+ {
+ history_strain_field[i][j].reinit(history_dof_handler.n_dofs());
+ }
+
history_strain_field_transfer0.interpolate(history_strain_field[0]);
if ( dim > 1)
{