Hello, I just stumbled over a bug in DofMap which pops up when using MONOMIAL with --node_major_dofs. To reproduce just add "--node_major_dofs" to the options of introduction/ex2. This results at my build in:
Running ./example-dbg eqn_sys.dat --node_major_dofs Assertion `dofid < this->end_dof(proc_id)' failed. dofid = 72 this->end_dof(proc_id) = 72 Stack frames: 10 0: libMesh::print_trace(std::ostream&) 1: libMesh::MacroFunctions::report_error(char const*, int, char const*, char const*) 2: libMesh::DofMap::distribute_dofs(libMesh::MeshBase&) 3: libMesh::System::init_data() 4: libMesh::ExplicitSystem::init_data() 5: libMesh::System::init() 6: libMesh::EquationSystems::init() 7: ./example-dbg() [0x417964] 8: __libc_start_main 9: ./example-dbg() [0x4170b1] I wanted to use a constant monomial to have two element based variables to implement a three-field-variational formulation for incompressibility. Is there any better type available in libmesh for this? Robert ------------------------------------------------------------------------------ October Webinars: Code for Performance Free Intel webinars can help you accelerate application performance. Explore tips for MPI, OpenMP, advanced profiling, and more. Get the most from the latest Intel processors and coprocessors. See abstracts and register > http://pubads.g.doubleclick.net/gampad/clk?id=60134071&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
