u1 = c[iii, j1] - v[j1] - umin