test negative epsilon

go.sh

@@ -1,12 +1,16 @@
 #!/bin/bash
-echo Compile with gcc -O2
+path=$PWD
+file=compare.cc
+
+echo Compile with gcc
 rm -f *.bin
 rm -f ../scattnlay
 #g++ -Ofast -std=c++11 compare.cc nmie.cc nmie-wrapper.cc -lm -lrt -o scattnlay.bin -static
 # g++ -Ofast -std=c++11 compare.cc nmie.cc nmie-wrapper.cc -lm -lrt -o scattnlay-pg.bin -static -pg
 
 #google profiler  ######## FAST!!!
- g++ -Ofast -std=c++11 compare.cc nmie.cc nmie-old.cc -lm -lrt -o scattnlay.bin /usr/lib/libtcmalloc.so.4 -fno-builtin-malloc -fno-builtin-calloc -fno-builtin-realloc -fno-builtin-free -march=native -mtune=native -msse4.2
+echo Uncomment next line to compile compare.cc
+# g++ -Ofast -std=c++11 $file nmie.cc nmie-old.cc -lm -lrt -o scattnlay.bin /usr/lib/libtcmalloc.so.4 -fno-builtin-malloc -fno-builtin-calloc -fno-builtin-realloc -fno-builtin-free -march=native -mtune=native -msse4.2
 
 #  g++ -Ofast -std=c++11 compare.cc nmie.cc nmie-wrapper.cc -lm -lrt -o scattnlay-g.bin -ltcmalloc -fno-builtin-malloc -fno-builtin-calloc -fno-builtin-realloc -fno-builtin-free -g
 
@@ -22,6 +26,7 @@ cd tests/shell
 # 	echo $file
 #     fi
 # done
+
 PROGRAM='../../../scattnlay'
 time  ASAN_SYMBOLIZER_PATH=/usr/bin/llvm-symbolizer-3.4  $PROGRAM -l 5 0.4642 1.8000 1.7000 0.7114 0.8000 0.7000 0.7393 1.2000 0.0900 0.9168 2.8000 0.2000 1.0000 1.5000 0.4000  -t 0.0 90.0 5 -c test01
 #result
@@ -78,3 +83,12 @@ echo
 # echo ~/../mmedia/soft/cov-analysis-linux64-7.6.0/bin/cov-build --dir cov-int make standalone
 # echo      tar -zcvf cov-int.tar.gz cov-int
 # echo to submit to coverity
+
+cd $path
+file=test-negative-epsilon.cc
+g++ -Ofast -std=c++11 $file nmie.cc -lm -lrt -o $file.bin /usr/lib/libtcmalloc.so.4 -fno-builtin-malloc -fno-builtin-calloc -fno-builtin-realloc -fno-builtin-free -march=native -mtune=native -msse4.2
+
+#DEBUG!
+#clang++ -g -O1 -fsanitize=address  -fno-optimize-sibling-calls -fno-omit-frame-pointer -std=c++11 $file nmie.cc -lm -lrt -o $file.bin
+
+./$file.bin

test-negative-epsilon.cc

