From: li pan <li76pan@ya...>  20060327 14:26:18

Deal all, I met a strang problem. Maybe you can also try it. Just write a simple code of building a mesh and defining an equation_system, adding some variables and so on. int main (int argc, char** argv) { // Initialize Petsc, like in example 2. libMesh::init (argc, argv); const unsigned int dim = 3; // Create a dimdimensional mesh. Mesh mesh (dim); std::cout<<DIM; MeshTools::Generation::build_cube(mesh, 1,2,2, 0.,1., 0.,1., 0.,1., HEX8); mesh.print_info(); EquationSystems equation_systems (mesh); TransientLinearImplicitSystem & system = equation_systems.add_system<TransientLinearImplicitSystem>("non_linear"); system.add_variable("u",FIRST); system.add_variable("v",FIRST); system.add_variable("w",FIRST); system.attach_assemble_function(assemble_nonlinear ); equation_systems.init(); equation_systems.parameters.set<unsigned int>("linear solver maximum iterations")=1000; equation_systems.parameters.set<Real>("linear solver tolerance")=1.e10; } That's it. Compile and run the program. You'll get error report: [0]PETSC ERROR: MatSeqAIJSetPreallocation_SeqAIJ() line 2631 in src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: Argument out of range! [0]PETSC ERROR: nnz cannot be greater than row length: local row 6 value 69 rowlength 54! [0]PETSC ERROR: MatCreateSeqAIJ() line 2544 in src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: User provided function() line 138 in unknowndirectory/src/numerics/petsc_matrix.C [unset]: aborting job: application called MPI_Abort(comm=0x84000000, 63)  process 0 But if I use MeshTools::Generation::build_cube(mesh, 1,2,1, 0.,1., 0.,1., 0.,1., HEX8); or some other grid dimensions, the problem disappear. Do you know why? I configured my petsc with downloadfblaslapack=1 downloadmpich=1 withshared=1 pan __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com 