NYs = self.iYobs.size
NY = Y.size
cols = self.iYobs
rows = np.arange(NYs)
vals = np.full(NYs, 1.0)
Hy = sps.coo_matrix((vals, (rows, cols)), shape=(NYs, NY))
return sps.vstack([-Hy, np.sqrt(self.gamma) * self.L])
if __name__ == '__main__':
L = np.array([1.0, 1.0])
N = np.array([32, 32])
g = Geom(L, N)
g.calculate()
ul = 2.0
ur = 1.0
bc = BC(g)
bc.dirichlet(g, "left", ul)
bc.dirichlet(g, "right", ur)
# Hydraulic conductivity
npr.seed(0)
se = SEKernel(std_dev=1.0, cor_len=0.15, std_dev_noise=0.0)
CY = se.covar(g.cells.centroids.T, g.cells.centroids.T)
Nc = np.prod(N)
Y = se.sample(CY, npr.randn(Nc))
K = np.exp(Y)
