Commit 67ec4b07 authored by jonas.kusch's avatar jonas.kusch
Browse files

no more debug output

parent 1cf542b0
Pipeline #117777 failed with stages
in 26 minutes and 48 seconds
......@@ -557,11 +557,8 @@ void ICRU::GetTransportCoefficients( Matrix& xi ) {
Vector mu( _XMU.size() );
std::vector<std::vector<double>> dxs;
// get scattering cross sections for all energies E
std::cout << "Computing dxs" << std::endl;
angdcs( 3u, dxs );
std::cout << "DONE." << std::endl;
for( unsigned n = 0; n < mu.size(); ++n ) mu[n] = 1.0 - 2.0 * _XMU[n]; // mu starts at 1 and goes to 0
std::cout << "Reset mu DONE." << std::endl;
#pragma omp parallel for
for( unsigned i = 0; i < _E.size(); ++i ) {
// compute moments with trapezoidal rule
......@@ -573,7 +570,6 @@ void ICRU::GetTransportCoefficients( Matrix& xi ) {
}
}
}
std::cout << "Computation xi DONE." << std::endl;
xi *= 2.0 * PI * H2OMolecularDensity;
}
......
......@@ -77,12 +77,9 @@ CSDSolverTrafoFP::CSDSolverTrafoFP( Config* settings ) : SNSolver( settings ) {
double exponent = minExp + ( maxExp - minExp ) / ( _nEnergies - 1 ) * n;
_energies[n] = pow(10.0,exponent);
}*/
std::cout << "Setting up transport coefficients..." << std::endl;
ICRU database( mu, _energies );
database.GetTransportCoefficients( _xi );
std::cout << "Setting up transport coefficients DONE." << std::endl;
database.GetStoppingPower( _s );
std::cout << "Setting up stopping powers DONE." << std::endl;
/*
// print coefficients
std::cout<<"E = [";
......
Supports Markdown
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