Hello everyone, I am trying to implement a continuous mechanics simulation and I am encountering troubles trying to compute the divergence of the deformation gradient inverse. I need it to assemble the variational form terms. My code reads something like the following:
- take the gradient of the previous solution *fe_values[components].get_function_gradients(old_sol, local_old_grad);* - for every quadrature point "q" compute the deformation gradient (F = Identity + grad(u)) *Tensor<2, dim> F = Identity + local_old_grad(q);* *Tensor<2, dim> F_inv = invert(F);* Now I'd like to compute the divergence of *F_inv*, is there a way to do that? I was thinking about automatic differentiation but maybe there is a more straightforward way Thank you -- 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 visit https://groups.google.com/d/msgid/dealii/5e5ed814-540d-4f05-93b9-09e087f7001en%40googlegroups.com.