Commit 8ae90316 authored by Steffen Schotthöfer's avatar Steffen Schotthöfer
Browse files

fixed checkerboard moment method source term prefactor

parent 5caf02f6
......@@ -95,7 +95,7 @@ std::vector<VectorVector> Checkerboard_PN::GetExternalSource( const Vector& ener
VectorVector Q( _mesh->GetNumCells(), Vector( 1u, 0.0 ) );
auto cellMids = _mesh->GetCellMidPoints();
for( unsigned j = 0; j < cellMids.size(); ++j ) {
if( isSource( cellMids[j] ) ) Q[j] = std::sqrt( 4 * M_PI ); // isotropic source
if( isSource( cellMids[j] ) ) Q[j] = 1.0 / std::sqrt( 4 * M_PI ); // isotropic source
}
return std::vector<VectorVector>( 1u, Q );
}
......
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