Commit 6655b524 authored by jannick.wolters's avatar jannick.wolters
Browse files

changed sigma_t to VectorVector _sigmaT

parent 220f0e18
Pipeline #142984 passed with stage
in 20 minutes and 18 seconds
...@@ -59,19 +59,13 @@ CSDPNSolver::CSDPNSolver( Config* settings ) : PNSolver( settings ) { ...@@ -59,19 +59,13 @@ CSDPNSolver::CSDPNSolver( Config* settings ) : PNSolver( settings ) {
// std::cout << sigma_ref.rows() << std::endl; // std::cout << sigma_ref.rows() << std::endl;
// std::cout << _energies.size() << std::endl; // std::cout << _energies.size() << std::endl;
Matrix sigma_t( _energies.size(), sigma_ref.rows() ); _sigmaT = VectorVector( _polyDegreeBasis, Vector( _energies.size() ) );
for( unsigned idx_degree = 0; idx_degree < _polyDegreeBasis; ++idx_degree ) { for( unsigned idx_degree = 0; idx_degree < _polyDegreeBasis; ++idx_degree ) {
Vector xs_m = blaze::column( sigma_ref, idx_degree ); // Scattering cross section Moments Vector xs_m = blaze::column( sigma_ref, idx_degree ); // Scattering cross section Moments
Interpolation interp( E_ref, xs_m ); Interpolation interp( E_ref, xs_m );
blaze::column( sigma_t, idx_degree ) = interp( _energies ); _sigmaT[idx_degree] = interp( _energies );
} }
// std::cout << "here\n";
// std::cout << size( sigma_t ) << std::endl;
_sigmaT = sigma_t;
TextProcessingToolbox::PrintMatrix( sigma_t );
Interpolation interpS( E_tab, S_tab ); Interpolation interpS( E_tab, S_tab );
_s = interpS( _energies ); _s = interpS( _energies );
} }
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment