`PRIMA_INTEGER_KIND` slightly affects the result when compiled with `-Ofast -fno-stack-arrays`

libprima / prima

PRIMA is a package for solving general nonlinear optimization problems without using derivatives. It provides the reference implementation for Powell's derivative-free optimization methods, i.e., COBYLA, UOBYQA, NEWUOA, BOBYQA, and LINCOA. PRIMA means Reference Implementation for Powell's methods with Modernization and Amelioration, P for Powell.

http://libprima.net

BSD 3-Clause "New" or "Revised" License

304 stars 40 forks source link

`PRIMA_INTEGER_KIND` slightly affects the result when compiled with `-Ofast -fno-stack-arrays` #186

Open zaikunzhang opened 5 months ago

zaikunzhang commented 5 months ago

According to https://github.com/primalib/prima/actions/workflows/profile_integer_kind.yml, PRIMA_INTEGER_KIND slightly affects the result of the following solvers when compiled with -Ofast -fno-stack-arrays.

PRIMA_INTEGER_KIND = 16: uobyqa, newuoa, bobyqa, cobyla
PRIMA_INTEGER_KIND = 64: cobyla

No visible difference is observed when compiled with -ffast-math.

This is strange.