[Getfem-users] Quadrature question.

[Ronan Perrussel]

Dear getfem users,

I am using getfem for solving an axisymmetric problem and I have to 
integrate a term of the following form :
\int_\Omega 1/r u v dr dz
where u and v are respectively solution and test function, r is the 
radial distance and z is the height.
Of course u and v are enforced to be zero on the axis. Consequently the 
integrand should be regular even for the elements touching the axis.
I wonder how to perform an accurate quadrature of this term particularly 
for the elements touching the axis because I do not know how to describe 
the function 1/r properly by a mesh_fem object (only in the area close 
to the axis).
I could do something by hand for the elements touching the axis but it 
would be helpful if someone had a better proposition.

Thanks in advance for your proposition(s).
Best regards,

Ronan

ronan_perrussel.vcf
Description: Vcard