Sinv&Mprd Procedures

[rus_nur]

I need to retranslate some source, which one was developed
for EC(which one very similar to IBM370)
I need sources for procedures: "SINV", "MPRD"

[durisinm]

I searched Google with the keywords Fortran, sinv, and mprd. There was one hit for http://pdp-10.trailing-edge.com/decuslib10-04/01/43,50342/progms.all.

The MPRD function seems to be for multiplying two matrices. The SINV function is for inverting a symmetric positive definite matrix by factorization.

If you have access to the IMSL routines that come with CVF Pro, then you can find equivalent functions there. The Compaq CXML library may also have equivalent functions, and I believe that Intel also has a math library.

Mike

[Jugoslav_Dujic]

As Mike said, these routines can be easily replaced using CXML or IMSL functions. CXML is based on Lapack, (which is the ancestor of Linpack) and BLAS routines, which are widely recognized as "standard" routines for linear algebra. Even if linking with CXML is not an option, you can find sources of Linpack/Lapack and underlying BLAS routines on Netlib (www.netlib.org/linpack or /lapack) -- I believe sources for matrix decomposition in Linpack and Lapack are identical.

In Linpack, look for sgefa (factorization) + sgedi (inversion), or, MUCH faster, replace inversion with linear system solution (sgesl) if appropriate. If the system is symmetric, appropriate routines are ssifa, ssidi, ssisl.

Matrix multiplication is probably part of BLAS (can't tell you offhand), but you can use intrinsic F90 MATMUL instead.

Jugoslav