It seems like there is a bug in description of RES and AYB in
$ less sla_lin_berr.f
*> RES is *DOUBLE PRECISION* array, dimension (N,NRHS)
*> The residual matrix, i.e., the matrix R in the relative backward
*> error formula above.
*> \param[in] AYB
*> AYB is*DOUBLE PRECISION *array, dimension (N, NRHS)
*> The denominator in the relative backward error formula above, i.e.,
*> the matrix abs(op(A_s))*abs(Y) + abs(B_s). The matrices A, Y, and B
*> are from iterative refinement (see cla_gerfsx_extended.f).
Both RES, and AYB have to be described as *REAL*.
-------------- next part --------------
An HTML attachment was scrubbed...