VTK output in kernel calulation as an option

examples/KernelV0-2.cpp

@@ -24,7 +24,7 @@ using namespace Lemma;
 int main(int argc, char** argv) {
 
     if (argc<3) {
-        std::cout << "./KernelV0-2 Kernel.yaml TxString RxString  \n";
+        std::cout << "./KernelV0-2 Kernel.yaml TxString RxString  vtkoutput<true/false> \n";
         return(EXIT_SUCCESS);
     }
 
@@ -33,7 +33,14 @@ int main(int argc, char** argv) {
 
     std::vector<std::string> tx = {std::string(argv[2])};
     std::vector<std::string> rx = {std::string(argv[3])};
-    Kern->CalculateK0( tx, rx, false ); // 3rd argument is vtk output
+
+    std::cout << "argv[4]\t" << argv[4] << std::endl;
+    if( std::string(argv[4]) == "true") {
+        Kern->CalculateK0( tx, rx, true ); // 3rd argument is vtk output
+    } else {
+        std::cout << "not using VTK" << std::endl;
+        Kern->CalculateK0( tx, rx, false ); // 3rd argument is vtk output
+    }
 
     // TODO fix python post-processing so this is not necessary
     // Save in simplified format for easy python plotting

src/KernelV0.cpp

@@ -373,6 +373,7 @@ namespace Lemma {
         // Compute Mn0
         Vector3r Mn0 = ComputeMn0(1.0, B0);
         Real Mn0Abs = Mn0.norm();
+        //std::cout << "Mn0\t" << Mn0.transpose() << std::endl;
 
         // Compute phase delay
         // TODO add transmiiter current phase and delay induced apparent time phase!