contrib/gcc/lambda-mat.c

1 /* Integer matrix math routines
32 /* Allocate a matrix of M rows x  N cols.  */
48 /* Copy the elements of M x N matrix MAT1 to MAT2.  */
60 /* Store the N x N identity matrix in MAT.  */
72 /* Return true if MAT is the identity matrix of SIZE */
95 /* Negate the elements of the M x N matrix MAT1 and store it in MAT2.  */
106 /* Take the transpose of matrix MAT1 and store it in MAT2.
107    MAT1 is an M x N matrix, so MAT2 must be N x M.  */
146    MAT1 is an M x R matrix, and MAT2 is R x N.  The resulting MAT2
167 /* Get column COL from the matrix MAT and store it in VEC.  MAT has
196 /* Swap rows R1 and R2 in matrix MAT.  */
208 /* Add a multiple of row R1 of matrix MAT with N columns to row R2:
223 /* Negate row R1 of matrix MAT which has N columns.  */
231 /* Multiply row R1 of matrix MAT with N columns by CONST1.  */
242 /* Exchange COL1 and COL2 in matrix MAT. M is the number of rows.  */
257 /* Add a multiple of column C1 of matrix MAT with M rows to column C2:
272 /* Negate column C1 of matrix MAT which has M rows.  */
283 /* Multiply column C1 of matrix MAT with M rows by CONST1.  */
294 /* Compute the inverse of the N x N matrix MAT and store it in INV.
302    and the identity matrix in parallel.  The inverse is the result of applying
303    the same operations on the identity matrix that reduce MAT to the identity
304    matrix.
306    When MAT is a 2 x 2 matrix, we don't go through the whole process, because
355      INV which starts as the identity matrix.  N is the number of rows and
391   /* Reduce TEMP from a lower triangular to the identity matrix.  Also compute
393      diagonal of TEMP.  If one of these elements is 0, the matrix has 0 as an
403       /* The matrix must not be singular.  */
437 /* Decompose a N x N matrix MAT to a product of a lower triangular H
438    and a unimodular U matrix such that MAT = H.U.  N is the size of
482 /* Given an M x N integer matrix A, this function determines an M x
483    M unimodular matrix U, and an M x N echelon matrix S such that
527 /* Given an M x N integer matrix A, this function determines an M x M
528    unimodular matrix V, and an M x N echelon matrix S such that "A =
633 /* Multiply a vector VEC by a matrix MAT.
634    MAT is an M*N matrix, and VEC is a vector with length N.  The result
638 lambda_matrix_vector_mult (lambda_matrix matrix, int m, int n,
646       dest[i] += matrix[i][j] * vec[j];
649 /* Print out an M x N matrix MAT to OUTFILE.  */
652 print_lambda_matrix (FILE * outfile, lambda_matrix matrix, int m, int n)
657     print_lambda_vector (outfile, matrix[i], n);