.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "gallery/comparisons/3D_triaxial_SimPEG.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_gallery_comparisons_3D_triaxial_SimPEG.py: 3. SimPEG: 3D with tri-axial anisotropy ======================================= `SimPEG `_ is an open source python package for simulation and gradient based parameter estimation in geophysical applications. Here we compare ``emg3d`` with ``SimPEG`` using the forward solver ``Pardiso``. .. GENERATED FROM PYTHON SOURCE LINES 9-21 .. code-block:: default import os import pooch import emg3d import numpy as np import matplotlib.pyplot as plt plt.style.use('bmh') # Adjust this path to a folder of your choice. data_path = os.path.join('..', 'download', '') .. GENERATED FROM PYTHON SOURCE LINES 23-25 Model parameters ---------------- .. GENERATED FROM PYTHON SOURCE LINES 25-41 .. code-block:: default # Depths (0 is sea-surface); # hence a deep sea case where we can ignore the air. water_depth = 3000 target_x = np.r_[-500, 500] target_y = target_x target_z = -water_depth + np.r_[-400, -100] # Resistivities res_sea = 0.33 res_back = [1., 2., 3.] # Background in x-, y-, and z-directions res_target = 100. # Acquisition frequency frequency = 1.0 .. GENERATED FROM PYTHON SOURCE LINES 42-44 Grid ---- .. GENERATED FROM PYTHON SOURCE LINES 44-56 .. code-block:: default vx, vz = np.arange(-20, 21)*100, np.arange(-34, -19)*100 grid = emg3d.meshes.construct_mesh( frequency=frequency, properties=1, center=(0, 0, 0), vector=(vx, vx, vz), center_on_edge=True, ) grid .. raw:: html
 MESH EXTENT CELL WIDTH FACTOR dir nC TensorMesh 131,072 cells x 64 -5,253.29 5,253.29 100.00 515.62 1.15 y 64 -5,253.29 5,253.29 100.00 515.62 1.15 z 32 -6,583.88 1,183.88 100.00 731.62 1.25