Run your code with option "-eps_view -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_package mumps -mat_mumps_icntl_28 2 -mat_mumps_icntl_29 2" and send us entire output of '-eps_view'.
Hong On Mon, Jun 2, 2014 at 2:13 AM, 张国熙 <[email protected]> wrote: > Hi, all > > Recently I'm trying to improve the performance of my program, so I add -help > option to see options available. I think I use MPIAIJ type of MAT, but I > receive information of SEQAIJ. > > Is it means actually I am using a sequential matrix type? > If not, how could I get information of MPIAIJ to improve performance? > > Code of generating mate > MatCreate(PETSC_COMM_WORLD,&A); > MatSetType(A,MATMPIAIJ); > MatSetSizes(A,PETSC_DECIDE, PETSC_DECIDE,size,size); > MatMPIAIJSetPreallocationCSR(A,Ap,Ai,temp); > > Code of getting information > char common_options[] = "-help -st_ksp_type preonly -st_pc_type lu > -st_pc_factor_mat_solver_package mumps -mat_mumps_icntl_28 2 > -mat_mumps_icntl_29 2"; > > ierr = PetscOptionsInsertString(common_options);CHKERRQ(ierr); > > Thank you very much > Guoxi
