Re: [Libmesh-users] example parallelmesh???

[Robert . Bodner]

Hello, 

I haven't tried the suggestion of Derek yet. Instead, I
wanted to try to split a smaller mesh with nemesis and see
if my code runs in parallel (configure with
--enable-parmesh)

I split the mesh with the following code: 
Mesh mesh;
                  
                TetGenIO TETGEN(mesh);
                TETGEN.read("meshfile.ele");

                
                set_subdomains(&TETGEN, mesh);   
            set_boundary_ids(mesh); 
                std::cout<<"tetgen finished reading"<<std::endl;
                
                mesh.find_neighbors();
                
                std::cout<<"neighbors found"<<std::endl;
                
                MetisPartitioner().partition(mesh,4);
                
                mesh.print_info();
                
                std::cout<<"partition done"<<std::endl;
        
        
                Nemesis_IO NEM ( mesh );
                
                std::cout<<"write"<<std::endl;
                
                {  
                OStringStream outmesh;
                outmesh << "mesh1/" <<
WriteMesh[actual_pulse]<<".nemesis";
                NEM.write(outmesh.str().c_str());
                } 

------------------------------------------------------------------------------
RSA&#174; Conference 2012
Save $700 by Nov 18
Register now&#33;
http://p.sf.net/sfu/rsa-sfdev2dev1
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users