Schuster_enforce_constraints_on_v function

Compute v matrix subject to constraints on rows 1..r-1.