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
|
/* Copyright (C) 2009-2011, 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"
#include <stdlib.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <string.h>
#include <errno.h>
#include <stdio.h>
#include <math.h>
int main(void)
{
const char *filedir = "dirfile";
const char *format = "dirfile/format";
int j, error, ge_error, r = 0;
gd_entry_t e;
#ifdef GD_NO_C99_API
const double a[4][2] = {{1, 29.03}, {0.3, 12.34}, {0.5, 99.55}, {1.8, 45.32}};
#else
const double _Complex a[4] = {
CMPLX(1, 29.03), CMPLX(0.3, 12.34),
CMPLX(0.5, 99.55), CMPLX(1.8, 45.32)};
#endif
DIRFILE *D;
rmdirfile();
D = gd_open(filedir, GD_RDWR | GD_CREAT | GD_VERBOSE);
gd_add_phase(D, "new", "in", 3, 0);
#ifdef GD_NO_C99_API
gd_madd_cpolynom(D, "new", "meta", 3, "in", (const double *)a);
#else
gd_madd_cpolynom(D, "new", "meta", 3, "in", a);
#endif
error = gd_error(D);
/* check */
gd_entry(D, "new/meta", &e);
ge_error = gd_error(D);
CHECKI(ge_error, 0);
if (!r) {
CHECKI(e.field_type, GD_POLYNOM_ENTRY);
CHECKS(e.in_fields[0], "in");
CHECKI(e.fragment_index, 0);
CHECKI(e.EN(polynom,poly_ord), 3);
CHECKX(e.flags, GD_EN_COMPSCAL | GD_EN_CALC);
for (j = 0; j < 4; ++j)
CHECKCi(j,e.EN(polynom,ca)[j], a[j]);
gd_free_entry_strings(&e);
}
gd_discard(D);
unlink(format);
rmdir(filedir);
CHECKI(error, GD_E_OK);
return r;
}
|
