I've attached a test code that freezes on the call
to mesh.sub_point_locator(). You can get the mesh from here
<https://dl.dropboxusercontent.com/u/10916994/test_mesh.exo> (I figured the
list would strip the mesh attachment).

The mesh is pretty messy (I didn't make it!) but still this seems like a
bug to me.

Let me know if anyone else can reproduce this?

Thanks,
David
// C++ include files that we need
#include <iostream>

// libMesh include
#include "libmesh/libmesh.h"
#include "libmesh/mesh.h"
#include "libmesh/equation_systems.h"

// Bring in everything from the libMesh namespace
using namespace libMesh;

// Begin the main program.
int main (int argc, char** argv)
{
  // Initialize libMesh and any dependent libraries
  LibMeshInit init (argc, argv);

  // Create a 3D mesh distributed across the default MPI communicator.
  Mesh mesh(init.comm());
  mesh.read("test_mesh.exo");

  // Print information about the mesh to the screen.
  mesh.print_info();

  AutoPtr<PointLocatorBase> point_locator = mesh.sub_point_locator();

  // All done.
  return 0;
}
------------------------------------------------------------------------------
Dive into the World of Parallel Programming The Go Parallel Website, sponsored
by Intel and developed in partnership with Slashdot Media, is your hub for all
things parallel software development, from weekly thought leadership blogs to
news, videos, case studies, tutorials and more. Take a look and join the 
conversation now. http://goparallel.sourceforge.net/
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to