mirror of
https://bitbucket.org/cosmicvoids/vide_public.git
synced 2025-07-04 07:11:12 +00:00
15 lines
606 B
C
15 lines
606 B
C
if(Work.iquad && Work.iquad_covar)
|
|
{
|
|
calc_rquad(Work.r_quad,rquad,Work.n_quad);
|
|
for(i=0;i<Work.n_quad;++i)
|
|
for(j=0;j<Work.n_quad;++j)
|
|
{
|
|
if(Work.r_quad[i]<Work.rhlo || Work.r_quad[j]<Work.rhlo)continue;
|
|
e=esys_quad(Work.r_quad[i])*rquad[i]*esys_quad(Work.r_quad[j]);
|
|
chi2c+=(rquad[i]-Work.data_h[i])*(rquad[j]-Work.data_h[j])*
|
|
(Work.covar_q[i][j]);
|
|
if(!ThisTask && OUTPUTZ)
|
|
printf("CHIQUAD%d %d %d %d %f %e %e %e %e %e %e\n",Work.imodel,iter,i,j,Work.r_quad[j],
|
|
rquad[i],Work.data_h[i],rquad[j],Work.data_h[j],Work.covar_h[i][j],e);
|
|
}
|
|
}
|