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
|
/* Copyright (C) 2013, 2017 D.V. Wiebe
*
***************************************************************************
*
* This file is part of the GetData project.
*
* GetData is free software; you can redistribute it and/or modify it under
* the terms of the GNU Lesser General Public License as published by the
* Free Software Foundation; either version 2.1 of the License, or (at your
* option) any later version.
*
* GetData 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 Lesser General Public
* License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with GetData; if not, write to the Free Software Foundation, Inc.,
* 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include "test.h"
int main(void)
{
const char *filedir = "dirfile";
const char *format = "dirfile/format";
const char *data = "dirfile/data";
int8_t c[8], d;
struct stat buf;
int fd, q, n, e1, e2, r = 0;
DIRFILE *D;
rmdirfile();
mkdir(filedir, 0700);
for (q = 0; q < 8; ++q)
c[q] = (int8_t)(39 + q * 2);
MAKEFORMATFILE(format, "lincom LINCOM 1 data 1;1 0;-1\ndata RAW INT8 8\n");
D = gd_open(filedir, GD_RDWR | GD_UNENCODED | GD_VERBOSE);
n = gd_putdata(D, "lincom", 5, 0, 1, 0, GD_INT8, c);
e1 = gd_error(D);
CHECKI(e1, 0);
CHECKI(n,8);
e2 = gd_close(D);
CHECKI(e2, 0);
if (stat(data, &buf)) {
perror("stat");
r = 1;
} else {
CHECKI(buf.st_size, 48 * sizeof(int8_t));
fd = open(data, O_RDONLY | O_BINARY);
q = 0;
while (read(fd, &d, sizeof(int8_t))) {
if (q < 40 || q > 48) {
CHECKIi(q,d,0);
} else {
/* if:
* L = (1+i) D - i
* then:
* D = (L + i) * (1 / (1+i))
* = (L + i) * ((1-i) / 2)
* = [ (L + 1) + (1 - L) i ] / 2
* so:
* Re(D) = (L + 1) / 2
* and, if:
* L = 39 + 2 (q - 40),
* then:
* Re(D) = (40 + 2(q-40)) / 2 = 20 + q - 40 = q - 20
*/
CHECKIi(q,d,q - 20);
}
q++;
}
close(fd);
}
unlink(data);
unlink(format);
rmdir(filedir);
return r;
}
|
