/orange/adamginsburg/repos/imf/imf/imf.py:921: RuntimeWarning: divide by zero encountered in log
  if not np.logical_or(np.isfinite(np.log(mmin)), np.isfinite(np.log(tolerance))):

M>8 siblings per primary (Kroupa ODF, closed form):
  m1=   8.0 Msun  N(M>8)=  1.00  N_sib=  0.00
  m1=  10.0 Msun  N(M>8)=  1.35  N_sib=  0.35
  m1=  12.0 Msun  N(M>8)=  1.72  N_sib=  0.72
  m1=  16.0 Msun  N(M>8)=  2.55  N_sib=  1.55
  m1=  20.0 Msun  N(M>8)=  3.47  N_sib=  2.47
  m1=  25.0 Msun  N(M>8)=  4.77  N_sib=  3.77
  m1=  32.0 Msun  N(M>8)=  6.85  N_sib=  5.85
  m1=  40.0 Msun  N(M>8)=  9.66  N_sib=  8.66
  m1=  60.0 Msun  N(M>8)= 19.28  N_sib= 18.28
  m1=  85.0 Msun  N(M>8)= 40.43  N_sib= 39.43
  m1= 120.0 Msun  N(M>8)=131.26  N_sib=130.26

RMS embedded SF sample in LF range: 1571 sources
  primaries with m1>=8 Msun: 1282 (82%)
  median primary mass: 12.4 Msun
  <N_sib (M>8)> per OB primary = 5.95
  <N_OB8 per OB primary>       = 6.95

  => total M>8 OB stars / primaries = 6.95
  => leading-order tau0 correction factor (1+<N_sib>) = 6.95
     tau0: 4.2 Myr -> 29.2 Myr  (leading order)
validation: imf _opt_sample vs closed form (real Kroupa)
  Mcl=   40  m1=   4.8  N(>8) sampled=  0  closed-form=  0.0
  Mcl=   80  m1=   7.8  N(>8) sampled=  0  closed-form=  0.0