@@ -0,0 +1,97 @@
+/**
+ * @file   test-negative-epsilon.cc
+ * @author Konstantin Ladutenko <kostyfisik at gmail (.) com>
+ * @date   Mon Mar  9 13:21:37 2015
+ * 
+ * @brief  test negative epsilon case
+ * 
+//    This file is part of scattnlay                                                //
+//                                                                                  //
+//    This program is free software: you can redistribute it and/or modify          //
+//    it under the terms of the GNU General Public License as published by          //
+//    the Free Software Foundation, either version 3 of the License, or             //
+//    (at your option) any later version.                                           //
+//                                                                                  //
+//    This program is distributed in the hope that it will be useful,               //
+//    but WITHOUT ANY WARRANTY; without even the implied warranty of                //
+//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the                 //
+//    GNU General Public License for more details.                                  //
+//                                                                                  //
+//    The only additional remark is that we expect that all publications            //
+//    describing work using this software, or all commercial products               //
+//    using it, cite the following reference:                                       //
+//    [1] O. Pena and U. Pal, "Scattering of electromagnetic radiation by           //
+//        a multilayered sphere," Computer Physics Communications,                  //
+//        vol. 180, Nov. 2009, pp. 2348-2354.                                       //
+//                                                                                  //
+//    You should have received a copy of the GNU General Public License             //
+//    along with this program.  If not, see <http://www.gnu.org/licenses/>.         //
+ * 
+ */
+#include "nmie.h"
+#include <stdio.h>
+template<class T> inline T pow2(const T value) {return value*value;};
+const double PI=3.14159265358979323846;  
+int main(int argc, char *argv[]) {
+  try {
+    double lambda_work_ = 3.75; // cm
+    double a_ = 0.75*lambda_work_;  // 2.8125 cm - size of PEC core
+    double min_index_ = 1e-11;
+    nmie::MultiLayerMie multi_layer_mie;  
+    // Only PEC target
+    multi_layer_mie.SetTargetPEC(a_);
+    multi_layer_mie.SetWavelength(lambda_work_);
+    multi_layer_mie.RunMieCalculations();
+    double PEC_Qsca = multi_layer_mie.GetQsca();
+    double PEC_r = multi_layer_mie.GetTotalRadius();
+    double PEC_RCS = PEC_Qsca*PI*pow2(PEC_r);
+    
+    // PEC target covered with with air layer
+    multi_layer_mie.SetCoatingWidth({0.1});
+    multi_layer_mie.SetCoatingIndex({{1.0,0.0}});
+    multi_layer_mie.RunMieCalculations();
+    double Qsca1 = multi_layer_mie.GetQsca();
+    double total_r1 = multi_layer_mie.GetTotalRadius();
+    double initial_RCS1 = Qsca1*PI*pow2(total_r1);
+    printf("RCS = %g cm^2 with (r=%g) and  RCS=%g cm^2 without (r=%g)air coating.\n",
+	   initial_RCS1, total_r1, 
+	   PEC_RCS, PEC_r);
+
+    //multi_layer_mie.SetMaxTermsNumber(150);
+
+    // Bi-layer, inner layer =  lossless metall.
+    std::vector<std::complex<double>> cindex;
+    cindex.clear();
+    double n=min_index_;
+    double k=std::sqrt(0.29);
+    cindex.push_back(std::complex<double>(n, k));
+    n=std::sqrt(24.6);
+    k=min_index_;
+    cindex.push_back(std::complex<double>(n, k));
+
+    multi_layer_mie.SetCoatingWidth({0.1,0.1});
+    multi_layer_mie.SetCoatingIndex(cindex);
+    multi_layer_mie.RunMieCalculations();
+    double Qsca = multi_layer_mie.GetQsca();
+    double total_r = multi_layer_mie.GetTotalRadius();
+    double initial_RCS = Qsca*PI*pow2(total_r);
+    printf("RCS=%g for bi-layer coating (total R=%g).\n", initial_RCS,total_r);
+
+    n=1.0;
+    k=min_index_;
+    cindex.push_back(std::complex<double>(n, k));
+    multi_layer_mie.SetCoatingWidth({0.1,0.1,0.1});
+    multi_layer_mie.SetCoatingIndex(cindex);
+    multi_layer_mie.RunMieCalculations();
+    Qsca = multi_layer_mie.GetQsca();
+    total_r = multi_layer_mie.GetTotalRadius();
+    initial_RCS = Qsca*PI*pow2(total_r);
+    printf("RCS=%g for bi-layer+air coating (total R=%g).\n", initial_RCS,total_r);
+    
+  } catch( const std::invalid_argument& ia ) {
+    // Will catch if  multi_layer_mie fails or other errors.
+    std::cerr << "Invalid argument: " << ia.what() << std::endl;
+    return -1;
+  }  
+    return 0;
+}