<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
  <title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
I guess you are using an old version of GetDP. The treatment of
attractors has changed.<br>
Try e.g. the nightly-build.<br>
Ruth<br>
<br>
Lei Zhang wrote:
<blockquote
 cite="mid:44d29c40811192341x5587c3c2qa6104183c93a5272@mail.gmail.com"
 type="cite">
  <pre wrap="">Hello,

When I saved your example in a cube.geo file, and tried gmsh cube.geo,
I got the following error,
Error   : 'cube.geo', line 25 : parse error (Field)
Error   : 'cube.geo', line 26 : parse error (Field)
Error   : 'cube.geo', line 27 : parse error (Field)
Error   : 'cube.geo', line 30 : parse error (Field)
Error   : 'cube.geo', line 31 : parse error (Field)
Error   : 'cube.geo', line 35 : parse error (Field)

What's wrong?

Thanks,
Lei


On Thu, Nov 20, 2008 at 3:03 PM, Christophe Geuzaine
<a class="moz-txt-link-rfc2396E" href="mailto:cgeuzaine@ulg.ac.be"><cgeuzaine@ulg.ac.be></a> wrote:
  </pre>
  <blockquote type="cite">
    <pre wrap="">Srinath Madhavan wrote:
    </pre>
    <blockquote type="cite">
      <pre wrap="">Hi all,

Is it possible to define a virtual line/curve in the center of a unit cube and get gmsh to refine the mesh all around the line/curve and coarsen the mesh as it moves away from the line toward the walls?

In gambit, I would simply create a virtual line/curve and apply a size function to it. How would I do this in gmsh?
      </pre>
    </blockquote>
    <pre wrap="">Here's an example:

// a cube
lc = 0.3;
Point(1) = {0, 0, 0, lc};
Point(2) = {1, 0,  0, lc} ;
Point(3) = {1, 1, 0, lc} ;
Point(4) = {0, 1, 0, lc} ;
Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;
Line Loop(5) = {4,1,-2,3} ;
Plane Surface(6) = {5} ;
Extrude {0,0,1} { Surface{6}; }

// the "virtual" curve
Point(15) = {0.1, 0.3, 0.5};
Point(16) = {0.2, 0.3, 0.3};
Point(17) = {0.3, 0.2, 0.3};
Point(18) = {0.4, 0.1, 0.4};
Point(19) = {0.5, 0.2, 0.4};
Spline(29) = {15, 16, 17, 18, 19};

// Attractor field returns the distance to the curve (actually, the
// distance to 100 equidistant points on the curve)
Field[1] = Attractor;
Field[1].NNodesByEdge = 100;
Field[1].EdgesList = {29};

// Matheval field returns "distance squared + lc/20"
Field[2] = MathEval;
Field[2].F = Sprintf("F1^2 + %g", lc / 20);

// Use the minimum of "distance squared + lc/20" and the char lengths
// defined at geometrical points as the new char length
Background Field = 2;




    </pre>
    <blockquote type="cite">
      <pre wrap="">Thanks & Best Regards,
Srinath Madhavan

Free as in Freedom. May the source be with you!




_______________________________________________
gmsh mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a>
<a class="moz-txt-link-freetext" href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a>


      </pre>
    </blockquote>
    <pre wrap="">
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
<a class="moz-txt-link-freetext" href="http://www.montefiore.ulg.ac.be/~geuzaine">http://www.montefiore.ulg.ac.be/~geuzaine</a>

_______________________________________________
gmsh mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a>
<a class="moz-txt-link-freetext" href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a>

    </pre>
  </blockquote>
  <pre wrap=""><!---->
_______________________________________________
gmsh mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gmsh@geuz.org">gmsh@geuz.org</a>
<a class="moz-txt-link-freetext" href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a>
  </pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">-- 
Dr. Ir. Ruth V. Sabariego
University of Liege, Institut Montefiore,
Dept. of Electrical Engineering & Computer Science, 
Applied & Computational Electromagnetics (ACE),
Sart Tilman Campus, Grande Traverse, 10 (B28), B-4000 LIEGE, Belgium
phone: +32-4-3663737 -- fax: +32-4-3662910 -- <a class="moz-txt-link-freetext" href="http://elap.montefiore.ulg.ac.be/">http://elap.montefiore.ulg.ac.be/</a></pre>
</body>
</html>