contrib/gcc/lambda-mat.c

576 lambda_matrix_first_nz_vec (lambda_matrix mat, int rowsize, int colsize,
585 	  && (lambda_vector_first_nz (mat[j], colsize, startrow) < colsize))
598 			       int colsize, int k, lambda_vector x)
606   M1 = lambda_matrix_new (colsize, colsize);
607   lambda_matrix_transpose (B, M1, rowsize, colsize);
610   M2 = lambda_matrix_new (colsize, colsize);
611   lambda_matrix_mult (B, M1, M2, rowsize, colsize, rowsize);
614   M3 = lambda_matrix_new (colsize, colsize);
618   lambda_matrix_mult (M1, M3, M2, colsize, rowsize, rowsize);
621   lambda_matrix_mult (M2, B, M1, colsize, rowsize, colsize);
622   lambda_matrix_negate (M1, M1, colsize, colsize);
624   I = lambda_matrix_new (colsize, colsize);
625   lambda_matrix_id (I, colsize);
627   lambda_matrix_add_mc (I, determinant, M1, 1, M2, colsize, colsize);
629   lambda_matrix_get_column (M2, colsize, k - 1, x);