> 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?
That's the way to go. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
