mirror of
https://bitbucket.org/cosmicvoids/vide_public.git
synced 2025-07-04 07:11:12 +00:00
60 lines
1.1 KiB
C
60 lines
1.1 KiB
C
#include <stdlib.h>
|
|
#include <math.h>
|
|
#include <stdio.h>
|
|
#ifdef PARALLEL
|
|
#include <mpi.h>
|
|
#endif
|
|
#include "header.h"
|
|
|
|
void linlin_bins()
|
|
{
|
|
FILE *fp,*fp2;
|
|
char fname[100];
|
|
int nsize_rpi,nsize_rsig,i,j;
|
|
double rp,rs,r,dx1,dx2,dx3,dx4,dx5,rslo,rshi,drs,binsize_rpi=2;
|
|
float x1;
|
|
|
|
|
|
nsize_rpi=15;
|
|
nsize_rsig=15;
|
|
binsize_rpi=2.0;
|
|
|
|
sprintf(fname,"%s.linlin",Task.root_filename);
|
|
fp=fopen(fname,"w");
|
|
|
|
rshi=0;
|
|
for(i=1;i<=nsize_rsig;++i)
|
|
{
|
|
rslo=rshi;
|
|
rshi=rslo+binsize_rpi;
|
|
drs=rshi-rslo;
|
|
rs=0.5*(rshi+rslo);
|
|
|
|
for(j=1;j<=nsize_rpi;++j)
|
|
{
|
|
rp=(j-0.5)*binsize_rpi;
|
|
r=sqrt(rs*rs+rp*rp);
|
|
if(!KAISER)
|
|
{
|
|
dx4=two_halo(rs,rp);
|
|
dx1=one_halo(rs,rp);
|
|
}
|
|
else
|
|
{
|
|
dx1=0;
|
|
dx4=0;
|
|
}
|
|
/*
|
|
dx5=integrated_bin(rslo,(j-1)*binsize_rpi,drs,binsize_rpi,10);
|
|
*/
|
|
dx5=xi2d_interp(rslo,(j-1)*binsize_rpi,rslo+drs,j*binsize_rpi);
|
|
fprintf(fp,"%e %e %e %e %e\n",rs,rp,dx5,dx1,dx4);
|
|
fprintf(stdout,"%e %e %e %e %e\n",rs,rp,dx5,dx1,dx4);
|
|
fflush(stdout);
|
|
fflush(fp);
|
|
}
|
|
}
|
|
fclose(fp);
|
|
|
|
|
|
}
|