[Getdp] solving trivial example using edge elements

[William McLean]

I am new to getdp and want to understand how to use 3D edge finite
elements.  As a first step, I tried solving the easiest problem I could
think of: find the vector field u satisfying

	(u,v) + (curl u,curl v) = (f,v)  for all v in H_0(curl,Omega),

        u = g on the Boundary of Omega,

taking a constant source term f=(1,0,0) and constant Dirichlet data
g=(1,0,0), so that the exact solution is the constant vector field
u=(1,0,0).  For the domain Omega, I chose just the unit cube (0,1)^3.

Unfortunately, my computed solution uh is nothing like the exact
solution.

The attached tar file holds the geometric description file cube.geo,
the getdp program edge_elts.pro, the solver parameter file solver.par
and a 3-line shell script run.sh.  Typing

> sh run.sh

will generate the mesh, perform the FE discretization and display
the computed solution uh.

I hope someone can point out where I am going wrong.

Bill McLean.





-- 
William McLean
School of Mathematical Sciences
University of Bath
Bath BA2 7AY, United Kingdom
Tel: +44 (0) 1225 386989
Fax: +44 (0) 1225 386492
-------------- next part --------------
A non-text attachment was scrubbed...
Name: edge_elts.tgz
Type: application/x-compressed-tar
Size: 36196 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20070514/e447da97/attachment.tgz>