branch: ignore_empty_integration_points
commit d5cebf00a9e3db80c66ab11b2249bc9ec11a859d
Author: lj <l...@plaxis.nl>
Date: Thu May 17 17:40:48 2018 +0200
Use consistent type.
---
src/getfem_integration.cc | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/getfem_integration.cc b/src/getfem_integration.cc
index 6e4b8f6..32a4826 100644
--- a/src/getfem_integration.cc
+++ b/src/getfem_integration.cc
@@ -253,7 +253,7 @@ namespace getfem {
GMM_ASSERT1(!valid, "Impossible to modify a valid integration method.");
if (gmm::abs(w) > 1.0E-15 || include_empty) {
++f;
- if (gmm::abs(w) <= 1.0E-15) w = 0.0;
+ if (gmm::abs(w) <= 1.0E-15) w = scalar_type(0);
GMM_ASSERT1(f <= cvr->structure()->nb_faces(), "Wrong argument.");
size_type i = pt_to_store[f].search_node(pt);
if (i == size_type(-1)) {
