fix: rbf interpolation functions should expect squared distance

ozwaldorf · Jul 31, 2023 · ce3b0f8 · ce3b0f8
1 parent 7d70d99
commit ce3b0f8
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/interpolation/rbf/gaussian.rs b/src/interpolation/rbf/gaussian.rs
@@ -6,7 +6,7 @@ pub struct GaussianFn {
 }
 impl RadialBasisFn for GaussianFn {
     fn radial_basis(&self, distance: f64) -> f64 {
-        (-self.shape * distance.powf(2.0)).exp()
+        (-self.shape * distance).exp()
     }
 }
 

diff --git a/src/interpolation/rbf/shepard.rs b/src/interpolation/rbf/shepard.rs
@@ -6,7 +6,7 @@ pub struct InverseDistanceFn {
 }
 impl RadialBasisFn for InverseDistanceFn {
     fn radial_basis(&self, distance: f64) -> f64 {
-        1.0 / distance.powf(self.power)
+        1.0 / distance.sqrt().powf(self.power)
     }
 }