I have a question regarding the fitting code to SDSS power spectrum of CosmoMC.
In the May 2006 version seems the following part has been revised again:
Code: Select all
w=1/(mset%mpk_sdev(:)*mset%mpk_sdev(:))
normV = sum(mpk_WPth*mpk_WPth*w)
tmp=sum(mpk_WPth*mset%mpk_P*w)/normV
chisq = sum(mset%mpk_P*(mset%mpk_P - mpk_WPth*tmp)*w + log(normV)
I am puzzuled about this.
Thank you for any suggestion.
Jun-Qing Xia