[Gmsh] gmsh fails to create points inside a loop

David Colignon David.Colignon at ulg.ac.be
Wed Jun 2 07:47:06 CEST 2010


On 02/06/10 07:42, jean pierre aubry wrote:

> hello
>
> well,  that's exactly  it
> thanks!!
>
> i will have to be a little bit more clever in the design of the loops to
> create the geometry in this case
>
> it would be much easier  if i wecould tell
>
> Line (10x) = Rotate [ .... {Duplicata { Line {100};}]
> enforcing the value 10x instead of letting Gmsh use newl


Hi Jean-Pierre,

it was possible before, but it created far more problems and had many side effects...

>
> a way of retrieving the number of the two ends of a line would be useful too


http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Transformations

Boundary { transform-list }
     (Not a transformation per-se.) Returns the boundary of the elementary entities in transform-list.

with

      transform-list:
        Point | Line | Surface | Volume { expression-list }; ... |
        Duplicata { Point | Line | Surface | Volume { expression-list }; ... } |
        transform

Regards,

Dave


-- 
David Colignon, Ph.D.
Collaborateur Logistique du F.R.S.-FNRS
CÉCI - Consortium des Équipements de Calcul Intensif
ACE - Applied & Computational Electromagnetics
Sart-Tilman B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
WWW:    http://hpc.montefiore.ulg.ac.be/
Agenda: http://www.google.com/calendar/embed?src=david.colignon%40gmail.com

>
> thanks again
>
> jean pierre aubry
>
> ----- Original Message -----
> From: "Eric Nutsch"<ericnutsch at gmail.com>
> To: "jean pierre aubry"<jeanpierre at lamachine.fr>
> Cc: "gmsh"<gmsh at geuz.org>
> Sent: Wednesday, June 02, 2010 5:07 AM
> Subject: Re: [Gmsh] gmsh fails to create points inside a loop
>
>
>> Though I havent had time to review your code explicitly, I know one
>> possible stumble...
>>
>> Gmsh can not create two points at the exact same point in space(in
>> addition to not having the same identifier). This is a likely mistake
>> when using loops, and I find myself running into that problem
>> frequently.
>>
>> Hope that helps,
>> Eric
>>
>>
>>
>>
>> On Tue, Jun 1, 2010 at 8:44 AM, jean pierre aubry
>> <jeanpierre at lamachine.fr>  wrote:
>>> hello
>>>
>>> i have the following loop to create an axis symetric struture
>>>
>>>
>>> For i In {0:num:pas}
>>>
>>> .........................
>>> Point(10000*i
>>>
> +100538)={-5800*Cos(a*i)-(0*Sin(a*i)),-5800*Sin(a*i)+(0*Cos(a*i)),6060,e};
>>>
>>> Point(10000*i
>>>
> +100701)={-4435*Cos(a*i)-(30*Sin(a*i)),-4435*Sin(a*i)+(30*Cos(a*i)),8766,e};
>>> Point(10000*i
>>>
> +100702)={-4435*Cos(a*i)-(-30*Sin(a*i)),-4435*Sin(a*i)+(-30*Cos(a*i)),8766,e
> };
>>>
>>> ..........................
>>>
>>> EndFor
>>>
>>> for i =0 all 3 points are created
>>>
>>> for the next st ps the point Point(10000*i+100701) is not created
>>>
>>> and this happens at a few other places within the loop which contains
>>> about 150 point to be repeated 16 times
>>>
>>> the point can however be created outside of the loop with the same
>>> parameter
>>>
>>> any idea of what's happening
>>>
>>> i wnat to create a structure similar to this one which was created with
>>> gmsh, many thanks!!
>>>
>>> http://www.code-aster.org/V2/spip.php?article298
>>>
>>> but much more complicated and largelly parametric
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> gmsh mailing list
>>> gmsh at geuz.org
>>> http://www.geuz.org/mailman/listinfo/gmsh
>>>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh