[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
s.popinet at niwa.co.nz
Fri May 15 02:52:55 UTC 2009
The following commit has been merged in the upstream branch:
commit f06afa12ae8353ff0900d76ea7ada327c796a6ae
Author: Stephane Popinet <s.popinet at niwa.co.nz>
Date: Thu Oct 13 07:18:33 2005 +1000
Thin-wall and waves test fixes
darcs-hash:20051012211833-fbd8f-e36c443b7883363506aa33f66258d221ba9c6a85.gz
diff --git a/test/circle/thin/error.ref b/test/circle/thin/error.ref
index 6ba8e5e..dab3967 100644
--- a/test/circle/thin/error.ref
+++ b/test/circle/thin/error.ref
@@ -2,5 +2,5 @@
4 6.417e-02 7.586e-02 2.048e-01
5 2.960e-02 3.379e-02 1.364e-01
6 1.426e-02 1.605e-02 9.463e-02
-7 7.001e-03 7.844e-03 6.646e-02
-8 3.468e-03 3.880e-03 4.687e-02
+7 7.001e-03 7.843e-03 6.646e-02
+8 3.468e-03 3.879e-03 4.686e-02
diff --git a/test/circle/thin/order.ref b/test/circle/thin/order.ref
index 32f6c49..499a967 100644
--- a/test/circle/thin/order.ref
+++ b/test/circle/thin/order.ref
@@ -2,4 +2,4 @@
5 1.1163 1.16674 0.586372
6 1.05362 1.07402 0.527474
7 1.02634 1.0331 0.509811
-8 1.01346 1.01535 0.504129
+8 1.01346 1.01572 0.504129
diff --git a/test/circle/thin/res-7.ref b/test/circle/thin/res-7.ref
index 7cbf3f4..7808fff 100644
--- a/test/circle/thin/res-7.ref
+++ b/test/circle/thin/res-7.ref
@@ -1,10 +1,10 @@
-0.38 1 2.824e+02 0.63
-0.71 2 7.886e+01 1.5
-1.04 3 2.151e+01 2
-1.37 4 5.840e+00 2.3
-1.69 5 1.585e+00 2.6
-2.03 6 4.299e-01 2.7
-2.36 7 1.166e-01 2.8
-2.69 8 3.165e-02 2.9
-3.02 9 8.587e-03 3
-3.33 10 2.330e-03 3.1
+0.4 1 2.895e+02 0.61
+0.73 2 8.132e+01 1.5
+1.07 3 2.256e+01 2
+1.39 4 6.258e+00 2.3
+1.7 5 1.736e+00 2.5
+2.03 6 4.814e-01 2.7
+2.36 7 1.335e-01 2.8
+2.68 8 3.703e-02 2.9
+3.05 9 1.027e-02 3
+3.36 10 2.848e-03 3
diff --git a/test/waves/waves.gfs b/test/waves/waves.gfs
index 5aae817..cae7fb8 100644
--- a/test/waves/waves.gfs
+++ b/test/waves/waves.gfs
@@ -55,8 +55,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
U = {
@include <gsl/gsl_sf_bessel.h>
@@ -68,10 +68,10 @@
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) -
+ 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) -
+ 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);
}
@@ -85,10 +85,10 @@
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) -
+ 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) -
+ 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);
}
@@ -110,8 +110,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -124,8 +124,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -138,8 +138,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -152,8 +152,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -166,8 +166,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -180,8 +180,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -194,8 +194,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -208,8 +208,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -222,8 +222,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -236,8 +236,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -250,8 +250,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -264,8 +264,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -278,8 +278,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -292,8 +292,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -306,8 +306,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -320,8 +320,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -334,8 +334,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -348,8 +348,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -362,8 +362,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -376,8 +376,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -390,8 +390,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -404,8 +404,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -418,8 +418,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -432,8 +432,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -446,8 +446,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -460,8 +460,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -474,8 +474,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -488,8 +488,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -502,8 +502,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -516,8 +516,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -530,8 +530,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -544,8 +544,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -558,8 +558,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -572,8 +572,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -586,8 +586,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -600,8 +600,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -614,8 +614,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -628,8 +628,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -642,8 +642,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -656,8 +656,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -670,8 +670,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -684,8 +684,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -698,8 +698,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -712,8 +712,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -726,8 +726,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -740,8 +740,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -754,8 +754,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -768,8 +768,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -782,8 +782,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -796,8 +796,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -810,8 +810,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -824,8 +824,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -838,8 +838,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -852,8 +852,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -866,8 +866,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -880,8 +880,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -894,8 +894,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -908,8 +908,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -922,8 +922,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -936,8 +936,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -950,8 +950,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -964,8 +964,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -978,8 +978,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -992,8 +992,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1006,8 +1006,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1020,8 +1020,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1034,8 +1034,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1048,8 +1048,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1062,8 +1062,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1076,8 +1076,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1090,8 +1090,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1104,8 +1104,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1118,8 +1118,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1132,8 +1132,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1146,8 +1146,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1160,8 +1160,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1174,8 +1174,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1188,8 +1188,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1202,8 +1202,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1216,8 +1216,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1230,8 +1230,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1244,8 +1244,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1258,8 +1258,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1272,8 +1272,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1286,8 +1286,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1300,8 +1300,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1314,8 +1314,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1328,8 +1328,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1342,8 +1342,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1356,8 +1356,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1370,8 +1370,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1384,8 +1384,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1398,8 +1398,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1412,8 +1412,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1426,8 +1426,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1440,8 +1440,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1454,8 +1454,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1468,8 +1468,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1482,8 +1482,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1496,8 +1496,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1510,8 +1510,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1524,8 +1524,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1538,8 +1538,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1552,8 +1552,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1566,8 +1566,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
@@ -1580,8 +1580,8 @@
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);
+ double a = 1./2555510.;
+ return a*cos (k*theta)*gsl_sf_bessel_Inu (k, r/D);
}
unbiased = 1
}
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list