LU factorization with maximum performances on FPS architectures 38/64 bit

A technique for dense linear system solution is presented which reaches the maximum performances on attached processors like FPS-120, 5000 and X64 using the Fortran language with calls to the vector routines. Starting from the Dongarra's LU factorization algorithm the key idea is to carry out a pseudo-transposition of the lower triangular matrix L (including the main diagonal) around the minor diagonal. The pseudo-transposition allows to carry out all the matrix vector operations involved in LU factorization with only stride 1 dot product operations which, using the TM Auxiliary Memory and the TMDOT routine, can be executed in the FPS processor obtaining the maximum speed. Since the algorithm uses only vector instructions it is fully portable on all the FPS 38/64 bit machines and in general on all the vector computers with a similar memory structure. Furthermore the algorithm can be easily translated into the new FORTRAN 8X, which will probably become the standard for future SIMD computers for numerical applications. The algorithm has been implemented on a FPS-100 yielding the asymptotic speed r?=8 MegaFLOPS (FPS-100 peak performances) and the half performances length N1/2 = 235. The N1/2 value could be lowered by using the APAL Assembly Language to code some critical parts, losing however the code portability.