geode-explicit

import os from geode_explicit import RegularGrid from geode_common import Point3D, BoundingBox # Define the bounding box for the grid min_point = Point3D([0.0, 0.0, 0.0]) max_point = Point3D([10.0, 10.0, 10.0]) grid_bbox = BoundingBox(min_point, max_point) # Create a 3D regular grid with 10x10x10 cells grid_dimensions = [10, 10, 10] grid = RegularGrid(grid_bbox, grid_dimensions) print(f"Created a RegularGrid with {grid.nb_cells()} cells.") print(f"Grid origin: {grid.grid_origin()}") print(f"Grid cell size: {grid.cell_size()}") # Save the grid to a file output_filename = "my_regular_grid.grd" grid.save_regular_grid(output_filename) print(f"Grid saved to {output_filename}") # Load the grid from the file loaded_grid = RegularGrid.load_regular_grid(output_filename) print(f"Loaded a RegularGrid with {loaded_grid.nb_cells()} cells.") # Clean up the created file os.remove(output_filename)