Rename sample_prior_perm

notebooks/ES_2D_Heat_Equation.py

@@ -69,15 +69,14 @@
 
 
 # %%
-def sample_prior_perm(N):
+def prior_conductivity(N):
     mesh = np.meshgrid(np.linspace(0, 1, nx), np.linspace(0, 1, nx))
-    lperms = np.exp(geostat.gaussian_fields(mesh, rng, N, r=0.8))
-    return lperms
+    return np.exp(geostat.gaussian_fields(mesh, rng, N, r=0.8))
 
 
 # Evensens' formulation of the Ensemble Smoother has the prior as
 # a (nx * nx, N) matrix, i.e (number of parameters, N).
-A = sample_prior_perm(N).T
+A = prior_conductivity(N).T
 
 # We'll also need a list of matrices to run simulations in parallel later on.
 # A list is also a bit easier to interactively visualize.
@@ -109,7 +108,7 @@ def interactive_prior_fields(n):
 
 # %%
 # Set the coefficient of heat transfer for each grid cell.
-# alpha_t = sample_prior_perm(1).T.reshape(nx, nx)
+# alpha_t = prior_conductivity(1).T.reshape(nx, nx)
 # Let's use as true parameter field one relization from the prior
 # to make it easier for the Ensemble Smoother to find the solution.
 alpha_t = alphas[0]