Commit 08250fe8 authored by Ruben Farinelli's avatar Ruben Farinelli
Browse files

Deleted a repeated function in MC routine and modified the number of angles 

for the iteration method
parent 580ff8c2

mc_slab.c

+3 −7
Original line number Diff line number Diff line
@@ -301,6 +301,7 @@ int slab_mc(int nphot, int seed) 
        Magfield_lab[0] = Magfield_lab_new[0];

        Magfield_lab[1] = Magfield_lab_new[1];

        Magfield_lab[2] = Magfield_lab_new[2];

        

      }

      else
@@ -516,6 +517,7 @@ void comptonization(double E_lab, 


  status = depolarization_probability(E_erf, E_new_erf, costheta_erf, phik, PolarFlag);

  

 

  if (status == 1)

  {

    exit(1);
@@ -537,12 +539,6 @@ void comptonization(double E_lab, 


  elmag_lorentztransform(ElectField_new_erf, MagField_Erf, el_vel_vect, Efield_lab_new, Magfield_lab_new, 2);



  // =========================================================================== */

  // Back convert from (K1)' in ERF to lab*/

  // =========================================================================== */



  versor_boost(k_new_erf, el_vel_vect, k_new_lab, 2);



  // ================================================================== */

  // Back convert electron energy from ERF to lab */

  // ================================================================== */

solve_rte.c

+1 −10
Original line number Diff line number Diff line
@@ -32,7 +32,7 @@ int solve_rte(int k_iter, int seed_distribution, int Nstep_tau) 
  printf("\n\nReceveing Nstep_tau %d\n", Nstep_tau);

  int ii, jj, kk, status;



  Nstep_mu = 20;

  Nstep_mu = 30;

  /* Nstep_tau = 100;*/



  tau0 = disktau / 2.;
@@ -191,15 +191,6 @@ int solve_rte(int k_iter, int seed_distribution, int Nstep_tau) 
    ptr_data->Spline_I2d_r_downstream = Spline_I2d_r_downstream;









   /* printf("CHECK HERE !\n");

    double zio=0;

        int status = gsl_spline2d_eval_extrap_e(Spline_I2d_l_downstream, 2*tau0, 0.5, xacc_2d, yacc_2d, &zio);



       printf("================> status %d val %lf\n", status, zio);

*/



    /*=======================================================================*/

    /*Solve the RTE at the quadrature points*/

    /*=======================================================================*/