simpeg.potential_fields.magnetics.Simulation3DIntegral#
- class simpeg.potential_fields.magnetics.Simulation3DIntegral(mesh, chi=None, chiMap=None, model_type='scalar', is_amplitude_data=False, engine='geoana', numba_parallel=True, **kwargs)[source]#
Bases:
BasePFSimulationMagnetic simulation in integral form.
- Parameters:
- mesh
discretize.TreeMeshordiscretize.TensorMesh Mesh use to run the magnetic simulation.
- survey
simpeg.potential_fields.magnetics.Survey Magnetic survey with information of the receivers.
- active_cells(
n_cells)numpy.ndarray,optional Array that indicates which cells in
meshare active cells.- chi
numpy.ndarray,optional Susceptibility array for the active cells in the mesh.
- chiMap
Mapping,optional Model mapping.
- model_type
str,optional Whether the model are susceptibilities of the cells (
"scalar"), or effective susceptibilities ("vector").- is_amplitude_databool,
optional If True, the returned fields will be the amplitude of the magnetic field. If False, the fields will be returned unmodified.
- sensitivity_dtype
numpy.dtype,optional Data type that will be used to build the sensitivity matrix.
- store_sensitivities{“ram”, “disk”, “forward_only”}
Options for storing sensitivity matrix. There are 3 options
‘ram’: sensitivities are stored in the computer’s RAM
‘disk’: sensitivities are written to a directory
‘forward_only’: you intend only do perform a forward simulation and sensitivities do not need to be stored
- sensitivity_path
str,optional Path to store the sensitivity matrix if
store_sensitivitiesis set to"disk". Default to “./sensitivities”.- engine{“geoana”, “choclo”},
optional Choose which engine should be used to run the forward model.
- numba_parallelbool,
optional If True, the simulation will run in parallel. If False, it will run in serial. If
engineis not"choclo"this argument will be ignored.- ind_active
np.ndarrayofintor bool Deprecated since version 0.23.0: Argument
ind_activeis deprecated in favor ofactive_cellsand will be removed in SimPEG v0.24.0.
- mesh
Attributes
The linear operator.
M: ndarray
Active cells in the mesh.
Magnetic susceptibility (si) physical property model.
Derivative of Magnetic Susceptibility (SI) wrt the model.
Mapping of the inversion model to Magnetic Susceptibility (SI).
A list of solver objects to clean when the model is updated
SimPEG
Counterobject to store iterations and run-times.HasModel.deleteTheseOnModelUpdate has been deprecated.
Engine that will be used to run the simulation.
active_cells.ind_active has been deprecated.
The model for a linear problem physical property model.
Mesh for the integral potential field simulations.
The inversion model.
model_type.modelType has been deprecated.
Derivative of The model for a linear problem wrt the model.
Mapping of the inversion model to The model for a linear problem.
Type of magnetization model
Number of data
Number of processes to use for forward modeling.
True if a model is necessary
Run simulation in parallel or single-threaded when using Numba.
dtype of the sensitivity matrix.
Path to directory where sensitivity file is stored.
Options for storing sensitivities.
The survey for the simulation.
Verbose progress printout.
ampDeriv
is_amplitude_data
tmi_projection
Methods
Jtvec(m, v[, f])Dot product between transposed sensitivity matrix and a vector.
Jtvec_approx(m, v[, f])Approximation of the Jacobian transpose times a vector for the model provided.
Jvec(m, v[, f])Dot product between sensitivity matrix and a vector.
Jvec_approx(m, v[, f])Approximation of the Jacobian times a vector for the model provided.
compute_amplitude(b_xyz)Compute amplitude of the magnetic field
dpred([m, f])Predicted data for the model provided.
evaluate_integral(receiver_location, components)Load in the active nodes of a tensor mesh and computes the magnetic forward relation between a cuboid and a given observation location outside the Earth [obsx, obsy, obsz]
fields(model)Return the computed geophysical fields for the model provided.
getJ(m[, f])Sensitivity matrix \(\mathbf{J}\).
getJtJdiag(m[, W, f])Compute diagonal of \(\mathbf{J}^T \mathbf{J}`\).
Return linear operator.
make_synthetic_data(m[, relative_error, ...])Make synthetic data for the model and Gaussian noise provided.
normalized_fields(fields)Return the normalized B fields
residual(m, dobs[, f])The data residual.
Galleries and Tutorials using simpeg.potential_fields.magnetics.Simulation3DIntegral#
Forward Simulation of Total Magnetic Intensity Data
Forward Simulation of Gradiometry Data for Magnetic Vector Models
Sparse Norm Inversion for Total Magnetic Intensity Data on a Tensor Mesh
Cross-gradient Joint Inversion of Gravity and Magnetic Anomaly Data
Joint PGI of Gravity + Magnetic on an Octree mesh using full petrophysical information
Joint PGI of Gravity + Magnetic on an Octree mesh without petrophysical information