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

eval loop bug corrected, poisson solver to fix

Browse Source
This commit is contained in:
Vasco C. B. Ferreira
2026-03-14 16:52:14 +01:00
parent f2ee1606b7
commit 9df98147b9
2 changed files with 536 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
nufi_solver.hpp
+24 -2
View File
@@ -184,13 +184,35 @@ inline void NuFISolver::run()
std::cout << "Timestep " << it << " / " << Nt << std::endl << std::endl;
// compute rho
for(size_t i = 0; i<Nx; i++)
double dx = Parameters::SPLINE_DX;
double x = Parameters::X_DOMAIN_LEFT;
for(size_t i = 0; i<Nx; i++, x+=dx)
{
double ith_rho = eval_rho(it, i, coeffs.get(), Parameters::NV);
double ith_rho = eval_rho(it, x, coeffs.get(), Parameters::NV);
AssertThrow(std::isfinite(ith_rho), ExcMessage("NaN detected in rho"));
rho.get()[i] = ith_rho;
}
// bool rho_written = false;
//
// if (!rho_written)
// {
// std::ofstream rho_file("rho_initial.dat");
//
// if (!rho_file)
// throw std::runtime_error("Could not open rho_initial.dat");
//
// rho_file.precision(16);
// rho_file << std::scientific;
//
// for (size_t i = 0; i < Nx; ++i)
// rho_file << rho.get()[i] << "\n";
//
// rho_written = true;
// }
poisson.set_rhs_function(std::make_unique<ChargeDensity_NuFI<1>>(rho.get(), Parameters::SPLINE_NX));
poisson.solve_step();