> I'm wondering if it's possible/straightforward to have per-cell quadrature > weights? I want to keep the quadrature positions the same but update the > quadrature weights.
If you only have a few different sets of quadratures, I would think having one fe_values object for each quadrature might be an easy solution. If not, here is something else to try: I think none of the data inside fe_values except .JxW() depends on the quadrature weights, just the locations. So to change the weights you could use fe_values.JxW(q)/old_weight[q]*new_weight[q] instead of fe_values.JxW(q) There are a few exceptions like get_function_values() that depend on the quadrature weights, but if you only want to do assembly, this hack should work. -- Timo Heister http://www.math.clemson.edu/~heister/ -- 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. For more options, visit https://groups.google.com/d/optout.