Hello,
I’m experiencing a runtime error when export the equation systems if I use more 
than one EquationSystems with different variables in each system.
In particular, I created 2 EquationSystems: systems1 with one P1 variable and 
systems2  with two P1 variables. After initializing, I write the 
EquationSystems on two separate exodus files
and here it breaks at runtime.
If I initialize the systems with the same number of variables , everything 
works fine.
Do you have any suggestions?
I attaching the output error and the simple code I used.
Thanks,
Simone

// Output

Creating Systems1:
Initializing s1
Creating Systems2:
Initializing s2
Exporting s1
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably 
memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see 
http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to 
find memory corruption errors
[0]PETSC ERROR: likely location of problem given in stack below
[0]PETSC ERROR: ---------------------  Stack Frames 
------------------------------------
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,
[0]PETSC ERROR:       INSTEAD the line number of the start of the function
[0]PETSC ERROR:       is given.
[0]PETSC ERROR: --------------------- Error Message 
--------------------------------------------------------------
[0]PETSC ERROR: Signal received
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for 
trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.7.2, Jun, 05, 2016
[0]PETSC ERROR: ./test_libmesh_example on a linux-opt named aorta by srossi Fri 
Aug 19 22:48:40 2016
[0]PETSC ERROR: Configure options 
--prefix=/not_backed_up/srossi/TPL/petsc/3.7.2 --PETSC_ARCH=linux-opt 
--with-default-arch=0 --with-c++-support --with-hypre=1 --download-hypre=1 
--with-x=0
[0]PETSC ERROR: #1 User provided function() line 0 in  unknown file



// main.cpp

#include "libmesh/libmesh.h"
#include "libmesh/mesh.h"
#include "libmesh/mesh_generation.h"
#include "libmesh/equation_systems.h"
#include "libmesh/explicit_system.h"
#include "libmesh/exodusII_io.h"

int main (int argc, char ** argv)
{
       using namespace libMesh;
LibMeshInit init (argc, argv);
Mesh mesh(init.comm());
std::cout << "Creating Mesh ..." << std::endl;
MeshTools::Generation::build_cube (mesh,
                                      2 , 2 , 2,
                                      0., 1.,
                                      0., 1.,
                                      0., 1.,
                                      TET4);

   std::cout << "Creating Systems1: " << std::endl;
   EquationSystems systems1(mesh);
   ExplicitSystem& s1 = systems1.add_system<ExplicitSystem>("s1");
   s1.add_variable("s1", FIRST);
   std::cout << "Initializing s1" << std::endl;
   systems1.init();

   std::cout << "Creating Systems2: " << std::endl;
   EquationSystems systems2(mesh);
   ExplicitSystem& s2 = systems2.add_system<ExplicitSystem>("s2");
   s2.add_variable("s2a", FIRST);
   s2.add_variable("s2b", FIRST); // Commenting this line below it runs ok
   std::cout << "Initializing s2" << std::endl;
   systems2.init();

   std::cout << "Exporting s1" << std::endl;
   ExodusII_IO(systems1.get_mesh()).write_equation_systems("s1.e", systems1);
   std::cout << "Exporting s2" << std::endl;
   ExodusII_IO(systems2.get_mesh()).write_equation_systems("s2.e", systems2);

   return 0;
}

------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to