Skip to content

[TEST] Add private test for Terranigma 2025_I #1029

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 14 commits into
base: main
Choose a base branch
from
Open

[TEST] Add private test for Terranigma 2025_I #1029

wants to merge 14 commits into from

Conversation

Leguark
Copy link
Member

@Leguark Leguark commented May 18, 2025

Description

Added a test for the 2025_I case in the private test suite. The test creates a geological model with multiple structural elements, sets up the model's extent, adds orientations, and computes the model using PyTorch backend. The test also includes visualization code using gempy_viewer to plot the 3D model with specific display settings.

Copy link
Member Author

Leguark commented May 18, 2025

This stack of pull requests is managed by Graphite. Learn more about stacking.

@Leguark Leguark changed the title [TEST] Added private test [TEST] Add private test for Terranigma 2025_I May 18, 2025
@Leguark Leguark marked this pull request as ready for review May 18, 2025 10:06
@graphite-app graphite-app bot added the gempy 3 Will come with the next major update label May 18, 2025
Copy link

graphite-app bot commented May 18, 2025

Graphite Automations

"Add gempy label" took an action on this PR • (05/18/25)

1 label was added to this PR based on Miguel de la Varga's automation.

@Leguark Leguark marked this pull request as draft May 18, 2025 12:28
@Leguark Leguark marked this pull request as ready for review May 19, 2025 11:19
@Leguark Leguark force-pushed the new_test branch 2 times, most recently from a533f60 to 7b82958 Compare May 19, 2025 11:40
@Leguark Leguark force-pushed the new_test branch 2 times, most recently from 74392a9 to cd25bd1 Compare June 13, 2025 08:59
)

# * The scalar fields can be found for dense and octree grids:
print(geo_model.solutions.raw_arrays.scalar_field_matrix)
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Here I leave the code to grab the scalar field


# * Here it is the way of rescaling one of the axis. Input transform
# * is used (by default) to rescale data into a unit cube but it accepts any transformation matrix.
geo_model.input_transform.scale[2] *= proper_rescale
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is the proper way to do what you were trying. In this model seems to work good

'show_xlabels': False,
'show_ylabels': False,
},
transformed_data=True # * This is interesting, transformed data shows the data as it goes to the interpolation (after applying the transform)
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This arg helps to see what you are interpolating exactly. I had to patch gempy_viewer so you will need to pull or update gempy_viewer

toc = time.perf_counter()
elapsed = toc - tic
print(f"Evaluate model on regular grid centers: {int(elapsed / 60)} minutes {int(elapsed % 60)} seconds")
print(f"Number of cells evaluated: {centers.shape[0]}")
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Added number of cells evaluation

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Blocked gempy 3 Will come with the next major update
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants