Modified solver routines cr.c, bcgs.c and bcgsl.c