[Mesa-users] STELLA Errors and Number of Frequency Bins

[Jake Nicoll]

Jared,

Thank you for looking into this. I will continue to run with as many frequency bins as I can.

I attached a comparison of the UVBRI quasi-bolometric light curves from the mesa.lbol files corresponding to runs with various frequencies. Ignoring the fact that I have not achieved a perfect fit or extended the light curves past the default 170 days yet, you can observe that the light curves seem to be "converging" as the number of frequency bins increases. However, there is some significant oscillation of the tails as one increases from 40 to 60 to 80 bins, and this oscillation makes the dependence of the light curve on the number of frequency bins unclear. It makes me wonder a little whether jumping to, say, 200 frequency bins would significantly impact the light curve. Just food for thought.

Thanks again,
-Jake

________________________________
From: Jared Goldberg <goldberg at physics.ucsb.edu>
Sent: Friday, May 15, 2020 4:39 PM
To: Jake Nicoll <jwnicoll at outlook.com>
Cc: mesa-users at lists.mesastar.org <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] STELLA Errors and Number of Frequency Bins

Hi Jake,
I am running STELLA in the 12778 version of MESA with the mesa.abn and mesa.hyd files created using the example_make_pre_ccsn and example_ccsn_IIP test suite cases (with certain parameters adjusted for specific cases). I am trying to use more frequency bins to better fit individual U, B, V, R, and I color curves to observational data. However, I sometimes encounter errors which disappear when using fewer bins.
For example, I successfully ran STELLA with the attached mesa.hyd and the mesa.abn files with 400 spatial zones and 120 frequency bins. However, when I tried to use 150 frequency bins, I received the error: "0Matrix is numerically singular-decomposition aborted", and STELLA failed to run.

I have seen the error myself, and can reproduce it with your input files. I also ran your mesa.hyd and mesa.abn files on 2 different computer architectures. On my Arch linux desktop I found the same singular matrix message you did, whereas on my mac laptop the calculation ran to ~160 days with 150 frequency bins, but crashed with a "file already open" warning (memory-related) without generating any mesa.lbol or useful output. We suspect a possible memory leak, such as in the opacity table calculations, which is being looked into further. Moreover, the numbers of grid points (for the opacity tables) and frequency bins are somehow capped for computational reasons, and these limits are hard-coded into the version of STELLA that ships with MESA. Thus for the memory issue with finer frequency binning, unfortunately there is simply an upper limit on the number of bins you can use. My solution has been "run with 100-120 bins and see if the results you care about are converged" and the occasional corollary "don't trust the colors much if they're not converged."

I am wondering if someone can help me understand these errors and let me know if there is a way to work around them that does not involve using fewer frequency bins.
I have also received error messages which read "Segmentation fault, invalid memory reference" for other sets of .abn and .hyd files.  I am not sure if these error messages are related, but I have worked around both of them by reducing the number of frequency bins I run with.

I believe the issues are related, and again are likely related to these memory limitations and the calculation of the opacity tables.  The "invalid memory reference" segfault, in fact, occurred for me after the opacity table calculation as STELLA started its hydrodynamical when I tried to run these same .abn/.hyd files with 200 frequency bins. Unfortunately, the best work-around is to use fewer frequency bins, and check convergence. (Say, run with 100, 120, and 130 bins, and see how different the colors look).

Cheers,
~Jared


On Fri, May 8, 2020 at 1:19 PM Jake Nicoll via Mesa-users <mesa-users at lists.mesastar.org<mailto:mesa-users at lists.mesastar.org>> wrote:
Hi all,

I am running STELLA in the 12778 version of MESA with the mesa.abn and mesa.hyd files created using the example_make_pre_ccsn and example_ccsn_IIP test suite cases (with certain parameters adjusted for specific cases). I am trying to use more frequency bins to better fit individual U, B, V, R, and I color curves to observational data. However, I sometimes encounter errors which disappear when using fewer bins.

For example, I successfully ran STELLA with the attached mesa.hyd and the mesa.abn files with 400 spatial zones and 120 frequency bins. However, when I tried to use 150 frequency bins, I received the error: "0Matrix is numerically singular-decomposition aborted", and STELLA failed to run.

I have also received error messages which read "Segmentation fault, invalid memory reference" for other sets of .abn and .hyd files. I am not sure if these error messages are related, but I have worked around both of them by reducing the number of frequency bins I run with.

I am wondering if someone can help me understand these errors and let me know if there is a way to work around them that does not involve using fewer frequency bins.

Best,
-Jake Nicoll
_______________________________________________
mesa-users at lists.mesastar.org<mailto:mesa-users at lists.mesastar.org>
https://lists.mesastar.org/mailman/listinfo/mesa-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200517/e52b4b73/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: LightCurvefreqcomp.pdf
Type: application/pdf
Size: 24190 bytes
Desc: LightCurvefreqcomp.pdf
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200517/e52b4b73/attachment.pdf>