Skip to content
Snippets Groups Projects
Commit f4bccbcc authored by Stéphane Del Pino's avatar Stéphane Del Pino
Browse files

Pass list of boundary conditions as reference (suggested by clazy)

parent 775a4038
No related branches found
No related tags found
1 merge request!37Feature/language
...@@ -26,7 +26,8 @@ struct GlaceScheme ...@@ -26,7 +26,8 @@ struct GlaceScheme
const MeshType& m_mesh; const MeshType& m_mesh;
GlaceScheme(const IMesh& mesh, std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>> bc_descriptor_list) GlaceScheme(const IMesh& mesh,
const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_descriptor_list)
: m_mesh{dynamic_cast<const MeshType&>(mesh)} : m_mesh{dynamic_cast<const MeshType&>(mesh)}
{ {
MeshDataType mesh_data(m_mesh); MeshDataType mesh_data(m_mesh);
...@@ -103,6 +104,12 @@ struct GlaceScheme ...@@ -103,6 +104,12 @@ struct GlaceScheme
if (t + dt > tmax) { if (t + dt > tmax) {
dt = tmax - t; dt = tmax - t;
} }
std::cout.setf(std::cout.scientific);
std::cout << "iteration " << rang::fg::cyan << std::setw(4) << iteration << rang::style::reset
<< " time=" << rang::fg::green << t << rang::style::reset << " dt=" << rang::fgB::blue << dt
<< rang::style::reset << '\n';
acoustic_solver.computeNextStep(t, dt, unknowns); acoustic_solver.computeNextStep(t, dt, unknowns);
block_eos.updatePandCFromRhoE(); block_eos.updatePandCFromRhoE();
...@@ -110,6 +117,9 @@ struct GlaceScheme ...@@ -110,6 +117,9 @@ struct GlaceScheme
t += dt; t += dt;
++iteration; ++iteration;
} }
std::cout << rang::style::bold << "Final time=" << rang::fgB::green << t << rang::style::reset << " reached after "
<< rang::fgB::cyan << iteration << rang::style::reset << rang::style::bold << " iterations"
<< rang::style::reset << '\n';
vtk_writer.write(m_mesh, vtk_writer.write(m_mesh,
{NamedItemValue{"density", rhoj}, NamedItemValue{"velocity", unknowns.uj()}, {NamedItemValue{"density", rhoj}, NamedItemValue{"velocity", unknowns.uj()},
NamedItemValue{"coords", m_mesh.xr()}, NamedItemValue{"coords", m_mesh.xr()},
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment