SimPEG.maps.ParametricBlock#
- class SimPEG.maps.ParametricBlock(mesh, epsilon=1e-06, p=10, **kwargs)[source]#
Bases:
SimPEG.maps.BaseParametricMapping for a rectangular block within a wholespace.
This mapping is used when the cells lying below the Earth’s surface can be parameterized by rectangular block within a homogeneous medium. The model is defined by the physical property value for the background (\(\sigma_0\)), the physical property value for the block (\(\sigma_b\)), parameters for the center of the block (\(x_b [,y_b, z_b]\)) and parameters for the dimensions along each Cartesian direction (\(dx [,dy, dz]\))
For this mapping, the set of input model parameters are organized:
\[\begin{split}\mathbf{m} = \begin{cases} 1D: \;\; [\sigma_0, \;\sigma_b,\; x_b , \; dx] \\ 2D: \;\; [\sigma_0, \;\sigma_b,\; x_b , \; dx,\; y_b , \; dy] \\ 3D: \;\; [\sigma_0, \;\sigma_b,\; x_b , \; dx,\; y_b , \; dy,\; z_b , \; dz] \end{cases}\end{split}\]The mapping \(\mathbf{u}(\mathbf{m})\) from the model to the mesh is given by:
\[\mathbf{u}(\mathbf{m}) = \sigma_0 + (\sigma_b - \sigma_0) \bigg [ \frac{1}{2} + \pi^{-1} \arctan \bigg ( a \, \boldsymbol{\eta} \big ( x_b, y_b, z_b, dx, dy, dz \big ) \bigg ) \bigg ]\]where a is a parameter that impacts the sharpness of the arctan function, and
\[\boldsymbol{\eta} \big ( x_b, y_b, z_b, dx, dy, dz \big ) = 1 - \sum_{\xi \in (x,y,z)} \bigg [ \bigg ( \frac{2(\boldsymbol{\xi_c} - \xi_b)}{d\xi} \bigg )^2 + \varepsilon^2 \bigg ]^{p/2}\]Parameters \(p\) and \(\varepsilon\) define the parameters of the Ekblom function. \(\boldsymbol{\xi_c}\) is a place holder for vectors containing the x, [y and z] cell center locations of the mesh, \(\xi_b\) is a placeholder for the x[, y and z] location for the center of the block, and \(d\xi\) is a placeholder for the x[, y and z] dimensions of the block.
- Parameters
- mesh
discretize.BaseMesh A discretize mesh
- indActive
numpy.ndarray Active cells array. Can be a boolean
numpy.ndarrayof length mesh.nC or anumpy.ndarrayofintcontaining the indices of the active cells.- slope
float Directly define the constant a in the mapping function which defines the sharpness of the boundaries.
- slopeFact
float Scaling factor for the sharpness of the boundaries based on cell size. Using this option, we set a = slopeFact / dh.
- epsilon
float Epsilon value used in the ekblom representation of the block
- p
float p-value used in the ekblom representation of the block.
- mesh
Examples
Attributes
epsilon value used in the ekblom representation of the block.
Number of parameters the mapping acts on.
p-value used in the ekblom representation of the block.
Dimensions of the mapping
Methods
deriv(m)Derivative of the mapping with respect to the input parameters.
mDict(m)Return model parameters as a dictionary.
