From 59b95dd0964b13324b26daf14fe8dc25724c423b Mon Sep 17 00:00:00 2001 From: "Vasco C. B. Ferreira" Date: Wed, 8 Jul 2026 20:39:36 +0200 Subject: [PATCH] . --- nufi/cells.h | 3 --- src/nufi_solver.cc | 12 +++++++++++- 2 files changed, 11 insertions(+), 4 deletions(-) diff --git a/nufi/cells.h b/nufi/cells.h index 3369adf..e1f74d4 100644 --- a/nufi/cells.h +++ b/nufi/cells.h @@ -46,9 +46,6 @@ void CellLocator::rebuild(const DoFHandler &dof_handler, info.lower = cell->vertex(0); info.upper = cell->vertex(GeometryInfo::vertices_per_cell - 1); - if (info.lower[0] > info.upper[0]) - std::swap(info.lower, info.upper); - info.h = info.upper[0] - info.lower[0]; cells.push_back(info); diff --git a/src/nufi_solver.cc b/src/nufi_solver.cc index ee95a31..7171552 100644 --- a/src/nufi_solver.cc +++ b/src/nufi_solver.cc @@ -13,6 +13,7 @@ #include #include #include +#include #include #include "nufi/fields.h" @@ -151,7 +152,10 @@ NuFISolver::eval_rho(unsigned int n, std::vector &X, } void NuFISolver::run() { - std::cout << "Building E_sline\n\n"; + + //====//====// + // Run prep // + //====//====// using std::abs; using std::max; @@ -185,6 +189,10 @@ void NuFISolver::run() { const double x_min = Parameters::X_DOMAIN_LEFT; double dx = Parameters::CALC_DX; + //====//====// + // Time loop// + //====//====// + for (unsigned int it = 0; it < Nt; ++it) { stopwatch timer; @@ -239,6 +247,8 @@ void NuFISolver::run() { double refine_start = timer.elapsed(); poisson.coarse_and_refine_grid(it, phi_history); refine_time = timer.elapsed() - refine_start; + std::cout << "Refinement step done in " + << std::to_string(std::floor(refine_time)) << "[s]"; } double timer_elapsed = timer.elapsed();