-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathMatrixproccesor.java
More file actions
340 lines (300 loc) · 13.3 KB
/
Matrixproccesor.java
File metadata and controls
340 lines (300 loc) · 13.3 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
package processor;
import java.util.Scanner;
public class Matrixproccesor {
public static void main(String[] args) {
Scanner scanner = new Scanner(System.in);
OperateMatrix obj = new OperateMatrix();
int rowNumber1,columnNumber1,rowNumber2,columnNumber2;
while(true) {
switchPrint();
System.out.print("Your choice: ");
int choice = scanner.nextInt();
switch (choice) {
case 1:
System.out.print("Enter the size of the first matrix: ");
rowNumber1 = scanner.nextInt();
columnNumber1 = scanner.nextInt();
System.out.println("Enter first matrix: ");
double[][] matrix1 = fill(rowNumber1, columnNumber1);
System.out.print("Enter the size of the second matrix: ");
rowNumber2 = scanner.nextInt();
columnNumber2 = scanner.nextInt();
System.out.println("Enter second matrix:");
double[][] matrix2 = fill(rowNumber2, columnNumber2);
double[][] addAns = obj.matrixAdd(matrix1, matrix2);
System.out.println("The addition result is:");
printMatrix(addAns);
break;
case 2:
System.out.print("Enter matrix size: ");
rowNumber1 = scanner.nextInt();
columnNumber1 = scanner.nextInt();
System.out.println("Enter matrix: ");
matrix1 = fill(rowNumber1, columnNumber1);
System.out.print("Enter the constant to be multiplied with the matrix: ");
int constant = scanner.nextInt();
double[][] answer = obj.matrixMultiply(matrix1, constant);
System.out.println("The multiplication result is:");
printMatrix(answer);
break;
case 3:
System.out.print("Enter the size of the first matrix: ");
rowNumber1 = scanner.nextInt();
columnNumber1 = scanner.nextInt();
System.out.println("Enter first matrix: ");
matrix1 = fill(rowNumber1, columnNumber1);
System.out.print("Enter the size of the second matrix: ");
rowNumber2 = scanner.nextInt();
columnNumber2 = scanner.nextInt();
System.out.println("Enter second matrix:");
matrix2 = fill(rowNumber2, columnNumber2);
double[][] answerMul = obj.multiply_two_matrices(matrix1, matrix2);
printMatrix(answerMul);
break;
case 4:
System.out.print("\n1. Main diagonal\n");
System.out.println("2. Side diagonal");
System.out.println("3. Vertical diagonal");
System.out.println("4. Horizontal line");
System.out.print("Your choice: ");
int type = scanner.nextInt();
System.out.print("Enter matrix size: ");
rowNumber1 = scanner.nextInt();
columnNumber1 = scanner.nextInt();
System.out.println("Enter matrix: ");
matrix1 = fill(rowNumber1, columnNumber1);
switch(type){
case 1:
System.out.println("The result is:");
double[][] transposed = obj.transposeMatrix(matrix1,"main_diagonal");
for (int i = 0; i < matrix1[0].length; i++) {
for (int j = 0; j < matrix1.length; j++) {
System.out.print(transposed[i][j] + " ");
}
System.out.print("\n");
}
break;
case 2:
System.out.println("The result is:");
transposed = obj.transposeMatrix(matrix1,"secondary_diagonal");
printMatrix(transposed);
break;
case 3:
System.out.println("The result is:");
transposed = obj.transposeMatrix(matrix1,"vertical");
printMatrix(transposed);
break;
case 4:
System.out.println("The result is:");
transposed = obj.transposeMatrix(matrix1,"horizontal");
printMatrix(transposed);
break;
}
break;
case 5:
System.out.print("Enter matrix size: ");
rowNumber1 = scanner.nextInt();
columnNumber1 = scanner.nextInt();
System.out.println("Enter matrix: ");
matrix1 = fill(rowNumber1, columnNumber1);
double number = Determinant.determinantOfMatrix(matrix1,matrix1.length);
System.out.println(format(number));
break;
case 6:
System.out.print("Enter matrix size: ");
rowNumber1 = scanner.nextInt();
columnNumber1 = scanner.nextInt();
System.out.println("Enter matrix: ");
matrix1 = fill(rowNumber1, columnNumber1);
System.out.println("The result is:");
number = Determinant.determinantOfMatrix(matrix1,matrix1.length);
if(number == 0){
System.out.println("Inverse cannot be found, determinant zero");
}
else{
double[][] ans = obj.matrixMultiply(Determinant.invert(matrix1),1 / number);
printMatrix(ans);
}
break;
case 0:
System.exit(0);
}
}
}
public static double[][] fill(int row, int column) {
double[][] matrix = new double[row][column];
Scanner scanner = new Scanner(System.in);
for (int i = 0; i < row; i++) {
for (int j = 0; j < column; j++) {
matrix[i][j] = Double.parseDouble(scanner.next());
}
}
return matrix;
}
public static void printMatrix(double[][] matrix) {
for (int i = 0 ; i < matrix.length; i++) {
for (int j = 0; j < matrix[i].length; j++) {
matrix[i][j] = Double.parseDouble(format(matrix[i][j]));
System.out.print(matrix[i][j] + " ");
}
System.out.print("\n");
}
}
public static String format(double digit) {
if (digit == (long) digit) {
return String.format("%d",(long)digit);
} else {
return String.format("%s" , digit);
}
}
public static void switchPrint() {
System.out.println("1. Add matrices");
System.out.println("2. Multiply matrix to a constant");
System.out.println("3. Multiply matrices");
System.out.println("4. Transpose matrix");
System.out.println("5. Calculate a determinant");
System.out.println("0. Exit");
}
}
class OperateMatrix {
//a function to add two matrices,returns the summed matrix
double[][] matrixAdd(double[][] matrix1, double[][] matrix2) {
double[][] matrixSum = new double[matrix1.length][matrix1[0].length]; //the sum matrix also has the same dimensions
for (int i = 0; i < matrix1.length; i++) { //considering only matrix1 because both matrices have same dimensions
for (int j = 0; j < matrix1[i].length; j++) { //the restraint condition,specifies column number
matrixSum[i][j] = matrix1[i][j] + matrix2[i][j];
}
}
return matrixSum;
}
// a function to just multiply a single number with a matrix
double[][] matrixMultiply(double[][] matrix, double constant) {
for (int i = 0; i < matrix.length; i++) {
for (int j = 0; j < matrix[i].length; j++) {
matrix[i][j] = matrix[i][j] * constant;
}
}
return matrix;
}
double[][] multiply_two_matrices(double[][] matrix1, double[][] matrix2) {
if (matrix1[0].length != matrix2.length) {
System.out.println("Error!");
}
double[][] multipliedAns = new double[matrix1.length][matrix2[0].length];
int i, j, k;
for (i = 0; i < matrix1.length; i++) {
for (j = 0; j < matrix2[0].length; j++) {
for (k = 0; k < matrix2.length; k++)
multipliedAns[i][j] += matrix1[i][k] * matrix2[k][j];
}
}
return multipliedAns;
}
double[][] transposeMatrix(double[][] matrix, String type) {
double[][] transposed = new double[matrix.length][matrix[0].length];
switch (type) {
case "main_diagonal":
for (int i = 0; i < matrix.length; i++) {
for (int j = 0; j < matrix[i].length; j++) {
transposed[i][j] = matrix[j][i];
}
}
return transposed;
case "secondary_diagonal":
for (int i = 0; i < (matrix.length - 1); i++) {
for (int j = 0; j < (matrix[0].length - 1) - i; j++) {
double tmp = matrix[i][j];
matrix[i][j] = matrix[(matrix[0].length - 1) - j][(matrix[0].length - 1) - i];
matrix[(matrix[0].length - 1) - j][(matrix[0].length - 1) - i] = tmp;
}
}
return matrix;
case "vertical":
for (int i = 0; i < matrix.length; i++) {
for (int j = 0; j < matrix[i].length / 2; j++) {
double temp = matrix[i][j];
matrix[i][j] = matrix[i][matrix.length - 1 - j];
matrix[i][matrix.length - 1 - j] = temp;
}
}
return matrix;
case "horizontal":
for (int i = 0; i < matrix.length / 2; i++) {
for (int j = 0; j < matrix[i].length; j++) {
double temp = matrix[i][j];
matrix[i][j] = matrix[matrix.length - 1 - i][j];
matrix[matrix.length - 1 - i][j] = temp;
}
}
return matrix;
}
return null;
}
}
class Determinant {
// Dimension of input square matrix
static final int N = 4;
// Function to get cofactor of
// mat[p][q] in temp[][]. n is
// current dimension of mat[][]
static void getCofactor(double[][] mat,
double[][] temp, int p, int q, int n) {
int i = 0, j = 0;
// Looping for each element of
// the matrix
for (int row = 0; row < n; row++) {
for (int col = 0; col < n; col++) {
// Copying into temporary matrix
// only those element which are
// not in given row and column
if (row != p && col != q) {
temp[i][j++] = mat[row][col];
// Row is filled, so increase
// row index and reset col
//index
if (j == n - 1) {
j = 0;
i++;
}
}
}
}
}
/* Recursive function for finding determinant
of matrix. n is current dimension of mat[][]. */
static double determinantOfMatrix(double[][] mat, int n) {
double D = 0; // Initialize result
// Base case : if matrix contains single
// element
if (n == 1)
return mat[0][0];
// To store cofactors
double[][] temp = new double[N][N];
// To store sign multiplier
int sign = 1;
// Iterate for each element of first row
for (int f = 0; f < n; f++) {
// Getting Cofactor of mat[0][f]
getCofactor(mat, temp, 0, f, n);
D += sign * mat[0][f]
* determinantOfMatrix(temp, n - 1);
// terms are to be added with
// alternate sign
sign = -sign;
}
return D;
}
static double[][] invert(double[][] matrix) {
OperateMatrix obj = new OperateMatrix();
matrix = obj.transposeMatrix(matrix,"main_diagonal");
double[][] cofactor = new double[matrix.length][matrix.length];
double[][] temp = new double[N][N];
for(int i = 0; i < matrix.length; i++){
for(int j = 0;j < matrix.length; j++){
getCofactor(matrix,temp,i,j,matrix.length);
cofactor[i][j] = Math.pow(-1,((i+1)+(j+1))) * determinantOfMatrix(temp, matrix.length - 1);
}
}
return cofactor;
}
}