-
Notifications
You must be signed in to change notification settings - Fork 6
/
hypisom.c
38 lines (31 loc) · 1011 Bytes
/
hypisom.c
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
#error THIS FILE IS NOT USED
#include"typedef.h"
main (int argc, char *argv[])
{
matrix_TYP **x1, **x2, **S, *X;
int anz;
read_header(argc, argv);
if(FILEANZ != 3)
{
printf("usage: hypstab 'file1' 'file2' 'file3' ,\n");
printf("where 'file1' contains an 1 by n matrix x.\n");
printf("where 'file2' contains an 1 by n matrix y.\n");
printf("and 'file3' contains a n by n matrix S.\n");
printf("of sylvester Typ (-1, n-1) s.t. xSx^{tr} < 0 and ySy^{tr} < 0\n");
if (is_option('h')){
exit(0);
}
else{
exit(31);
}
}
x1 = mget_mat(FILENAMES[0], &anz);
x2 = mget_mat(FILENAMES[1], &anz);
S = mget_mat(FILENAMES[2], &anz);
X = hyperbolic_isometry(x1[0], x2[0], S[0]);
if(X == NULL)
printf("There exists no isometry\n");
else
put_mat(X, NULL, "isometry", 0);
exit(3);
}