Skip to content

[ENH] Add optional Transform input to meshes_to_subsurface #19

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

Draft
wants to merge 1 commit into
base: mixing_keops_torch
Choose a base branch
from
Draft
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 7 additions & 2 deletions gempy_engine/core/data/raw_arrays_solution.py
Original file line number Diff line number Diff line change
Expand Up @@ -4,6 +4,7 @@

import numpy as np

from .transforms import Transform
from ..backend_tensor import BackendTensor
from .interp_output import InterpOutput
from .dual_contouring_mesh import DualContouringMesh
Expand Down Expand Up @@ -112,7 +113,7 @@ def _set_scalar_field_at_surface_points(self, octree_output: OctreeLevel):

self.scalar_field_at_surface_points = temp_list

def meshes_to_subsurface(self):
def meshes_to_subsurface(self, input_transform: Transform | None = None):
ss = require_subsurface()
pd = require_pandas()

Expand All @@ -129,9 +130,13 @@ def meshes_to_subsurface(self):

concatenated_id_array = np.concatenate(vertex_id_array)
concatenated_cell_id_array = np.concatenate(cell_id_array)

all_vertex = np.concatenate(vertex)
if input_transform:
all_vertex = input_transform.apply_inverse(all_vertex)

meshes: ss.UnstructuredData = ss.UnstructuredData.from_array(
vertex=np.concatenate(vertex),
vertex=all_vertex,
cells=np.concatenate(simplex_list),
vertex_attr=pd.DataFrame({'id': concatenated_id_array}),
cells_attr=pd.DataFrame({'id': concatenated_cell_id_array})
Expand Down