From e96d0f1273608666e6022703462df724d4c28a9c Mon Sep 17 00:00:00 2001 From: Santiago Ospina De Los Rios <sospinar@gmail.com> Date: Fri, 16 Oct 2020 18:49:38 +0200 Subject: [PATCH] Do not clear vector but set it to zero --- dune/copasi/local_operator/diffusion_reaction/multidomain.hh | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/dune/copasi/local_operator/diffusion_reaction/multidomain.hh b/dune/copasi/local_operator/diffusion_reaction/multidomain.hh index 6710f77cf..577d97ff6 100644 --- a/dune/copasi/local_operator/diffusion_reaction/multidomain.hh +++ b/dune/copasi/local_operator/diffusion_reaction/multidomain.hh @@ -522,7 +522,8 @@ public: // evaluate concentrations at quadrature point - u_i.clear(), u_o.clear(); + std::fill(u_i.begin(),u_i.end(),0.); + std::fill(u_o.begin(),u_o.end(),0.); for (std::size_t comp = 0; comp < components_i; comp++) for (std::size_t dof = 0; dof < local_basis_i.size(); dof++) u_i[comp] += x_coeff_local_i(comp, dof) * phiu_i[dof]; -- GitLab