Mesh: Tensor: CreationΒΆ

For tensor meshes, there are some functions that can come in handy. For example, creating mesh tensors can be a bit time consuming, these can be created speedily by just giving numbers and sizes of padding. See the example below, that follows this notation:

h1 = (
       (cellSize, numPad, [, increaseFactor]),
       (cellSize, numCore),
       (cellSize, numPad, [, increaseFactor])
     )

Note

You can center your mesh by passing a ‘C’ for the x0[i] position. A ‘N’ will make the entire mesh negative, and a ‘0’ (or a 0) will make the mesh start at zero.

(Source code, png, hires.png, pdf)

../../_images/Mesh_Tensor_Creation-1.png
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
from SimPEG import Mesh
import matplotlib.pyplot as plt


def run(plotIt=True):
    """

        Mesh: Tensor: Creation
        ======================

        For tensor meshes, there are some functions that can come
        in handy. For example, creating mesh tensors can be a bit time
        consuming, these can be created speedily by just giving numbers
        and sizes of padding. See the example below, that follows this
        notation::

            h1 = (
                   (cellSize, numPad, [, increaseFactor]),
                   (cellSize, numCore),
                   (cellSize, numPad, [, increaseFactor])
                 )

        .. note::

            You can center your mesh by passing a 'C' for the x0[i] position.
            A 'N' will make the entire mesh negative, and a '0' (or a 0) will
            make the mesh start at zero.

    """
    h1 = [(10, 5, -1.3), (5, 20), (10, 3, 1.3)]
    M = Mesh.TensorMesh([h1, h1], x0='CN')
    if plotIt:
        M.plotGrid()

if __name__ == '__main__':
    run()
    plt.show()