Add _xi1&_xi2 in .h file

code/include/solvers/csdsolvertrafofp.h

@@ -26,7 +26,9 @@ class CSDSolverTrafoFP : public SNSolver
     double _beta; /*!  @brief  Coefficient of GFP operators (see Olbrant 2010, Appendix B)*/
 
     Matrix _xi;  /*!  @brief matrix of transport coefficients */
-
+    Vector _xi1;
+    Vector _xi2;
+    
     unsigned _FPMethod; /*!  @brief Encodes different ways of computing coefficients alpha, alpha2 & beta, _FPMethod == 1, 2 ,3 stand for methods with increasing accuracy (see Olbrant 2010, Appendix B)*/
 
     bool _RT; /*!  @brief radiotherapy application (on/off), if true use crosssections + stopping powers from database  */

code/src/solvers/snsolver.cpp

@@ -9,6 +9,7 @@
 // externals
 #include "spdlog/spdlog.h"
 #include <mpi.h>
+#include <iostream>
 
 SNSolver::SNSolver( Config* settings ) : Solver( settings ) {