C code 
home


Introduction: Literate programming is a way to enhance documentation of code, to help the reader understand what is being computed. Code for scientific computing should be written as a technical report, with the code itself embedded in the text. The tools we use are C and LaTeX, which are beautifully entwined in CWEB. As R is a powerful and popular system for statistical computing, an R interface will accompany most of our CWEB code. Our code is distributed under the GNU General Public License. We gratefully acknowledge the support of DFG through the SPP Priority Programme 1590, and the Priority Programme SPP 1819. Simulation code for rapid adaptation of a diploid population to a fixed optimum. A CWEB file (and the sha512sum) for simulating the time until the population has reached a fixed optimum. The documented code is here (and the sha512sum). Seedbank model with unbounded jumps. A CWEB file (and the sha512sum) for computing the time until most recent common ancestor of two lines (manuscript) with J. Blath, A.G. Casanova, and N. Kurt. Also the R interface (and the accompanying sha512sum). Expected normalised branch lengths associated with Lambdacoalescents. A CWEB tarball (and the sha512224sum) for computing the expected normalised branch lengths of a genealogy associated with certain Lambdacoalescents (manuscript) with M. Birkner and J. Blath. Expected normalised branch lengths associated with population growth. A CWEB tarball (and the sha512sum) for computing (in C via R) normalised expected branch lengths associated with population growth. And CWEB code for computing in C (and the sha512sum). A manuscript with M. Birkner, J. Blath, and F. Freund has been accepted by Genetics. Drawing branch lengths under population growth. A CWEB tarball (and the sha512224sum) for drawing (in C via R) branch lengths under either algebraic or exponential population growth. A manuscript with M. Birkner, J. Blath, and F. Freund is in preparation. Expected branch lengths under population growth models. A CWEB tarball (and the sha512224sum) for computing expected branch lengths under either algebraic or exponential population growth. Also computes distances between normalised expected branch lengths. Branch lengths under Lambda coalescents. A CWEB tarball (and the sha512224sum) for drawing branch lengths associated with Lambda coalescents. Coalescent simulator for a seed bank coalescent. A CWEB tarball (and the sha512224sum) with C code for generating random data under a seed bank coalescent (manuscript with J. Blath, AG. Casanova, N. Kurt, and M. WilkeBerenguer). Random branch lengths under Xicoalescents.A CWEB tarball (and the sha512224sum) for drawing branch lengths of a genealogy associated with Xicoalescents (manuscript with J Blath, MC Cronjager, and M Hammer). Exact normalised expected branch lengths under Xicoalescents. A CWEB tarball (and the sha512224sum) for computing exact (normalised) expected branch lengths of a genealogy associated with Xicoalescents (manuscript with J Blath, MC Cronjager, and M Hammer). Exact normalised expected branch lengths under Xicoalescents: R interface A CWEB tarball (and the sha512224sum) for computing exact (normalised) expected branch lengths of a genealogy associated with Xicoalescents in R with a C code as backbone; (manuscript with J Blath, MC Cronjager, and M Hammer). 