38   int i;  in lambda_matrix_new()  local
42   for (i = 0; i < m; i++)  in lambda_matrix_new()
43     mat[i] = lambda_vector_new (n);  in lambda_matrix_new()
54   int i;  in lambda_matrix_copy()  local
56   for (i = 0; i < m; i++)  in lambda_matrix_copy()
57     lambda_vector_copy (mat1[i], mat2[i], n);  in lambda_matrix_copy()
65   int i, j;  in lambda_matrix_id()  local
67   for (i = 0; i < size; i++)  in lambda_matrix_id()
69       mat[i][j] = (i == j) ? 1 : 0;  in lambda_matrix_id()
77   int i, j;  in lambda_matrix_id_p()  local
78   for (i = 0; i < size; i++)  in lambda_matrix_id_p()
81 	if (i == j)  in lambda_matrix_id_p()
83 	    if (mat[i][j] != 1)  in lambda_matrix_id_p()
88 	    if (mat[i][j] != 0)  in lambda_matrix_id_p()
100   int i;  in lambda_matrix_negate()  local
102   for (i = 0; i < m; i++)  in lambda_matrix_negate()
103     lambda_vector_negate (mat1[i], mat2[i], n);  in lambda_matrix_negate()
112   int i, j;  in lambda_matrix_transpose()  local
114   for (i = 0; i < n; i++)  in lambda_matrix_transpose()
116       mat2[i][j] = mat1[j][i];  in lambda_matrix_transpose()
126   int i;  in lambda_matrix_add()  local
128   for (i = 0; i < m; i++)  in lambda_matrix_add()
129     lambda_vector_add (mat1[i], mat2[i], mat3[i], n);  in lambda_matrix_add()
139   int i;  in lambda_matrix_add_mc()  local
141   for (i = 0; i < m; i++)  in lambda_matrix_add_mc()
142     lambda_vector_add_mc (mat1[i], const1, mat2[i], const2, mat3[i], n);  in lambda_matrix_add_mc()
154   int i, j, k;  in lambda_matrix_mult()  local
156   for (i = 0; i < m; i++)  in lambda_matrix_mult()
160 	  mat3[i][j] = 0;  in lambda_matrix_mult()
162 	    mat3[i][j] += mat1[i][k] * mat2[k][j];  in lambda_matrix_mult()
174   int i;  in lambda_matrix_get_column()  local
176   for (i = 0; i < n; i++)  in lambda_matrix_get_column()
177     vec[i] = mat[i][col];  in lambda_matrix_get_column()
185   int i;  in lambda_matrix_delete_rows()  local
189   for (i = to; i < rows; i++)  in lambda_matrix_delete_rows()
190     mat[i - dist] = mat[i];  in lambda_matrix_delete_rows()
192   for (i = rows - dist; i < rows; i++)  in lambda_matrix_delete_rows()
193     mat[i] = NULL;  in lambda_matrix_delete_rows()
214   int i;  in lambda_matrix_row_add()  local
219   for (i = 0; i < n; i++)  in lambda_matrix_row_add()
220     mat[r2][i] += const1 * mat[r1][i];  in lambda_matrix_row_add()
236   int i;  in lambda_matrix_row_mc()  local
238   for (i = 0; i < n; i++)  in lambda_matrix_row_mc()
239     mat[r1][i] *= const1;  in lambda_matrix_row_mc()
247   int i;  in lambda_matrix_col_exchange()  local
249   for (i = 0; i < m; i++)  in lambda_matrix_col_exchange()
251       tmp = mat[i][col1];  in lambda_matrix_col_exchange()
252       mat[i][col1] = mat[i][col2];  in lambda_matrix_col_exchange()
253       mat[i][col2] = tmp;  in lambda_matrix_col_exchange()
263   int i;  in lambda_matrix_col_add()  local
268   for (i = 0; i < m; i++)  in lambda_matrix_col_add()
269     mat[i][c2] += const1 * mat[i][c1];  in lambda_matrix_col_add()
277   int i;  in lambda_matrix_col_negate()  local
279   for (i = 0; i < m; i++)  in lambda_matrix_col_negate()
280     mat[i][c1] *= -1;  in lambda_matrix_col_negate()
288   int i;  in lambda_matrix_col_mc()  local
290   for (i = 0; i < m; i++)  in lambda_matrix_col_mc()
291     mat[i][c1] *= const1;  in lambda_matrix_col_mc()
347   int i, j;  in lambda_matrix_inverse_hard()  local
362       for (i = j; i < n; i++)  in lambda_matrix_inverse_hard()
363 	if (row[i] < 0)  in lambda_matrix_inverse_hard()
365 	    lambda_matrix_col_negate (temp, n, i);  in lambda_matrix_inverse_hard()
366 	    lambda_matrix_col_negate (inv, n, i);  in lambda_matrix_inverse_hard()
377 	  for (i = j + 1; i < n; i++)  in lambda_matrix_inverse_hard()
381 	      factor = -1 * row[i];  in lambda_matrix_inverse_hard()
385 	      lambda_matrix_col_add (temp, n, j, i, factor);  in lambda_matrix_inverse_hard()
386 	      lambda_matrix_col_add (inv, n, j, i, factor);  in lambda_matrix_inverse_hard()
413 	  for (i = 0; i < j; i++)  in lambda_matrix_inverse_hard()
414 	    lambda_matrix_col_mc (inv, n, i, diagonal);  in lambda_matrix_inverse_hard()
415 	  for (i = j + 1; i < n; i++)  in lambda_matrix_inverse_hard()
416 	    lambda_matrix_col_mc (inv, n, i, diagonal);  in lambda_matrix_inverse_hard()
422       for (i = j - 1; i >= 0; i--)  in lambda_matrix_inverse_hard()
424 	  if (row[i])  in lambda_matrix_inverse_hard()
426 	      int factor = -row[i];  in lambda_matrix_inverse_hard()
427 	      lambda_matrix_col_add (temp, n, j, i, factor);  in lambda_matrix_inverse_hard()
428 	      lambda_matrix_col_add (inv, n, j, i, factor);  in lambda_matrix_inverse_hard()
446   int i, j, factor, minimum_col;  in lambda_matrix_hermite()  local
456       for (i = j; i < n; i++)  in lambda_matrix_hermite()
458 	  if (row[i] < 0)  in lambda_matrix_hermite()
460 	      lambda_matrix_col_negate (H, n, i);  in lambda_matrix_hermite()
461 	      lambda_vector_negate (U[i], U[i], n);  in lambda_matrix_hermite()
472 	  for (i = j + 1; i < n; i++)  in lambda_matrix_hermite()
474 	      factor = row[i] / row[j];  in lambda_matrix_hermite()
475 	      lambda_matrix_col_add (H, n, j, i, -1 * factor);  in lambda_matrix_hermite()
476 	      lambda_matrix_row_add (U, n, i, j, factor);  in lambda_matrix_hermite()
493   int i, j, i0 = 0;  in lambda_matrix_right_hermite()  local
503 	  for (i = m - 1; i >= i0; i--)  in lambda_matrix_right_hermite()
505 	      while (S[i][j] != 0)  in lambda_matrix_right_hermite()
509 		  a = S[i-1][j];  in lambda_matrix_right_hermite()
510 		  b = S[i][j];  in lambda_matrix_right_hermite()
516 		  lambda_matrix_row_add (S, n, i, i-1, -factor);  in lambda_matrix_right_hermite()
517 		  lambda_matrix_row_exchange (S, i, i-1);  in lambda_matrix_right_hermite()
519 		  lambda_matrix_row_add (U, m, i, i-1, -factor);  in lambda_matrix_right_hermite()
520 		  lambda_matrix_row_exchange (U, i, i-1);  in lambda_matrix_right_hermite()
538   int i, j, i0 = 0;  in lambda_matrix_left_hermite()  local
548 	  for (i = m - 1; i >= i0; i--)  in lambda_matrix_left_hermite()
550 	      while (S[i][j] != 0)  in lambda_matrix_left_hermite()
554 		  a = S[i-1][j];  in lambda_matrix_left_hermite()
555 		  b = S[i][j];  in lambda_matrix_left_hermite()
561 		  lambda_matrix_row_add (S, n, i, i-1, -factor);  in lambda_matrix_left_hermite()
562 		  lambda_matrix_row_exchange (S, i, i-1);  in lambda_matrix_left_hermite()
564 		  lambda_matrix_col_add (V, m, i-1, i, factor);  in lambda_matrix_left_hermite()
565 		  lambda_matrix_col_exchange (V, m, i, i-1);  in lambda_matrix_left_hermite()
641   int i, j;  in lambda_matrix_vector_mult()  local
644   for (i = 0; i < m; i++)  in lambda_matrix_vector_mult()
646       dest[i] += matrix[i][j] * vec[j];  in lambda_matrix_vector_mult()
654   int i;  in print_lambda_matrix()  local
656   for (i = 0; i < m; i++)  in print_lambda_matrix()
657     print_lambda_vector (outfile, matrix[i], n);  in print_lambda_matrix()