Skip to content

Commit a158d1b

Browse files
committed
[WIP] Adding basic implementation of gravity
1 parent fd8dc51 commit a158d1b

File tree

3 files changed

+13
-1
lines changed

3 files changed

+13
-1
lines changed
Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1 +1,2 @@
1-
from _apparent_thickness import apparent_thickness_likelihood
1+
from _apparent_thickness import apparent_thickness_likelihood
2+
from _gravity_inv import gravity_inversion_likelihood

gempy_probability/modules/likelihoods/_apparent_thickness.py

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -5,6 +5,8 @@
55
def apparent_thickness_likelihood(model_solutions: gp.data.Solutions) -> torch.Tensor:
66
"""
77
This function computes the thickness of the geological layer.
8+
9+
Notes: This is not completed
810
"""
911
simulated_well = model_solutions.octrees_output[0].last_output_center.custom_grid_values
1012
thickness = simulated_well.sum()
Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,9 @@
1+
import pyro
2+
import torch
3+
import gempy as gp
4+
5+
def gravity_inversion_likelihood(model_solutions: gp.data.Solutions) -> torch.Tensor:
6+
simulated_geophysics = model_solutions.gravity
7+
pyro.deterministic(r'$\mu_{gravity}$', simulated_geophysics)
8+
9+
return simulated_geophysics

0 commit comments

Comments
 (0)