Added storage of full 2x2 form factor calculations to IInterferenceFunctionStrategy

Core/Algorithms/inc/IInterferenceFunctionStrategy.h

@@ -105,6 +105,9 @@ private:
 
     //! cached form factor evaluations
     mutable std::vector<complex_t> m_ff00, m_ff01, m_ff10, m_ff11;
+
+    //! cached polarized form factors
+    mutable std::vector<Eigen::Matrix2cd> m_ff_pol;
 };
 
 inline cvector_t IInterferenceFunctionStrategy::getQ(const cvector_t& k_i,

Core/Algorithms/src/IInterferenceFunctionStrategy.cpp

@@ -104,6 +104,7 @@ void IInterferenceFunctionStrategy::calculateFormFactorLists(
         (*it)->mp_ff->setSpecularInfo(p_in_coeffs, P_out_coeffs.get());
         Eigen::Matrix2cd ff_mat = (*it)->mp_ff->evaluatePol(k_i, k_f_bin,
                 alpha_f_bin, phi_f_bin);
+        m_ff_pol.push_back(ff_mat);
         m_ff00.push_back((complex_t)ff_mat(0,0));
         m_ff01.push_back((complex_t)ff_mat(0,1));
         m_ff10.push_back((complex_t)ff_mat(1,0));