Re: [Getfem-users] Coordinates of gauss points

[Yves Renard]

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