hiop icon indicating copy to clipboard operation
hiop copied to clipboard

Ma57 parameter fix

Open nychiang opened this issue 1 year ago • 4 comments

This PR changes some parameter settings in Ma57:

  1. 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.
  2. allocate more spaces to some work arrary, if icntl[9] > 1 (required by ma57)
  3. allow the combination of icntl[9]>1 and job = 0 or job=2 (but user needs to modify the source code to test these features --- per @pelesh request)

@pelesh @cnpetra

nychiang avatar Jul 25 '23 19:07 nychiang

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 avatar Jul 27 '23 00:07 pelesh

@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.

nychiang avatar Jul 27 '23 07:07 nychiang

@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.

I need to do some more testing. Preliminary tests did not work out of the box.

pelesh avatar Aug 11 '23 16:08 pelesh

@pelesh have you tried the changes in this PR?

nychiang avatar Nov 27 '23 22:11 nychiang