example-get-Mie.cc 2.7 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344
  1. //**********************************************************************************//
  2. // Copyright (C) 2009-2015 Ovidio Pena <ovidio@bytesfall.com> //
  3. // Copyright (C) 2013-2015 Konstantin Ladutenko <kostyfisik@gmail.com> //
  4. // //
  5. // This file is part of scattnlay //
  6. // //
  7. // This program is free software: you can redistribute it and/or modify //
  8. // it under the terms of the GNU General Public License as published by //
  9. // the Free Software Foundation, either version 3 of the License, or //
  10. // (at your option) any later version. //
  11. // //
  12. // This program is distributed in the hope that it will be useful, //
  13. // but WITHOUT ANY WARRANTY; without even the implied warranty of //
  14. // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the //
  15. // GNU General Public License for more details. //
  16. // //
  17. // The only additional remark is that we expect that all publications //
  18. // describing work using this software, or all commercial products //
  19. // using it, cite the following reference: //
  20. // [1] O. Pena and U. Pal, "Scattering of electromagnetic radiation by //
  21. // a multilayered sphere," Computer Physics Communications, //
  22. // vol. 180, Nov. 2009, pp. 2348-2354. //
  23. // //
  24. // You should have received a copy of the GNU General Public License //
  25. // along with this program. If not, see <http://www.gnu.org/licenses/>. //
  26. //**********************************************************************************//
  27. // This program returns expansion coefficents of Mie series
  28. #include "../../src/nmie-applied.h"
  29. int main(int argc, char *argv[]) {
  30. try {
  31. nmie::MultiLayerMie multi_layer_mie_;
  32. const std::complex<double> epsilon_Si(13.64, 0.047);
  33. const std::complex<double> epsilon_Ag(-28.05, 1.525);
  34. } catch( const std::invalid_argument& ia ) {
  35. // Will catch if multi_layer_mie fails or other errors.
  36. std::cerr << "Invalid argument: " << ia.what() << std::endl;
  37. return -1;
  38. }
  39. return 0;
  40. }