nmie.h 4.0 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758
  1. //**********************************************************************************//
  2. // Copyright (C) 2009-2015 Ovidio Pena <ovidio@bytesfall.com> //
  3. // //
  4. // This file is part of scattnlay //
  5. // //
  6. // This program is free software: you can redistribute it and/or modify //
  7. // it under the terms of the GNU General Public License as published by //
  8. // the Free Software Foundation, either version 3 of the License, or //
  9. // (at your option) any later version. //
  10. // //
  11. // This program is distributed in the hope that it will be useful, //
  12. // but WITHOUT ANY WARRANTY; without even the implied warranty of //
  13. // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the //
  14. // GNU General Public License for more details. //
  15. // //
  16. // The only additional remark is that we expect that all publications //
  17. // describing work using this software, or all commercial products //
  18. // using it, cite the following reference: //
  19. // [1] O. Pena and U. Pal, "Scattering of electromagnetic radiation by //
  20. // a multilayered sphere," Computer Physics Communications, //
  21. // vol. 180, Nov. 2009, pp. 2348-2354. //
  22. // //
  23. // You should have received a copy of the GNU General Public License //
  24. // along with this program. If not, see <http://www.gnu.org/licenses/>. //
  25. //**********************************************************************************//
  26. #define VERSION "0.3.1"
  27. #include <complex>
  28. #include <vector>
  29. int ScattCoeffs(int L, int pl, std::vector<double> x, std::vector<std::complex<double> > m, int nmax,
  30. std::vector<std::complex<double> > &an, std::vector<std::complex<double> > &bn);
  31. int nMie(int L, std::vector<double> x, std::vector<std::complex<double> > m,
  32. int nTheta, std::vector<double> Theta,
  33. double *Qext, double *Qsca, double *Qabs, double *Qbk, double *Qpr, double *g, double *Albedo,
  34. std::vector<std::complex<double> > &S1, std::vector<std::complex<double> > &S2);
  35. int nMie(int L, int pl, std::vector<double> x, std::vector<std::complex<double> > m,
  36. int nTheta, std::vector<double> Theta,
  37. double *Qext, double *Qsca, double *Qabs, double *Qbk, double *Qpr, double *g, double *Albedo,
  38. std::vector<std::complex<double> > &S1, std::vector<std::complex<double> > &S2);
  39. int nMie(int L, std::vector<double> x, std::vector<std::complex<double> > m,
  40. int nTheta, std::vector<double> Theta, int nmax,
  41. double *Qext, double *Qsca, double *Qabs, double *Qbk, double *Qpr, double *g, double *Albedo,
  42. std::vector<std::complex<double> > &S1, std::vector<std::complex<double> > &S2);
  43. int nMie(int L, int pl, std::vector<double> x, std::vector<std::complex<double> > m,
  44. int nTheta, std::vector<double> Theta, int nmax,
  45. double *Qext, double *Qsca, double *Qabs, double *Qbk, double *Qpr, double *g, double *Albedo,
  46. std::vector<std::complex<double> > &S1, std::vector<std::complex<double> > &S2);
  47. int nField(int L, int pl, std::vector<double> x, std::vector<std::complex<double> > m, int nmax,
  48. int ncoord, std::vector<double> Xp, std::vector<double> Yp, std::vector<double> Zp,
  49. std::vector<std::vector<std::complex<double> > >& E, std::vector<std::vector<std::complex<double> > >& H);