adjust_dual
adjust_dual(w + i, nr, w + i + nrc, nr <= nc, mu);
sum += adjust_dual(w + i, 1, w + i + nr, nc <= nr, mu);