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