Skip to content

access violation reading with refinement after deletion #242

@veenstrajelmer

Description

@veenstrajelmer

Describe the bug
The code below gives an access violation: OSError: exception: access violation reading 0x00000246736EA038

To Reproduce

import numpy as np
from meshkernel import GeometryList, DeleteMeshOption, ProjectionType, MakeGridParameters, MeshKernel, RefinementType, GriddedSamples, MeshRefinementParameters

lon_min, lon_max, lat_min, lat_max = 99, 114, -7, 13.5
dx = dy = 0.1
crs = 'EPSG:4326'

# create base grid
make_grid_parameters = MakeGridParameters(
    origin_x=lon_min,
    origin_y=lat_min,
    upper_right_x=lon_max,
    upper_right_y=lat_max,
    block_size_x=dx,
    block_size_y=dy,
    )

mk = MeshKernel(projection=ProjectionType.SPHERICAL)
mk.curvilinear_compute_rectangular_grid_on_extension(make_grid_parameters)
mk.curvilinear_convert_to_mesh2d() #convert to ugrid/mesh2d

xx = np.array([98.0, 107.5, 98.0, 98])
yy = np.array([2.5, -7.5, -7.5, 2.5])

delete_pol_geom = GeometryList(x_coordinates=xx, y_coordinates=yy)
mk.mesh2d_delete(geometry_list=delete_pol_geom, 
                        delete_option=DeleteMeshOption.INSIDE_NOT_INTERSECTED,
                        invert_deletion=False)

min_edge_size = 1250 # in meters

lon_np = np.array([lon_min, lon_max])
lat_np = np.array([lat_min, lat_max])
values_np = np.array([-9999, -9999, -9999, -9999])
gridded_samples = GriddedSamples(
    x_coordinates=lon_np,
    y_coordinates=lat_np,
    values=values_np)

#refinement
mesh_refinement_parameters = MeshRefinementParameters(
    min_edge_size=min_edge_size,
    refinement_type=RefinementType.WAVE_COURANT,
    smoothing_iterations=2,
    )

mk.mesh2d_refine_based_on_gridded_samples(
    gridded_samples=gridded_samples,
    mesh_refinement_params=mesh_refinement_parameters,
    use_nodal_refinement=True,
    )

Version info (please complete the following information):

  • OS: Windows
  • Version 8.1.3

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions