[mesa-users] Other_Energy_Implicit: Extra Heat L and R Partial Derivatives

[Thaddeus Komacek]

Hi all,

I am working on a project using other_energy_implicit, which currently 
takes as input partial derivatives of extra_heat with respect to density 
and temperature for the Henyey solver. It would be useful to also have 
the partials with respect to radius and luminosity, for improved 
convergence when considering heating that is tied to the star or planet 
structure.

It appears that in release 6794 the partials of extra_heat with respect 
to R are already coded ("d_extra_heat_dlnR00", "d_extra_heat_dlnRp1") in 
hydro_eqns and hydro_vars but commented out and not part of the 
star_info structure. Would there be a simple way to incorporate these 
into the solver used by MESA, and likewise include the luminosity partials?

Thank you,
Tad

-- 
Thaddeus Komacek
Graduate Student
Lunar and Planetary Laboratory
University of Arizona
Office Phone: 520-626-3814
http://www.lpl.arizona.edu/~tkomacek/
--