  In this note we describe two modifications of the ScaLAPACK subroutines
PxGEQPF for computing the QR factorization with the Businger-Golub column
pivoting. First, we resolve a subtle numerical instability in the same way as
we have done it for the LAPACK subroutines xGEQPF, xGEQP3 in 2006. [LAPACK
Working Note 176 (2006); ACM Trans. Math. Softw. 2008]. The problem originates
in the first release of LINPACK in the 1970's: due to severe cancellations in
the down-dating of partial column norms, the pivoting procedure may be in the
dark completely about the true norms of the pivot column candidates. This may
cause miss-pivoting, and as a result loss of the important rank revealing
structure of the computed triangular factor, with severe consequences on other
solvers that rely on the rank revealing pivoting. The instability is so subtle
that, e.g., inserting a WRITE statement or changing the process topology can
drastically change the result. Secondly, we also correct a programming error in
the complex subroutines PCGEQPF, PZGEQPF, which also causes wrong pivoting
because of erroneous use of PSCNRM2, PDZNRM2 for the explicit norm computation.
