summaryrefslogtreecommitdiff
path: root/ce
diff options
context:
space:
mode:
authorlacour2008-02-15 13:26:11 +0000
committerlacour2008-02-15 13:26:11 +0000
commit1a1e87e2d181ce8615e4703f690078277feb3d1f (patch)
tree4202e70e2e5d1efb2be5d7d7093acaa550a5e506 /ce
parent672cd423868bf7a7a12d2d9b275304a9b469494c (diff)
Remove test on sparc from ce/test/rx/general.
Add spoc project in CP. Add spoc coefficient check for matrix M. git-svn-id: svn+ssh://pessac/svn/cesar/trunk@1419 017c9cb6-072f-447c-8318-d5b54f68fe89
Diffstat (limited to 'ce')
-rwxr-xr-xce/test/rx/general/Makefile_HOST_LINUX_SPARC25
-rw-r--r--ce/test/rx/general/src/test_speed.c193
-rw-r--r--ce/test/rx/general/src/test_spoc.c97
-rw-r--r--ce/test/rx/general/src/test_spoc_opt.c346
4 files changed, 0 insertions, 661 deletions
diff --git a/ce/test/rx/general/Makefile_HOST_LINUX_SPARC b/ce/test/rx/general/Makefile_HOST_LINUX_SPARC
deleted file mode 100755
index 8af6f9fec5..0000000000
--- a/ce/test/rx/general/Makefile_HOST_LINUX_SPARC
+++ /dev/null
@@ -1,25 +0,0 @@
-BASE = ../../..
-
-CC=sparc-linux-gcc
-INCLUDES=ce/test/rx
-
-#EXTRA_HOST_CFLAGS=-O3
-EXTRA_HOST_LDFLAGS=-lm
-HOST_PROGRAMS = test_speed_host_linux_sparc
-test_speed_host_linux_sparc_SOURCES = test_speed.c
-test_speed_host_linux_sparc_MODULES = lib mac/common ce ce/test/common
-
-HOST_PROGRAMS += test_spoc_host_linux_sparc
-test_spoc_host_linux_sparc_SOURCES = test_spoc.c
-HOST_PROGRAMS += test_spoc_opt_host_linux_sparc
-test_spoc_opt_host_linux_sparc_SOURCES = test_spoc_opt.c
-
-#HOST_PROGRAMS += test_cei_host_linux_sparc
-#test_cei_host_linux_sparc_SOURCES = test_cei.c
-#test_cei_host_linux_sparc_MODULES = lib mac/common ce ce/test/common
-
-ce_MODULE_SOURCES = cei_param.c frame_measurement.c
-ce_test_common_MODULE_SOURCES = print_utils.c
-
-include $(BASE)/common/make/top.mk
-
diff --git a/ce/test/rx/general/src/test_speed.c b/ce/test/rx/general/src/test_speed.c
deleted file mode 100644
index 210b082cc3..0000000000
--- a/ce/test/rx/general/src/test_speed.c
+++ /dev/null
@@ -1,193 +0,0 @@
-/* Cesar project {{{
- *
- * Copyright (C) 2007 Spidcom
- *
- * <<<Licence>>>
- *
- * }}} */
-/**
- * \file ./src/test_speed.c
- * \brief « brief description »
- * \ingroup « module »
- *
- * « long description »
- */
-
-#include "common/std.h"
-#include "sys/time.h"
-#include "stdio.h"
-#include "string.h"
-#include "ce/inc/frame_measurement.h"
-#include "mac/common/tonemap.h"
-#include "ce/test/common/print_utils.h"
-#include "lib/bitstream.h"
-#define RXCE
-#include "ce/inc/cei.h"
-#include "math.h"
-
-
-int
-main (int argc, char **argv)
-{
- struct timeval tv0, tv1;
- double test_d;
- int r;
- double flottant = 0.123456789;
- double res = flottant;
- double x=0;
-
- gettimeofday (&tv0, NULL);
- for (r=-10000; r<10000; r++)
- {
- x = (double)(M_PI*r)/ 10000.0;
- for (test_d=0; test_d<100; test_d++)
- {
- res = test_d * x;
- }
- }
- gettimeofday (&tv1, NULL);
- printf ("test duration [nothing] %d \n", (int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
- gettimeofday (&tv0, NULL);
- for (r=-10000; r<10000; r++)
- {
- x = (double)(M_PI*r)/ 10000.0;
- for (test_d=0; test_d<100; test_d++)
- {
- res = test_d / x;
- }
- }
- gettimeofday (&tv1, NULL);
- printf ("test duration [div : dble/dble] %d \n", (int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
- gettimeofday (&tv0, NULL);
- for (r=-10000; r<10000; r++)
- {
- x = (double)(M_PI*r)/ 10000.0;
- for (test_d=0; test_d<100; test_d++)
- {
- res = sin (x);
- }
- }
- gettimeofday (&tv1, NULL);
- printf ("test duration [sin (dble)] %d \n", (int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
-
- uint word = 0x76543210;
- uint i, test;
- int nibble;
- bitstream_t stream_writer;
- gettimeofday (&tv0, NULL);
- gettimeofday (&tv1, NULL);
- printf ("getitmeofday duration %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
- gettimeofday (&tv0, NULL);
- for (test=0; test<1000; test++)
- {
- word = 0;
- for (i=0; i<8; i++)
- {
- nibble = i+1;
- word = (word >> 4) | (nibble << 28);
- }
- }
- gettimeofday (&tv1, NULL);
- printf ("nibble to u32 duration test1 %d /1000\n",(int) (tv1.tv_usec - tv0.tv_usec));
-
- gettimeofday (&tv0, NULL);
- for (test=0; test<1000; test++)
- {
- word = 0;
- for (i=0; i<8; i++)
- {
- nibble = i+1;
- word = word | (nibble << (i<<2));
- }
- }
- gettimeofday (&tv1, NULL);
- printf ("nibble to u32 duration test2 %d /1000\n",(int) (tv1.tv_usec - tv0.tv_usec));
-
- frame_measurement_init ();
-
- uint noise = 0;
- gettimeofday (&tv0, NULL);
- update_threshold (1, 1);
- for (i=0; i<1536; i++)
- {
- noise2mod (noise);
- }
- gettimeofday (&tv1, NULL);
- printf ("noise2mod duration --worst : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
-
- noise = 10000;
- gettimeofday (&tv0, NULL);
- update_threshold (1, 1);
- for (i=0; i<1536; i++)
- {
- noise2mod (noise);
-
- }
- gettimeofday (&tv1, NULL);
- printf ("noise2mod duration --best : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
-
-
- gettimeofday (&tv0, NULL);
- blk_t *blk = blk_alloc_desc ();
- gettimeofday (&tv1, NULL);
- printf ("first blk allocation duration : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
- blk_release_desc (blk);
- gettimeofday (&tv0, NULL);
- blk = blk_alloc_desc ();
- gettimeofday (&tv1, NULL);
- printf ("second blk allocation duration : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
- blk_release_desc (blk);
-
- blk_t *blk_last;
- gettimeofday (&tv0, NULL);
- blk_t *blk_first = blk_alloc_desc_range (10, &blk_last);
- gettimeofday (&tv1, NULL);
- printf ("multi-blk allocation duration (10) : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
- blk_release_desc_range (blk_first, blk_last);
-
-
- phy_chandata_t *last, *browser;
- phy_chandata_t *first = (phy_chandata_t *) blk_alloc_desc_range (PHY_CHANDATA_NRJ_BLK_NB, (blk_t **) &last);
- last->last = 1;
- last->blk.next = NULL;
- browser = first;
- noise = 0;
- int j, k=0;
- u32 * data = (u32 * ) browser->blk.data;
- for (i=0; i<PHY_CHANDATA_NRJ_BLK_NB; i++)
- {
- for (j=0; j<BLK_SIZE/4; j++)
- {
- data[j] = (THR3[1][k%7]) | ((THR3[1][(k+1)%7])<<16);
- k+=2;
- }
- browser = (phy_chandata_t *) (browser->blk.next);
- if (browser) data = (u32 * ) browser->blk.data;
- }
- tonemap_t *worst_tonemap = NULL;
- gettimeofday (&tv0, NULL);
- compute_worst_tonemap (&worst_tonemap, first);
- gettimeofday (&tv1, NULL);
- printf ("default tonemap computation duration : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
- ce_print_tonemap (worst_tonemap);
- blk_release_desc_range ((blk_t *) first, (blk_t *)last);
-
- u8 buf[1536];
- for (i=0; i<1536; i++) buf[i] = 0x00;
- u8 mask[(PHY_CARRIER_NB+7)/8];
- for (i=0; i<(PHY_CARRIER_NB+7)/8; i++) mask[i] = 0xFF;
- mask[i-1] = ((1<<(PHY_CARRIER_NB%8))-1);
- uint entry_nb;
- gettimeofday (&tv0, NULL);
- for (test=0; test<100; test++)
- {
- bitstream_init (&stream_writer, buf, 1518, BITSTREAM_WRITE);
- entry_nb = cei_encode_tm (&stream_writer, worst_tonemap->tmdma_desc_head, mask);
- }
- gettimeofday (&tv1, NULL);
- printf ("cei building duration : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
-
- tonemap_free (worst_tonemap);
- return 0;
-}
-
diff --git a/ce/test/rx/general/src/test_spoc.c b/ce/test/rx/general/src/test_spoc.c
deleted file mode 100644
index 44790ec09e..0000000000
--- a/ce/test/rx/general/src/test_spoc.c
+++ /dev/null
@@ -1,97 +0,0 @@
-/* Cesar project {{{
- *
- * Copyright (C) 2007 Spidcom
- *
- * <<<Licence>>>
- *
- * }}} */
-/**
- * \file ./src/test_speed.c
- * \brief « brief description »
- * \ingroup « module »
- *
- * « long description »
- */
-
-#include "common/std.h"
-#include "sys/time.h"
-#include "math.h"
-#include "stdio.h"
-
-#define Pmin 74
-#define N0 1155
-#define Pmax (Pmin + N0 - 1)
-
-double reg_alpha[41];
-double reg_beta[41];
-double reg_bloc2[21];
-double reg_a0, reg_c0;
-
-void
-spoc_update (double rho)
-{
- int k;
- double pirho = M_PI*rho;
- double x = pirho * Pmin;
- double xp = pirho * Pmax;
- double s, sp;
- double den = (xp - x);
-
- double x0 = 2*M_PI*1536*rho;
-
- for (k=-10; k<=10; k++)
- {
- reg_bloc2[k+10] = sin(x0+k*M_PI)/(x0+k*M_PI);
- }
-
- double tmp1 = sin (x)/x;
- double tmp2 = x*x;
- reg_a0 = (sin(xp)/xp - tmp1)/(xp*xp - tmp2);
- reg_c0 = tmp1 - reg_a0*tmp2;
-
- for (k=1; k<=20; k++)
- {
- x += (M_PI + pirho);
- xp += M_PI;
- s = sin (x) /x;
- sp = sin (xp) /xp;
- den -= pirho;
- reg_beta[k+20] = (sp-s)/(den);
- reg_alpha[k+20] = s - reg_beta[k+20]*x;
- }
- x = pirho * Pmin;
- xp = pirho * Pmax;
- den = (xp -x);
- for (k=1; k<=20; k++)
- {
- xp -= (M_PI + pirho);
- x -= M_PI;
- s = sin (x) /x;
- sp = sin (xp) /xp;
- den -= pirho;
- reg_beta[20-k] = (sp-s)/(den);
- reg_alpha[20-k] = s - reg_beta[20-k]*x;
- }
-}
-
-
-int
-main (int argc, char **argv)
-{
- struct timeval tv0, tv1;
- double rho = 300e-6;
- int t;
-
- gettimeofday (&tv0, NULL);
- for (t=0; t<1; t++) spoc_update (rho);
- gettimeofday (&tv1, NULL);
- printf ("duration spoc2 : %d\n",(int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
-
- gettimeofday (&tv0, NULL);
- for (t=0; t<1000; t++) spoc_update (rho);
- gettimeofday (&tv1, NULL);
- printf ("duration 1000*spoc2 : %d\n",(int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
-
- return 0;
-}
-
diff --git a/ce/test/rx/general/src/test_spoc_opt.c b/ce/test/rx/general/src/test_spoc_opt.c
deleted file mode 100644
index 31bef0ef87..0000000000
--- a/ce/test/rx/general/src/test_spoc_opt.c
+++ /dev/null
@@ -1,346 +0,0 @@
-/* Cesar project {{{
- *
- * Copyright (C) 2007 Spidcom
- *
- * <<<Licence>>>
- *
- * }}} */
-
-
-#include "common/std.h"
-#include "sys/time.h"
-#include "math.h"
-#include "stdio.h"
-#include <asm/leon.h>
-#include "ce/inc/frame_measurement.h"
-#include "mac/common/tonemap.h"
-
-//#define SINC_STORAGE
-//#define COMPARE
-
-double kaiser[10] = { 0.9905, 0.9622, 0.9164, 0.8546, 0.7792, 0.6929, 0.5989, 0.5005, 0.4011, 0.304 };
-
-#define Pmin 74
-#define N0 1155
-#define Pmax (Pmin + N0 - 1)
-
-//typedef int spoc_prec_t;
-typedef float spoc_prec_t;
-//typedef double spoc_prec_t; // en flottant sinx/x avec x=0.1ppm retourne 1
-
-spoc_prec_t reg_alpha[41];
-spoc_prec_t reg_beta[41];
-spoc_prec_t reg_a0, reg_c0;
-
-int alpha[3]={372, 1140, -3072};
-spoc_prec_t reg_bloc2[21][3];
-
-#ifdef COMPARE
-double reg_alpha_2[41];
-double reg_beta_2[41];
-double reg_a0_2, reg_c0_2;
-double reg_bloc2_2[21][3];
-#endif
-
-#ifdef SINC_STORAGE
-spoc_prec_t tmp_sinc[41];
-spoc_prec_t tmp_sinc_2[41];
-spoc_prec_t tmp_sincp[41];
-spoc_prec_t tmp_sincp_2[41];
-#endif
-
-inline void
-spoc_bloc2 (int n, spoc_prec_t pirho)
-{
- spoc_prec_t xp = alpha[n]*pirho;
- spoc_prec_t xm = xp;
- spoc_prec_t s = sin (xp);
- int sign = -1;
- int k;
-
- for (k=1; k<=10; k++)
- {
- xp += M_PI;
- xm -= M_PI;
- reg_bloc2[k+10][n] = kaiser[k-1] * sign*s/xp;
- reg_bloc2[10-k][n] = kaiser[k-1] * sign*s/xm;
- sign = -sign;
- }
-}
-
-inline void
-spoc_update (spoc_prec_t rho)
-{
- int k;
- spoc_prec_t x;
- spoc_prec_t xp;
- spoc_prec_t s;
- spoc_prec_t pirho;
- pirho = M_PI * rho;
-
- spoc_prec_t abs_rho = rho;
- if (abs_rho < 0 ) abs_rho = -abs_rho;
- if (abs_rho>=0.01*1e-6)
- {
- spoc_prec_t spr;
- spoc_prec_t cpr;
- spoc_prec_t sp;
- spoc_prec_t c;
- spoc_prec_t zs;
- spoc_prec_t den;
- spoc_prec_t den0, sin0, sinp0, xp0, x0;
-
-
- spr = sin (pirho);
- cpr = cos (pirho);
- x0 = x = pirho * Pmin;
- xp0 = xp = pirho * Pmax;
- //s = sin(x);
- sinp0 = sp = sin(xp0);
- c = cos(x);
- sin0 = zs = sin(x);
- den0 = den = (xp - x);
-
- reg_bloc2[10][0] = sin (alpha[0]*pirho)/(alpha[0]*pirho);
- reg_bloc2[10][1] = sin (alpha[1]*pirho)/(alpha[1]*pirho);
- reg_bloc2[10][2] = sin (alpha[2]*pirho)/(alpha[2]*pirho);
-
- //spoc_prec_t tmp1 = zs/x;
- //spoc_prec_t tmp2 = x*x;
- //reg_a0 = (sp/xp - tmp1)/(xp*xp - tmp2);
- //reg_c0 = tmp1 - reg_a0*tmp2;
- reg_a0 = (sp/xp - zs/x)/(xp*xp - x*x);
- //reg_a0 = rho;
- reg_c0 = zs/x - reg_a0*x*x;
-
- for (k=1; k<=20; k++)
- {
- x += (M_PI + pirho);
- xp += M_PI;
- s = -zs * cpr - c * spr;
- c = zs * spr - c * cpr;
- sp = -sp;
- den -= pirho;
-
- zs = s/x;
-#ifdef SINC_STORAGE
- tmp_sinc[k+20] = zs;
- tmp_sincp[k+20] = sp/xp;
-#endif
- reg_beta[k+20] = (sp/xp - zs)/den;
- reg_alpha[k+20] = zs - reg_beta[k+20]*x;
-
- zs = s;
- }
- x = x0;
- xp = xp0;
- s = sin0;
- c = cos(xp);
- zs = sinp0;
- den = den0;
- for (k=1; k<=20; k++)
- {
- xp -= (M_PI + pirho);
- x -= M_PI;
- sp = -zs * cpr + c * spr;
- c = -zs * spr - c * cpr;
- s = -s;
- den -= pirho;
-
- zs = sp/xp;
-#ifdef SINC_STORAGE
- tmp_sinc[20-k] = s/x;
- tmp_sincp[20-k] = zs;
-#endif
- //printf ("%05.10g-%05.10g=%05.10g\n", zs, s/x, zs - s/x);
- reg_beta[20-k] = (zs - s/x)/den;
- reg_alpha[20-k] = zs - reg_beta[20-k]*xp;
-
- zs = sp;
- }
- }
- else
- {
- reg_bloc2[10][0] = 1-((alpha[0]*pirho)*(alpha[0]*pirho))/6;
- reg_bloc2[10][1] = 1-((alpha[1]*pirho)*(alpha[1]*pirho))/6;
- reg_bloc2[10][2] = 1-((alpha[2]*pirho)*(alpha[2]*pirho))/6;
- reg_a0 = -1.0/6;
- reg_c0 = 1;
- int sign = 1;
- x = 0;
- for (k=1; k<=20; k++)
- {
- x += M_PI;
- reg_beta[k+20] = -sign/x;
- reg_alpha[k+20] = sign;
- reg_beta[20-k] = sign/x;
- reg_alpha[20-k] = sign;
- sign = -sign;
- }
- }
- /** SPOC BLOC_2 */
- spoc_bloc2 (0, pirho);
- spoc_bloc2 (1, pirho);
- spoc_bloc2 (2, pirho);
- }
-
-#ifdef COMPARE
-void
-spoc_update2 (double rho)
-{
- int k,i;
- double pirho = M_PI*rho;
- double x = pirho * Pmin;
- double xp = pirho * Pmax;
- double s, sp;
- double den = (xp - x);
-
- for (i=0; i<3; i++)
- {
- double x0 = alpha[i]*pirho;
- for (k=-10; k<=10; k++)
- {
- double kk = 1;
- if (k<0) kk = kaiser[-k-1];
- if (k>0) kk = kaiser[k-1];
- reg_bloc2_2[k+10][i] = kk * sin(x0+k*M_PI)/(x0+k*M_PI);
- }
- }
-
- double tmp1 = sin (x)/x;
- double tmp2 = x*x;
- reg_a0_2 = (sin(xp)/xp - tmp1)/(xp*xp - tmp2);
- reg_c0_2 = tmp1 - reg_a0_2*tmp2;
-
- for (k=1; k<=20; k++)
- {
- x += (M_PI + pirho);
- xp += M_PI;
- s = sin (x) /x;
- sp = sin (xp) /xp;
- den -= pirho;
-#ifdef SINC_STORAGE
- tmp_sinc_2[k+20] = s;
- tmp_sincp_2[k+20] = sp;
-#endif
- reg_beta_2[k+20] = (sp-s)/(den);
- reg_alpha_2[k+20] = s - reg_beta_2[k+20]*x;
- }
- x = pirho * Pmin;
- xp = pirho * Pmax;
- den = (xp -x);
- for (k=1; k<=20; k++)
- {
- xp -= (M_PI + pirho);
- x -= M_PI;
- s = sin (x) /x;
- sp = sin (xp) /xp;
-#ifdef SINC_STORAGE
- tmp_sinc_2[20-k] = s;
- tmp_sincp_2[20-k] = sp;
-#endif
- den -= pirho;
- reg_beta_2[20-k] = (sp-s)/(den);
- reg_alpha_2[20-k] = s - reg_beta_2[20-k]*x;
- }
-}
-#endif
-
-void
-get_error (double a, double b, double *error)
-{
- if (a > b)
- {
- if ((a-b)>*error) *error = (a-b);
- }
- else
- {
- if ((b-a)>*error) *error = (b-a);
- }
-}
-
-int
-main (int argc, char **argv)
-{
- struct timeval tv0, tv1;
- int test_rho = 0;
- int t;
- //int regv = sparc_leon3_get_cache_ctrl_reg();
- //printf("cache ctrl reg = 0x%x\n", regv);
-
- gettimeofday (&tv0, NULL);
- gettimeofday (&tv1, NULL);
- //printf ("first blk allocation duration : %d\n",(int) (tv1.tv_usec - tv0.tv_usec));
-
- double error = 0;
-
- for (test_rho=-30000; test_rho <= 30000; test_rho++)
- //for (test_rho=-500; test_rho < 500; test_rho++)
- //for (test_rho=22540; test_rho < 22541; test_rho++)
- {
- double rho = test_rho * 1e-8;
-
- //gettimeofday (&tv0, NULL);
- //for (t=0; t<1; t++)
- // sin (rho)/rho;
- //gettimeofday (&tv1, NULL);
- //printf ("duration first sin : %d\n",(int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
- //get_error (0,0, &error);
- //gettimeofday (&tv0, NULL);
- //for (t=0; t<1; t++)
- spoc_update (rho);
- //gettimeofday (&tv1, NULL);
- //printf ("duration 1000 spoc1 : %d\n",(int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
-
- //rho = 00.0*1e-6;
- //gettimeofday (&tv0, NULL);
- //for (t=0; t<1; t++)
- // spoc_update (rho);
- //gettimeofday (&tv1, NULL);
- //printf ("duration 1*spoc1 : %d\n",(int) ((tv1.tv_sec*1000000+tv1.tv_usec)-(tv0.tv_sec*1000000 + tv0.tv_usec)));
-#ifdef COMPARE
- //for (t=0; t<1; t++)
- spoc_update2 (rho);
-
- int k;
- for (k=-10; k<=10; k++)
- {
- //printf ("bloc2[%d] = %05.10g <==> %05.10g, diff=%05.10g\n", k, reg_bloc2[k+10][0], reg_bloc2_2[k+10][0], reg_bloc2[k+10][0] - reg_bloc2_2[k+10][0]);
- get_error (reg_bloc2[k+10][0], reg_bloc2_2[k+10][0], &error);
- //printf ("bloc2[%d] = %05.10g <==> %05.10g, diff=%05.10g\n", k, reg_bloc2[k+10][1], reg_bloc2_2[k+10][1], reg_bloc2[k+10][1] - reg_bloc2_2[k+10][1]);
- get_error (reg_bloc2[k+10][1], reg_bloc2_2[k+10][1], &error);
- //printf ("bloc2[%d] = %05.10g <==> %05.10g, diff=%05.10g\n", k, reg_bloc2[k+10][2], reg_bloc2_2[k+10][2], reg_bloc2[k+10][2] - reg_bloc2_2[k+10][2]);
- get_error (reg_bloc2[k+10][2], reg_bloc2_2[k+10][2], &error);
- }
- for (k=-20; k<20; k++)
- {
- if (k==0)
- {
- get_error( reg_a0, reg_a0_2, &error);
- get_error( reg_c0, reg_c0_2, &error);
- //printf ("a0 = %05.10g <==> %05.10g, diff=%05.10g\n", reg_a0, reg_a0_2, reg_a0 - reg_a0_2);
- //printf ("c0 = %05.10g <==> %05.10g, diff=%05.10g\n", reg_c0, reg_c0_2, reg_c0 - reg_c0_2);
-
- }
- else
- {
-#ifdef SINC_STORAGE
- printf ("sinc[%d]=%05.10g <==> %05.10g, diff=%05.10g\n", k, tmp_sinc[k+20], tmp_sinc_2[k+20], tmp_sinc[k+20] - tmp_sinc_2[k+20]);
- printf ("sincp[%d]=%05.10g <==> %05.10g, diff=%05.10g\n", k, tmp_sincp[k+20], tmp_sincp_2[k+20], tmp_sincp[k+20] - tmp_sincp_2[k+20]);
-#endif
- get_error (reg_beta[k+20], reg_beta_2[k+20], &error);
- //printf ("pente[%d]=%05.10g <==> %05.10g, diff=%05.10g\n", k, reg_beta[k+20], reg_beta_2[k+20], reg_beta[k+20] - reg_beta_2[k+20]);
- get_error (reg_alpha[k+20], reg_alpha_2[k+20], &error);
- //printf ("ordo[%d]=%05.10g <==> %05.10g, diff=%05.10g\n\n", k, reg_alpha[k+20], reg_alpha_2[k+20], reg_alpha[k+20] - reg_alpha_2[k+20]);
- }
-
- }
-#endif
- printf ("%d %05.18g %05.18g\n", test_rho, reg_a0, reg_c0);
- error = 0;
- }
- //printf ("worst of worst error = %05.12g\n", error);
-
- return 0;
-}
-