Hello, I need to calculate an integral of the squared gradient magnitude of the scalar I calculate using FEM.
Right now I iterate over all active cells and use fe_values class to get the gradients. After that I use a quadrature formula and for ever quadrature point add gradients[q_point].norm_square() * fe_values.JxW(q_point) to the total value of integral. Is this the correct way of calculate this integral? Is there a better (simpler code/more efficient) way of doing this? SY, Rostislav.
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
