max/TugJulia
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

feat: print sums of input and output field

Browse Source
This commit is contained in:
Max Lübke 2023-09-07 11:10:56 +02:00
parent 42ad07c252
commit 5490216bc0
1 changed files with 10 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
naaice/BTCS_2D_NAAICE.cpp
View File

@ -112,6 +112,8 @@ int main(int argc, char *argv[]) {
MatrixXd concentrations = rmVecTocmMatrix(init_values_vec, row, col); MatrixXd concentrations = rmVecTocmMatrix(init_values_vec, row, col);
grid.setConcentrations(concentrations); grid.setConcentrations(concentrations);
const double sum_init = concentrations.sum();
// // (optional) set alphas of the grid, e.g.: // // (optional) set alphas of the grid, e.g.:
const auto alphax_vec = CSVToVector<double>(INPUT_ALPHAX_FILE); const auto alphax_vec = CSVToVector<double>(INPUT_ALPHAX_FILE);
MatrixXd alphax = rmVecTocmMatrix(alphax_vec, row, col); MatrixXd alphax = rmVecTocmMatrix(alphax_vec, row, col);
@ -154,4 +156,12 @@ int main(int argc, char *argv[]) {
// // **** RUN SIMULATION **** // // **** RUN SIMULATION ****
simulation.run(); simulation.run();
const double sum_after = grid.getConcentrations().sum();
std::cout << "Sum of init field: " << std::to_string(sum_init)
<< "\nSum after iteration: " << std::to_string(sum_after)
<< std::endl;
return 0;
} }