Re: [deal.II] Re: Computing the solution gradient at the quadrature point on a face

2021-08-01 Thread Wolfgang Bangerth
On 8/1/21 12:05 PM, Michael Lee wrote: fe_face_values.*get_function_gradients *returns a rank 2 tensor as the solution is the displacement vector. No, it returns a vector of Tensor<1,dim> objects. If you want the gradients of the displacement vector, you need

Re: [deal.II] Computing the solution gradient at the quadrature point on a face

2021-08-01 Thread Jean-Paul Pelteret
Hi Michael, The problem here is that “scratch.solution_grads_u_total” is initialised with the cell FEValues, and so has n_cell_quadrature_points entries. You need to pass "scratch.fe_face_values_ref[u_fe].get_function_gradients(scratch.solution_total, scratch.solution_grads_u_total);”

[deal.II] Re: Computing the solution gradient at the quadrature point on a face

2021-08-01 Thread Michael Lee
fe_face_values.*get_function_gradients *returns a rank 2 tensor as the solution is the displacement vector. Let me be more specific. In the code gallery " Quasi-Static Finite-Strain Compressible Elasticity (The deal.II Library: The 'Quasi-Static Finite-Strain Compressible Elasticity' code

[deal.II] Re: Computing the solution gradient at the quadrature point on a face

2021-08-01 Thread Peter Munch
The return type should be `*Tensor<2,dim,NumberType>*` shouldn't it? Hope this helps. If not, could describe how the code does not work? It it not compiling? Is the result wrong? PM On Sunday, 1 August 2021 at 04:09:45 UTC+2 lian...@gmail.com wrote: > Dear All, > > I want to do a surface