hiop
hiop copied to clipboard
Ma57 parameter fix
This PR changes some parameter settings in Ma57:
- Set icntl[9] to 1 --- default number of iterative refinement to 1. Right now we initialize it as 10 and replace it to 1 before using it.
- allocate more spaces to some
work
arrary, if icntl[9] > 1 (required by ma57) - allow the combination of
icntl[9]>1
andjob = 0
orjob=2
(but user needs to modify the source code to test these features --- per @pelesh request)
@pelesh @cnpetra
I ran some tests with large ACOPF models with MA57 in this branch and convergence was much worse than when run with develop
branch. I ran it as is without any changes to MA57 parameters. This is my hiop.options
file:
KKTLinsys xdycyd
print_options yes
verbosity_level 3
tolerance 1e-6
#fixed_var relax
duals_init zero
fact_acceptor inertia_free
compute_mode cpu
linear_solver_sparse ma57
@pelesh did you use the same option file in develop
branch? I didn't change any setting in this PR (see the change log). I tried the latest develop
branch with your option file, and all the unit tests work fine.
@pelesh did you use the same option file in
develop
branch? I didn't change any setting in this PR (see the change log). I tried the latestdevelop
branch with your option file, and all the unit tests work fine.
I need to do some more testing. Preliminary tests did not work out of the box.
@pelesh have you tried the changes in this PR?