Maxim Yurkin

Here are two other projects, which struggle with the same problem with some success. So we may take some code or ideas from that: - [SMUTHI](https://pypi.org/project/SMUTHI/) - [IF-DDAM](https://github.com/pchaumet/IF-DDAM) But further...

Another relevant project (suggested by @AmosEgel) is https://github.com/modelics/strata . Probably we can just link to it for computation of the Green's tensor for multi-layered substrates.

r1304 - a85f1ce3c2a271722c78140a7476af40b4b547c0

/cc @fabiods

Gaussian densities have been used in Lim J., Kang S., Kim J., Kim W.Y., and Ryu S. Non-empirical atomistic dipole-interaction-model for quantum plasmon simulation of nanoparticles, [_Scientific Reports_ **7**, 15775](http://doi.org/10.1038/s41598-017-16053-6)...

This will probably require marking some command line options as disabled at compile time, with some reasonable processing of this flag afterwards.

/cc @palatni @dsmunev

/cc @inzhevatkin

The next step would be to allow the rotation of the bounding box to get the minimal-volume box. There are algorithms for that, including based on voxel representation - "oriented...

When this is implemented, it makes sense to produce warning if the real grid is smaller than given in the command line, at least, when `-save_geom` is used. The reason...