ceil_at function

Rounding with predefined thresholds