Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code

2024-05-07 Thread Matthew Knepley
On Tue, May 7, 2024 at 2:23 PM Shatanawi, Sawsan Muhammad via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hello everyone, I hope this email finds you well. My Name is Sawsan > Shatanawi, and I was developing a Fortran code for simulating groundwater > flow in a 3D system with nonlinear behavio

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Mark Adams
Yes, easy to solve. It is a kind of Lapalcian that comes up in graph algorithms. You mentioned problems with pinned points, with zero displacement. Conceptually you want to create A'A without any modification for pinned points. Then remove the rows and columns associated with pinned points. A'A ha

Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code

2024-05-07 Thread Barry Smith
Have you considered https://urldefense.us/v3/__https://www.pflotran.org/documentation/user_guide/how_to/installation/installation.html__;!!G_uCfscf7eWS!eZkBlApyezAyEBaDxGemEyUmTe87omvkv59rrf0Mq8L4bOTqEL9wynuMJ1ci9kRDOqucZBiYdGT2tRTVuCSJdSs$ > On May 7, 2024, at 2:22 PM, Shatanawi, Sawsan M

Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code

2024-05-07 Thread Shatanawi, Sawsan Muhammad via petsc-users
Hello everyone, I hope this email finds you well. My Name is Sawsan Shatanawi, and I was developing a Fortran code for simulating groundwater flow in a 3D system with nonlinear behavior. I solved the nonlinear system using the PCG solver and Picard iteration, but I did not get good resul

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Pierre Jolivet
OK, it’s very trivial to solve with algebraic solvers. If you are willing to share larger test cases, maybe then some issue will arise (please do not attach them to your mail, send a URL, if you cannot and must attach it to your mail, switch to the following mailing list: petsc-ma...@mcs.anl.gov

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Seiz,Marco
The system basically describes the displacements of particles relative to each other due to processes occurring on their interfaces, i. e. vacancy absorption. The idea is to say there's a displacement jump u_a - u_b = du_{a,b}, for every interface ZjQcmQRYFpfptBannerStart

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Seiz,Marco
Pierre, I've attached the dumps of the matrix + RHS for something of about 3k x 1k. Regarding the weird divergence behaviour, I tried again at home but I still get the same results. I am running a rolling release distribution on both machines, ZjQcmQRYFpfptBannerStart T

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Matthew Knepley
On Tue, May 7, 2024 at 8:28 AM Mark Adams wrote: > "A^T A being similar to a finite difference Poisson matrix if the rows > were permuted randomly." > Normal eqs are a good option in general for rectangular systems and we > have Poisson solvers. > I missed that. Why not first permute back to the

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Mark Adams
"A^T A being similar to a finite difference Poisson matrix if the rows were permuted randomly." Normal eqs are a good option in general for rectangular systems and we have Poisson solvers. I'm not sure what you mean by "permuted randomly." A random permutation of the matrix can kill performance bu

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Matthew Knepley
On Tue, May 7, 2024 at 5:12 AM Pierre Jolivet wrote: > On 7 May 2024, at 9: 10 AM, Marco Seiz wrote: Thanks > for the quick response! On 07. 05. 24 14: 24, Pierre Jolivet wrote: On 7 > May 2024, at 7: 04 AM, Marco Seiz wrote: This > Message Is From an External > ZjQcmQRYFpfptBannerStart > This

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Pierre Jolivet
> On 7 May 2024, at 9:10 AM, Marco Seiz wrote: > > Thanks for the quick response! > > On 07.05.24 14:24, Pierre Jolivet wrote: >> >> >>> On 7 May 2024, at 7:04 AM, Marco Seiz wrote: >>> >>> This Message Is From an External Sender >>> This message came from outside your organization. >>> He

Re: [petsc-users] Reasons for breakdown in preconditioned LSQR

2024-05-07 Thread Marco Seiz
Thanks for the quick response! On 07. 05. 24 14: 24, Pierre Jolivet wrote: > > >> On 7 May 2024, at 7: 04 AM, Marco Seiz wrote: >> >> This Message Is From an External Sender >> This message ZjQcmQRYFpfptBannerStart This Message Is From an External Sender