[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
popinet at users.sf.net
Fri May 15 02:54:08 UTC 2009
The following commit has been merged in the upstream branch:
commit b5a1bebfe031774f50f874e63020e3be698c3ee2
Author: Stephane Popinet <popinet at users.sf.net>
Date: Sun Mar 18 11:18:19 2007 +1100
Simplified circular wave test case (using GfsGlobal)
darcs-hash:20070318001819-d4795-60929c4b62cd8983b2a221f1036dd0d47ba47e6c.gz
diff --git a/test/waves/correlation.ref b/test/waves/correlation.ref
index 5062942..0661cf5 100644
--- a/test/waves/correlation.ref
+++ b/test/waves/correlation.ref
@@ -1,3 +1,3 @@
-5 2.4 0.9259
-6 0.4 1
-7 0 1
+5 2.3 0.9252
+6 0.3 0.9995
+7 0.1 1
diff --git a/test/waves/waves.gfs b/test/waves/waves.gfs
index 9361244..0803511 100644
--- a/test/waves/waves.gfs
+++ b/test/waves/waves.gfs
@@ -38,58 +38,57 @@
#
# Author: St\'ephane Popinet
# Command: sh waves.sh waves.gfs
-# Version: 0.8.0
+# Version: 1.0.0
# Required files: waves.sh solution.gfv correlation.ref
# Generated files: correlation correlation.tex solution.eps
#
1 0 GfsOcean GfsBox GfsGEdge {} {
Time { end = 37.80501984 dtmax = 0.1 }
PhysicalParams { g = 5.87060327757e-3 }
- Init {} {
- P = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx);
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
+ Global {
+ #include <gsl/gsl_sf_bessel.h>
+ @link -lgsl -lgslcblas
+
+ #define Ik(k,r,D) (gsl_sf_bessel_Inu ((k) - 1., (r)/(D))/(D)\
+ - (k)/(r)*gsl_sf_bessel_Inu ((k), (r)/(D)))
+
+ static double D = 8.83906519983e-2;
+ static double k = 3.;
+ static double sigma = 0.4986;
+ static double a = 1./2555510.;
+
+ static double pwave (double x, double y, double angle) {
+ double theta = atan2 (y, x) + angle*M_PI/180.;
+ double r = sqrt (x*x + y*y);
return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
- U = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- @define Ik(k,r,D) (gsl_sf_bessel_Inu ((k) - 1., (r)/(D))/(D)\
- - (k)/(r)*gsl_sf_bessel_Inu ((k), (r)/(D)))
- double theta = atan2(y,x);
+
+ static double ur (double theta, double r) {
+ return -a*D*D/5.87060327757e-3*sin (k*theta)*(sigma*Ik (k, r, D) -
+ k/r*gsl_sf_bessel_Inu (k, r/D));
+ }
+
+ static double vt (double theta, double r) {
+ return a*D*D/5.87060327757e-3*cos (k*theta)*(Ik (k, r, D) -
+ k*sigma/r*gsl_sf_bessel_Inu (k, r/D));
+ }
+
+ static double uwave (double x, double y) {
+ double theta = atan2 (y, x);
double r = sqrt (x*x + y*y);
- double D = 8.83906519983e-2;
- double k = 3.;
- double sigma = 0.4986;
- double a = 1./2555510.;
- double ur = -a*D*D/5.87060327757e-3*sin (k*theta)*(sigma*Ik (k, r, D) -
- k/r*gsl_sf_bessel_Inu (k, r/D));
- double vt = a*D*D/5.87060327757e-3*cos (k*theta)*(Ik (k, r, D) -
- k*sigma/r*gsl_sf_bessel_Inu (k, r/D));
- return ur*cos (theta) - vt*sin (theta);
+ return ur (theta, r)*cos (theta) - vt (theta, r)*sin (theta);
}
- V = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- @define Ik(k,r,D) (gsl_sf_bessel_Inu ((k) - 1., (r)/(D))/(D)\
- - (k)/(r)*gsl_sf_bessel_Inu ((k), (r)/(D)))
- double theta = atan2(y,x);
+
+ static double vwave (double x, double y) {
+ double theta = atan2 (y, x);
double r = sqrt (x*x + y*y);
- double D = 8.83906519983e-2;
- double k = 3.;
- double sigma = 0.4986;
- double a = 1./2555510.;
- double ur = -a*D*D/5.87060327757e-3*sin (k*theta)*(sigma*Ik (k, r, D) -
- k/r*gsl_sf_bessel_Inu (k, r/D));
- double vt = a*D*D/5.87060327757e-3*cos (k*theta)*(Ik (k, r, D) -
- k*sigma/r*gsl_sf_bessel_Inu (k, r/D));
- return ur*sin (theta) + vt*cos (theta);
+ return ur (theta, r)*sin (theta) + vt (theta, r)*cos (theta);
}
+ }
+ Init {} {
+ P = pwave(cx, cy, 0)
+ U = uwave(x, y)
+ V = vwave(x, y)
H = 1
}
Refine LEVEL
@@ -97,1490 +96,536 @@
AdvectionParams { scheme = none }
ApproxProjectionParams { tolerance = 1e-9 weighted = 0 }
SourceCoriolis {} 1.
- EventScript { start = end } { echo -n "-30 " }
EventFilter { istep = 1 } U 4
EventFilter { istep = 1 } V 4
+ EventScript { start = end } { echo -n "-30 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -30*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -30)
unbiased = 1
}
EventScript { start = end } { echo -n "-29 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -29*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -29)
unbiased = 1
}
EventScript { start = end } { echo -n "-28 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -28*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -28)
unbiased = 1
}
EventScript { start = end } { echo -n "-27 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -27*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -27)
unbiased = 1
}
EventScript { start = end } { echo -n "-26 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -26*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -26)
unbiased = 1
}
EventScript { start = end } { echo -n "-25 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -25*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -25)
unbiased = 1
}
EventScript { start = end } { echo -n "-24 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -24*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -24)
unbiased = 1
}
EventScript { start = end } { echo -n "-23 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -23*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -23)
unbiased = 1
}
EventScript { start = end } { echo -n "-22 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -22*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -22)
unbiased = 1
}
EventScript { start = end } { echo -n "-21 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -21*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -21)
unbiased = 1
}
EventScript { start = end } { echo -n "-20 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -20*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -20)
unbiased = 1
}
EventScript { start = end } { echo -n "-19 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -19*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -19)
unbiased = 1
}
EventScript { start = end } { echo -n "-18 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -18*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -18)
unbiased = 1
}
EventScript { start = end } { echo -n "-17 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -17*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -17)
unbiased = 1
}
EventScript { start = end } { echo -n "-16 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -16*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -16)
unbiased = 1
}
EventScript { start = end } { echo -n "-15 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -15*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -15)
unbiased = 1
}
EventScript { start = end } { echo -n "-14 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -14*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -14)
unbiased = 1
}
EventScript { start = end } { echo -n "-13 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -13*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -13)
unbiased = 1
}
EventScript { start = end } { echo -n "-12 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -12*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -12)
unbiased = 1
}
EventScript { start = end } { echo -n "-11 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -11*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -11)
unbiased = 1
}
EventScript { start = end } { echo -n "-10 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -10*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -10)
unbiased = 1
}
EventScript { start = end } { echo -n "-9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -9)
unbiased = 1
}
EventScript { start = end } { echo -n "-8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -8)
unbiased = 1
}
EventScript { start = end } { echo -n "-7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -7)
unbiased = 1
}
EventScript { start = end } { echo -n "-6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -6)
unbiased = 1
}
EventScript { start = end } { echo -n "-5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -5)
unbiased = 1
}
EventScript { start = end } { echo -n "-4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -4)
unbiased = 1
}
EventScript { start = end } { echo -n "-3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -3)
unbiased = 1
}
EventScript { start = end } { echo -n "-2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -2)
unbiased = 1
}
EventScript { start = end } { echo -n "-1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + -1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, -1)
unbiased = 1
}
EventScript { start = end } { echo -n "0 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0)
unbiased = 1
}
EventScript { start = end } { echo -n "0.1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.1)
unbiased = 1
}
EventScript { start = end } { echo -n "0.2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.2)
unbiased = 1
}
EventScript { start = end } { echo -n "0.3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.3)
unbiased = 1
}
EventScript { start = end } { echo -n "0.4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.4)
unbiased = 1
}
EventScript { start = end } { echo -n "0.5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.5)
unbiased = 1
}
EventScript { start = end } { echo -n "0.6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.6)
unbiased = 1
}
EventScript { start = end } { echo -n "0.7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.7)
unbiased = 1
}
EventScript { start = end } { echo -n "0.8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.8)
unbiased = 1
}
EventScript { start = end } { echo -n "0.9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 0.9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 0.9)
unbiased = 1
}
EventScript { start = end } { echo -n "1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1)
unbiased = 1
}
EventScript { start = end } { echo -n "1.1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.1)
unbiased = 1
}
EventScript { start = end } { echo -n "1.2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.2)
unbiased = 1
}
EventScript { start = end } { echo -n "1.3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.3)
unbiased = 1
}
EventScript { start = end } { echo -n "1.4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.4)
unbiased = 1
}
EventScript { start = end } { echo -n "1.5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.5)
unbiased = 1
}
EventScript { start = end } { echo -n "1.6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.6)
unbiased = 1
}
EventScript { start = end } { echo -n "1.7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.7)
unbiased = 1
}
EventScript { start = end } { echo -n "1.8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.8)
unbiased = 1
}
EventScript { start = end } { echo -n "1.9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 1.9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 1.9)
unbiased = 1
}
EventScript { start = end } { echo -n "2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2)
unbiased = 1
}
EventScript { start = end } { echo -n "2.1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.1)
unbiased = 1
}
EventScript { start = end } { echo -n "2.2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.2)
unbiased = 1
}
EventScript { start = end } { echo -n "2.3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.3)
unbiased = 1
}
EventScript { start = end } { echo -n "2.4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.4)
unbiased = 1
}
EventScript { start = end } { echo -n "2.5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.5)
unbiased = 1
}
EventScript { start = end } { echo -n "2.6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.6)
unbiased = 1
}
EventScript { start = end } { echo -n "2.7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.7)
unbiased = 1
}
EventScript { start = end } { echo -n "2.8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.8)
unbiased = 1
}
EventScript { start = end } { echo -n "2.9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 2.9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 2.9)
unbiased = 1
}
EventScript { start = end } { echo -n "3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3)
unbiased = 1
}
EventScript { start = end } { echo -n "3.1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.1)
unbiased = 1
}
EventScript { start = end } { echo -n "3.2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.2)
unbiased = 1
}
EventScript { start = end } { echo -n "3.3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.3)
unbiased = 1
}
EventScript { start = end } { echo -n "3.4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.4)
unbiased = 1
}
EventScript { start = end } { echo -n "3.5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.5)
unbiased = 1
}
EventScript { start = end } { echo -n "3.6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.6)
unbiased = 1
}
EventScript { start = end } { echo -n "3.7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.7)
unbiased = 1
}
EventScript { start = end } { echo -n "3.8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.8)
unbiased = 1
}
EventScript { start = end } { echo -n "3.9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 3.9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 3.9)
unbiased = 1
}
EventScript { start = end } { echo -n "4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4)
unbiased = 1
}
EventScript { start = end } { echo -n "4.1 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.1*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.1)
unbiased = 1
}
EventScript { start = end } { echo -n "4.2 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.2*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.2)
unbiased = 1
}
EventScript { start = end } { echo -n "4.3 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.3*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.3)
unbiased = 1
}
EventScript { start = end } { echo -n "4.4 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.4*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.4)
unbiased = 1
}
EventScript { start = end } { echo -n "4.5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.5)
unbiased = 1
}
EventScript { start = end } { echo -n "4.6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.6)
unbiased = 1
}
EventScript { start = end } { echo -n "4.7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.7)
unbiased = 1
}
EventScript { start = end } { echo -n "4.8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.8)
unbiased = 1
}
EventScript { start = end } { echo -n "4.9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 4.9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 4.9)
unbiased = 1
}
EventScript { start = end } { echo -n "5 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 5*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 5)
unbiased = 1
}
EventScript { start = end } { echo -n "6 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 6*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 6)
unbiased = 1
}
EventScript { start = end } { echo -n "7 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 7*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 7)
unbiased = 1
}
EventScript { start = end } { echo -n "8 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 8*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 8)
unbiased = 1
}
EventScript { start = end } { echo -n "9 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 9*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 9)
unbiased = 1
}
EventScript { start = end } { echo -n "10 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 10*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 10)
unbiased = 1
}
EventScript { start = end } { echo -n "11 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 11*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 11)
unbiased = 1
}
EventScript { start = end } { echo -n "12 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 12*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 12)
unbiased = 1
}
EventScript { start = end } { echo -n "13 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 13*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 13)
unbiased = 1
}
EventScript { start = end } { echo -n "14 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 14*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 14)
unbiased = 1
}
EventScript { start = end } { echo -n "15 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 15*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 15)
unbiased = 1
}
EventScript { start = end } { echo -n "16 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 16*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 16)
unbiased = 1
}
EventScript { start = end } { echo -n "17 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 17*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 17)
unbiased = 1
}
EventScript { start = end } { echo -n "18 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 18*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 18)
unbiased = 1
}
EventScript { start = end } { echo -n "19 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 19*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 19)
unbiased = 1
}
EventScript { start = end } { echo -n "20 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 20*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 20)
unbiased = 1
}
EventScript { start = end } { echo -n "21 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 21*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 21)
unbiased = 1
}
EventScript { start = end } { echo -n "22 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 22*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 22)
unbiased = 1
}
EventScript { start = end } { echo -n "23 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 23*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 23)
unbiased = 1
}
EventScript { start = end } { echo -n "24 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 24*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 24)
unbiased = 1
}
EventScript { start = end } { echo -n "25 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 25*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 25)
unbiased = 1
}
EventScript { start = end } { echo -n "26 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 26*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 26)
unbiased = 1
}
EventScript { start = end } { echo -n "27 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 27*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 27)
unbiased = 1
}
EventScript { start = end } { echo -n "28 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 28*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 28)
unbiased = 1
}
EventScript { start = end } { echo -n "29 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 29*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 29)
unbiased = 1
}
EventScript { start = end } { echo -n "30 " }
OutputCorrelation { start = end } stdout { v = P } {
- s = {
- @include <gsl/gsl_sf_bessel.h>
- @link -lgsl -lgslcblas
- double theta = atan2(cy,cx) + 30*M_PI/180.;
- double r = sqrt (cx*cx + cy*cy);
- double D = 8.83906519983e-2;
- double k = 3.;
- double a = 1./2555510.;
- return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
- }
+ s = pwave(cx, cy, 30)
unbiased = 1
}
OutputSimulation { start = end } sim-LEVEL {}
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list