search

inducer / pytential

Evaluate layer and volume potentials accurately. Solve integral equations.
https://pypi.python.org/pypi/pytential
25 stars 15 forks source link

Update Yukawa-Beltrami tests #207

Closed alexfikl closed 1 year ago

alexfikl commented 1 year ago
Managed to get the failing tests passing for me (ran it on koelsch), but not very satisfied with the results. Before, it was getting a small error but not converging h Error Running EOC
9.89001398e-01 6.15050357e-04
4.95775984e-01 1.32145031e-03 -1.11
Overall -1.11

and now it's doing

h Error Running EOC
9.89001398e-01 2.14125798e-02
4.95775984e-01 6.24579805e-03 1.78
Overall 1.78

which is a bit more converging, but not great. What do you think? I'm inclined to blame the myriad of hypersingular kernels in there.

Fixes #205

alexfikl commented 1 year ago

Might be worth nothing that the right preconditioned Yukawa-Beltrami operator worked with the previous settings, it's just the left-preconditioned one that was failing. This also adds the right-preconditioned ones as a slow test, since it takes about 10 minutes to run.

This is how the two versions look (the S0 IntG and friends use a Laplace kernel and the ones without a 0 use the Yukawa kernel)

https://github.com/inducer/pytential/blob/3a7807beec59f50c9d6046cf9b200d756915279d/pytential/symbolic/pde/beltrami.py#L180-L198

inducer commented 1 year ago

LGTM! Thanks for investigating.