Mercurial > fem-fenics-eugenio

annotate inst/example/NavierStokes/lshape-domain.m @ 174:e85390ed620d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add copyright.
author gedeone-octave <marcovass89@hotmail.it>
date Sat, 26 Oct 2013 17:24:56 +0100
parents d03d998e64af
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
174
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
1 ## Copyright (C) 2013 Marco Vassallo <gedeone-octave@users.sourceforge.net>
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
2 ##
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
3 ## This program is free software; you can redistribute it and/or modify it under
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
4 ## the terms of the GNU General Public License as published by the Free Software
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
5 ## Foundation; either version 3 of the License, or (at your option) any later
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
6 ## version.
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
7 ##
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
8 ## This program is distributed in the hope that it will be useful, but WITHOUT
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
9 ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
10 ## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
11 ## details.
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
12 ##
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
13 ## You should have received a copy of the GNU General Public License along with
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
14 ## this program; if not, see <http://www.gnu.org/licenses/>.
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
15
e85390ed620d Add copyright.
gedeone-octave <marcovass89@hotmail.it>
parents: 159
diff changeset
16
76
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
17 name = [tmpnam ".geo"];
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
18 fid = fopen (name, "w");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
19 fputs (fid,"Point (1) = {0, 0, 0, 0.1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
20 fputs (fid,"Point (2) = {1, 0, 0, 0.1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
21 fputs (fid,"Point (3) = {1, 0.5, 0, 0.1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
22 fputs (fid,"Point (4) = {0.5, 0.5, 0, 0.1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
23 fputs (fid,"Point (5) = {0.5, 1, 0, 0.1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
24 fputs (fid,"Point (6) = {0, 1, 0,0.1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
25
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
26 fputs (fid,"Line (1) = {5, 6};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
27 fputs (fid,"Line (2) = {2, 3};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
28
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
29 fputs (fid,"Line(3) = {6,1,2};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
30 fputs (fid,"Line(4) = {5,4,3};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
31 fputs (fid,"Line Loop(7) = {3,2,-4,1};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
32 fputs (fid,"Plane Surface(8) = {7};\n");
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
33 fclose (fid);
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
34 msho = msh2m_gmsh (canonicalize_file_name (name)(1:end-4),"scale", 1,"clscale", .2);
102
d880c19f691b Add comment to the script and pot the mesh.
gedeone-octave <marcovass89@hotmail.it>
parents: 76
diff changeset
35 trimesh (msho.t(1:3,:)', msho.p(1,:)', msho.p(2,:)');
76
6cb2a3e3be68 New examples which solve vector field problems
gedeone-octave <marcovass89@hotmail.it>
parents:
diff changeset
36 unlink (canonicalize_file_name (name));