Remove unnecessary exp/log calls in EOS_Omni

Create issue
Issue #1172 new
Erik Schnetter created an issue

EOS_Omni seems to call exp/log more often than necessary in the nuc_eos table lookup routines. Use algebraic identities to remove them.


Comments (2)

  1. Frank Löffler
    • removed comment

    Can you please provide specific examples, or maybe even a patch? I fear otherwise this ticket is just going to rot here.

  2. anonymous
    • removed comment

    Here is one example:

      barotropiceos_logrho(:) = log10(10.0d0**barotropiceos_logrho(:) * rho_gf)

    This code first exponentiates a quantity, and then takes a logarithm. The code contains many other similar expressions. I believe this code was implemented as a literal translation from a publication, without optimizing its numerics.

  3. Log in to comment