Symmetry for Helmholtz problems

[luiz....@eac.ufsm.br]

Hello all,

I am trying to solve Helmholtz problems that would really benefit from a symmetry condition to establish an infinite plane. Reading the documentation I didn't seem to find any native function to do so. I am not familiar with C++ to go inside the source code and alter the Helmholtz kernel, but if such solutions are said to be possible, I can attempt it.

Are there any known strategies to implement symmetry in Bempp without having to actually mesh the symmetrical geometry?

I appreciate the help!

Luiz Augusto Alvim

[Matthew Scroggs]

Hello,

I don't think this is something that Bempp can do. If you want to add a custom kernel to do this, I'd recommend using the new bempp-cl (https://github.com/bempp/bempp-cl). With this you could use a numba-compiled python based function to define your kernel which may be easier than diving into C++. Although there's currently no documentation on how to add kernels so it may be easier for now to just mesh the geometry.

All the best,

Matthew

[luiz....@eac.ufsm.br]

Hello Matthew,

Thanks for your response. I have finally managed to update all my code to run with numba kernels successfully and am currently attempting to define the custom kernels for the symmetry problem. The computational cost of doubling the mesh size is too great and in the long run, surely the symmetry condition could be of great use.

Thank you all for making such a powerful tool open for us to tweak and hopefully the custom kernels will be successfully implemented.

Best wishes,

Luiz Augusto Alvim