vcbferreira/NuFI_deal.ii
1
0
Fork 0
mirror of https://codeberg.org/vcbferreira/NuFI_deal.ii synced 2026-08-12 14:33:18 +02:00
Code Issues Packages Projects Releases Wiki Activity

bug corection, removal of bloat

Browse Source
This commit is contained in:
Vasco C. B. Ferreira
2026-07-27 20:28:33 +02:00
parent 73fa2fed30
commit fbe57b2a12
2 changed files with 6 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files
nufi/poisson_problem.h
-7
View File
@@ -434,13 +434,6 @@ template <int dim> void PoissonProblem<dim>::coarse_and_refine_grid(size_t it) {
std::string grid_file_name = std::string grid_file_name =
Parameters::PLOT_DIR + "grid_" + std::to_string(it); Parameters::PLOT_DIR + "grid_" + std::to_string(it);
save_grid_to_file(grid_file_name); save_grid_to_file(grid_file_name);
// std::vector<double> Ex =
// sample_electric_potential(Parameters::X_DOMAIN_LEFT,
// Parameters::X_DOMAIN_RIGHT,
// Parameters::PLOT_NX);
//
// save_space_vector(Ex, "Ex_after_coarsed", it);
} }
template <int dim> void PoissonProblem<dim>::estimate_error() { template <int dim> void PoissonProblem<dim>::estimate_error() {
src/nufi_solver.cc
+6 -4
View File
@@ -193,7 +193,7 @@ void NuFISolver::run() {
std::vector<double> int_E_squared; std::vector<double> int_E_squared;
int_E_squared.reserve(Nt); int_E_squared.reserve(Nt);
std::vector<double> int_E_squared_times; // <-- add this std::vector<double> int_E_squared_times;
int_E_squared_times.reserve(Nt); int_E_squared_times.reserve(Nt);
std::vector<GridStructure<1>> grid_versions; std::vector<GridStructure<1>> grid_versions;
@@ -215,7 +215,8 @@ void NuFISolver::run() {
<< "\n"; << "\n";
std::ofstream error_file(Parameters::PLOT_DIR + "error_estimate.dat"); std::ofstream error_file(Parameters::PLOT_DIR + "error_estimate.dat");
error_file << "it error_estimate\n"; error_file << "# nufi Kelly l2 error estimate\n";
error_file << "# it l2_error_estimate\n";
[[maybe_unused]] const double x_min = Parameters::X_DOMAIN_LEFT; [[maybe_unused]] const double x_min = Parameters::X_DOMAIN_LEFT;
[[maybe_unused]] double dx = Parameters::CALC_DX; [[maybe_unused]] double dx = Parameters::CALC_DX;
@@ -237,8 +238,9 @@ void NuFISolver::run() {
<< std::endl; << std::endl;
// START: diagnostics // START: diagnostics
// std::cout << "cells = " << poisson.triangulation.n_active_cells() << "\n" std::cout << "cells = " << poisson.get_triangulation().n_active_cells()
// << " dofs = " << poisson.dof_handler.n_dofs() << "\n"; << "\n"
<< " dofs = " << poisson.get_dof_handler().n_dofs() << "\n";
// double min_h = 1e100; // double min_h = 1e100;
// double max_h = 0; // double max_h = 0;
// //