[petsc-users] Use MPIAIJ, but receive help information of seqaij

[Hong Zhang]

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, 张国熙 <altriaex86 at gmail.com> 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