[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
popinet at users.sf.net
Fri May 15 02:55:37 UTC 2009
The following commit has been merged in the upstream branch:
commit db89e08b3f57eddbe66769e7623845a0542e0af6
Author: Stephane Popinet <popinet at users.sf.net>
Date: Wed Jul 23 14:10:41 2008 +1000
Cleaner axi test case parameter file
darcs-hash:20080723041041-d4795-11848af15297c3bb412acade92b160157de1e5c8.gz
diff --git a/test/axi/axi.gfs b/test/axi/axi.gfs
index 902eacf..17d5052 100644
--- a/test/axi/axi.gfs
+++ b/test/axi/axi.gfs
@@ -39,13 +39,14 @@
# Generated files: error.eps order.eps isolines.eps
#
1 0 GfsAxi GfsBox GfsGEdge {} {
- Time { iend = 0 end = 1 }
+ Time { end = 0 }
+ PhysicalParams { L = 50 }
AdvectionParams { scheme = none }
ApproxProjectionParams { tolerance = 1e-10 }
Refine 4
- Refine (LEVEL + 50.*(x*x + y*y)*(4. - LEVEL))
+ Refine (LEVEL + 1./50.*(x*x + y*y)*(4. - LEVEL))
Global {
- #define A0 0.01
+ #define A0 0.5
#define U0 1.
}
Solid (ellipse (0., 0., A0, A0))
@@ -53,8 +54,7 @@
U = U0
Phi = {
double r = sqrt (cx*cx + cy*cy);
- double theta = atan2 (cy, cx);
- return U0*A0*A0*A0*cos(theta)/(2.*r*r);
+ return U0*A0*A0*A0*cx/(2.*r*r*r);
}
}
OutputErrorNorm { start = end } { awk '{ print LEVEL " " $7 " " $9}' } { v = U } {
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list