[Mesa-users] using FreeEOS_builder

[Ebraheem Farag]

Hello all,

I'm wondering if anyone has used the FreeEOS builder package in MESA
recently to build new tables or run the test case.

>From $MESA_DIR/eos/eosFreeEOS_builder. I have version 2.2.1 downloaded and
installed following https://cococubed.com/mesa_market/free_eos.html.

I'm in MESA r15140 for this, but it's essentially identical to r22.11.1,
which I have also tested with identical results. I checked the mesa github
to see what changed between versions, and it only seems to be related to
the location of the FreeEOS package.
------------------
The primary issue I am running into is with the test case, where I get an
error "STOP bfgs_linesearch: internal bracketing logic error" before the
 './free_eos_table   0.02   0.70   mesa-FreeEOS_02z70x.data' is able to
write that last few table lines.

I've attached the build and install log, and a text file of test case
output ./free_eos_table 0.02 0.70 mesa-FreeEOS_02z70x.data
<https://drive.google.com/file/d/1dmaVTlvS9JtzkcaT3eezpLzxJwAsgxT6/view?usp=gmail>
.

I know this error is coming from line 715 in 'free_eos-2.2.1/src/bfgs.f',
but I have little idea what is happening:

     *do* *while*(need345.or.*.not.*acceptable)

        *if*(not_need345) *then*

          *if*(debug) *then*

            *write*(0,*) 'ai, fai, fprimeai = ', ai, fai, fprimeai

            *write*(0,*) 'bi, fbi = ', bi, fbi

            *write*(0,*) 'fai.eq.fmin', fai*.eq.*fmin

            *write*(0,*) 'fai.le.f0+ai*rho*fprime0',

     &        fai*.le.*f0+ai*rho*fprime0

            *write*(0,*) '(bi-ai)*fprimeai.lt.0.d0.and.'//

     &        'abs(fprimeai).gt.-sigma*fprime0 = ',

     &        (bi-ai)*fprimeai*.lt.*0.d0.and.

     &        abs(fprimeai)*.gt.*-sigma*fprime0

            *write*(0,*) 'fbi.gt.f0+bi*rho*fprime0.or.fbi.ge.fai = ',

     &        fbi*.gt.*f0+bi*rho*fprime0.or.fbi*.ge.*fai

          *endif*

C          test assertion given by PMOO, 2.6.3

          *if*(*.not.*(

     &      (fai*.eq.*fmin)*.and.*

     &      (fai*.le.*f0+ai*rho*fprime0)*.and.*

     &      ((bi-ai)*fprimeai*.lt.*0.d0)*.and.*

     &      (abs(fprimeai)*.gt.*-sigma*fprime0)*.and.*

     &      (fbi*.gt.*f0+bi*rho*fprime0.or.fbi*.ge.*fai)))

     &      *stop* 'bfgs_linesearch: internal bracketing logic error'

p.s. I'm on Mac OS 12.6, but I get the same result on Ubuntu 20LTS. I
rebuilt and ran in multiple MESASDKs on Mac os and Ubuntu from r15140 to
date, to ensure I get the same result regardless. (The build files look
identical between each and the test case fails at the same location.)

Any help/advice would be much appreciated, thank you!

-EbF
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230505/6759a7ff/attachment.htm>
-------------- next part --------------
/Applications/mesasdk/share/doc/free_eos/AUTHORS
/Applications/mesasdk/share/doc/free_eos/COPYING
/Applications/mesasdk/share/doc/free_eos/ChangeLog
/Applications/mesasdk/share/doc/free_eos/ChangeLog.prior_SF
/Applications/mesasdk/share/doc/free_eos/NEWS
/Applications/mesasdk/share/doc/free_eos/README
/Applications/mesasdk/lib/libfree_eos.0.0.8.dylib
/Applications/mesasdk/lib/libfree_eos.0.dylib
/Applications/mesasdk/lib/libfree_eos.dylib
/Applications/mesasdk/share/doc/free_eos/README.developers
/Applications/mesasdk/include/free_eos/bfgs.h
/Applications/mesasdk/include/free_eos/constants.h
/Applications/mesasdk/include/free_eos/diagnostics.h
/Applications/mesasdk/include/free_eos/eff_fit.h
/Applications/mesasdk/include/free_eos/eos_free_calc_argument.h
/Applications/mesasdk/include/free_eos/excitation_block.h
/Applications/mesasdk/include/free_eos/helium1.h
/Applications/mesasdk/include/free_eos/ionization.h
/Applications/mesasdk/include/free_eos/nuvar.h
/Applications/mesasdk/include/free_eos/quad4block.h
/Applications/mesasdk/include/free_eos/statistical_weights.h
/Applications/mesasdk/share/doc/free_eos/README.free_eos_test
/Applications/mesasdk/share/doc/free_eos/README.test_rosenbrock
-------------- next part --------------
A non-text attachment was scrubbed...
Name: make_install.out
Type: application/octet-stream
Size: 7484 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230505/6759a7ff/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: make.out
Type: application/octet-stream
Size: 49855 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230505/6759a7ff/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cmake.out
Type: application/octet-stream
Size: 2128 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230505/6759a7ff/attachment-0002.obj>