<div dir="ltr"><div>Hello everyone,</div><div>I am using slepc4py solving a generalized eigenvalue problem, with a large spase matrix.  The following is my python code. The matrix is first created using PETSc.Mat().createAIJ() and saved to file matrix-A.dat . </div><div><br></div><div><br></div><div><div>   opts = PETSc.Options()</div><div>   opts.setValue('-st_pc_factor_mat_solver_package','mumps')</div></div><div><br></div><div><div>    viewer = PETSc.Viewer().createBinary('matrix-A.dat', 'r',comm=PETSc.COMM_SELF)</div><div>    A = PETSc.Mat().create(comm=PETSc.COMM_SELF)</div><div>    A.load(viewer)</div><div>    A.setFromOptions()</div><div>    A.setUp()</div><div>       </div><div>    Print('load petsc binary finished')</div><div>    </div><div>    t0 = time()</div><div>    </div><div>    E = SLEPc.EPS().create(comm=PETSc.COMM_SELF)</div><div>    E.setOperators(A,None)<br></div><div>    E.setDimensions(noEigs2solv,PETSc.DECIDE)           # set number of  eigenvalues to compute</div><div>    E.setWhichEigenpairs(E.Which.TARGET_MAGNITUDE)</div><div>    E.setProblemType(SLEPc.EPS.ProblemType.GNHEP)    # generalized non-Hermitian eigenvalue problem</div><div>    </div><div>    E.setTolerances(arnoldiTol,1000)     # set tolerance </div><div>    if SIGMA is not None:</div><div>        E.setTarget(SIGMA)               # set the desired eigenvalue </div><div>    </div><div>    if options_v0 is not None:</div><div>        E.setInitialSpace(options_v0)    # set the initial vector</div><div>       </div><div>    st = E.getST()</div><div>    st.setType('sinvert')</div><div>           </div><div>    E.setFromOptions()</div><div>    Print('Start solving ...')</div><div>    E.solve()</div></div><div><br></div><div>My problem here is, when I run this code on cluster, I would end up with the following error. Strange thing is that I used the save library on my local PC, and it works without any error.  If any one of you could give me hint what is wrong?  Thanks</div><br><div><br></div><div><div> ** ERROR RETURN ** FROM ZMUMPS INFO(1)=   -1</div><div> ** INFO(2)=               0</div><div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range</div><div>[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger</div><div>[0]PETSC ERROR: or see <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind</a></div><div>[0]PETSC ERROR: or try <a href="http://valgrind.org">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors</div><div>[0]PETSC ERROR: likely location of problem given in stack below</div><div>[0]PETSC ERROR: ---------------------  Stack Frames ------------------------------------</div><div>[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,</div><div>[0]PETSC ERROR:       INSTEAD the line number of the start of the function</div><div>[0]PETSC ERROR:       is given.</div><div>[0]PETSC ERROR: [0] PetscInitializeMUMPS line 1188 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/petsc/src/mat/impls/aij/mpi/mumps/mumps.c</div><div>[0]PETSC ERROR: [0] MatGetFactor_aij_mumps line 2125 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/petsc/src/mat/impls/aij/mpi/mumps/mumps.c</div><div>[0]PETSC ERROR: [0] MatGetFactor line 4341 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/petsc/src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: [0] PCSetUp_LU line 59 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/petsc/src/ksp/pc/impls/factor/lu/lu.c</div><div>[0]PETSC ERROR: [0] PCSetUp line 886 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/petsc/src/ksp/pc/interface/precon.c</div><div>[0]PETSC ERROR: [0] KSPSetUp line 294 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/petsc/src/ksp/ksp/interface/itfunc.c</div><div>[0]PETSC ERROR: [0] STSetUp_Sinvert line 96 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/slepc-3.8.2/src/sys/classes/st/impls/sinvert/sinvert.c</div><div>[0]PETSC ERROR: [0] STSetUp line 233 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/slepc-3.8.2/src/sys/classes/st/interface/stsolve.c</div><div>[0]PETSC ERROR: [0] EPSSetUp line 104 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/slepc-3.8.2/src/eps/interface/epssetup.c</div><div>[0]PETSC ERROR: [0] EPSSolve line 129 /zhome/academic/HLRS/iag/iagyonwu/.local/lib/slepc-3.8.2/src/eps/interface/epssolve.c</div><div>application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0</div><div>[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=59</div></div><div><br></div><div>with regards</div><div>Yong</div></div>