Good day,
I have a program, that uses getfem::new_interpolated_fem. The call to
this
functions occurs many times in the program. The program, actually, simulates
a movement of conductor in a electrostatic field and I use a call to this
function to enforce Dirichlet conditions with largrange multiplier along the
boundary of the conductor. So, for some meshes (not all and not always) when
a conductor reaches a certain position the call to
getfem::pfem ifem=getfem::new_interpolated_fem(mf_electric,mim_line);
hangs. Hangs for a very long time (for at least a night, while it shouldn't
take even a second for other positions of the same conductor). It doesn't
produce an error and uses all the computing power, like it entered an
infinite loop.
in the above case, for instance,
mf_electric is FEM_QK(2,1), built on the mesh with GT_QK(2,1), using
getfem::regular_unit_mesh (with noise ON) and mim_line is IM_GAUSS1D(10).
Do you, gentlemen, see any reason why this function would hang?
Regards,
Andriy
--
__________________________________________
Andriy Andreykiv (PhD, MSc)
Delft University of Technology
Faculty of Mechanical Engineering
Fundamentals of Microsystems
Mekelweg 2
2628 CD Delft
The Netherlands
E-mail : [EMAIL PROTECTED]
Tel. : (31) 15 2786818
Fax. : (31) 15 2789475
www : http://www-tm.wbmt.tudelft.nl/~andrico
private: (31) 6 47376804
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
