Retain mean mode from white noise

Ceyron · Oct 22, 2024 · 6122236 · 6122236
1 parent 47107ee
commit 6122236
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/exponax/ic/_gaussian_random_field.py b/exponax/ic/_gaussian_random_field.py
@@ -66,7 +66,7 @@ def __call__(
         wavenumer_norm_grid = jnp.linalg.norm(wavenumber_grid, axis=0, keepdims=True)
         amplitude = jnp.power(wavenumer_norm_grid, -self.powerlaw_exponent / 2.0)
         amplitude = (
-            amplitude.flatten().at[0].set(0.0).reshape(wavenumer_norm_grid.shape)
+            amplitude.flatten().at[0].set(1.0).reshape(wavenumer_norm_grid.shape)
         )
 
         real_key, imag_key = jr.split(key, 2)