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
|
/* IGraph library.
Copyright (C) 2022 The igraph development team <igraph@igraph.org>
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <https://www.gnu.org/licenses/>.
*/
#include <igraph.h>
#include "core/math.h"
#include "test_utilities.h"
int main(void) {
igraph_matrix_complex_t c;
igraph_matrix_t real;
igraph_matrix_t imag;
int e[] = {1, 2, 3, 4};
int e2[] = {5, 6, 7, 8};
matrix_init_int_row_major(&real, 2, 2, e);
matrix_init_int_row_major(&imag, 2, 2, e2);
printf("Complex matrix:\n");
igraph_matrix_complex_create(&c, &real, &imag);
igraph_matrix_complex_fprint(&c, stdout);
printf("Real part:\n");
igraph_matrix_resize(&real, 0, 0);
igraph_matrix_complex_real(&c, &real);
igraph_matrix_print(&real);
printf("Imaginary part:\n");
igraph_matrix_resize(&imag, 0, 0);
igraph_matrix_complex_imag(&c, &imag);
igraph_matrix_print(&imag);
printf("Real and imaginary part:\n");
igraph_matrix_resize(&real, 0, 0);
igraph_matrix_resize(&imag, 0, 0);
igraph_matrix_complex_realimag(&c, &real, &imag);
igraph_matrix_print(&real);
igraph_matrix_print(&imag);
igraph_matrix_complex_destroy(&c);
igraph_matrix_destroy(&real);
igraph_matrix_destroy(&imag);
{
printf("Complex matrix from polar:\n");
igraph_matrix_complex_t p;
igraph_real_t r_e[] = {1, 2, 3, 4};
igraph_real_t theta_e[] = {0, .5 * M_PI, M_PI, 1.5 * M_PI};
igraph_matrix_t r;
igraph_matrix_t theta;
igraph_matrix_view(&r, r_e, 2, 2);
igraph_matrix_view(&theta, theta_e, 2, 2);
igraph_matrix_complex_create_polar(&p, &r, &theta);
print_matrix_complex_round(&p);
igraph_matrix_complex_destroy(&p);
}
VERIFY_FINALLY_STACK();
{
igraph_real_t e3[] = {1, 2, 3};
igraph_real_t e4[] = {5, 6, 7, 8};
printf("Check if unequal number of imaginary and real rows is handled correctly.\n");
igraph_matrix_view(&real, e3, 1, 2);
igraph_matrix_view(&imag, e4, 2, 2);
CHECK_ERROR(igraph_matrix_complex_create(&c, &real, &imag), IGRAPH_EINVAL);
CHECK_ERROR(igraph_matrix_complex_create_polar(&c, &real, &imag), IGRAPH_EINVAL);
}
{
igraph_real_t e3[] = {1, 2, 3};
igraph_real_t e4[] = {5, 6, 7, 8};
printf("Check if unequal number of imaginary and real columns is handled correctly.\n");
igraph_matrix_view(&real, e3, 2, 1);
igraph_matrix_view(&imag, e4, 2, 2);
CHECK_ERROR(igraph_matrix_complex_create(&c, &real, &imag), IGRAPH_EINVAL);
CHECK_ERROR(igraph_matrix_complex_create_polar(&c, &real, &imag), IGRAPH_EINVAL);
}
VERIFY_FINALLY_STACK();
return 0;
}
|
