Embedding rivers in triangulated irregular networks with linear programming