def _get_minimal_smoothing(N):
    """An heuristic function for the "most optimal"
    (in some sense) smoothing parameter for inverting
    an integral operator of quadrature dimension *N*."""
    
    return 9.5e-11*N**1.71