@@ -134,14 +134,11 @@ SHARE
134134 }
135135
136136 *my = ((packing_factor / unit_cell_volume) * 100 * sigma);
137- if (data_transfer.Inhomogenous_incoherent_struct->verbosity){
138- printf("\nDEBUG STATEMENT FOR INHOMOGENOUS INCOHERENT PROCESS:\n"
139- "mu=%g,\t sigma = %g,\t unit cell = %g,\t packing factor = %g\n"
140- "Neutron ray x,y,z = %g, %g, %g,\t Neutron speed = %g, %g, %g\t Neutron time = %g\n",
141- *my, sigma, unit_cell_volume, packing_factor,
142- _particle->x, _particle->y, _particle->z,
143- _particle->vx, _particle->vy, _particle->vz, _particle->t
144- );
137+ if (data_transfer.Inhomogenous_incoherent_struct->verbosity) {
138+ printf ("\nDEBUG STATEMENT FOR INHOMOGENOUS INCOHERENT PROCESS:\n"
139+ "mu=%g,\t sigma = %g,\t unit cell = %g,\t packing factor = %g\n"
140+ "Neutron ray x,y,z = %g, %g, %g,\t Neutron speed = %g, %g, %g\t Neutron time = %g\n",
141+ *my, sigma, unit_cell_volume, packing_factor, _particle->x, _particle->y, _particle->z, _particle->vx, _particle->vy, _particle->vz, _particle->t);
145142 }
146143 // printf("\nMu=%g\tsig=%g\tunit=%g\tpack=%g\n", *my, sigma, unit_cell_volume, packing_factor);
147144 return 1;
@@ -174,7 +171,7 @@ SHARE
174171 if (rand01 () < QE_sampling_frequency) {
175172 v_i = k_length * K2V;
176173 E_i = VS2E * v_i * v_i;
177- dE = lorentzian_width* tan (PI / 2 * randpm1 ());
174+ dE = lorentzian_width * tan (PI / 2 * randpm1 ());
178175 E_f = E_i + dE;
179176 if (E_f <= 0)
180177 return 0;
0 commit comments