[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Coordinates of gauss points
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] Coordinates of gauss points |
Date: |
Wed, 11 Feb 2015 12:24:36 +0100 (CET) |
Dear Eleonora,
If you just want to print the coordinate of the Gauss point inside de generic
assembly you can just add a
Print(X) inside your assembly string.
If you do not want to change the result of your computation your can insert it
inside an instuction having no effect such as
... *(1+Norm(Print(X)-X))
or something similar.
Yves.
----- Original Message -----
From: "Eleonora Piersanti" <address@hidden>
To: "Yves Renard" <address@hidden>
Cc: address@hidden
Sent: Wednesday, February 11, 2015 11:49:18 AM
Subject: Re: [Getfem-users] Coordinates of gauss points
Thank you for the answer.
Since I need to assamble a complicate term I am using the high level generic
assembly. For debugging purpose, I am using the Print command, as suggested in
the guide, that print on screen the value of a quantity in each Gauss Point. I
would recover the coordinates of those points so that I can understand if what
I am doing is correct.
I tried with fem_interpolation_context but I have problems to understand how to
call the constructor properly. Is there a way to save X variable in a vector
with high level features?
Sorry, maybe the question is trivial.
Kind regards,
Eleonora
Il Mercoledì 11 Febbraio 2015 8:16, Yves Renard <address@hidden> ha
scritto:
Dear Eleonora,
It is of course possible. It depends in which context you need them.
In the high level generic assembly, the current position (i.e. the coordinates
of the current Gauss point) is simply denoted by X.
In te internal of getfem, when elementary computations occurs, there is a
special object named "fem_interpolation_context"
. If ctx is such an object, ctx.xref() is the coordinate of the Gauss point on
the reference element and ctx.xreal() on the real one.
If you need the Gauss points on each element, you can loop over the elements
and the Gauss points inside an element and using a fem_interpolation_context
object have the coordinate of each Gauss point on the real element.
Don't forget also that there is a specific structure named "im_data" which is
specially dedicated to store some data on each Gauss points and which can be
used in generic assembly.
Yves.
----- Original Message -----
From: "Eleonora Piersanti" <address@hidden>
To: address@hidden
Sent: Tuesday, February 10, 2015 7:49:00 PM
Subject: [Getfem-users] Coordinates of gauss points
Hi to everyone,
is it possible to obtain the coordinates of the Gauss points?
Thanks,
Eleonora
_______________________________________________
Getfem-users mailing list
address@hidden
https://mail.gna.org/listinfo/getfem-users