~ubuntu-branches/ubuntu/saucy/rheolef/saucy-proposed

« back to all changes in this revision

Viewing changes to nfem/ptst/my_cube_PH_gmsh.mshcad

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito
  • Date: 2012-04-06 09:12:21 UTC
  • mfrom: (1.1.5)
  • Revision ID: package-import@ubuntu.com-20120406091221-m58me99p1nxqui49
Tags: 6.0-1
* New upstream release 6.0 (major changes):
  - massively distributed and parallel support
  - full FEM characteristic method (Lagrange-Gakerkin method) support
  - enhanced users documentation 
  - source code supports g++-4.7 (closes: #667356)
* debian/control: dependencies for MPI distributed solvers added
* debian/rules: build commands simplified
* debian/librheolef-dev.install: man1/* to man9/* added
* debian/changelog: package description rewritted (closes: #661689)

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
//  gmsh -2 my_cube_PH_gmsh.mshcad -o my_cube_PH-10.msh
 
2
//  ../pbin/msh2geo < my_cube_PH-10.msh | ../pbin/geo -upgrade -geo - > my_cube_PH-10.geo
 
3
//  Omega = ]a,b[x]c,d[x]f,g[
 
4
    a = 0; c = 0; f = 0;
 
5
    b = 1; d = 1; g = 1;
 
6
    n = 5;
 
7
    hloc = 1.0/n;
 
8
    Point(1) = {a, c, f, hloc};
 
9
    Point(2) = {b, c, f, hloc};
 
10
    Point(6) = {b, c, g, hloc};
 
11
    Point(5) = {a, c, g, hloc};
 
12
    Line(1) = {1,2};
 
13
    Line(6) = {2,6};
 
14
    Line(9) = {5,6};
 
15
    Line(5) = {1,5};
 
16
    Line Loop(300) = {1,6,-9,-5};
 
17
    Mesh.Algorithm = 1;
 
18
    Plane Surface(3) = {300} ;
 
19
    angle = 90.0;
 
20
    Mesh.RecombinationAlgorithm = 0;
 
21
    Recombine Surface {3} = angle; //  mixt qua & tri
 
22
    extr[] = Extrude{0,d-c,0}{ Surface{3}; Layers{n}; Recombine;};
 
23
    Physical Surface("boundary")= {-extr[2], -extr[4], 3, -extr[3], -extr[0], -extr[5]};
 
24
    Physical Surface("bottom")  = {-extr[2]};
 
25
    Physical Surface("top")     = {-extr[4]};
 
26
    Physical Surface("left")    = {3};
 
27
    Physical Surface("front")   = {-extr[3]};
 
28
    Physical Surface("right")   = {-extr[0]};
 
29
    Physical Surface("back")    = {-extr[5]};
 
30
    Physical Volume("interior") = {extr[1]};