getfem-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Getfem-users] Help regarding calculation of (n-1) dimensional integ


From: Arvind Ajoy
Subject: Re: [Getfem-users] Help regarding calculation of (n-1) dimensional integrals in an (n) dim problem
Date: Thu, 7 Aug 2008 17:07:43 +0530

Dear Ronan

Thanks for the suggestion.

I use the matlab interface to getfem. When I try the equivalent interface command to the one
you suggested,
s1 = gf_slice({'planar',+1,[0;0;0],[1;0;0]},m,0)
I get an error that   "Argument 3 is out of bounds : 0 not in [1...1000]". It looks like the minimum refinement parameter is 1, not 0.

Thanks again for you help.

Regards
Arvind

On Wed, Aug 6, 2008 at 8:29 PM, Ronan Perrussel <address@hidden> wrote:
Dear Arvind,

Maybe it will not solve your problem but do you try to adjust the number of refinement (I suppose if you set this value to zero it will be ok)...
See here :
http://download.gna.org/getfem/doc/getfemuser/getfemuser_18.html#id2
with mesh_slicer::exec(int nrefine)

I hope it helps...

Ronan

Arvind Ajoy a écrit :
Dear Getfem users

I am solving Poisson's equation self-consistently with Schrodinger's equation in a cylindrical
geometry. I use Gmsh to create the mesh, which generates an extruded prismatic mesh. The
axis of the cylinder is oriented along the x-direction.

>From a solution U(x,y,z) obtained on this mesh, I need to compute the following quantities
(in TeX notation)
1. \int \int _{yz} U(x_o,y,z) dy dz     , at different locations x_o
2. \int _{x}  U(x,y,z) dx

I tried using  mesh_slice; however, I am unable to generate a slice without mesh refinement!

Does anyone have any ideas as to how I should go about solving this problem?

Look forward to your replies

Regards
Arvind
------------------------------------------------------------------------

_______________________________________________
Getfem-users mailing list
address@hidden
https://mail.gna.org/listinfo/getfem-users
 



reply via email to

[Prev in Thread] Current Thread [Next in Thread]