For functions depending on two variables, we automatically construct tri6 angulations subject to the condition that the continuous, piecewise linear approxi7 mation, underor overestimation never deviates more than a given δ -tolerance from 8 the original function over a given domain. This tolerance is ensured by solving sub9 problems over each triangle to global optimality. The continuous, piec...