[Libmesh-users] 1-D mesh

[John P. <pet...@cf...>]

That is a good question.  In developing the library, we
often skipped over the implementation of algorithms for
1D meshes.  I believe at the time we wanted to do "real"
problems as quickly as possible, while handling 1D problems
in Matlab if they came along.


As a first cut to why the find_neighbors() routine won't
currently work in 1D, the Elem::side() member, which returns a
pointer to an Elem, throws an error for 1D elements,
see include/geom/edge.h

  /**
   * The \p Elem::side() member makes no sense for edges.
   */
  AutoPtr<Elem> side (const unsigned int) const
  { error(); AutoPtr<Elem> ap(NULL); return ap; }


I have no doubt that find_neighbors() can and should be
implemented for 1D meshes, and I have a feeling it should
be pretty easy.  I'll try to work up something in the near
future.

-John

Manav Bhatia writes:
 > Please ignore my previous mail with subject line "1-D
 > mesh".
 > 
 > Here is an addendum to my question on the
 > find_neighbours. 
 > 
 > Basically I am trying to make a 1-D mesh, but the code
 > aborts with an error due to this condition in
 > MeshBase::find_neighbours().
 > 
 >   if (_dim == 1)
 >     error();
 > 
 > Is there a way to make a 1-D mesh that I may be
 > missing? 
 > 
 > Thanks
 > Manav