Scattering_BIE_QPAX

Quadrature by Parity eXpansions method for computing the scattering by a high aspect ratio ellipse using boundary integral equation methods.

Reference

Codes are associated with the article:

• C. Carvalho, A. D. Kim, L. Lewis, and Z. Moitier, Quadrature by Parity Asymptotic eXpansions (QPAX) for scattering by high aspect ratio particles, SIAM Multiscale Modeling & Simulation (2021). [arXiv, DOI]

and the proceeding:

• C. Carvalho, A. D. Kim, and Z. Moitier, Quadrature by Parity Asymptotic eXpansions (QPAX) for light scattering by high aspect ratio plasmonic particle, SIAM Multiscale Modeling & Simulation (2021). [Submited]

The bibtex entry can be access by the command.

python3 -c "from src import __bibtex__; print(__bibtex__)"

Requirements

• Python version:

  • Tested on Python 3.9;
  • Should works on Python ≥3.7 but not tested.

• Require the following libraries:

  • For numerical computations: NumPy and SciPy;

  • For symbolic computations: SymPy;

  • For visualization: Matplotlib and seaborn;

  • For Jupyter notebook: IPython and JupyterLab.

Might works with previous versions of the libraries but if it does not works try to update the libraries for example through pip

python3 -m pip install --user --upgrade -r requirements.txt

and for the dev dependencies

python3 -m pip install --user --upgrade -r requirements-dev.txt

Install

Clone from GitHub repository:

git clone https://github.com/zmoitier/Scattering_BIE_QPAX.git

Instructions for usage

Run the Jupyter notebook notebook_article/run_fig_n.ipynb to get the corresponding figure n in the article and run the Jupyter notebook notebook_proceeding/run_fig_n.ipynb to get the corresponding figure n in the proceeding.

Symbolic expressions

The Symbolic_expansions folder provides inner expansions for Laplace and Helmholtz double-layer potentials using Mathematica notebook or SymPy (see section 3 and 4 or the associated manuscript).

Contact

If you have any questions or suggestions please feel free to create an issue in this repository.