1 /* Lambda matrix transformations.
2 Copyright (C) 2003, 2004 Free Software Foundation, Inc.
3 Contributed by Daniel Berlin <dberlin@dberlin.org>.
4
5 This file is part of GCC.
6
7 GCC is free software; you can redistribute it and/or modify it under
8 the terms of the GNU General Public License as published by the Free
9 Software Foundation; either version 2, or (at your option) any later
10 version.
11
12 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
13 WARRANTY; without even the implied warranty of MERCHANTABILITY or
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
15 for more details.
16
17 You should have received a copy of the GNU General Public License
18 along with GCC; see the file COPYING. If not, write to the Free
19 Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA
20 02110-1301, USA. */
21
22 #include "config.h"
23 #include "system.h"
24 #include "coretypes.h"
25 #include "tm.h"
26 #include "ggc.h"
27 #include "tree.h"
28 #include "target.h"
29 #include "varray.h"
30 #include "lambda.h"
31
32 /* Allocate a new transformation matrix. */
33
34 lambda_trans_matrix
lambda_trans_matrix_new(int colsize,int rowsize)35 lambda_trans_matrix_new (int colsize, int rowsize)
36 {
37 lambda_trans_matrix ret;
38
39 ret = ggc_alloc (sizeof (*ret));
40 LTM_MATRIX (ret) = lambda_matrix_new (rowsize, colsize);
41 LTM_ROWSIZE (ret) = rowsize;
42 LTM_COLSIZE (ret) = colsize;
43 LTM_DENOMINATOR (ret) = 1;
44 return ret;
45 }
46
47 /* Return true if MAT is an identity matrix. */
48
49 bool
lambda_trans_matrix_id_p(lambda_trans_matrix mat)50 lambda_trans_matrix_id_p (lambda_trans_matrix mat)
51 {
52 if (LTM_ROWSIZE (mat) != LTM_COLSIZE (mat))
53 return false;
54 return lambda_matrix_id_p (LTM_MATRIX (mat), LTM_ROWSIZE (mat));
55 }
56
57
58 /* Compute the inverse of the transformation matrix MAT. */
59
60 lambda_trans_matrix
lambda_trans_matrix_inverse(lambda_trans_matrix mat)61 lambda_trans_matrix_inverse (lambda_trans_matrix mat)
62 {
63 lambda_trans_matrix inverse;
64 int determinant;
65
66 inverse = lambda_trans_matrix_new (LTM_ROWSIZE (mat), LTM_COLSIZE (mat));
67 determinant = lambda_matrix_inverse (LTM_MATRIX (mat), LTM_MATRIX (inverse),
68 LTM_ROWSIZE (mat));
69 LTM_DENOMINATOR (inverse) = determinant;
70 return inverse;
71 }
72
73
74 /* Print out a transformation matrix. */
75
76 void
print_lambda_trans_matrix(FILE * outfile,lambda_trans_matrix mat)77 print_lambda_trans_matrix (FILE *outfile, lambda_trans_matrix mat)
78 {
79 print_lambda_matrix (outfile, LTM_MATRIX (mat), LTM_ROWSIZE (mat),
80 LTM_COLSIZE (mat));
81 }
82