Re: [petsc-users] MPI barrier issue using MatZeroRows

2023-11-28 Thread Amneet Bhalla
I added that option but the code still gets stuck at the same call MatZeroRows with 3 processors. On Tue, Nov 28, 2023 at 7:23 PM Amneet Bhalla wrote: > > > On Tue, Nov 28, 2023 at 6:42 PM Barry Smith wrote: > >> >> for (int comp = 0; comp < 2; ++comp) >> { >>

Re: [petsc-users] MPI barrier issue using MatZeroRows

2023-11-28 Thread Amneet Bhalla
On Tue, Nov 28, 2023 at 6:42 PM Barry Smith wrote: > > for (int comp = 0; comp < 2; ++comp) > { > ... > for (Box::Iterator bc(bc_coef_box); bc; bc++) > { >.. > if

Re: [petsc-users] MPI barrier issue using MatZeroRows

2023-11-28 Thread Barry Smith
for (int comp = 0; comp < 2; ++comp) { ... for (Box::Iterator bc(bc_coef_box); bc; bc++) { .. if (IBTK::abs_equal_eps(b, 0.0)) {

Re: [petsc-users] [EXTERNAL] Re: Unexpected performance losses switching to COO interface

2023-11-28 Thread Jed Brown
"Fackler, Philip via petsc-users" writes: > That makes sense. Here are the arguments that I think are relevant: > > -fieldsplit_1_pc_type redundant -fieldsplit_0_pc_type sor -pc_type fieldsplit > -pc_fieldsplit_detect_coupling​ What sort of physics are in splits 0 and 1? SOR is not a good GPU

[petsc-users] MPI barrier issue using MatZeroRows

2023-11-28 Thread Amneet Bhalla
Hi Folks, I am using MatZeroRows() to set Dirichlet boundary conditions. This works fine for the serial run and the solver produces correct results (verified through analytical solution). However, when I run the case in parallel, the simulation gets stuck at MatZeroRows(). My understanding is

Re: [petsc-users] Fortran problem MatGetValuesLocal

2023-11-28 Thread Mark Adams
It looks like we don't have a Fortran interface for this. In the source code: 2366: Developer Note:2367: This is labeled with C so does not automatically generate Fortran stubs and interfaces2368: because it requires multiple Fortran interfaces depending on which arguments are scalar or

Re: [petsc-users] Fortran problem MatGetValuesLocal

2023-11-28 Thread Barry Smith
This is fixed in branch barry/2023-11-28/add-matsetvalueslocal-fortran/release see also https://gitlab.com/petsc/petsc/-/merge_requests/7065 > On Nov 28, 2023, at 8:46 AM, Frank Bramkamp wrote: > > Dear PETSc team, > > > We are using the latest petsc version 3.20.1, intel compiler 2023,

[petsc-users] Fortran problem MatGetValuesLocal

2023-11-28 Thread Frank Bramkamp
Dear PETSc team, We are using the latest petsc version 3.20.1, intel compiler 2023, and we found the following problem: We want to call the function MatGetValuesLocal to extract a block sub-matrix from an assembled matrix (e.g. a 5x5 blocked sub matrix). We use the matrix format MatCreateBAIJ

[petsc-users] msh4 geenerated from neper

2023-11-28 Thread Sharan Roongta
Hello, I have problem loading the msh4 file generated from neper. The mesh format is 4.1 The error is: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Unexpected data in file [0]PETSC ERROR: File is not a valid