[Mesa-users] Question about RSP and equilibrium quantities

[Adam Jermyn]

If you don’t want to push dt way down you can interpolate in time and integrate over one period.

A neat example of this is using linear interpolation, so that the integral gives a trapezoid rule. For genuinely periodic functions the trapezoidal rule converges exponentially in the number of samples. The output from RSP may not be precisely periodic, so this won’t push the variance  down to machine precision, but it should be enough that what’s left is “real” (in the sense of being variation from period to period in what RSP outputs) rather than an aliasing effect.

-Adam
On Jan 21, 2020, 7:11 PM -0500, Bill Paxton <paxton at kitp.ucsb.edu>, wrote:
>
> > On Jan 21, 2020, at 4:08 PM, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> >
> > I increased the resolution by turning down the history interval to 1. Is there a way to force more steps per cycle?
>
> RSP adjusts the timestep at each cycle according to the recent period and the control parameter
>
> RSP_target_steps_per_cycle
>
> The default value for that is 600.     Try 1200 or 1800.
>
> -b
>
>
> _______________________________________________
> 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/20200121/784f66e2/attachment.htm>