Mercurial > gemma

diff pkg/mesh/raster.go @ 5703:d2ccf6bb6940 sr-v2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make plane eval for z-values of triangles numerial more robust.
author Sascha L. Teichmann <sascha.teichmann@intevation.de>
date Mon, 19 Feb 2024 17:48:13 +0100
parents 1222b777f51f
children
line wrap: on
line diff
--- a/pkg/mesh/raster.go	Wed Feb 14 22:38:14 2024 +0100
+++ b/pkg/mesh/raster.go	Mon Feb 19 17:48:13 2024 +0100
@@ -104,8 +104,10 @@
 					n++
 				}
 
-				if n > 0 {
+				if n > 0 && !math.IsNaN(sum) {
 					row[j] = sum / float64(n)
+				} else {
+					row[j] = noData
 				}
 				px += r.CellSize
 			}
@@ -168,7 +170,7 @@
 					n++
 				}
 
-				if n > 0 {
+				if n > 0 && !math.IsNaN(sum) {
 					row[j] -= sum / float64(n)
 				} else {
 					row[j] = noData