[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