Hi, I am trying to use PCBDDC for the vector based FEM. (Complex system, double precision ) My code can work well with *asm*, petsc-3.21.1/petsc/arch-linux-c-opt/bin/mpirun -n 8 ./app -pc_type asm -pc_asm_overlap 6 -ksp_converged_reason -ksp_view -ksp_gmres_modifiedgramschmidt -ksp_gmres_restart 1500 -ksp_rtol 1e-8 -ksp_monitor -ksp_max_it 100000
When I tried BDDC, it was stuck for solving the linear system (it can not print anything for ksp_monitor). I did the conversion for matrix, * Mat J;* * MatConvert(A, MATIS, MAT_INITIAL_MATRIX, &J);* * KSPSetOperators(ksp, A, J);* * MatDestroy(&J);* * KSPSetInitialGuessNonzero(ksp, PETSC_TRUE);* * KSPSetFromOptions(ksp);* petsc-3.21.1/petsc/arch-linux-c-debug/bin/mpirun -n 2 ./app -ksp_type cg -pc_type bddc -ksp_monitor -mat_type is Do you have any suggestions? Thanks , Xiaodong On Mon, Jul 29, 2024 at 6:19 PM neil liu <liufi...@gmail.com> wrote: > When I compile with real data, > it shows the attached error. > > The data file is in binary format, right? > > > > On Mon, Jul 29, 2024 at 5:36 PM Stefano Zampini <stefano.zamp...@gmail.com> > wrote: > >> Your PETSc installation is for complex, data is for real >> >> On Mon, Jul 29, 2024, 23:14 neil liu <liufi...@gmail.com> wrote: >> >>> I compiled Petsc with single precision. However, it is not converged >>> with the data. Please see the attached file. On Mon, Jul 29, 2024 at 4: 25 >>> PM Barry Smith <bsmith@ petsc. dev> wrote: This can happen if the data >>> was stored in single precision >>> ZjQcmQRYFpfptBannerStart >>> This Message Is From an External Sender >>> This message came from outside your organization. >>> >>> ZjQcmQRYFpfptBannerEnd >>> I compiled Petsc with single precision. However, it is not converged >>> with the data. >>> >>> Please see the attached file. >>> >>> On Mon, Jul 29, 2024 at 4:25 PM Barry Smith <bsm...@petsc.dev> wrote: >>> >>>> >>>> This can happen if the data was stored in single precision and PETSc >>>> was built for double. >>>> >>>> >>>> On Jul 29, 2024, at 3:55 PM, neil liu <liufi...@gmail.com> wrote: >>>> >>>> This Message Is From an External Sender >>>> This message came from outside your organization. >>>> Dear Petsc developers,, >>>> >>>> I am trying to run >>>> https://urldefense.us/v3/__https://petsc.org/release/src/ksp/ksp/tutorials/ex72.c.html__;!!G_uCfscf7eWS!f-0Lu-su_W4yzz3WY5WMTUfWrVt3g96PdhasAO6-CtYzQzHqw-HviyclvMgCjclVqn6LQ6joNJmKan04GFddJg$ >>>> >>>> <https://urldefense.us/v3/__https://petsc.org/release/src/ksp/ksp/tutorials/ex72.c.html__;!!G_uCfscf7eWS!ZG4gvmS6hQD8ymbvCUDfAatzRUJHzmWO-hOgp9m0xXuAXgIB-fxe_xspYs3WEPi_Ed0UFLMHKanYuYWrTlQGrA$> >>>> with >>>> >>>> petsc-3.21.1/petsc/arch-linux-c-opt/bin/mpirun -n 2 ./ex72 -f >>>> /Documents/PetscData/poisson_DMPLEX_32x32_16.dat -pc_type bddc -ksp_type cg >>>> -ksp_norm_type natural -ksp_error_if_not_converged -mat_type is >>>> >>>> The file was downloaded and put in the directory PetscData. >>>> >>>> The error is shown as follows, >>>> >>>> 0]PETSC ERROR: --------------------- Error Message >>>> -------------------------------------------------------------- >>>> [0]PETSC ERROR: Read from file failed >>>> [0]PETSC ERROR: Read past end of file >>>> [0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>>> program crashed before usage or a spelling mistake, etc! >>>> [0]PETSC ERROR: Option left: name:-ksp_error_if_not_converged (no >>>> value) source: command line >>>> [0]PETSC ERROR: Option left: name:-ksp_norm_type value: natural >>>> source: command line >>>> [0]PETSC ERROR: Option left: name:-ksp_type value: cg source: command >>>> line >>>> [0]PETSC ERROR: Option left: name:-pc_type value: bddc source: >>>> command line >>>> [0]PETSC ERROR: See >>>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!f-0Lu-su_W4yzz3WY5WMTUfWrVt3g96PdhasAO6-CtYzQzHqw-HviyclvMgCjclVqn6LQ6joNJmKan22rsyOwQ$ >>>> >>>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!ZG4gvmS6hQD8ymbvCUDfAatzRUJHzmWO-hOgp9m0xXuAXgIB-fxe_xspYs3WEPi_Ed0UFLMHKanYuYVR0x14Xg$> >>>> for trouble shooting. >>>> [0]PETSC ERROR: Petsc Release Version 3.21.1, unknown >>>> [0]PETSC ERROR: ./ex72 on a arch-linux-c-opt named >>>> Mon Jul 29 15:50:04 2024 >>>> [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran >>>> --with-cxx=g++ --download-fblaslapack --download-mpich >>>> --with-scalar-type=complex --download-triangle --with-debugging=no >>>> [0]PETSC ERROR: #1 PetscBinaryRead() at >>>> /home/xxxxxx/Documents/petsc-3.21.1/petsc/src/sys/fileio/sysio.c:327 >>>> [0]PETSC ERROR: #2 PetscViewerBinaryWriteReadAll() at >>>> /home/xiaodong.liu/Documents/petsc-3.21.1/petsc/src/sys/classes/viewer/impls/binary/binv.c:1077 >>>> [0]PETSC ERROR: #3 PetscViewerBinaryReadAll() at >>>> /home/xiaodong.liu/Documents/petsc-3.21.1/petsc/src/sys/classes/viewer/impls/binary/binv.c:1119 >>>> [0]PETSC ERROR: #4 MatLoad_MPIAIJ_Binary() at >>>> Documents/petsc-3.21.1/petsc/src/mat/impls/aij/mpi/mpiaij.c:3093 >>>> [0]PETSC ERROR: #5 MatLoad_MPIAIJ() at >>>> /Documents/petsc-3.21.1/petsc/src/mat/impls/aij/mpi/mpiaij.c:3035 >>>> [0]PETSC ERROR: #6 MatLoad() at >>>> /Documents/petsc-3.21.1/petsc/src/mat/interface/matrix.c:1344 >>>> [0]PETSC ERROR: #7 MatLoad_IS() at >>>> /Documents/petsc-3.21.1/petsc/src/mat/impls/is/matis.c:2575 >>>> [0]PETSC ERROR: #8 MatLoad() at >>>> /home/Documents/petsc-3.21.1/petsc/src/mat/interface/matrix.c:1344 >>>> [0]PETSC ERROR: #9 main() at ex72.c:105 >>>> [0]PETSC ERROR: PETSc Option Table entries: >>>> [0]PETSC ERROR: -f >>>> /Documents/PetscData/poisson_DMPLEX_32x32_16.dat (source: command line) >>>> [0]PETSC ERROR: -ksp_error_if_not_converged (source: command line) >>>> [0]PETSC ERROR: -ksp_norm_type natural (source: command line) >>>> [0]PETSC ERROR: -ksp_type cg (source: command line) >>>> [0]PETSC ERROR: -mat_type is (source: command line) >>>> [0]PETSC ERROR: -pc_type bddc (source: command line) >>>> [0]PETSC ERROR: ----------------End of Error Message -------send entire >>>> error message to petsc-ma...@mcs.anl.gov---------- >>>> application called MPI_Abort(MPI_COMM_SELF, 66) - process 0 >>>> >>>> >>>>
