Commit 85aa2001 authored by Steffen Schotthöfer's avatar Steffen Schotthöfer
Browse files

check if we need area cleaned values

parent 8c075531
Pipeline #101184 passed with stages
in 30 minutes and 33 seconds
#include <Python.h>
#include <mpi.h>
#include "io.h"
......@@ -6,27 +5,14 @@
#include "settings/config.h"
#include "solvers/sphericalharmonics.h"
#include <fstream>
#include <iostream>
#include <string>
// ----
int main( int argc, char** argv ) {
MPI_Init( &argc, &argv );
// wchar_t* program = Py_DecodeLocale( argv[0], NULL );
// Py_SetProgramName( program );
std::string filename = ParseArguments( argc, argv );
// CD Load Settings from File
Config* config = new Config( filename );
// Test the physics reader
// Physics testPhysic();
// testPhysic.ReadENDL( "ENDL_H.txt" );
// Print input file and run info to file
PrintLogHeader( filename );
......@@ -37,7 +23,6 @@ int main( int argc, char** argv ) {
solver->Solve();
solver->Save();
// if( Py_IsInitialized() ) Py_Finalize();
MPI_Finalize();
return EXIT_SUCCESS;
}
......@@ -135,6 +135,7 @@ void MNSolver::Solve() {
double dFlux = 1e10;
Vector fluxNew( _nCells, 0.0 );
Vector fluxOld( _nCells, 0.0 );
// Vector solArea( _nTotalEntries, 0.0 );
double mass1 = 0;
for( unsigned i = 0; i < _nCells; ++i ) {
......@@ -150,10 +151,6 @@ void MNSolver::Solve() {
if( rank == 0 ) log->info( "{:10} {:10}", "t", "dFlux" );
if( rank == 0 ) log->info( "{:03.8f} {:01.5e} {:01.5e}", -1.0, dFlux, mass1 );
// Time measurement
// auto start = chrono::steady_clock::now();
// auto end = chrono::steady_clock::now();
// Loop over energies (pseudo-time of continuous slowing down approach)
for( unsigned idx_energy = 0; idx_energy < _nEnergies; idx_energy++ ) {
......@@ -163,7 +160,9 @@ void MNSolver::Solve() {
// ------- Reconstruction Step ----------------
_optimizer->Solve( _alpha[idx_cell], _sol[idx_cell], _moments, idx_cell );
// solArea = _sol[idx_cell] * _areas[idx_cell];
_optimizer->Solve( _alpha[idx_cell], _sol[idx_cell], _moments, idx_cell ); // check if we need volume cleaned values
// ------- Relizablity Reconstruction Step ----
......
......@@ -6,10 +6,12 @@
int main( int argc, char** argv ) {
MPI_Init( &argc, &argv );
wchar_t* program = Py_DecodeLocale( argv[0], NULL );
Py_SetProgramName( program );
const int result = Catch::Session().run( argc, argv );
if( Py_IsInitialized() ) Py_Finalize();
MPI_Finalize();
return result;
}
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