Re: [petsc-users] Filling non-zero values of a Petsc matrix using numpy arrays with non-zero indices and values (petsc4py)

2023-08-18 Thread Zhang, Hong via petsc-users
You can use this to build a PETSc matrix with the index arrays ai,aj and the value array aa: PETSc.Mat().createAIJ(size=(nrows,ncols), csr=(ai,aj,aa)) Hong (Mr.) > On Aug 17, 2023, at 7:37 AM, Erik Kneller via petsc-users > wrote: > > Hi All, > > I need to fill non-zero values of a

Re: [petsc-users] Configure AMGX

2023-08-18 Thread Zisheng Ye
Thanks, that branch works. From: Satish Balay Sent: Friday, August 18, 2023 11:01 AM To: Zisheng Ye Cc: PETSc users list Subject: Re: [petsc-users] Configure AMGX [External Sender] Can you try the update in branch "balay/amgx-cuda-12"? Satish On Fri, 18 Aug

Re: [petsc-users] Configure AMGX

2023-08-18 Thread Satish Balay via petsc-users
Can you try the update in branch "balay/amgx-cuda-12"? Satish On Fri, 18 Aug 2023, Zisheng Ye wrote: > Dear PETSc team > > I am configuring AMGX package under the main branch with CUDA 12.1. But it > can't get through. Can you help to solve the problem? I have attached the > configure.log to

Re: [petsc-users] PetscCall( ) in fortran

2023-08-18 Thread Satish Balay via petsc-users
I think gfortran defaults to fixed form for .F and free-form for .F90 This can be changed with FFLAGS=-ffree-form - but yeah - switching the suffix might be more suitable.. In addition - PETSc configure attempts to add in "-ffree-line-length-none -ffree-line-length-0" options - so that extra