diff options
author | xleroy <xleroy@fca1b0fc-160b-0410-b1d3-a4f43f01ea2e> | 2006-09-17 12:04:56 +0000 |
---|---|---|
committer | xleroy <xleroy@fca1b0fc-160b-0410-b1d3-a4f43f01ea2e> | 2006-09-17 12:04:56 +0000 |
commit | 2ec5b3fb2ccb0120be641e077089f3da5e53d8a3 (patch) | |
tree | d1e6f6a9a3d99f0095dc96fe320214de68918158 /test/c/spectral.c | |
parent | e37d620f5b9b05e16563545cba9c538f8d31c746 (diff) | |
download | compcert-2ec5b3fb2ccb0120be641e077089f3da5e53d8a3.tar.gz compcert-2ec5b3fb2ccb0120be641e077089f3da5e53d8a3.zip |
Davantage de tests
git-svn-id: https://yquem.inria.fr/compcert/svn/compcert/trunk@104 fca1b0fc-160b-0410-b1d3-a4f43f01ea2e
Diffstat (limited to 'test/c/spectral.c')
-rw-r--r-- | test/c/spectral.c | 81 |
1 files changed, 81 insertions, 0 deletions
diff --git a/test/c/spectral.c b/test/c/spectral.c new file mode 100644 index 00000000..68579d9e --- /dev/null +++ b/test/c/spectral.c @@ -0,0 +1,81 @@ +/* -*- mode: c -*- + * + * The Great Computer Language Shootout + * http://shootout.alioth.debian.org/ + * + * Contributed by Sebastien Loisel + */ + +#include <stdio.h> +#include <stdlib.h> +#include <math.h> + +double eval_A(int i, int j) { return 1.0/((i+j)*(i+j+1)/2+i+1); } + +void eval_A_times_u(int N, const double u[], double Au[]) +{ + int i,j; + for(i=0;i<N;i++) + { + Au[i]=0; + for(j=0;j<N;j++) Au[i]+=eval_A(i,j)*u[j]; + } +} + +void eval_At_times_u(int N, const double u[], double Au[]) +{ + int i,j; + for(i=0;i<N;i++) + { + Au[i]=0; + for(j=0;j<N;j++) Au[i]+=eval_A(j,i)*u[j]; + } +} + +void eval_AtA_times_u(int N, const double u[], double AtAu[]) +{ + double *v = malloc(N * sizeof(double)); + eval_A_times_u(N,u,v); + eval_At_times_u(N,v,AtAu); + free(v); +} + +int main(int argc, char *argv[]) +{ + int i; + int N = ((argc == 2) ? atoi(argv[1]) : 2500); + double * u, * v, vBv, vv; + u = malloc(N * sizeof(double)); + v = malloc(N * sizeof(double)); + for(i=0;i<N;i++) u[i]=1; + for(i=0;i<10;i++) + { + eval_AtA_times_u(N,u,v); + eval_AtA_times_u(N,v,u); + } + vBv=vv=0; + for(i=0;i<N;i++) { vBv+=u[i]*v[i]; vv+=v[i]*v[i]; } + printf("%0.9f\n",sqrt(vBv/vv)); + return 0; +} + +/******** + build & benchmark results + +BUILD COMMANDS FOR: spectralnorm.gcc + +Fri Sep 15 20:48:08 PDT 2006 + +/usr/bin/gcc -pipe -Wall -O3 -fomit-frame-pointer -funroll-loops -march=pentium4 -lm spectralnorm.c -o spectralnorm.gcc_run + +================================================================= +COMMAND LINE (%A is single numeric argument): + +spectralnorm.gcc_run %A + +N=2500 + +PROGRAM OUTPUT +============== +1.274224153 +********/ |