Add reference samples in divergence outputs

[TuomasBorman]

Include reference samples to output of divergence functions. That requires new parameter to handle the name. I am thinking if we just convert name to vector. That must be length of three: c("divergence", "time_diff", "ref_samples")