From a7f2bf39fbd12c162aab932a1825dcf794917493 Mon Sep 17 00:00:00 2001 From: Eivind Jahren Date: Wed, 13 Sep 2023 09:11:21 +0200 Subject: [PATCH] Rename sample_prior_perm --- notebooks/ES_2D_Heat_Equation.py | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/notebooks/ES_2D_Heat_Equation.py b/notebooks/ES_2D_Heat_Equation.py index b25ca13..082711b 100644 --- a/notebooks/ES_2D_Heat_Equation.py +++ b/notebooks/ES_2D_Heat_Equation.py @@ -69,15 +69,14 @@ # %% -def sample_prior_perm(N): +def sample_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 = sample_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 = sample_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]