diff --git a/src/lbm/lattice_model/SmagorinskyLES.h b/src/lbm/lattice_model/SmagorinskyLES.h
index 215b9652e1fc189eca63df0aefacfabbf0aad925..169439e77c9d4341c86654f70756ee193797ec16 100644
--- a/src/lbm/lattice_model/SmagorinskyLES.h
+++ b/src/lbm/lattice_model/SmagorinskyLES.h
@@ -144,7 +144,7 @@ void SmagorinskyLES< LatticeModel_T, Filter_T >::operator()( IBlock * block, con
             {
                real_t qij = real_t(0);
                for( auto d = Stencil_T::begin(); d != Stencil_T::end(); ++d )
-                  qij += nonEquilibrium[ d.toIdx() ] * stencil::c[alpha][*d] * stencil::c[beta][*d];
+                  qij += nonEquilibrium[ d.toIdx() ] * real_c(stencil::c[alpha][*d]) * real_c(stencil::c[beta][*d]);
                filteredMeanMomentum += qij * qij;
             }
          }