485   vector<double> e_b_sq;  in collision_stopping_power()  local
504       e_b_sq.push_back(elm.ionization_energy_[j] * elm.ionization_energy_[j]);  in collision_stopping_power()
521   double rho = sternheimer_adjustment(f, e_b_sq, e_p_sq, n_conduction, log_I,  in collision_stopping_power()
539     double delta = density_effect(f, e_b_sq, e_p_sq, n_conduction, rho, E,  in collision_stopping_power()
1103   const vector<double>& e_b_sq, double e_p_sq, double n_conduction,  in sternheimer_adjustment()  argument
1121       double e_r_sq = e_b_sq[i] * rho * rho + 2.0 / 3.0 * f[i] * e_p_sq;  in sternheimer_adjustment()
1123       gp += e_b_sq[i] * f[i] * rho / e_r_sq;  in sternheimer_adjustment()
1147 double density_effect(const vector<double>& f, const vector<double>& e_b_sq,  in density_effect()  argument
1164       beta_0_sq += f[i] * e_p_sq / (e_b_sq[i] * rho * rho);  in density_effect()
1183       double c = e_b_sq[i] * rho * rho / e_p_sq + w_sq;  in density_effect()
1209     double l_sq = e_b_sq[i] * rho * rho / e_p_sq + 2.0 / 3.0 * f[i];  in density_effect()