test/tf.h

   1 #define MAIN_(x) A_CAST_2(x, _tf)
   2 #include "test.h"
   3 #include "a/tf.h"
   4
   5 static a_float u[] = {
   6     A_FLOAT_C(2000.0),
   7     A_FLOAT_C(-652.99418),
   8     A_FLOAT_C(-344.66975),
   9     A_FLOAT_C(-168.84826),
  10     A_FLOAT_C(-133.5101),
  11     A_FLOAT_C(-109.5296),
  12     A_FLOAT_C(-86.454203),
  13     A_FLOAT_C(-67.783481),
  14     A_FLOAT_C(-53.22429),
  15     A_FLOAT_C(-41.808457),
  16     A_FLOAT_C(-32.835939),
  17     A_FLOAT_C(-25.786913),
  18     A_FLOAT_C(-20.2509),
  19     A_FLOAT_C(-15.903432),
  20     A_FLOAT_C(-12.489385),
  21     A_FLOAT_C(-9.8083792),
  22     A_FLOAT_C(-7.703033),
  23     A_FLOAT_C(-6.0497455),
  24     A_FLOAT_C(-4.7514505),
  25     A_FLOAT_C(-3.731923),
  26     A_FLOAT_C(-2.9313047),
  27     A_FLOAT_C(-2.3025903),
  28     A_FLOAT_C(-1.8088678),
  29     A_FLOAT_C(-1.4211508),
  30     A_FLOAT_C(-1.1166775),
  31     A_FLOAT_C(-0.87757332),
  32     A_FLOAT_C(-0.68980192),
  33     A_FLOAT_C(-0.54234094),
  34     A_FLOAT_C(-0.42653479),
  35     A_FLOAT_C(-0.33558647),
  36     A_FLOAT_C(-0.26415848),
  37     A_FLOAT_C(-0.2080594),
  38     A_FLOAT_C(-0.16399783),
  39     A_FLOAT_C(-0.12938911),
  40     A_FLOAT_C(-0.10220358),
  41     A_FLOAT_C(-0.080847398),
  42     A_FLOAT_C(-0.064068962),
  43     A_FLOAT_C(-0.050885415),
  44     A_FLOAT_C(-0.040524947),
  45     A_FLOAT_C(-0.032381475),
  46     A_FLOAT_C(-0.025979056),
  47     A_FLOAT_C(-0.020943941),
  48     A_FLOAT_C(-0.016982637),
  49     A_FLOAT_C(-0.013864667),
  50     A_FLOAT_C(-0.011409044),
  51     A_FLOAT_C(-0.0094736412),
  52     A_FLOAT_C(-0.00794685),
  53     A_FLOAT_C(-0.0067410251),
  54     A_FLOAT_C(-0.0057873396),
  55     A_FLOAT_C(-0.0050317445),
  56     A_FLOAT_C(-0.0044317954),
  57     A_FLOAT_C(-0.0039541603),
  58     A_FLOAT_C(-0.0035726634),
  59     A_FLOAT_C(-0.0032667484),
  60     A_FLOAT_C(-0.0030202717),
  61     A_FLOAT_C(-0.002820555),
  62     A_FLOAT_C(-0.0026576412),
  63     A_FLOAT_C(-0.0025237103),
  64     A_FLOAT_C(-0.00241262),
  65     A_FLOAT_C(-0.0023195462),
  66     A_FLOAT_C(-0.0022406992),
  67     A_FLOAT_C(-0.0021731023),
  68     A_FLOAT_C(-0.0021144168),
  69     A_FLOAT_C(-0.002062805),
  70     A_FLOAT_C(-0.0020168229),
  71     A_FLOAT_C(-0.0019753356),
  72     A_FLOAT_C(-0.0019374509),
  73     A_FLOAT_C(-0.0019024669),
  74     A_FLOAT_C(-0.0018698319),
  75     A_FLOAT_C(-0.0018391113),
  76     A_FLOAT_C(-0.0018099631),
  77     A_FLOAT_C(-0.0017821178),
  78     A_FLOAT_C(-0.0017553627),
  79     A_FLOAT_C(-0.00172953),
  80     A_FLOAT_C(-0.001704487),
  81     A_FLOAT_C(-0.0016801285),
  82     A_FLOAT_C(-0.0016563712),
  83     A_FLOAT_C(-0.0016331487),
  84     A_FLOAT_C(-0.0016104081),
  85     A_FLOAT_C(-0.0015881069),
  86     A_FLOAT_C(-0.0015662109),
  87     A_FLOAT_C(-0.0015446926),
  88     A_FLOAT_C(-0.0015235294),
  89     A_FLOAT_C(-0.0015027029),
  90     A_FLOAT_C(-0.0014821978),
  91     A_FLOAT_C(-0.0014620014),
  92     A_FLOAT_C(-0.0014421028),
  93     A_FLOAT_C(-0.0014224928),
  94     A_FLOAT_C(-0.0014031634),
  95     A_FLOAT_C(-0.0013841076),
  96     A_FLOAT_C(-0.0013653192),
  97     A_FLOAT_C(-0.0013467927),
  98     A_FLOAT_C(-0.0013285228),
  99     A_FLOAT_C(-0.001310505),
 100     A_FLOAT_C(-0.0012927348),
 101     A_FLOAT_C(-0.0012752081),
 102     A_FLOAT_C(-0.001257921),
 103     A_FLOAT_C(-0.0012408699),
 104     A_FLOAT_C(-0.0012240512),
 105     A_FLOAT_C(-0.0012074614),
 106 };
 107
 108 int MAIN(int argc, char *argv[]) // NOLINT(misc-definitions-in-headers)
 109 {
 110     main_init(argc, argv, 1);
 111
 112     a_float num[] = {A_FLOAT_C(6.59492796e-05), A_FLOAT_C(6.54019884e-05)};
 113     a_float den[] = {A_FLOAT_C(-1.97530991), A_FLOAT_C(0.97530991)};
 114
 115     a_tf ctx;
 116     a_float input[A_LEN(num)];
 117     a_float output[A_LEN(den)];
 118     a_tf_init(&ctx, A_LEN(num), num, input, A_LEN(den), den, output);
 119
 120     for (unsigned int i = 0; i < A_LEN(u); ++i)
 121     {
 122         a_tf_iter(&ctx, u[i]);
 123         debug(A_FLOAT_PRI("+", "f,") A_FLOAT_PRI("+", "f,") A_FLOAT_PRI("+", "f\n"),
 124               A_FLOAT_C(0.001) * a_float_c(i), A_FLOAT_C(1.0), *ctx.output);
 125     }
 126
 127 #if defined(__cplusplus) && (__cplusplus > 201100L)
 128     A_ASSERT_BUILD(std::is_pod<a_tf>::value);
 129 #endif /* __cplusplus */
 130
 131     return 0;
 132 }