search

inducer / modepy

Modes and nodes for high-order discretizations
Other
16 stars 7 forks source link

Mass matrix quad #83

Closed inducer closed 5 months ago

inducer commented 5 months ago

This adds two new mass matrix functions, plus a slew of more-or-less related cleanups. Best read commit-by-commit.

inducer commented 5 months ago

Thanks for taking a look!

Also, regarding that lumping in grudge. Is it possible with this interface? If I understand correctly, it would just mean passing a Gauss-Legendre quadrature or a Lobatto-G-L to integrate in nodal/modal_quad_mass_matrix?

I'm planning to add a separate function that returns just the diagonal, after appropriate checking. But this seemed like a nice, self-contained set of changes, so I figured I'd do that separately. I'm also still puzzling over the interface.