[Getfem-users] slicer_half_space problems (bug or feature?)

[Roman Putanowicz]

Dear All,

I have a problem with slicer_half_space : 
I would like produce a line slice over a 2D mesh.
Everything works, fields are interpolated, etc.
but the actual number of edges and nodes is
twice as it should be. The edges are just doubled
and have the same orientation.

I have called slicer.merge_nodes() but without
any result.

I would be grateful if you could tell me if
this is GetFEM++ feature (and why so) or a simple bug?

In the attachment there is a simple program which 
produces  quad mesh in a rectangle and then a 
slice of the mesh.

Thank you in advance,

Regards

Roman 

PS: In the program I write out mesh quality because
otherwise the file is empty. Is there any other way
export a mesh without a field defined over it?
-- 
Roman Putanowicz, PhD  < address@hidden  >
Institute for Computational Civil Engng (L-5)
Dept. of Civil Engng, Cracow Univ. of Technology
www.l5.pk.edu.pl, tel. +48 12 628 2569, fax 2034

line_slice.cxx
Description: Text Data