Working and crosschecked:

  - Computation of Gamma(t), drho(t) without replica
  - Indexing and error propagation