Working on the C_xi.c implementation of the DK14 profile.

src/C_xi.c

@@ -1,4 +1,5 @@
 #include "C_xi.h"
+#include "C_peak_height.h"
 #include "C_power.h"
 
 #include "gsl/gsl_integration.h"
@@ -255,15 +256,16 @@ double xi_mm_at_R_exact(double R, double*k, double*P, int Nk){
  */
 
 int calc_xi_DK(double*R, int NR, double M, double rhos, double rs, double alpha, double be, double se, double beta, double gamma, int delta, double om, double*xi){
+  double rhom = rhomconst*om; //SM h^2/Mpc^3
   xi[0] = 0;
   double*rho_ein = (double*)malloc(NR*sizeof(double));
   double*f_trans = (double*)malloc(NR*sizeof(double));
   double*rho_outer = (double*)malloc(NR*sizeof(double));
   int i;
   double g_b = gamma/beta;
-  double r_t = 0; //NEED nu for this
-  double rhom = rhomconst*om; //SM h^2/Mpc^3
+  //Compute R200m
   double Rdelta = pow(M/(1.33333333333*M_PI*rhom*delta), 0.33333333333);
+  double r_t = 0; //r_t = (1.9-0.18*nu)*R200m; //NEED nu for this
   if (rhos < 0){ //means it wasn't passed in
     rhos = rhos_einasto_at_M(M, rs, alpha, delta, om);
   }