Introduced a version of the tron solver compatible with Total Quasi-Newton approach.
Two minor bug fixes: always call the constraint function before the nonsquare quasi-Newton method and fix the counting error in the number of BQP iterations.