-
-
Notifications
You must be signed in to change notification settings - Fork 253
[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
base: main
Are you sure you want to change the base?
Conversation
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. |
a533f60
to
7b82958
Compare
74392a9
to
cd25bd1
Compare
# Conflicts: # test/test_private/test_terranigma/test_2025_I.py
Added logic to compute and print the number of cells evaluated during octree processing. Also included output for evaluated cell count on the regular grid centers.
Modified the `ve` value to a fixed scale of 1. and enabled `transformed_data` in the 3D plot test for better coverage.
) | ||
|
||
# * The scalar fields can be found for dense and octree grids: | ||
print(geo_model.solutions.raw_arrays.scalar_field_matrix) |
There was a problem hiding this comment.
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 |
There was a problem hiding this comment.
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) |
There was a problem hiding this comment.
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]}") |
There was a problem hiding this comment.
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
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.