Beispiel #1
0
extern inline void montecarlo_compute_distances(rpacket_t *packet, storage_model_t *storage)
{
  // Check if the last line was the same nu as the current line.
  if (rpacket_get_close_line(packet))
    {
      // If so set the distance to the line to 0.0
      rpacket_set_d_line(packet, 0.0);
      // Reset close_line.
      rpacket_set_close_line(packet, false);
    }
  else
    {
      rpacket_set_d_boundary(packet, compute_distance2boundary(packet, storage));
      double d_line;
      compute_distance2line(packet, storage, &d_line);
      rpacket_set_d_line(packet, d_line);
      rpacket_set_d_electron(packet, compute_distance2electron(packet, storage));
    }
}
Beispiel #2
0
bool
test_rpacket_get_d_electron(double value) {
    rpacket_t rp;
    rpacket_set_d_electron(&rp, value);
    return value==rpacket_get_d_electron(&rp);
}