docs: def MultimodeCombine(pupils):
      """
      Return the instantaneous coherent fluxes and photometric fluxes for a
      multiway multimode combiner (no spatial filtering)
      """
      fluxes=[np.vdot(pupils[i],pupils[i]).real for i in range(len(pupils))]
      coherentFluxes=[np.vdot(pupils[i],pupils[j])
                      for i in range(1,len(pupils))
                      for j in range(i)]
      return fluxes,coherentFluxes