coaddition, still misses error propagation
adapt to multi wave format
add cholesky_solve routine