From 08250fe8af770e5536d5542377099f31a40b80cd Mon Sep 17 00:00:00 2001
From: Ruben Farinelli <ruben@MacBook-Air-di-Ruben.local>
Date: Wed, 27 Sep 2023 20:25:05 +0200
Subject: [PATCH] Deleted a repeated function in MC routine and modified the
number of angles for the iteration method
---
mc_slab.c | 10 +++-------
solve_rte.c | 11 +----------
2 files changed, 4 insertions(+), 17 deletions(-)
diff --git a/mc_slab.c b/mc_slab.c
index a127c0c..9929141 100644
--- a/mc_slab.c
+++ b/mc_slab.c
@@ -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
@@ -515,7 +516,8 @@ 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 */
// ================================================================== */
diff --git a/solve_rte.c b/solve_rte.c
index 44139b8..9d7ee65 100644
--- a/solve_rte.c
+++ b/solve_rte.c
@@ -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*/
/*=======================================================================*/
--
GitLab