Re: [Libmesh-devel] GmshIO::read

2008-01-28 Thread Roy Stogner
On Mon, 28 Jan 2008, David Knezevic wrote: > I came across a problem in GmshIO::read. The issue is that in some cases a > gmsh mesh contains lower dimensional elements which are just used for > specifying boundary info so they aren't added to the libMesh mesh during > mesh.read() and therefore

[Libmesh-devel] GmshIO::read

2008-01-28 Thread David Knezevic
I came across a problem in GmshIO::read. The issue is that in some cases a gmsh mesh contains lower dimensional elements which are just used for specifying boundary info so they aren't added to the libMesh mesh during mesh.read() and therefore assert( ... || e->id() == _elements.size()); in Ser

Re: [Libmesh-devel] ParallelMesh example

2008-01-28 Thread Roy Stogner
On Mon, 28 Jan 2008, Umar Farooq wrote: > The two matrices have a different non-zero structure...entries in diferent > indicies. Oh, I understand - yes, this should be different depending on how many processors you use. libMesh numbers degrees of freedom based on its element numbering, which in

Re: [Libmesh-devel] ParallelMesh example

2008-01-28 Thread Roy Stogner
On Mon, 28 Jan 2008, Umar Farooq wrote: > I was wondering if anyone has an example of how to use the ParallelMesh Ideally, all the examples will eventually work with ParallelMesh, and the only concern for the user will be to decide whether their code is better run with the memory overhead of Ser

[Libmesh-devel] ParallelMesh example

2008-01-28 Thread Umar Farooq
Hey, I was wondering if anyone has an example of how to use the ParallelMesh with an assemble function that loops through element nodes for imposing boundary conditions, after the element matrix has been added to the system matrix. Ex9 was a good template for my problem a while ago, but Im not su