ivovellekoop/wavesim

High-accuracy simulation of light propagation

Wavesim - A fast and accurate method for solving the Helmholtz and time-independent Maxwell's equation

When using this code, please refer to:

G. Osnabrugge, S. Leedumrongwatthanakun, I.M. Vellekoop - A convergent Born series for solving the inhomogeneous Helmholtz equation in arbitrarily large media, Journal of Computational Physics Volume 322, 1 October 2016, Pages 113–124, doi:10.1016/j.jcp.2016.06.034 Freely available at: http://www.sciencedirect.com/science/article/pii/S0021999116302595

G. Osnabrugge, M. Benedictus, I.M. Vellekoop - An ultra-thin boundary layer for high-accuracy simulations, Optics Express 29 (2), 11 January 2021, Pages 1649-1658, doi:10.1364/OE.412833

wavesim.org

We are working to improve and accelerate Wavesim further. On 27 February 2024, we released a new version that uses a CUDA-based acceleration module (cumex) to provide around 2x speed up. Want to find out more? Want to participate in the forum for discussions, queries, and requests? Then please visit www.wavesim.org.

wavesim's Issues

ARL, PBL is bugged

Excellent work. However I tried to run the examples are given, and there are problems with nonperiodic conditions crash (ARL,PBL).
I would love it if you could fix it.

"Index exceeds the number of array elements (128).

Error in Simulation.default_callback (line 303)
sroi = ax(obj.roi(:,2));

Error in Simulation/next (line 212)
obj.callback(obj, state);

Error in WaveSimBase/run_algorithm (line 160)
state = next(obj, state, can_terminate);

Error in Simulation/exec (line 178)
state = run_algorithm(obj, state);

Error in pulse2d (line 36)
E_set(:,:,f) = exec(sim, source);"