Making enhanced, interpolated Nanjing lambdas the default

[ilyamandel]

The goal of this PR is to make enhanced, interpolated Nanjing lambdas the default when these are used.

The PR also adds a long-promised Matlab post-processing script for comparing runs. (In order to ease run comparison, several outputs are added to the switch log by default.)

This PR fixes the documentation bug in #1377 .

Attached are some of the outputs of the run comparison script between two runs:

./COMPAS -n 10000 --switch-log
vs.
./COMPAS -n 10000 --switch-log --common-envelope-lambda-nanjing-enhanced --common-envelope-lambda-nanjing-interpolate-in-mass --common-envelope-lambda-nanjing-interpolate-in-metallicity

The differences are sufficiently small (for these solar metallicity runs) that this change appears to be safe. Note the more detailed investigation in, e.g., #1343 (reply in thread) .

DCOs: #Merging DNS #Merging NSBH #Merging BBH % BNS via CE % NSBH via CE % BBH via CE
Default: 3 22 29 100 95 86
Enhanced: 3 22 26 100 95 85

LMXBs: #LMXB #NS LMXB
Default: 19 10
Enhanced: 19 10

Supernova statistics: Default Enhanced Interpolated \lambda

Total number of binaries simulated: 10000 10000
Total number of supernovae: 5852 5904
Of these, complete collapse to BH: 809 825
Number of binaries with two SNe: 1854 1867
Number of binaries with one SN: 2144 2170
Number of binaries unbound by SN: 3100 3115
Unbound by first SN: 2712 2726
Unbound by second SN: 388 389

Total number of USSN: 180 186
Total number of ECSN: 176 176
Total number of AIC: 0 0
Total number of PISN: 0 0
Total number of PPISN: 0 0

#SNe from mass-transfering progenitors: 1497 1538
#Stripped-envelope SNe: 2054 2098
Of these, not AIC or ECSN: 1880 1924
Of these, stripped by winds, no RLOF: 562 565
Of these, # stripped by RLOF: 1318 1359
Of these, previous CE: 242 269
Or double-core CE simultaneous with SN: 52 60

[jeffriley]

Thanks @ilyamandel . Just a quick review:

(a) in changelog.h you have "Added several outputs to the Switch log..."

"the Switch log" -> BSE switch log

which leads me to,

(b) why not add stellar type, mass, and radius to the SSE switch log? (and if you do, you'll need to make the appropriate changes to the documentation)

EDIT: And as long as you are changing program-options-list-defaults.rst can you fix "multicative" in the description of --timestep-multipliers?

[ilyamandel]

@jeffriley : I hadn't added it to the SSE Switch log because I didn't need it for my particular test, but I agree it's a good idea. So added.

[jeffriley]

Thanks @ilyamandel !