Re: [Getfem-users] Coordinates of gauss points

[Eleonora Piersanti]

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