An implementation of Helmholtz equation, for acoustics

Here is a code we have been working on.

Since you guys have much more experience in the appropriate ways to write weak forms (for performance, correctness, coding standards, or something else), I would appreciate a lot if you have a look at my implementation of the weak form of the Helmholtz equation, and comment on that, specially if I’m missing some naive points, which is entirely possible.

The file where you can find this weak form is this: hef/solver.py · master · CFD hunks / HEF-Acoustics · GitLab

Thank you so much!