Nav apraksta

Konstantin Ladutenko 9f9c1d38f9 Makefile changes 10 gadi atpakaļ
debian b911a457c8 Deleted some garbage files generated during compilation. 10 gadi atpakaļ
doc c5f9db4d87 Small changes and new example image 10 gadi atpakaļ
tests fca8e2c3fc removed dependance from bessel.cc 10 gadi atpakaļ
.gitignore c78969a511 Re-added scattnlay.cpp, it is needed to compile Python extension without cython installed. 10 gadi atpakaļ
CHANGES 1b7786ad5a Initial commit based on Ovidio version of python-scattnlay 0.3.0 10 gadi atpakaļ
COPYING 1b7786ad5a Initial commit based on Ovidio version of python-scattnlay 0.3.0 10 gadi atpakaļ
LICENSE 1b7786ad5a Initial commit based on Ovidio version of python-scattnlay 0.3.0 10 gadi atpakaļ
MANIFEST.in 8d6d84aa20 Calculation of electric field outside the particle is mostly working. Added a python test for this calculation. 10 gadi atpakaļ
Makefile 9f9c1d38f9 Makefile changes 10 gadi atpakaļ
PKG-INFO a02a844cdc Modifications to complete the porting to C++. Most of the work was done by 'kostyfisik'. 10 gadi atpakaļ
README.md 2700b4a599 small changes 10 gadi atpakaļ
compare.cc bf1cc1ce47 Small changes 10 gadi atpakaļ
go.sh fca8e2c3fc removed dependance from bessel.cc 10 gadi atpakaļ
nmie-old.cc 3a6320a005 Revised calculation of electric field. Everything seems right now for the calculation outside the particle. Also did several small format changes. 10 gadi atpakaļ
nmie-old.h a098084807 Core calculations ported to C++ class. 10 gadi atpakaļ
nmie-wrapper.cc 553a8d05d1 Reversed changes to scattnlay.pyx because the python extension would compile well but it was not returning the complex vectors (S1 and S2). If changed again we must verify that tes04.py works!!! 10 gadi atpakaļ
nmie-wrapper.h fd8d92ade7 Separating the calculation core in a new class. It is not yet ready. 10 gadi atpakaļ
nmie.cc fca8e2c3fc removed dependance from bessel.cc 10 gadi atpakaļ
nmie.h 9b6ce87c48 Cleaned code and removed unused functions 10 gadi atpakaļ
push-to-github.sh eea51ce5ca Changes to push script 10 gadi atpakaļ
py_nmie.cc 8dafd465ef Python wrapper working with new class. 10 gadi atpakaļ
py_nmie.h 8dafd465ef Python wrapper working with new class. 10 gadi atpakaļ
scattnlay.cpp 03967872ac Fixed compilation of python extension. Calculation of an_bulk_, bn_bulk_ seems broken. Field calculation using these values gives an error. 10 gadi atpakaļ
scattnlay.pyx 8dafd465ef Python wrapper working with new class. 10 gadi atpakaļ
setup.py 9f9c1d38f9 Makefile changes 10 gadi atpakaļ
setup_cython.py 9f9c1d38f9 Makefile changes 10 gadi atpakaļ
standalone.cc 5ce07e9840 removed warnings about unsigned int comparison 10 gadi atpakaļ
test-negative-epsilon.cc 097982a8b2 seem to work ok with optimizer 10 gadi atpakaļ

README.md

output example Output example: Field distribution inside layered Si\Ag\Si sphere and Poynting vector distribution in Ag sphere with poweflow lines calculated with Scattnlay.

How to use scattnlay

Table of contents:

Compile Code:

Compilation options

  • make source - Create source package (python library)
  • make install - Install on local system (python library)
  • make buildrpm - Generate a rpm package (python library)
  • make builddeb - Generate a deb package (python library)
  • make standalone - Create a standalone program
  • make clean - Delete temporal files

Use:

  1. Python library
    • Use scattnlay directly
from scattnlay import scattnlay
...
x = ...
m = ...
terms, Qext, Qsca, Qabs, Qbk, Qpr, g, Albedo, S1, S2 = scattnlay(x, m)
...
  • Execute some of the test scripts (located in the folder 'tests/python') Example:
./test01.py
  1. Standalone program
    • Execute scattnlay directly Usage:
scattnlay -l Layers x1 m1.r m1.i [x2 m2.r m2.i ...] [-c comment]
  • Execute some of the test scripts (located in the folder 'tests/shell') Example:
./test01.sh > test01.csv
  1. C++ library
    try {
      MultiLayerMie multi_layer_mie;
      multi_layer_mie.SetLayersSize(x);
      multi_layer_mie.SetLayersIndex(m);

      multi_layer_mie.RunMieCalculation();

      *Qsca = multi_layer_mie.GetQsca();
      *Qabs = multi_layer_mie.GetQabs();
    } catch(const std::invalid_argument& ia) {
      // Will catch if  multi_layer_mie fails or other errors.
      std::cerr << "Invalid argument: " << ia.what() << std::endl;
      throw std::invalid_argument(ia);
      return -1;
    }

Papers

  1. "Scattering of electromagnetic radiation by a multilayered sphere" O. Pena and U. Pal, Computer Physics Communications, vol. 180, Nov. 2009, pp. 2348-2354. http://dx.doi.org/10.1016/j.cpc.2009.07.010

  2. "Reduction of scattering using thin all-dielectric shells designed by stochastic optimizer" Konstantin Ladutenko, Ovidio Peña-Rodríguez, Irina Melchakova, Ilya Yagupov, and Pavel Belov J. Appl. Phys., vol. 116, pp. 184508, 2014 http://dx.doi.org/10.1063/1.4900529

Acknowledgment

We expect that all publications describing work using this software, or all commercial products using it, cite the following reference:

O. Pena and U. Pal, "Scattering of electromagnetic radiation by a multilayered sphere," Computer Physics Communications, vol. 180, Nov. 2009, pp. 2348-2354.

License

GPL v3+