Page 1 of 1

CAMB: question on extra evolution variable

Posted: June 13 2019
by Angela Chen
Hi Anthony,

I added a set of extra evolution variables for dark radiation, and I have some concerns about my practice here.
What I did was to add EV%dr_ix as the equation index of dark radiation, and add the same number of equations as massless neutrino in subroutine 'SetupScalarArrayIndices', then did adjustments to subroutine 'initial' and derivs. If my understanding is right, now when I'm doing GaugeInterface_EvolveScal(EV,tau,y,State%Transfer_Times(itf),tol1,ind,c,w), y variable has clx,q,pi for all components including my extension for dr, but since I did nothing to EV%OutputTransfer, I won't get the transfer data of dr in the end.
My question is, will this practice mess up the output of other transfer functions? Is not getting dr perturbation evolution the only aftermath that I didn't do OutputTransfer and only adapting 'y' structure in the code? Specifically, I was worried about that EV%OutputTransfer components mismatching the differential equation components would cause some transfer data to be put into wrong boxes.

Thanks,
Angela

Re: CAMB: question on extra evolution variable

Posted: June 13 2019
by Antony Lewis
I think not setting OutputTransfer should be fine if you don't want to look at/calculate anything with the transfer functions for your new components.

Re: CAMB: question on extra evolution variable

Posted: June 13 2019
by Angela Chen
Antony Lewis wrote:
June 13 2019
I think not setting OutputTransfer should be fine if you don't want to look at/calculate anything with the transfer functions for your new components.
Thanks Antony, sorry for typing your name wrong! Another further concern is, would lacking this part of dr transfer function information cause error in the integration of Cls? and later in the calculation of lensed Cl?

Re: CAMB: question on extra evolution variable

Posted: June 14 2019
by Antony Lewis
No, the sources for the CL are calculated as separate things.

Re: CAMB: question on extra evolution variable

Posted: June 18 2019
by Angela Chen
Antony Lewis wrote:
June 14 2019
No, the sources for the CL are calculated as separate things.
Ok, thanks!