Added support for quad precision, small hack for single precision in pipeCG: same problem as in double precision with wrong vector updates. doesn't help much for convergence though.
Forgot to add riemannsolvers_qp.f