def calcPv(D, H):
    """
    Calculate geometric albedo given asteroid diameter and H magnitude.

    Parameters
    ----------
    D : float or `~numpy.ndarray` (N)
        Asteroid diameter in meters.
    H : float or `~numpy.ndarray` (N)
        Absolute H magnitude.

    Returns
    -------
    float or `~numpy.ndarray` (N)
        Geometric albedo.
    """
    return (1.329e6 / (D * 10**(H / 5)))**2