salvus.mesh.algorithms.optimize_dt

def dgamma(
    m: UnstructuredMesh, ij: npt.NDArray, rho: float = 1.0
) -> npt.NDArray: ...

def dgamma_h(
    m: UnstructuredMesh,
    h_ek_ref: float | npt.NDArray,
    g_exponent: int = 3,
    rho: float = 35.0,
) -> npt.NDArray: ...

def gamma(
    m: UnstructuredMesh,
    ij: npt.NDArray,
    rho: float = 1.0,
    return_gamma_max: bool = False,
    return_gamma_array: bool = False,
) -> npt.NDArray | float: ...

def gamma_h(
    m: UnstructuredMesh,
    h_ek_ref: float | npt.NDArray,
    g_exponent: int = 3,
    rho: float = 35.0,
    return_gamma_h_max: bool = False,
) -> float | npt.NDArray: ...

def optimize_dt(
    m: UnstructuredMesh,
    fixed_side_sets: list[list[set[tuple[int, int]]]] | None = None,
    maxiter: int = 50,
    rho: float = 35.0,
    rho_h: float = 35.0,
    weight_h: float = 0.0,
    h_ek_ref: npt.NDArray | float = 0.0,
    g_exponent: int = 3,
    scale_mesh: bool = True,
    scale_factor: float = 10000.0,
    verbose: bool = False,
) -> UnstructuredMesh: ...

def optimize_dt_locally(
    m: UnstructuredMesh,
    fixed_side_set_names: list[str],
    maxiter: int = 50,
    rho: float = 35.0,
    rho_h: float = 35.0,
    weight_h: float = 0.0,
    h_ek_ref: npt.NDArray[np.floating] | float = 0.0,
    g_exponent: int = 3,
    element_count_opt: int = 100,
    halo_width: int = 3,
    scale_factor: float = 10000.0,
    verbose: bool = False,
) -> UnstructuredMesh: ...