rnatm

commit 18e67850a59852942ec3a9834a15d8612bb3bb32
parent f1adca8ea68930072871f4674d4eb5319b89e031
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Wed, 27 Jul 2022 11:29:20 +0200

Updating voxel values

Empty voxels store degenerate ranges rather than 0 to prevent empty
space from forcing octree refinement

Diffstat:
M
src/rnatm_octree.c
 | 
11
+----------

1 file changed, 1 insertion(+), 10 deletions(-)
diff --git a/src/rnatm_octree.c b/src/rnatm_octree.c
@@ -486,7 +486,6 @@ vx_get(const size_t xyz[3], const uint64_t mcode, void* dst, void* context)
   struct build_octrees_context* ctx = context;
   float* vx = NULL;
   uint64_t ivx, ipart;
-  float kext_min, kext_max;
   int log2_part_def;
   ASSERT(xyz && dst && ctx);
   (void)xyz;
@@ -513,15 +512,7 @@ vx_get(const size_t xyz[3], const uint64_t mcode, void* dst, void* context)
   }
 
   vx = partition_get_voxel(ctx->part, 0, 0, ivx);
-  kext_min = vx[voxel_idata(RNATM_RADCOEF_Kext, RNATM_SVX_OP_MIN)];
-  kext_max = vx[voxel_idata(RNATM_RADCOEF_Kext, RNATM_SVX_OP_MAX)];
-
-  if(kext_min > kext_max) {
-    /* The voxel is not covered by any tetrahedron. Its coefficients are zero */
-    memset(dst, 0, NFLOATS_PER_VOXEL * sizeof(float));
-  } else {
-    memcpy(dst, vx, NFLOATS_PER_VOXEL * sizeof(float));
-  }
+  memcpy(dst, vx, NFLOATS_PER_VOXEL * sizeof(float));
 }
 
 static void