On 10/25/19 11:06 AM, 'Andreas Hegendörfer' via deal.II User Group wrote: > > for (unsigned int i=0; i<cell->get_fe().dofs_per_cell; ++i) > { > strain=fe_values[displacement].symmetric_gradient (i, q_point) > *solution [local_dof_indices[i]]+strain;
You can do the first part simpler by using fe_values[displacement].get_function_symmetric_gradients(...) > } > strain_solution[local_dof_indices[q_point]]=strain; So here you put 'strain' into a solution vector. This solution vector must be associated with a finite element. And if it is associated... > In the programm I do have only one hp::DoFHandler<dim> dof_handler > overall. I am using the two elements: substructure_fe (FE_Q<dim>(2), > dim, FE_Nothing<dim>(), 1 ) and piezo_fe (FE_Q<dim>(2), dim, > FE_Q<dim>(2), 1). ...with this element, then you will necessarily have a continuous space into which you put the strains. If you don't like that, you'll have to create a different (discontinuous) FE, DoFHandler, and associated vector. > I save the result of the last calculation at one support point in > strain_solution[local_dof_indices[q_point]]. So the ohter calculation > results are lost but it should be no effort to save them, too. Correct. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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/1e523ccf-c490-a493-0599-d5b90e5c476d%40colostate.edu.