/* Last edited on 2023-09-07 19:03:37 by stolfi */

void ralent_choose_arad(int32_t ni, r2_t arad[]);
  /* Stores into {arad[0..ni-1]} a random alignment vector
    meant to be the radius of the basic ellipsoid {\RE}.
    Some of the coordinates will be zero, the others
    will be in the range {[1.0 _ 2.0)}. */
    
void ralent_choose_arad(int32_t ni, r2_t arad[])
  { 
    double rmax = 4.999;
    double rmin = 1.500;
    r2_t zfrac = ...;
    r2_align_throw_arad(ni, zfrac, rmin, rmax, arad, TRUE);
    return;
  }