diff --git a/src/lbm/SuViscoelasticity.h b/src/lbm/SuViscoelasticity.h
index 80f0abf15be03c1f401be351b1106e925b581c45..d7a2644467f999598a69536a58718e067d81b41c 100644
--- a/src/lbm/SuViscoelasticity.h
+++ b/src/lbm/SuViscoelasticity.h
@@ -223,7 +223,7 @@ public:
                   if (boundaryHandling->isDomain(cell + *d)) {
                      f[i] += stress->get(cell + *d)(k, i) * d.c(k);
                   } else if(boundaryHandling->isDomain(cell - *d)){
-                     f[i] += -stress->get(cell - *d)(k, i) * d.c(k) + stress->get(cell)(k, i) * d.c(k) * real_c(2);
+                     f[i] += -stress->get(cell - *d)(k, i) * d.c(k) + stress->get(cell)(k, i) * d.c(k) * real_c(2.0);
                   } else {
                      f[i] += real_t(0);
                   }