> On 30-Dec-2017, at 11:40 PM, luca.heltai <[email protected]> wrote: > > I’m thinking of a matrix of size (n_quadrature_points x n_active_cells) x > n_dofs, and then you slice the results cellwise instead of repeatedly > calling get_function_values. > > once: > > M[q+active_cell_index*n_dofs_per_cell, i] = fe_values.shape_value(i,q);
This looks like a Vandermonde matrix which would be identical on every cell. In that case, it is not necessary to store it for every cell. Thanks praveen -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
