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
|
typedef struct
{
int size1;
int size2;
int** array;
}
matrix;
matrix* matrix_alloc(int size1, int size2)
{
matrix* m = malloc(sizeof(matrix));
m->size1 = size1;
m->size2 = size2;
m->array = malloc(size1*size2*sizeof(int));
int i;
for (i = 0; i < size1; i++)
m->array[i] = malloc(size2*sizeof(int));
return m;
}
matrix* cpy_matrix(matrix* m)
{
matrix* cpy = matrix_alloc(m->size1, m->size2);
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
cpy->array[i][j] = m->array[i][j];
return cpy;
}
matrix* sin_matrix(matrix* m)
{
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
m->array[i][j] = 255*sin((double)m->array[i][j]);
return m;
}
matrix* cos_matrix(matrix* m)
{
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
m->array[i][j] = 255*cos((double)m->array[i][j]);
return m;
}
matrix* sqrt_matrix(matrix* m)
{
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
m->array[i][j] = (int)sqrt((double)m->array[i][j]);
return m;
}
void matrix_add(matrix* m1, matrix* m2)
{
int i, j;
for (i = 0; i < m1->size1; i++)
for (j = 0; j < m1->size2; j++)
m1->array[i][j] += m2->array[i][j];
}
void matrix_sub(matrix* m1, matrix* m2)
{
int i, j;
for (i = 0; i < m1->size1; i++)
for (j = 0; j < m1->size2; j++)
m1->array[i][j] -= m2->array[i][j];
}
void matrix_add_constant(matrix* m, int c)
{
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
m->array[i][j] += c;
}
void matrix_mul_elements(matrix* m1, matrix* m2)
{
int i, j;
for (i = 0; i < m1->size1; i++)
for (j = 0; j < m1->size2; j++)
m1->array[i][j] *= m2->array[i][j];
}
void matrix_scale(matrix* m, int c)
{
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
m->array[i][j] *= c;
}
matrix* mod256(matrix* m)
{
int i, j;
for (i = 0; i < m->size1; i++)
for (j = 0; j < m->size2; j++)
m->array[i][j] %= 256;
return m;
}
|