def METHOD_NAME(par):
    """Create small dataset with several hyperbolic events and check output
    contains the events at correct time and correct amplitude
    """
    # Data creation
    t0 = 50
    vrms_x = 1.0
    vrms_y = 1.0
    amp = 0.6
    # Create axes
    t, _, x, y = makeaxis(par)
    # Create data
    d, dwav = hyperbolic3d(x, y, t, t0, vrms_x, vrms_y, amp, wav)
    # Assert shape
    assert d.shape[0] == par["ny"]
    assert d.shape[1] == par["nx"]
    assert d.shape[2] == par["nt"]
    assert dwav.shape[0] == par["ny"]
    assert dwav.shape[1] == par["nx"]
    assert dwav.shape[2] == par["nt"]
    # Assert correct position of event
    assert_array_equal(d[par["ny"] // 2, par["nx"] // 2, t0], amp)