A different forward dense matvec to test in MPI environment.
Try an Allgatherv() operation in the rmatvec since Allreduce() appears to be slower.
(Based on previous testing.)
Added the specified 'data_suffix' argument to the stored Jacobian approximation.
Implemented distributed storage of approximate Jacobian when in MPI environment.
Fixed a typo in specifying output file names.
A cost-saving measure: put a cap on the number of BQP (not CG) iterations in TR subproblem.
Model change: specify indices in numpy arrays instead of lists for speed.
Fixed array type for empty arrays.
Add a function to compute the infeasibility term explicitly to observe breakdown of aug. Lag. function for debugging.
Remove iterations limit on TR solver for testing.
Typo fix in log output.
Fixed an important indexing error.
A fix for the MPI case where the number of processors is less than the number of constraints.
Modify total quasi-Newton classes to restart Jacobian approximation after every nonlinear minimization.
Small fix to last commit.
Bring back cap on BQP iterations.
A different approach to restarting the approximate Jacobian, only restart after an excessive number of SBMIN iterations.
Introduce option for restarting the Jacobian every so often during SBMIN solve.
In Total QN version of SBMIN, include an update if the step is rejected and the option is set to true.
Merge branch 'mpi-matvecs-dev' into auglag-QN