[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
stephane.popinet at paradise.net.nz
Fri May 15 02:51:58 UTC 2009
The following commit has been merged in the upstream branch:
commit 8b719076027351a41c291900a627a8d8791ea7e1
Author: Stephane Popinet <stephane.popinet at paradise.net.nz>
Date: Mon Jun 13 20:01:28 2005 +1000
New 'Reynolds' test case, replacement for 'stationary'
Using auto-documenting with gfs2doc.
darcs-hash:20050613100128-fa380-d57485c9c501e7ec2709c356d0f09faf278bd941.gz
diff --git a/test/euler/reynolds/check.py b/test/euler/reynolds/check.py
new file mode 100644
index 0000000..e106665
--- /dev/null
+++ b/test/euler/reynolds/check.py
@@ -0,0 +1,66 @@
+from sys import *
+from math import *
+
+class Curve:
+ def read(self,file,x,y):
+ line = file.readline()
+ while line:
+ record = line.split()
+ self.l.append((float(record[x-1]),float(record[y-1])))
+ line = file.readline()
+ def __init__(self,f=None,x=1,y=2):
+ self.l = []
+ if f:
+ self.read(open(f,'r'),x,y)
+ def value(self,x):
+ l,u = 0,len(self.l)-1
+ while u - l > 1:
+ i = (u + l)/2
+ if self.l[i][0] > x:
+ u = i
+ else:
+ l = i
+ x0,x1,y0,y1 = self.l[l][0],self.l[u][0],self.l[l][1],self.l[u][1]
+ return y0 + (x - x0)/(x1 - x0)*(y1 - y0)
+ def __sub__(self,other):
+ c = Curve()
+ if len(self.l) < len(other.l):
+ for p in self.l:
+ c.l.append((p[0], p[1] - other.value(p[0])))
+ else:
+ for p in other.l:
+ c.l.append((p[0], self.value(p[0]) - p[1]))
+ return c
+ def sum(self,f=lambda x: x):
+ s = 0.
+ p1 = None
+ for p in self.l:
+ if p1:
+ s += (p[0] - p1[0])*f((p[1] + p1[1])/2.)
+ p1 = p
+ return s
+ def max(self,f=lambda x: x):
+ m = None
+ for p in self.l:
+ v = f(p[1])
+ if not m or v > m:
+ m = v
+ return m
+ def min(self,f=lambda x: x):
+ m = None
+ for p in self.l:
+ v = f(p[1])
+ if not m or v < m:
+ m = v
+ return m
+ def mean(self):
+ return self.sum()/self.sum(lambda x: 1.)
+ def norm1(self):
+ return self.sum(lambda x: abs(x))/self.sum(lambda x: 1.)
+ def norm2(self):
+ return sqrt(self.sum(lambda x: x*x)/self.sum(lambda x: 1.))
+ def normi(self):
+ return self.max(lambda x: abs(x))
+ def write(self):
+ for p in self.l:
+ print p[0],p[1]
diff --git a/test/euler/reynolds/div5.ref b/test/euler/reynolds/div5.ref
new file mode 100644
index 0000000..762147c
--- /dev/null
+++ b/test/euler/reynolds/div5.ref
@@ -0,0 +1,81 @@
+Divergence time: 0 first: 1.051e-15 second: 1.330e-15 infty: 4.441e-15
+Divergence time: 0.0252425 first: 3.092e-02 second: 3.805e-02 infty: 9.125e-02
+Divergence time: 0.0504502 first: 3.245e-02 second: 4.028e-02 infty: 1.005e-01
+Divergence time: 0.0756655 first: 3.264e-02 second: 4.086e-02 infty: 1.062e-01
+Divergence time: 0.100885 first: 3.282e-02 second: 4.180e-02 infty: 1.096e-01
+Divergence time: 0.126098 first: 3.324e-02 second: 4.226e-02 infty: 1.111e-01
+Divergence time: 0.151315 first: 3.305e-02 second: 4.237e-02 infty: 1.129e-01
+Divergence time: 0.176554 first: 3.271e-02 second: 4.224e-02 infty: 1.148e-01
+Divergence time: 0.201797 first: 3.318e-02 second: 4.281e-02 infty: 1.150e-01
+Divergence time: 0.22703 first: 3.294e-02 second: 4.249e-02 infty: 1.155e-01
+Divergence time: 0.252277 first: 3.349e-02 second: 4.305e-02 infty: 1.166e-01
+Divergence time: 0.277513 first: 3.335e-02 second: 4.273e-02 infty: 1.170e-01
+Divergence time: 0.302764 first: 3.272e-02 second: 4.229e-02 infty: 1.169e-01
+Divergence time: 0.327997 first: 3.324e-02 second: 4.255e-02 infty: 1.168e-01
+Divergence time: 0.353214 first: 3.325e-02 second: 4.239e-02 infty: 1.165e-01
+Divergence time: 0.378429 first: 3.284e-02 second: 4.235e-02 infty: 1.147e-01
+Divergence time: 0.403651 first: 3.307e-02 second: 4.250e-02 infty: 1.159e-01
+Divergence time: 0.428874 first: 3.299e-02 second: 4.236e-02 infty: 1.166e-01
+Divergence time: 0.4541 first: 3.292e-02 second: 4.224e-02 infty: 1.180e-01
+Divergence time: 0.479329 first: 3.286e-02 second: 4.220e-02 infty: 1.175e-01
+Divergence time: 0.504559 first: 3.280e-02 second: 4.209e-02 infty: 1.175e-01
+Divergence time: 0.529791 first: 3.295e-02 second: 4.214e-02 infty: 1.170e-01
+Divergence time: 0.555025 first: 3.326e-02 second: 4.235e-02 infty: 1.163e-01
+Divergence time: 0.580263 first: 3.330e-02 second: 4.242e-02 infty: 1.150e-01
+Divergence time: 0.605505 first: 3.330e-02 second: 4.239e-02 infty: 1.147e-01
+Divergence time: 0.630748 first: 3.328e-02 second: 4.233e-02 infty: 1.159e-01
+Divergence time: 0.655995 first: 3.325e-02 second: 4.228e-02 infty: 1.168e-01
+Divergence time: 0.681246 first: 3.321e-02 second: 4.224e-02 infty: 1.176e-01
+Divergence time: 0.706502 first: 3.319e-02 second: 4.222e-02 infty: 1.181e-01
+Divergence time: 0.731764 first: 3.318e-02 second: 4.221e-02 infty: 1.186e-01
+Divergence time: 0.757032 first: 3.317e-02 second: 4.219e-02 infty: 1.188e-01
+Divergence time: 0.782306 first: 3.317e-02 second: 4.219e-02 infty: 1.190e-01
+Divergence time: 0.807586 first: 3.317e-02 second: 4.218e-02 infty: 1.190e-01
+Divergence time: 0.832872 first: 3.317e-02 second: 4.217e-02 infty: 1.190e-01
+Divergence time: 0.858163 first: 3.346e-02 second: 4.235e-02 infty: 1.191e-01
+Divergence time: 0.883458 first: 3.378e-02 second: 4.262e-02 infty: 1.192e-01
+Divergence time: 0.908753 first: 3.383e-02 second: 4.269e-02 infty: 1.191e-01
+Divergence time: 0.934045 first: 3.384e-02 second: 4.264e-02 infty: 1.190e-01
+Divergence time: 0.959337 first: 3.380e-02 second: 4.257e-02 infty: 1.188e-01
+Divergence time: 0.984633 first: 3.378e-02 second: 4.252e-02 infty: 1.186e-01
+Divergence time: 1.00993 first: 3.377e-02 second: 4.249e-02 infty: 1.185e-01
+Divergence time: 1.03524 first: 3.375e-02 second: 4.247e-02 infty: 1.184e-01
+Divergence time: 1.06055 first: 3.373e-02 second: 4.245e-02 infty: 1.183e-01
+Divergence time: 1.08587 first: 3.391e-02 second: 4.264e-02 infty: 1.183e-01
+Divergence time: 1.11117 first: 3.404e-02 second: 4.282e-02 infty: 1.184e-01
+Divergence time: 1.13647 first: 3.407e-02 second: 4.285e-02 infty: 1.184e-01
+Divergence time: 1.16177 first: 3.407e-02 second: 4.277e-02 infty: 1.184e-01
+Divergence time: 1.18708 first: 3.407e-02 second: 4.267e-02 infty: 1.185e-01
+Divergence time: 1.21239 first: 3.407e-02 second: 4.259e-02 infty: 1.186e-01
+Divergence time: 1.23771 first: 3.406e-02 second: 4.253e-02 infty: 1.187e-01
+Divergence time: 1.26303 first: 3.404e-02 second: 4.249e-02 infty: 1.188e-01
+Divergence time: 1.28835 first: 3.412e-02 second: 4.263e-02 infty: 1.190e-01
+Divergence time: 1.31368 first: 3.417e-02 second: 4.264e-02 infty: 1.192e-01
+Divergence time: 1.339 first: 3.422e-02 second: 4.257e-02 infty: 1.193e-01
+Divergence time: 1.36433 first: 3.428e-02 second: 4.262e-02 infty: 1.195e-01
+Divergence time: 1.38967 first: 3.429e-02 second: 4.259e-02 infty: 1.197e-01
+Divergence time: 1.415 first: 3.431e-02 second: 4.251e-02 infty: 1.198e-01
+Divergence time: 1.44034 first: 3.428e-02 second: 4.244e-02 infty: 1.199e-01
+Divergence time: 1.46568 first: 3.423e-02 second: 4.239e-02 infty: 1.200e-01
+Divergence time: 1.49103 first: 3.419e-02 second: 4.235e-02 infty: 1.200e-01
+Divergence time: 1.51637 first: 3.416e-02 second: 4.233e-02 infty: 1.201e-01
+Divergence time: 1.54172 first: 3.414e-02 second: 4.232e-02 infty: 1.200e-01
+Divergence time: 1.56708 first: 3.413e-02 second: 4.231e-02 infty: 1.201e-01
+Divergence time: 1.59243 first: 3.406e-02 second: 4.223e-02 infty: 1.200e-01
+Divergence time: 1.6178 first: 3.396e-02 second: 4.211e-02 infty: 1.199e-01
+Divergence time: 1.64316 first: 3.399e-02 second: 4.218e-02 infty: 1.197e-01
+Divergence time: 1.66853 first: 3.398e-02 second: 4.223e-02 infty: 1.194e-01
+Divergence time: 1.6939 first: 3.396e-02 second: 4.226e-02 infty: 1.194e-01
+Divergence time: 1.71927 first: 3.396e-02 second: 4.227e-02 infty: 1.194e-01
+Divergence time: 1.74464 first: 3.388e-02 second: 4.219e-02 infty: 1.193e-01
+Divergence time: 1.77002 first: 3.374e-02 second: 4.205e-02 infty: 1.192e-01
+Divergence time: 1.79541 first: 3.378e-02 second: 4.211e-02 infty: 1.192e-01
+Divergence time: 1.82082 first: 3.367e-02 second: 4.202e-02 infty: 1.192e-01
+Divergence time: 1.84621 first: 3.371e-02 second: 4.208e-02 infty: 1.192e-01
+Divergence time: 1.87161 first: 3.356e-02 second: 4.197e-02 infty: 1.189e-01
+Divergence time: 1.89701 first: 3.353e-02 second: 4.202e-02 infty: 1.188e-01
+Divergence time: 1.9224 first: 3.345e-02 second: 4.198e-02 infty: 1.189e-01
+Divergence time: 1.9478 first: 3.342e-02 second: 4.201e-02 infty: 1.191e-01
+Divergence time: 1.9732 first: 3.343e-02 second: 4.204e-02 infty: 1.192e-01
+Divergence time: 1.9986 first: 3.329e-02 second: 4.188e-02 infty: 1.191e-01
+Divergence time: 2 first: 1.042e-02 second: 1.399e-02 infty: 5.882e-02
diff --git a/test/euler/reynolds/div6.ref b/test/euler/reynolds/div6.ref
new file mode 100644
index 0000000..a114d32
--- /dev/null
+++ b/test/euler/reynolds/div6.ref
@@ -0,0 +1,161 @@
+Divergence time: 0 first: 2.185e-15 second: 3.067e-15 infty: 1.332e-14
+Divergence time: 0.0125302 first: 4.099e-03 second: 5.305e-03 infty: 3.318e-02
+Divergence time: 0.025057 first: 4.457e-03 second: 5.817e-03 infty: 3.658e-02
+Divergence time: 0.0375852 first: 4.532e-03 second: 5.980e-03 infty: 3.035e-02
+Divergence time: 0.0501124 first: 4.560e-03 second: 6.147e-03 infty: 4.063e-02
+Divergence time: 0.062639 first: 4.638e-03 second: 6.283e-03 infty: 3.565e-02
+Divergence time: 0.0751658 first: 4.730e-03 second: 6.440e-03 infty: 3.617e-02
+Divergence time: 0.0876912 first: 4.845e-03 second: 6.625e-03 infty: 4.406e-02
+Divergence time: 0.100217 first: 4.867e-03 second: 6.742e-03 infty: 4.734e-02
+Divergence time: 0.112742 first: 4.895e-03 second: 6.821e-03 infty: 4.418e-02
+Divergence time: 0.125269 first: 4.962e-03 second: 6.917e-03 infty: 4.458e-02
+Divergence time: 0.137795 first: 4.974e-03 second: 6.987e-03 infty: 4.415e-02
+Divergence time: 0.150321 first: 5.057e-03 second: 7.192e-03 infty: 3.966e-02
+Divergence time: 0.162847 first: 5.090e-03 second: 7.260e-03 infty: 4.421e-02
+Divergence time: 0.175375 first: 5.054e-03 second: 7.076e-03 infty: 3.751e-02
+Divergence time: 0.187902 first: 4.977e-03 second: 6.936e-03 infty: 4.334e-02
+Divergence time: 0.200429 first: 5.010e-03 second: 7.025e-03 infty: 4.447e-02
+Divergence time: 0.212956 first: 5.043e-03 second: 7.068e-03 infty: 4.101e-02
+Divergence time: 0.225484 first: 4.997e-03 second: 6.960e-03 infty: 3.819e-02
+Divergence time: 0.238012 first: 5.028e-03 second: 7.100e-03 infty: 4.240e-02
+Divergence time: 0.250539 first: 4.997e-03 second: 7.000e-03 infty: 4.064e-02
+Divergence time: 0.263067 first: 4.975e-03 second: 6.976e-03 infty: 4.058e-02
+Divergence time: 0.275595 first: 5.006e-03 second: 7.022e-03 infty: 3.961e-02
+Divergence time: 0.288121 first: 4.957e-03 second: 6.969e-03 infty: 4.193e-02
+Divergence time: 0.300647 first: 4.998e-03 second: 7.168e-03 infty: 4.829e-02
+Divergence time: 0.313172 first: 4.918e-03 second: 6.965e-03 infty: 4.450e-02
+Divergence time: 0.3257 first: 4.889e-03 second: 6.933e-03 infty: 5.175e-02
+Divergence time: 0.338226 first: 4.934e-03 second: 7.037e-03 infty: 5.656e-02
+Divergence time: 0.350754 first: 4.951e-03 second: 7.058e-03 infty: 5.226e-02
+Divergence time: 0.36328 first: 4.947e-03 second: 6.993e-03 infty: 5.058e-02
+Divergence time: 0.375806 first: 4.970e-03 second: 6.969e-03 infty: 4.899e-02
+Divergence time: 0.388332 first: 4.996e-03 second: 6.966e-03 infty: 3.967e-02
+Divergence time: 0.400857 first: 5.018e-03 second: 6.976e-03 infty: 4.367e-02
+Divergence time: 0.413382 first: 5.008e-03 second: 6.926e-03 infty: 4.087e-02
+Divergence time: 0.425907 first: 4.982e-03 second: 6.806e-03 infty: 4.141e-02
+Divergence time: 0.438433 first: 4.978e-03 second: 6.815e-03 infty: 4.048e-02
+Divergence time: 0.450959 first: 4.963e-03 second: 6.804e-03 infty: 3.950e-02
+Divergence time: 0.463485 first: 4.920e-03 second: 6.761e-03 infty: 3.881e-02
+Divergence time: 0.476011 first: 4.894e-03 second: 6.752e-03 infty: 3.841e-02
+Divergence time: 0.488537 first: 4.854e-03 second: 6.713e-03 infty: 3.830e-02
+Divergence time: 0.501064 first: 4.832e-03 second: 6.685e-03 infty: 3.832e-02
+Divergence time: 0.51359 first: 4.829e-03 second: 6.665e-03 infty: 3.828e-02
+Divergence time: 0.526117 first: 4.827e-03 second: 6.647e-03 infty: 3.821e-02
+Divergence time: 0.538644 first: 4.826e-03 second: 6.634e-03 infty: 3.814e-02
+Divergence time: 0.551171 first: 4.825e-03 second: 6.623e-03 infty: 3.812e-02
+Divergence time: 0.563699 first: 4.825e-03 second: 6.617e-03 infty: 3.810e-02
+Divergence time: 0.576227 first: 4.824e-03 second: 6.612e-03 infty: 3.809e-02
+Divergence time: 0.588755 first: 4.824e-03 second: 6.608e-03 infty: 3.808e-02
+Divergence time: 0.601284 first: 4.824e-03 second: 6.605e-03 infty: 3.808e-02
+Divergence time: 0.613813 first: 4.823e-03 second: 6.588e-03 infty: 3.807e-02
+Divergence time: 0.626341 first: 4.826e-03 second: 6.599e-03 infty: 3.806e-02
+Divergence time: 0.63887 first: 4.827e-03 second: 6.597e-03 infty: 3.806e-02
+Divergence time: 0.651399 first: 4.828e-03 second: 6.603e-03 infty: 3.806e-02
+Divergence time: 0.663928 first: 4.829e-03 second: 6.592e-03 infty: 3.805e-02
+Divergence time: 0.676458 first: 4.832e-03 second: 6.605e-03 infty: 3.813e-02
+Divergence time: 0.688987 first: 4.834e-03 second: 6.604e-03 infty: 3.819e-02
+Divergence time: 0.701516 first: 4.835e-03 second: 6.611e-03 infty: 3.812e-02
+Divergence time: 0.714046 first: 4.837e-03 second: 6.614e-03 infty: 3.806e-02
+Divergence time: 0.726576 first: 4.838e-03 second: 6.599e-03 infty: 3.819e-02
+Divergence time: 0.739105 first: 4.842e-03 second: 6.620e-03 infty: 3.824e-02
+Divergence time: 0.751636 first: 4.846e-03 second: 6.610e-03 infty: 3.817e-02
+Divergence time: 0.764166 first: 4.848e-03 second: 6.619e-03 infty: 3.809e-02
+Divergence time: 0.776696 first: 4.849e-03 second: 6.623e-03 infty: 3.814e-02
+Divergence time: 0.789227 first: 4.851e-03 second: 6.612e-03 infty: 3.854e-02
+Divergence time: 0.801758 first: 4.855e-03 second: 6.627e-03 infty: 3.855e-02
+Divergence time: 0.814289 first: 4.857e-03 second: 6.611e-03 infty: 3.863e-02
+Divergence time: 0.826821 first: 4.861e-03 second: 6.630e-03 infty: 3.863e-02
+Divergence time: 0.839353 first: 4.863e-03 second: 6.617e-03 infty: 3.856e-02
+Divergence time: 0.851884 first: 4.867e-03 second: 6.638e-03 infty: 3.848e-02
+Divergence time: 0.864416 first: 4.868e-03 second: 6.641e-03 infty: 3.844e-02
+Divergence time: 0.876948 first: 4.869e-03 second: 6.633e-03 infty: 3.850e-02
+Divergence time: 0.889479 first: 4.873e-03 second: 6.633e-03 infty: 3.850e-02
+Divergence time: 0.902011 first: 4.878e-03 second: 6.647e-03 infty: 3.853e-02
+Divergence time: 0.914543 first: 4.879e-03 second: 6.636e-03 infty: 3.854e-02
+Divergence time: 0.927074 first: 4.880e-03 second: 6.610e-03 infty: 3.846e-02
+Divergence time: 0.939606 first: 4.890e-03 second: 6.638e-03 infty: 3.857e-02
+Divergence time: 0.952137 first: 4.898e-03 second: 6.660e-03 infty: 3.858e-02
+Divergence time: 0.964668 first: 4.903e-03 second: 6.668e-03 infty: 3.849e-02
+Divergence time: 0.9772 first: 4.906e-03 second: 6.684e-03 infty: 3.852e-02
+Divergence time: 0.989731 first: 4.913e-03 second: 6.704e-03 infty: 3.865e-02
+Divergence time: 1.00226 first: 4.919e-03 second: 6.713e-03 infty: 3.866e-02
+Divergence time: 1.01479 first: 4.922e-03 second: 6.723e-03 infty: 3.861e-02
+Divergence time: 1.02732 first: 4.922e-03 second: 6.710e-03 infty: 3.866e-02
+Divergence time: 1.03986 first: 4.922e-03 second: 6.674e-03 infty: 3.864e-02
+Divergence time: 1.05239 first: 4.931e-03 second: 6.716e-03 infty: 3.857e-02
+Divergence time: 1.06492 first: 4.934e-03 second: 6.669e-03 infty: 3.839e-02
+Divergence time: 1.07745 first: 4.936e-03 second: 6.722e-03 infty: 3.728e-02
+Divergence time: 1.08998 first: 4.950e-03 second: 6.748e-03 infty: 3.840e-02
+Divergence time: 1.10252 first: 4.957e-03 second: 6.773e-03 infty: 3.868e-02
+Divergence time: 1.11505 first: 4.958e-03 second: 6.773e-03 infty: 3.863e-02
+Divergence time: 1.12758 first: 4.959e-03 second: 6.795e-03 infty: 3.854e-02
+Divergence time: 1.14011 first: 4.967e-03 second: 6.792e-03 infty: 3.802e-02
+Divergence time: 1.15265 first: 4.955e-03 second: 6.744e-03 infty: 4.116e-02
+Divergence time: 1.16518 first: 4.938e-03 second: 6.754e-03 infty: 4.138e-02
+Divergence time: 1.17771 first: 4.949e-03 second: 6.799e-03 infty: 4.151e-02
+Divergence time: 1.19024 first: 4.960e-03 second: 6.804e-03 infty: 4.141e-02
+Divergence time: 1.20278 first: 4.949e-03 second: 6.795e-03 infty: 4.182e-02
+Divergence time: 1.21531 first: 4.946e-03 second: 6.808e-03 infty: 4.859e-02
+Divergence time: 1.22784 first: 4.946e-03 second: 6.837e-03 infty: 4.806e-02
+Divergence time: 1.24037 first: 4.948e-03 second: 6.838e-03 infty: 4.244e-02
+Divergence time: 1.25291 first: 4.963e-03 second: 6.860e-03 infty: 4.354e-02
+Divergence time: 1.26544 first: 4.955e-03 second: 6.840e-03 infty: 4.386e-02
+Divergence time: 1.27797 first: 4.945e-03 second: 6.874e-03 infty: 4.859e-02
+Divergence time: 1.2905 first: 4.949e-03 second: 6.870e-03 infty: 4.375e-02
+Divergence time: 1.30303 first: 4.966e-03 second: 6.889e-03 infty: 4.362e-02
+Divergence time: 1.31557 first: 4.962e-03 second: 6.880e-03 infty: 4.731e-02
+Divergence time: 1.3281 first: 4.918e-03 second: 6.808e-03 infty: 4.370e-02
+Divergence time: 1.34063 first: 4.898e-03 second: 6.829e-03 infty: 4.848e-02
+Divergence time: 1.35317 first: 4.900e-03 second: 6.832e-03 infty: 5.229e-02
+Divergence time: 1.3657 first: 4.929e-03 second: 6.887e-03 infty: 5.230e-02
+Divergence time: 1.37823 first: 4.943e-03 second: 6.866e-03 infty: 4.734e-02
+Divergence time: 1.39077 first: 4.910e-03 second: 6.803e-03 infty: 4.754e-02
+Divergence time: 1.4033 first: 4.917e-03 second: 6.825e-03 infty: 5.314e-02
+Divergence time: 1.41584 first: 4.904e-03 second: 6.848e-03 infty: 5.795e-02
+Divergence time: 1.42837 first: 4.912e-03 second: 6.829e-03 infty: 5.359e-02
+Divergence time: 1.4409 first: 4.918e-03 second: 6.795e-03 infty: 5.446e-02
+Divergence time: 1.45344 first: 4.927e-03 second: 6.780e-03 infty: 5.355e-02
+Divergence time: 1.46597 first: 4.920e-03 second: 6.811e-03 infty: 5.524e-02
+Divergence time: 1.47851 first: 4.911e-03 second: 6.778e-03 infty: 5.286e-02
+Divergence time: 1.49104 first: 4.898e-03 second: 6.711e-03 infty: 5.779e-02
+Divergence time: 1.50357 first: 4.909e-03 second: 6.759e-03 infty: 5.867e-02
+Divergence time: 1.51611 first: 4.926e-03 second: 6.761e-03 infty: 5.850e-02
+Divergence time: 1.52864 first: 4.929e-03 second: 6.773e-03 infty: 5.671e-02
+Divergence time: 1.54118 first: 4.899e-03 second: 6.644e-03 infty: 4.476e-02
+Divergence time: 1.55371 first: 4.900e-03 second: 6.641e-03 infty: 4.351e-02
+Divergence time: 1.56625 first: 4.884e-03 second: 6.622e-03 infty: 4.800e-02
+Divergence time: 1.57878 first: 4.887e-03 second: 6.630e-03 infty: 4.910e-02
+Divergence time: 1.59132 first: 4.883e-03 second: 6.596e-03 infty: 5.304e-02
+Divergence time: 1.60385 first: 4.909e-03 second: 6.647e-03 infty: 4.095e-02
+Divergence time: 1.61639 first: 4.893e-03 second: 6.606e-03 infty: 5.055e-02
+Divergence time: 1.62892 first: 4.872e-03 second: 6.561e-03 infty: 4.546e-02
+Divergence time: 1.64146 first: 4.869e-03 second: 6.554e-03 infty: 4.611e-02
+Divergence time: 1.65399 first: 4.879e-03 second: 6.565e-03 infty: 5.310e-02
+Divergence time: 1.66653 first: 4.856e-03 second: 6.559e-03 infty: 4.624e-02
+Divergence time: 1.67907 first: 4.852e-03 second: 6.538e-03 infty: 4.603e-02
+Divergence time: 1.6916 first: 4.837e-03 second: 6.547e-03 infty: 4.548e-02
+Divergence time: 1.70414 first: 4.837e-03 second: 6.552e-03 infty: 4.626e-02
+Divergence time: 1.71668 first: 4.847e-03 second: 6.540e-03 infty: 5.308e-02
+Divergence time: 1.72921 first: 4.863e-03 second: 6.585e-03 infty: 4.438e-02
+Divergence time: 1.74175 first: 4.857e-03 second: 6.589e-03 infty: 4.653e-02
+Divergence time: 1.75428 first: 4.843e-03 second: 6.556e-03 infty: 4.545e-02
+Divergence time: 1.76682 first: 4.864e-03 second: 6.600e-03 infty: 4.608e-02
+Divergence time: 1.77936 first: 4.869e-03 second: 6.644e-03 infty: 5.305e-02
+Divergence time: 1.79189 first: 4.873e-03 second: 6.630e-03 infty: 4.361e-02
+Divergence time: 1.80443 first: 4.869e-03 second: 6.595e-03 infty: 4.365e-02
+Divergence time: 1.81697 first: 4.873e-03 second: 6.658e-03 infty: 4.542e-02
+Divergence time: 1.8295 first: 4.855e-03 second: 6.625e-03 infty: 4.606e-02
+Divergence time: 1.84204 first: 4.870e-03 second: 6.623e-03 infty: 5.304e-02
+Divergence time: 1.85458 first: 4.841e-03 second: 6.544e-03 infty: 3.999e-02
+Divergence time: 1.86711 first: 4.814e-03 second: 6.501e-03 infty: 4.411e-02
+Divergence time: 1.87965 first: 4.827e-03 second: 6.576e-03 infty: 4.933e-02
+Divergence time: 1.89219 first: 4.848e-03 second: 6.603e-03 infty: 5.234e-02
+Divergence time: 1.90473 first: 4.875e-03 second: 6.640e-03 infty: 5.302e-02
+Divergence time: 1.91726 first: 4.876e-03 second: 6.631e-03 infty: 4.033e-02
+Divergence time: 1.9298 first: 4.863e-03 second: 6.613e-03 infty: 4.685e-02
+Divergence time: 1.94234 first: 4.856e-03 second: 6.634e-03 infty: 5.298e-02
+Divergence time: 1.95487 first: 4.884e-03 second: 6.710e-03 infty: 5.791e-02
+Divergence time: 1.96741 first: 4.871e-03 second: 6.650e-03 infty: 5.302e-02
+Divergence time: 1.97995 first: 4.857e-03 second: 6.588e-03 infty: 3.997e-02
+Divergence time: 1.99249 first: 4.859e-03 second: 6.579e-03 infty: 4.587e-02
+Divergence time: 2 first: 3.352e-03 second: 4.571e-03 infty: 2.964e-02
diff --git a/test/euler/reynolds/div7.ref b/test/euler/reynolds/div7.ref
new file mode 100644
index 0000000..1c9b834
--- /dev/null
+++ b/test/euler/reynolds/div7.ref
@@ -0,0 +1,321 @@
+Divergence time: 0 first: 4.516e-15 second: 6.442e-15 infty: 2.842e-14
+Divergence time: 0.00625377 first: 5.257e-04 second: 8.349e-04 infty: 1.666e-02
+Divergence time: 0.0125071 first: 5.953e-04 second: 1.016e-03 infty: 1.737e-02
+Divergence time: 0.0187607 first: 6.167e-04 second: 1.039e-03 infty: 1.542e-02
+Divergence time: 0.0250141 first: 6.453e-04 second: 1.119e-03 infty: 1.960e-02
+Divergence time: 0.0312674 first: 6.796e-04 second: 1.164e-03 infty: 1.574e-02
+Divergence time: 0.0375208 first: 6.960e-04 second: 1.186e-03 infty: 1.597e-02
+Divergence time: 0.0437743 first: 7.080e-04 second: 1.228e-03 infty: 2.090e-02
+Divergence time: 0.0500275 first: 7.190e-04 second: 1.265e-03 infty: 2.266e-02
+Divergence time: 0.0562806 first: 7.442e-04 second: 1.302e-03 infty: 1.846e-02
+Divergence time: 0.0625338 first: 7.638e-04 second: 1.335e-03 infty: 2.122e-02
+Divergence time: 0.0687869 first: 7.765e-04 second: 1.365e-03 infty: 2.058e-02
+Divergence time: 0.0750401 first: 7.825e-04 second: 1.379e-03 infty: 2.266e-02
+Divergence time: 0.0812931 first: 8.038e-04 second: 1.413e-03 infty: 2.118e-02
+Divergence time: 0.0875462 first: 7.931e-04 second: 1.388e-03 infty: 1.784e-02
+Divergence time: 0.0937993 first: 8.249e-04 second: 1.445e-03 infty: 1.957e-02
+Divergence time: 0.100052 first: 8.023e-04 second: 1.402e-03 infty: 2.134e-02
+Divergence time: 0.106305 first: 8.167e-04 second: 1.462e-03 infty: 2.299e-02
+Divergence time: 0.112559 first: 8.260e-04 second: 1.494e-03 infty: 2.184e-02
+Divergence time: 0.118812 first: 8.352e-04 second: 1.504e-03 infty: 2.263e-02
+Divergence time: 0.125065 first: 8.475e-04 second: 1.538e-03 infty: 2.248e-02
+Divergence time: 0.131318 first: 8.570e-04 second: 1.563e-03 infty: 2.326e-02
+Divergence time: 0.137571 first: 8.706e-04 second: 1.573e-03 infty: 2.146e-02
+Divergence time: 0.143824 first: 8.529e-04 second: 1.533e-03 infty: 2.254e-02
+Divergence time: 0.150077 first: 8.654e-04 second: 1.586e-03 infty: 2.258e-02
+Divergence time: 0.15633 first: 8.736e-04 second: 1.590e-03 infty: 2.298e-02
+Divergence time: 0.162583 first: 8.713e-04 second: 1.566e-03 infty: 2.276e-02
+Divergence time: 0.168836 first: 8.434e-04 second: 1.513e-03 infty: 2.244e-02
+Divergence time: 0.175089 first: 8.665e-04 second: 1.602e-03 infty: 2.232e-02
+Divergence time: 0.181343 first: 8.769e-04 second: 1.591e-03 infty: 2.202e-02
+Divergence time: 0.187596 first: 8.853e-04 second: 1.591e-03 infty: 2.196e-02
+Divergence time: 0.193849 first: 8.134e-04 second: 1.428e-03 infty: 2.254e-02
+Divergence time: 0.200102 first: 8.671e-04 second: 1.572e-03 infty: 2.066e-02
+Divergence time: 0.206355 first: 8.804e-04 second: 1.521e-03 infty: 2.210e-02
+Divergence time: 0.212608 first: 8.066e-04 second: 1.419e-03 infty: 2.191e-02
+Divergence time: 0.218862 first: 8.534e-04 second: 1.579e-03 infty: 2.221e-02
+Divergence time: 0.225115 first: 8.630e-04 second: 1.575e-03 infty: 2.176e-02
+Divergence time: 0.231368 first: 8.617e-04 second: 1.543e-03 infty: 1.920e-02
+Divergence time: 0.237621 first: 8.379e-04 second: 1.505e-03 infty: 2.171e-02
+Divergence time: 0.243875 first: 8.636e-04 second: 1.603e-03 infty: 2.276e-02
+Divergence time: 0.250128 first: 8.671e-04 second: 1.599e-03 infty: 2.178e-02
+Divergence time: 0.256381 first: 8.655e-04 second: 1.593e-03 infty: 2.197e-02
+Divergence time: 0.262634 first: 8.711e-04 second: 1.605e-03 infty: 2.160e-02
+Divergence time: 0.268888 first: 8.901e-04 second: 1.648e-03 infty: 2.170e-02
+Divergence time: 0.275141 first: 8.513e-04 second: 1.557e-03 infty: 2.195e-02
+Divergence time: 0.281394 first: 8.556e-04 second: 1.562e-03 infty: 2.282e-02
+Divergence time: 0.287647 first: 8.579e-04 second: 1.571e-03 infty: 2.193e-02
+Divergence time: 0.2939 first: 8.618e-04 second: 1.580e-03 infty: 2.201e-02
+Divergence time: 0.300154 first: 8.586e-04 second: 1.565e-03 infty: 2.191e-02
+Divergence time: 0.306407 first: 8.649e-04 second: 1.614e-03 infty: 2.190e-02
+Divergence time: 0.31266 first: 8.605e-04 second: 1.597e-03 infty: 2.328e-02
+Divergence time: 0.318913 first: 8.517e-04 second: 1.553e-03 infty: 2.196e-02
+Divergence time: 0.325166 first: 8.568e-04 second: 1.585e-03 infty: 2.134e-02
+Divergence time: 0.331419 first: 8.503e-04 second: 1.558e-03 infty: 2.148e-02
+Divergence time: 0.337673 first: 8.507e-04 second: 1.562e-03 infty: 2.221e-02
+Divergence time: 0.343926 first: 8.548e-04 second: 1.580e-03 infty: 2.113e-02
+Divergence time: 0.350179 first: 8.473e-04 second: 1.553e-03 infty: 2.293e-02
+Divergence time: 0.356432 first: 8.459e-04 second: 1.569e-03 infty: 2.396e-02
+Divergence time: 0.362685 first: 8.520e-04 second: 1.573e-03 infty: 2.128e-02
+Divergence time: 0.368938 first: 8.492e-04 second: 1.553e-03 infty: 2.255e-02
+Divergence time: 0.375191 first: 8.490e-04 second: 1.549e-03 infty: 2.423e-02
+Divergence time: 0.381444 first: 8.450e-04 second: 1.548e-03 infty: 2.368e-02
+Divergence time: 0.387697 first: 8.450e-04 second: 1.581e-03 infty: 2.891e-02
+Divergence time: 0.39395 first: 8.461e-04 second: 1.570e-03 infty: 2.543e-02
+Divergence time: 0.400203 first: 8.447e-04 second: 1.550e-03 infty: 2.624e-02
+Divergence time: 0.406456 first: 8.428e-04 second: 1.538e-03 infty: 2.469e-02
+Divergence time: 0.412709 first: 8.432e-04 second: 1.548e-03 infty: 2.482e-02
+Divergence time: 0.418962 first: 8.333e-04 second: 1.497e-03 infty: 2.134e-02
+Divergence time: 0.425215 first: 8.379e-04 second: 1.531e-03 infty: 2.185e-02
+Divergence time: 0.431468 first: 8.336e-04 second: 1.504e-03 infty: 2.318e-02
+Divergence time: 0.437721 first: 8.308e-04 second: 1.501e-03 infty: 2.571e-02
+Divergence time: 0.443975 first: 8.382e-04 second: 1.522e-03 infty: 2.615e-02
+Divergence time: 0.450228 first: 8.374e-04 second: 1.506e-03 infty: 2.121e-02
+Divergence time: 0.456481 first: 8.283e-04 second: 1.459e-03 infty: 2.143e-02
+Divergence time: 0.462734 first: 8.275e-04 second: 1.483e-03 infty: 2.307e-02
+Divergence time: 0.468987 first: 8.303e-04 second: 1.488e-03 infty: 2.228e-02
+Divergence time: 0.47524 first: 8.289e-04 second: 1.464e-03 infty: 2.129e-02
+Divergence time: 0.481493 first: 8.291e-04 second: 1.468e-03 infty: 1.842e-02
+Divergence time: 0.487746 first: 8.324e-04 second: 1.467e-03 infty: 2.156e-02
+Divergence time: 0.493999 first: 8.336e-04 second: 1.471e-03 infty: 1.789e-02
+Divergence time: 0.500252 first: 8.324e-04 second: 1.471e-03 infty: 2.015e-02
+Divergence time: 0.506505 first: 8.325e-04 second: 1.476e-03 infty: 1.842e-02
+Divergence time: 0.512758 first: 8.310e-04 second: 1.469e-03 infty: 1.717e-02
+Divergence time: 0.519011 first: 8.264e-04 second: 1.458e-03 infty: 1.809e-02
+Divergence time: 0.525264 first: 8.206e-04 second: 1.442e-03 infty: 2.153e-02
+Divergence time: 0.531517 first: 8.114e-04 second: 1.442e-03 infty: 2.351e-02
+Divergence time: 0.53777 first: 8.093e-04 second: 1.436e-03 infty: 2.500e-02
+Divergence time: 0.544023 first: 8.107e-04 second: 1.458e-03 infty: 2.614e-02
+Divergence time: 0.550277 first: 8.118e-04 second: 1.458e-03 infty: 2.613e-02
+Divergence time: 0.55653 first: 8.155e-04 second: 1.459e-03 infty: 2.630e-02
+Divergence time: 0.562783 first: 8.186e-04 second: 1.446e-03 infty: 2.035e-02
+Divergence time: 0.569036 first: 8.179e-04 second: 1.456e-03 infty: 2.409e-02
+Divergence time: 0.575289 first: 8.196e-04 second: 1.461e-03 infty: 2.438e-02
+Divergence time: 0.581543 first: 8.258e-04 second: 1.490e-03 infty: 2.131e-02
+Divergence time: 0.587796 first: 8.260e-04 second: 1.487e-03 infty: 2.290e-02
+Divergence time: 0.594049 first: 8.315e-04 second: 1.483e-03 infty: 1.938e-02
+Divergence time: 0.600302 first: 8.320e-04 second: 1.489e-03 infty: 2.161e-02
+Divergence time: 0.606556 first: 8.404e-04 second: 1.516e-03 infty: 2.254e-02
+Divergence time: 0.612809 first: 8.332e-04 second: 1.478e-03 infty: 2.113e-02
+Divergence time: 0.619062 first: 8.366e-04 second: 1.481e-03 infty: 2.292e-02
+Divergence time: 0.625315 first: 8.335e-04 second: 1.465e-03 infty: 2.116e-02
+Divergence time: 0.631569 first: 8.313e-04 second: 1.468e-03 infty: 2.275e-02
+Divergence time: 0.637822 first: 8.309e-04 second: 1.474e-03 infty: 2.342e-02
+Divergence time: 0.644075 first: 8.295e-04 second: 1.470e-03 infty: 2.113e-02
+Divergence time: 0.650329 first: 8.257e-04 second: 1.468e-03 infty: 2.288e-02
+Divergence time: 0.656582 first: 8.230e-04 second: 1.461e-03 infty: 2.312e-02
+Divergence time: 0.662835 first: 8.225e-04 second: 1.459e-03 infty: 2.058e-02
+Divergence time: 0.669088 first: 8.176e-04 second: 1.440e-03 infty: 2.082e-02
+Divergence time: 0.675342 first: 8.178e-04 second: 1.430e-03 infty: 1.979e-02
+Divergence time: 0.681595 first: 8.182e-04 second: 1.423e-03 infty: 1.587e-02
+Divergence time: 0.687849 first: 8.136e-04 second: 1.415e-03 infty: 2.142e-02
+Divergence time: 0.694102 first: 8.100e-04 second: 1.411e-03 infty: 2.292e-02
+Divergence time: 0.700355 first: 8.053e-04 second: 1.408e-03 infty: 2.298e-02
+Divergence time: 0.706608 first: 8.052e-04 second: 1.400e-03 infty: 1.945e-02
+Divergence time: 0.712862 first: 8.014e-04 second: 1.395e-03 infty: 2.154e-02
+Divergence time: 0.719115 first: 7.949e-04 second: 1.379e-03 infty: 2.296e-02
+Divergence time: 0.725368 first: 7.907e-04 second: 1.362e-03 infty: 1.955e-02
+Divergence time: 0.731622 first: 7.856e-04 second: 1.345e-03 infty: 1.930e-02
+Divergence time: 0.737875 first: 7.806e-04 second: 1.345e-03 infty: 2.116e-02
+Divergence time: 0.744128 first: 7.804e-04 second: 1.341e-03 infty: 1.897e-02
+Divergence time: 0.750382 first: 7.783e-04 second: 1.330e-03 infty: 1.569e-02
+Divergence time: 0.756635 first: 7.806e-04 second: 1.346e-03 infty: 1.678e-02
+Divergence time: 0.762889 first: 7.773e-04 second: 1.343e-03 infty: 1.893e-02
+Divergence time: 0.769142 first: 7.748e-04 second: 1.341e-03 infty: 1.968e-02
+Divergence time: 0.775396 first: 7.750e-04 second: 1.339e-03 infty: 1.727e-02
+Divergence time: 0.781649 first: 7.728e-04 second: 1.339e-03 infty: 1.775e-02
+Divergence time: 0.787902 first: 7.713e-04 second: 1.341e-03 infty: 2.023e-02
+Divergence time: 0.794156 first: 7.714e-04 second: 1.346e-03 infty: 2.068e-02
+Divergence time: 0.800409 first: 7.731e-04 second: 1.346e-03 infty: 1.905e-02
+Divergence time: 0.806663 first: 7.723e-04 second: 1.340e-03 infty: 1.852e-02
+Divergence time: 0.812916 first: 7.756e-04 second: 1.362e-03 infty: 1.955e-02
+Divergence time: 0.81917 first: 7.739e-04 second: 1.352e-03 infty: 2.072e-02
+Divergence time: 0.825423 first: 7.747e-04 second: 1.341e-03 infty: 1.965e-02
+Divergence time: 0.831677 first: 7.784e-04 second: 1.355e-03 infty: 1.694e-02
+Divergence time: 0.83793 first: 7.808e-04 second: 1.378e-03 infty: 1.821e-02
+Divergence time: 0.844183 first: 7.803e-04 second: 1.379e-03 infty: 2.019e-02
+Divergence time: 0.850437 first: 7.804e-04 second: 1.362e-03 infty: 2.065e-02
+Divergence time: 0.85669 first: 7.780e-04 second: 1.343e-03 infty: 1.772e-02
+Divergence time: 0.862943 first: 7.849e-04 second: 1.401e-03 infty: 1.919e-02
+Divergence time: 0.869197 first: 7.869e-04 second: 1.406e-03 infty: 1.908e-02
+Divergence time: 0.87545 first: 7.847e-04 second: 1.392e-03 infty: 2.155e-02
+Divergence time: 0.881703 first: 7.833e-04 second: 1.363e-03 infty: 1.885e-02
+Divergence time: 0.887956 first: 7.850e-04 second: 1.380e-03 infty: 1.870e-02
+Divergence time: 0.894209 first: 7.790e-04 second: 1.354e-03 infty: 1.965e-02
+Divergence time: 0.900463 first: 7.837e-04 second: 1.372e-03 infty: 1.961e-02
+Divergence time: 0.906716 first: 7.858e-04 second: 1.393e-03 infty: 1.808e-02
+Divergence time: 0.912969 first: 7.829e-04 second: 1.365e-03 infty: 1.810e-02
+Divergence time: 0.919223 first: 7.793e-04 second: 1.341e-03 infty: 1.992e-02
+Divergence time: 0.925476 first: 7.868e-04 second: 1.387e-03 infty: 2.534e-02
+Divergence time: 0.931729 first: 7.893e-04 second: 1.410e-03 infty: 2.502e-02
+Divergence time: 0.937983 first: 7.837e-04 second: 1.391e-03 infty: 2.130e-02
+Divergence time: 0.944236 first: 7.812e-04 second: 1.366e-03 infty: 1.873e-02
+Divergence time: 0.950489 first: 7.823e-04 second: 1.383e-03 infty: 2.170e-02
+Divergence time: 0.956743 first: 7.828e-04 second: 1.392e-03 infty: 2.589e-02
+Divergence time: 0.962996 first: 7.830e-04 second: 1.391e-03 infty: 2.483e-02
+Divergence time: 0.969249 first: 7.841e-04 second: 1.392e-03 infty: 2.687e-02
+Divergence time: 0.975502 first: 7.854e-04 second: 1.402e-03 infty: 2.341e-02
+Divergence time: 0.981755 first: 7.861e-04 second: 1.396e-03 infty: 2.175e-02
+Divergence time: 0.988008 first: 7.882e-04 second: 1.400e-03 infty: 2.175e-02
+Divergence time: 0.994262 first: 7.895e-04 second: 1.410e-03 infty: 2.246e-02
+Divergence time: 1.00051 first: 7.877e-04 second: 1.381e-03 infty: 1.883e-02
+Divergence time: 1.00677 first: 7.875e-04 second: 1.386e-03 infty: 2.432e-02
+Divergence time: 1.01302 first: 7.836e-04 second: 1.370e-03 infty: 1.920e-02
+Divergence time: 1.01927 first: 7.844e-04 second: 1.392e-03 infty: 2.647e-02
+Divergence time: 1.02552 first: 7.852e-04 second: 1.389e-03 infty: 2.725e-02
+Divergence time: 1.03178 first: 7.843e-04 second: 1.384e-03 infty: 2.519e-02
+Divergence time: 1.03803 first: 7.828e-04 second: 1.387e-03 infty: 2.514e-02
+Divergence time: 1.04429 first: 7.826e-04 second: 1.385e-03 infty: 2.745e-02
+Divergence time: 1.05054 first: 7.891e-04 second: 1.406e-03 infty: 2.198e-02
+Divergence time: 1.05679 first: 7.894e-04 second: 1.403e-03 infty: 2.082e-02
+Divergence time: 1.06305 first: 7.854e-04 second: 1.379e-03 infty: 2.111e-02
+Divergence time: 1.0693 first: 7.881e-04 second: 1.387e-03 infty: 2.233e-02
+Divergence time: 1.07555 first: 7.887e-04 second: 1.395e-03 infty: 2.668e-02
+Divergence time: 1.0818 first: 7.870e-04 second: 1.394e-03 infty: 2.200e-02
+Divergence time: 1.08806 first: 7.805e-04 second: 1.338e-03 infty: 1.939e-02
+Divergence time: 1.09431 first: 7.805e-04 second: 1.342e-03 infty: 2.601e-02
+Divergence time: 1.10056 first: 7.777e-04 second: 1.336e-03 infty: 2.212e-02
+Divergence time: 1.10682 first: 7.789e-04 second: 1.335e-03 infty: 2.380e-02
+Divergence time: 1.11307 first: 7.799e-04 second: 1.336e-03 infty: 1.901e-02
+Divergence time: 1.11932 first: 7.784e-04 second: 1.336e-03 infty: 2.615e-02
+Divergence time: 1.12557 first: 7.734e-04 second: 1.318e-03 infty: 2.946e-02
+Divergence time: 1.13183 first: 7.766e-04 second: 1.320e-03 infty: 2.010e-02
+Divergence time: 1.13808 first: 7.758e-04 second: 1.319e-03 infty: 2.211e-02
+Divergence time: 1.14433 first: 7.745e-04 second: 1.318e-03 infty: 2.389e-02
+Divergence time: 1.15059 first: 7.670e-04 second: 1.279e-03 infty: 1.878e-02
+Divergence time: 1.15684 first: 7.645e-04 second: 1.267e-03 infty: 1.666e-02
+Divergence time: 1.16309 first: 7.647e-04 second: 1.272e-03 infty: 1.732e-02
+Divergence time: 1.16934 first: 7.651e-04 second: 1.291e-03 infty: 2.419e-02
+Divergence time: 1.1756 first: 7.637e-04 second: 1.277e-03 infty: 2.363e-02
+Divergence time: 1.18185 first: 7.621e-04 second: 1.274e-03 infty: 1.846e-02
+Divergence time: 1.1881 first: 7.636e-04 second: 1.279e-03 infty: 1.693e-02
+Divergence time: 1.19436 first: 7.603e-04 second: 1.254e-03 infty: 1.638e-02
+Divergence time: 1.20061 first: 7.563e-04 second: 1.249e-03 infty: 2.454e-02
+Divergence time: 1.20686 first: 7.552e-04 second: 1.244e-03 infty: 2.442e-02
+Divergence time: 1.21312 first: 7.554e-04 second: 1.254e-03 infty: 2.214e-02
+Divergence time: 1.21937 first: 7.556e-04 second: 1.249e-03 infty: 1.779e-02
+Divergence time: 1.22562 first: 7.568e-04 second: 1.256e-03 infty: 1.928e-02
+Divergence time: 1.23188 first: 7.575e-04 second: 1.251e-03 infty: 1.590e-02
+Divergence time: 1.23813 first: 7.559e-04 second: 1.246e-03 infty: 1.731e-02
+Divergence time: 1.24438 first: 7.543e-04 second: 1.239e-03 infty: 1.813e-02
+Divergence time: 1.25064 first: 7.533e-04 second: 1.246e-03 infty: 2.342e-02
+Divergence time: 1.25689 first: 7.514e-04 second: 1.248e-03 infty: 2.484e-02
+Divergence time: 1.26314 first: 7.534e-04 second: 1.250e-03 infty: 2.073e-02
+Divergence time: 1.2694 first: 7.512e-04 second: 1.232e-03 infty: 1.818e-02
+Divergence time: 1.27565 first: 7.509e-04 second: 1.224e-03 infty: 1.753e-02
+Divergence time: 1.2819 first: 7.493e-04 second: 1.211e-03 infty: 2.169e-02
+Divergence time: 1.28816 first: 7.475e-04 second: 1.214e-03 infty: 2.475e-02
+Divergence time: 1.29441 first: 7.489e-04 second: 1.224e-03 infty: 2.256e-02
+Divergence time: 1.30066 first: 7.488e-04 second: 1.224e-03 infty: 2.375e-02
+Divergence time: 1.30692 first: 7.481e-04 second: 1.227e-03 infty: 1.900e-02
+Divergence time: 1.31317 first: 7.488e-04 second: 1.217e-03 infty: 1.840e-02
+Divergence time: 1.31942 first: 7.489e-04 second: 1.217e-03 infty: 1.837e-02
+Divergence time: 1.32568 first: 7.502e-04 second: 1.225e-03 infty: 1.807e-02
+Divergence time: 1.33193 first: 7.511e-04 second: 1.227e-03 infty: 1.764e-02
+Divergence time: 1.33819 first: 7.474e-04 second: 1.202e-03 infty: 1.520e-02
+Divergence time: 1.34444 first: 7.449e-04 second: 1.201e-03 infty: 2.344e-02
+Divergence time: 1.35069 first: 7.467e-04 second: 1.214e-03 infty: 2.422e-02
+Divergence time: 1.35695 first: 7.470e-04 second: 1.214e-03 infty: 1.835e-02
+Divergence time: 1.3632 first: 7.482e-04 second: 1.216e-03 infty: 1.847e-02
+Divergence time: 1.36945 first: 7.463e-04 second: 1.208e-03 infty: 1.849e-02
+Divergence time: 1.37571 first: 7.429e-04 second: 1.197e-03 infty: 1.844e-02
+Divergence time: 1.38196 first: 7.426e-04 second: 1.200e-03 infty: 1.770e-02
+Divergence time: 1.38821 first: 7.419e-04 second: 1.189e-03 infty: 1.697e-02
+Divergence time: 1.39447 first: 7.395e-04 second: 1.177e-03 infty: 1.497e-02
+Divergence time: 1.40072 first: 7.375e-04 second: 1.176e-03 infty: 1.763e-02
+Divergence time: 1.40697 first: 7.392e-04 second: 1.186e-03 infty: 1.838e-02
+Divergence time: 1.41323 first: 7.396e-04 second: 1.193e-03 infty: 1.823e-02
+Divergence time: 1.41948 first: 7.374e-04 second: 1.181e-03 infty: 1.833e-02
+Divergence time: 1.42573 first: 7.358e-04 second: 1.172e-03 infty: 1.828e-02
+Divergence time: 1.43199 first: 7.344e-04 second: 1.172e-03 infty: 1.753e-02
+Divergence time: 1.43824 first: 7.338e-04 second: 1.172e-03 infty: 1.784e-02
+Divergence time: 1.44449 first: 7.324e-04 second: 1.166e-03 infty: 1.827e-02
+Divergence time: 1.45075 first: 7.333e-04 second: 1.174e-03 infty: 1.767e-02
+Divergence time: 1.457 first: 7.307e-04 second: 1.165e-03 infty: 1.768e-02
+Divergence time: 1.46325 first: 7.290e-04 second: 1.157e-03 infty: 1.838e-02
+Divergence time: 1.46951 first: 7.294e-04 second: 1.161e-03 infty: 1.848e-02
+Divergence time: 1.47576 first: 7.292e-04 second: 1.161e-03 infty: 1.843e-02
+Divergence time: 1.48201 first: 7.300e-04 second: 1.160e-03 infty: 1.837e-02
+Divergence time: 1.48827 first: 7.293e-04 second: 1.165e-03 infty: 1.833e-02
+Divergence time: 1.49452 first: 7.298e-04 second: 1.166e-03 infty: 1.830e-02
+Divergence time: 1.50077 first: 7.284e-04 second: 1.168e-03 infty: 2.116e-02
+Divergence time: 1.50703 first: 7.274e-04 second: 1.163e-03 infty: 2.465e-02
+Divergence time: 1.51328 first: 7.263e-04 second: 1.156e-03 infty: 2.021e-02
+Divergence time: 1.51953 first: 7.284e-04 second: 1.163e-03 infty: 2.010e-02
+Divergence time: 1.52579 first: 7.257e-04 second: 1.161e-03 infty: 2.502e-02
+Divergence time: 1.53204 first: 7.252e-04 second: 1.157e-03 infty: 2.168e-02
+Divergence time: 1.53829 first: 7.264e-04 second: 1.162e-03 infty: 2.009e-02
+Divergence time: 1.54455 first: 7.278e-04 second: 1.166e-03 infty: 2.025e-02
+Divergence time: 1.5508 first: 7.268e-04 second: 1.161e-03 infty: 2.169e-02
+Divergence time: 1.55705 first: 7.261e-04 second: 1.166e-03 infty: 2.524e-02
+Divergence time: 1.56331 first: 7.259e-04 second: 1.163e-03 infty: 2.082e-02
+Divergence time: 1.56956 first: 7.241e-04 second: 1.158e-03 infty: 2.276e-02
+Divergence time: 1.57581 first: 7.254e-04 second: 1.163e-03 infty: 2.433e-02
+Divergence time: 1.58207 first: 7.274e-04 second: 1.172e-03 infty: 2.676e-02
+Divergence time: 1.58832 first: 7.285e-04 second: 1.168e-03 infty: 2.463e-02
+Divergence time: 1.59457 first: 7.258e-04 second: 1.160e-03 infty: 2.453e-02
+Divergence time: 1.60083 first: 7.288e-04 second: 1.172e-03 infty: 2.723e-02
+Divergence time: 1.60708 first: 7.294e-04 second: 1.179e-03 infty: 2.951e-02
+Divergence time: 1.61333 first: 7.301e-04 second: 1.172e-03 infty: 2.116e-02
+Divergence time: 1.61959 first: 7.310e-04 second: 1.171e-03 infty: 2.663e-02
+Divergence time: 1.62584 first: 7.306e-04 second: 1.171e-03 infty: 2.778e-02
+Divergence time: 1.63209 first: 7.298e-04 second: 1.173e-03 infty: 2.966e-02
+Divergence time: 1.63835 first: 7.317e-04 second: 1.169e-03 infty: 2.280e-02
+Divergence time: 1.6446 first: 7.337e-04 second: 1.173e-03 infty: 1.916e-02
+Divergence time: 1.65085 first: 7.338e-04 second: 1.170e-03 infty: 2.355e-02
+Divergence time: 1.65711 first: 7.341e-04 second: 1.178e-03 infty: 2.375e-02
+Divergence time: 1.66336 first: 7.343e-04 second: 1.165e-03 infty: 2.281e-02
+Divergence time: 1.66961 first: 7.337e-04 second: 1.166e-03 infty: 1.645e-02
+Divergence time: 1.67587 first: 7.353e-04 second: 1.168e-03 infty: 1.723e-02
+Divergence time: 1.68212 first: 7.358e-04 second: 1.172e-03 infty: 1.901e-02
+Divergence time: 1.68837 first: 7.360e-04 second: 1.171e-03 infty: 1.689e-02
+Divergence time: 1.69463 first: 7.350e-04 second: 1.167e-03 infty: 1.552e-02
+Divergence time: 1.70088 first: 7.347e-04 second: 1.161e-03 infty: 1.468e-02
+Divergence time: 1.70713 first: 7.360e-04 second: 1.166e-03 infty: 1.646e-02
+Divergence time: 1.71339 first: 7.356e-04 second: 1.164e-03 infty: 1.645e-02
+Divergence time: 1.71964 first: 7.345e-04 second: 1.157e-03 infty: 1.444e-02
+Divergence time: 1.72589 first: 7.342e-04 second: 1.161e-03 infty: 1.481e-02
+Divergence time: 1.73215 first: 7.356e-04 second: 1.164e-03 infty: 1.654e-02
+Divergence time: 1.7384 first: 7.347e-04 second: 1.159e-03 infty: 1.563e-02
+Divergence time: 1.74465 first: 7.330e-04 second: 1.152e-03 infty: 1.424e-02
+Divergence time: 1.75091 first: 7.328e-04 second: 1.155e-03 infty: 1.391e-02
+Divergence time: 1.75716 first: 7.338e-04 second: 1.158e-03 infty: 1.559e-02
+Divergence time: 1.76342 first: 7.337e-04 second: 1.160e-03 infty: 1.597e-02
+Divergence time: 1.76967 first: 7.321e-04 second: 1.150e-03 infty: 1.581e-02
+Divergence time: 1.77592 first: 7.301e-04 second: 1.145e-03 infty: 1.358e-02
+Divergence time: 1.78218 first: 7.297e-04 second: 1.143e-03 infty: 1.378e-02
+Divergence time: 1.78843 first: 7.302e-04 second: 1.148e-03 infty: 1.541e-02
+Divergence time: 1.79468 first: 7.297e-04 second: 1.144e-03 infty: 1.466e-02
+Divergence time: 1.80094 first: 7.282e-04 second: 1.138e-03 infty: 1.318e-02
+Divergence time: 1.80719 first: 7.272e-04 second: 1.139e-03 infty: 1.326e-02
+Divergence time: 1.81344 first: 7.279e-04 second: 1.143e-03 infty: 1.480e-02
+Divergence time: 1.8197 first: 7.275e-04 second: 1.139e-03 infty: 1.430e-02
+Divergence time: 1.82595 first: 7.258e-04 second: 1.138e-03 infty: 1.468e-02
+Divergence time: 1.8322 first: 7.259e-04 second: 1.133e-03 infty: 1.320e-02
+Divergence time: 1.83846 first: 7.254e-04 second: 1.136e-03 infty: 1.426e-02
+Divergence time: 1.84471 first: 7.246e-04 second: 1.130e-03 infty: 1.321e-02
+Divergence time: 1.85096 first: 7.241e-04 second: 1.133e-03 infty: 1.374e-02
+Divergence time: 1.85722 first: 7.243e-04 second: 1.134e-03 infty: 1.413e-02
+Divergence time: 1.86347 first: 7.234e-04 second: 1.128e-03 infty: 1.371e-02
+Divergence time: 1.86973 first: 7.226e-04 second: 1.129e-03 infty: 1.409e-02
+Divergence time: 1.87598 first: 7.222e-04 second: 1.124e-03 infty: 1.322e-02
+Divergence time: 1.88223 first: 7.215e-04 second: 1.126e-03 infty: 1.322e-02
+Divergence time: 1.88849 first: 7.218e-04 second: 1.127e-03 infty: 1.349e-02
+Divergence time: 1.89474 first: 7.214e-04 second: 1.124e-03 infty: 1.322e-02
+Divergence time: 1.90099 first: 7.203e-04 second: 1.120e-03 infty: 1.322e-02
+Divergence time: 1.90725 first: 7.195e-04 second: 1.119e-03 infty: 1.322e-02
+Divergence time: 1.9135 first: 7.197e-04 second: 1.121e-03 infty: 1.322e-02
+Divergence time: 1.91975 first: 7.195e-04 second: 1.123e-03 infty: 1.322e-02
+Divergence time: 1.92601 first: 7.186e-04 second: 1.114e-03 infty: 1.322e-02
+Divergence time: 1.93226 first: 7.177e-04 second: 1.115e-03 infty: 1.322e-02
+Divergence time: 1.93851 first: 7.181e-04 second: 1.117e-03 infty: 1.322e-02
+Divergence time: 1.94477 first: 7.179e-04 second: 1.119e-03 infty: 1.322e-02
+Divergence time: 1.95102 first: 7.178e-04 second: 1.119e-03 infty: 1.323e-02
+Divergence time: 1.95728 first: 7.166e-04 second: 1.110e-03 infty: 1.322e-02
+Divergence time: 1.96353 first: 7.153e-04 second: 1.109e-03 infty: 1.322e-02
+Divergence time: 1.96978 first: 7.156e-04 second: 1.110e-03 infty: 1.322e-02
+Divergence time: 1.97604 first: 7.159e-04 second: 1.113e-03 infty: 1.323e-02
+Divergence time: 1.98229 first: 7.161e-04 second: 1.115e-03 infty: 1.323e-02
+Divergence time: 1.98854 first: 7.152e-04 second: 1.109e-03 infty: 1.324e-02
+Divergence time: 1.9948 first: 7.148e-04 second: 1.110e-03 infty: 1.324e-02
+Divergence time: 2 first: 6.304e-04 second: 9.776e-04 infty: 1.130e-02
diff --git a/test/euler/stationary1/stationary1.gfs b/test/euler/reynolds/reynolds.gfs
similarity index 90%
rename from test/euler/stationary1/stationary1.gfs
rename to test/euler/reynolds/reynolds.gfs
index 9ad508d..8574161 100644
--- a/test/euler/stationary1/stationary1.gfs
+++ b/test/euler/reynolds/reynolds.gfs
@@ -27,9 +27,9 @@
# \end{figure}
#
# Author: St\'ephane Popinet
-# Command: sh stationary1.sh
+# Command: sh reynolds.sh reynolds.gfs
# Version: 0.6.4
-# Required files: stationary1.sh
+# Required files: reynolds.sh div5.ref div6.ref div7.ref reynolds.ref check.py
# Running time: 3 minutes
# Generated files: divmax.eps reynolds.eps divL2.eps kinetic.eps
#
diff --git a/test/euler/reynolds/reynolds.ref b/test/euler/reynolds/reynolds.ref
new file mode 100644
index 0000000..72eb102
--- /dev/null
+++ b/test/euler/reynolds/reynolds.ref
@@ -0,0 +1,3 @@
+5 9908.61
+6 79507.8
+7 630992
diff --git a/test/euler/reynolds/reynolds.sh b/test/euler/reynolds/reynolds.sh
new file mode 100644
index 0000000..0096413
--- /dev/null
+++ b/test/euler/reynolds/reynolds.sh
@@ -0,0 +1,49 @@
+rm -f reynolds
+
+for level in 5 6 7; do
+ sed "s/LEVEL/$level/" < $1 | gerris2D - | awk -v m=1 -v level=$level '{
+ time = $3
+ ke = $5
+ if (time == 0)
+ ke0 = ke;
+ }END{
+ a = -log(ke/ke0)/time
+ nu = a/(4.*(2.*m*3.14159265359)^2)
+ print level " " 1./nu
+ }' >> reynolds
+done
+
+cat <<EOF | gnuplot
+ set term postscript eps color solid 20
+ set output 'divmax.eps'
+ set xlabel 'Time'
+ set ylabel 'Divergence Max'
+ plot [0:2]'div5' u 3:9 t "5" w l lw 2, 'div6' u 3:9 t "6" w l lw 2, 'div7' u 3:9 t "7" w l lw 2, 'div5.ref' u 3:9 t "5 (ref)" w l lw 2, 'div6.ref' u 3:9 t "6 (ref)" w l lw 2, 'div7.ref' u 3:9 t "7 (ref)" w l lw 2
+ set output 'divL2.eps'
+ set ylabel 'Divergence L2'
+ plot [0:2]'div5' u 3:7 t "5" w l lw 2, 'div6' u 3:7 t "6" w l lw 2, 'div7' u 3:7 t "7" w l lw 2, 'div5.ref' u 3:7 t "5 (ref)" w l lw 2, 'div6.ref' u 3:7 t "6 (ref)" w l lw 2, 'div7.ref' u 3:7 t "7 (ref)" w l lw 2
+ set output 'kinetic.eps'
+ set ylabel 'Kinetic energy'
+ plot [0:2]'kinetic5' u 3:5 t "5" w l lw 2, 'kinetic6' u 3:5 t "6" w l lw 2, 'kinetic7' u 3:5 t "7" w l lw 2
+ set output 'reynolds.eps'
+ set xlabel 'Level'
+ set ylabel 'Effective Reynolds number'
+ set logscale y
+ plot 'reynolds' u 1:2 t "" w lp lw 2, 'reynolds.ref' u 1:2 t "ref" w lp lw 2
+EOF
+
+if cat <<EOF | python ; then :
+from check import *
+from sys import *
+if (Curve('div5',3,9) - Curve('div5.ref',3,9)).max() > 0.05*Curve('div5.ref',3,9).mean() or\
+ (Curve('div6',3,9) - Curve('div6.ref',3,9)).max() > 0.05*Curve('div6.ref',3,9).mean() or\
+ (Curve('div7',3,9) - Curve('div7.ref',3,9)).max() > 0.05*Curve('div7.ref',3,9).mean() or\
+ (Curve('div5',3,7) - Curve('div5.ref',3,7)).max() > 0.05*Curve('div5.ref',3,7).mean() or\
+ (Curve('div6',3,7) - Curve('div6.ref',3,7)).max() > 0.05*Curve('div6.ref',3,7).mean() or\
+ (Curve('div7',3,7) - Curve('div7.ref',3,7)).max() > 0.05*Curve('div7.ref',3,7).mean() or\
+ (Curve('reynolds',1,2) - Curve('reynolds.ref',1,2)).max() > 0.0 :
+ exit(1)
+EOF
+else
+ exit 1
+fi
diff --git a/test/euler/stationary1/stationary1.sh b/test/euler/stationary1/stationary1.sh
deleted file mode 100644
index 75f8c63..0000000
--- a/test/euler/stationary1/stationary1.sh
+++ /dev/null
@@ -1,36 +0,0 @@
-rm -f reynolds
-
-for level in 5 6 7; do
- sed "s/LEVEL/$level/" < stationary1.gfs | gerris2D - | awk -v m=1 -v level=$level '{
- time = $3
- ke = $5
- if (time == 0)
- ke0 = ke;
- }END{
- a = -log(ke/ke0)/time
- nu = a/(4.*(2.*m*3.14159265359)^2)
- print level " " 1./nu
- }' >> reynolds
-done
-
-cat <<EOF | gnuplot
- set term postscript eps
- set output 'divmax.eps'
- set xlabel 'Time'
- set ylabel 'Divergence Max'
- plot [0:2]'div5' u 3:9 t "5" w l, 'div6' u 3:9 t "6" w l, 'div7' u 3:9 t "7" w l
- set output 'divL2.eps'
- set ylabel 'Divergence L2'
- plot [0:2]'div5' u 3:7 t "5" w l, 'div6' u 3:7 t "6" w l, 'div7' u 3:7 t "7" w l
- set output 'kinetic.eps'
- set ylabel 'Kinetic energy'
- plot [0:2]'kinetic5' u 3:5 t "5" w l, 'kinetic6' u 3:5 t "6" w l, 'kinetic7' u 3:5 t "7" w l
- set output 'reynolds.eps'
- set xlabel 'Level'
- set ylabel 'Effective Reynolds number'
- set nokey
- set logscale y
- plot 'reynolds' u 1:2 w lp
-EOF
-
-rm -f div? kinetic? reynolds
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list