DM scattering terms in the nu sectro are zero

nuDM.ini

@@ -23,7 +23,7 @@ N_ur = 3.046
 -------------------
 ---> g-cdm coupling
 -------------------
-u_nuDM_0   = 0.
+u_nuDM_0   = 1.e-5
 n_nuDM     = 0.
 alpha_nuDM = 0., 1.0, 1.5 # the firsts to values are for detla and theta and not used!
 			  # the final value defines all alpha_l for all higher values of l
@@ -36,7 +36,7 @@ sBBN file  = ./bbn/sBBN.dat
 ---> OUTPUT
 -----------
 output               = mPk
-root                 = output/2018-12-07_mixed-damping-study/lcdm_extra_
+root                 = output/2018-12-07_mixed-damping-study/n0_u1em5_freeNU_
 #format               = class
 k_per_decade_for_pk  = 100
 P_k_max_h/Mpc        = 100

source/perturbations.c

@@ -7386,7 +7386,7 @@ int perturb_derivs(double tau,
           -(1.-ppt->three_ceff2_ur)*a_prime_over_a*y[pv->index_pt_theta_ur];
 
 	if(pth->has_coupling_nuDM==_TRUE_)
-	  dy[pv->index_pt_theta_ur] += pvecthermo[pth->index_th_dmu_nuDM]*(y[pv->index_pt_theta_cdm]-y[pv->index_pt_theta_ur]);
+	  dy[pv->index_pt_theta_ur] += 0.*pvecthermo[pth->index_th_dmu_nuDM]*(y[pv->index_pt_theta_cdm]-y[pv->index_pt_theta_ur]);
 
         if(ppw->approx[ppw->index_ap_ufa] == (int)ufa_off) {
 
@@ -7398,21 +7398,21 @@ int perturb_derivs(double tau,
                  // non-standard term, non-zero if cvis2_ur not 1/3
                  -(1.-ppt->three_cvis2_ur)*(8./15.*(y[pv->index_pt_theta_ur]+metric_shear)));
 	  if(pth->has_coupling_nuDM==_TRUE_)
-	    dy[pv->index_pt_shear_ur] += -0.5*ppt->alpha_nuDM[2]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_shear_ur];
+	    dy[pv->index_pt_shear_ur] += 0.*-0.5*ppt->alpha_nuDM[2]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_shear_ur];
 
           /** - -----> Exact ur l=3 */
           l = 3;
           dy[pv->index_pt_l3_ur] = k/(2.*l+1.)*
             (l*2.*s_l[l]*s_l[2]*y[pv->index_pt_shear_ur]-(l+1.)*s_l[l+1]*y[pv->index_pt_l3_ur+1]);
 	  if(pth->has_coupling_nuDM==_TRUE_)
-	    dy[pv->index_pt_l3_ur] -= ppt->alpha_nuDM[l]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_l3_ur];
+	    dy[pv->index_pt_l3_ur] -= 0.*ppt->alpha_nuDM[l]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_l3_ur];
 
           /** - -----> Exactc ur l>3 */
           for (l = 4; l < pv->l_max_ur; l++) {
             dy[pv->index_pt_delta_ur+l] = k/(2.*l+1)*
               (l*s_l[l]*y[pv->index_pt_delta_ur+l-1]-(l+1.)*s_l[l+1]*y[pv->index_pt_delta_ur+l+1]);
 	    if(pth->has_coupling_nuDM==_TRUE_)
-	      dy[pv->index_pt_delta_ur+l] -= ppt->alpha_nuDM[l]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_delta_ur+l];
+	      dy[pv->index_pt_delta_ur+l] -= 0.*ppt->alpha_nuDM[l]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_delta_ur+l];
 	  }
 
           /** - -----> Exact ur lmax_ur */
@@ -7420,7 +7420,7 @@ int perturb_derivs(double tau,
           dy[pv->index_pt_delta_ur+l] =
             k*(s_l[l]*y[pv->index_pt_delta_ur+l-1]-(1.+l)*cotKgen*y[pv->index_pt_delta_ur+l]);
 	  if(pth->has_coupling_nuDM==_TRUE_)
-	    dy[pv->index_pt_delta_ur+l] -= ppt->alpha_nuDM[l]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_delta_ur+l];
+	    dy[pv->index_pt_delta_ur+l] -= 0.*ppt->alpha_nuDM[l]*pvecthermo[pth->index_th_dmu_nuDM]*y[pv->index_pt_delta_ur+l];
         }
 
         else {