From 1145d8a3eb335bc72beff3a810ec7aabc9fd9263 Mon Sep 17 00:00:00 2001 From: Yi Zhang Date: Thu, 21 Nov 2024 16:39:14 +0800 Subject: [PATCH] tmp --- lib/mesh/regular_grid.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/lib/mesh/regular_grid.cpp b/lib/mesh/regular_grid.cpp index b30462c..28bc240 100644 --- a/lib/mesh/regular_grid.cpp +++ b/lib/mesh/regular_grid.cpp @@ -1082,7 +1082,7 @@ void gctl::regular_grid::boolean(std::string newname, std::string datname, std:: { for (size_t i = 0; i < val_ptr->size(); i++) { - if (fabs(val_ptr2->at(i)) < 1e-8) val_ptr->at(i) = GCTL_BDL_MAX; + if (fabs(val_ptr2->at(i)) < 1e-8 || isnan(val_ptr2->at(i))) val_ptr->at(i) = GCTL_BDL_MAX; else val_ptr->at(i) = val_ptr1->at(i); } } @@ -1090,7 +1090,7 @@ void gctl::regular_grid::boolean(std::string newname, std::string datname, std:: { for (size_t i = 0; i < val_ptr->size(); i++) { - if (fabs(val_ptr2->at(i)) < 1e-8) val_ptr->at(i) = val_ptr1->at(i); + if (fabs(val_ptr2->at(i)) < 1e-8 || isnan(val_ptr2->at(i))) val_ptr->at(i) = val_ptr1->at(i); else val_ptr->at(i) = GCTL_BDL_MAX; } }