<!-- BaNnErBlUrFlE-BoDy-start -->
<!-- Preheader Text : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">
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
</div>
<!-- Preheader Text : END -->
<!-- Email Banner : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerStart</div>
<!--[if ((ie)|(mso))]>
<table border="0" cellspacing="0" cellpadding="0" width="100%" style="padding: 16px 0px 16px 0px; direction: ltr" ><tr><td>
<table border="0" cellspacing="0" cellpadding="0" style="padding: 0px 10px 5px 6px; width: 100%; border-radius:4px; border-top:4px solid #90a4ae;background-color:#D0D8DC;"><tr><td valign="top">
<table align="left" border="0" cellspacing="0" cellpadding="0" style="padding: 4px 8px 4px 8px">
<tr><td style="color:#000000; font-family: 'Arial', sans-serif; font-weight:bold; font-size:14px; direction: ltr">
This Message Is From an External Sender
</td></tr>
<tr><td style="color:#000000; font-weight:normal; font-family: 'Arial', sans-serif; font-size:12px; direction: ltr">
This message came from outside your organization.
</td></tr>
</table>
</td></tr></table>
</td></tr></table>
<![endif]-->
<![if !((ie)|(mso))]>
<div dir="ltr" id="pfptBannery32hwk8" style="all: revert !important; display:block !important; text-align: left !important; margin:16px 0px 16px 0px !important; padding:8px 16px 8px 16px !important; border-radius: 4px !important; min-width: 200px !important; background-color: #D0D8DC !important; background-color: #D0D8DC; border-top: 4px solid #90a4ae !important; border-top: 4px solid #90a4ae;">
<div id="pfptBannery32hwk8" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBannery32hwk8" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-weight:bold !important; font-weight:bold; font-size:14px !important; line-height:18px !important; line-height:18px">
This Message Is From an External Sender
</div>
<div id="pfptBannery32hwk8" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-weight:normal; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-size:12px !important; line-height:18px !important; line-height:18px; margin-top:2px !important;">
This message came from outside your organization.
</div>
</div>
<div style="clear: both !important; display: block !important; visibility: hidden !important; line-height: 0 !important; font-size: 0.01px !important; height: 0px"> </div>
</div>
<![endif]>
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerEnd</div>
<!-- Email Banner : END -->
<!-- BaNnErBlUrFlE-BoDy-end -->
<html>
<head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBannery32hwk8 { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButtony32hwk8:hover, .pfptPrimaryButtony32hwk8:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButtony32hwk8:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<p class="elementToProof" style="line-height: normal; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36); background-color: white;">Hello everyone,</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36);"> </span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36);">I hope this email finds you well.</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36);"><br>
<br>
</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36);"> My Name is Sawsan Shatanawi,</span><span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);"> 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 results although I checked my matrix and RHS and everything, I decided to change my solver
to Newton Rapson method.<br>
I checked PETSc documents but I have a few questions:<br>
1) My groundwater system is time-dependent, so should I use TS only instead of SNES?<br>
2) My system has its deltaT, would using deltaT as dt affect my solver, or is it better to use TS-PETSc dt? Also, would using PETSc dt affect the simulation of the groundwater system</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);">3) I want my Jacobian matrix to be calculated by PETSc automatically
</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);">4) Do I need to define and calculate the residual vector?</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);"> </span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);">My A-Matrix contains coefficients and external sources and my RHS vector includes the boundary conditions <br>
</span><span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36);"><br>
<br>
</span></p>
<p style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(36, 36, 36);">Please find the attached file contains a draft of my code</span><span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);"><br>
<br>
</span></p>
<p class="elementToProof" style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);">Thank you in advance for your time and help.</span></p>
<p class="elementToProof" style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);"><br>
</span></p>
<div style="background-color: white; padding: 0in; border-width: 1pt; border-style: solid; border-color: rgb(217, 217, 227);">
<p style="line-height: normal; background-color: white; margin: 15pt 0in; padding: 0in; font-family: Calibri, sans-serif; font-size: 11pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);">Best regards,</span></p>
<p style="line-height: normal; background-color: white; margin: 15pt 0in; padding: 0in; font-family: Calibri, sans-serif; font-size: 11pt;">
<span style="font-family: "Segoe UI", sans-serif; font-size: 11.5pt; color: rgb(55, 65, 81);"> Sawsan</span></p>
</div>
<div class="elementToProof" style="line-height: normal; background-color: white; margin: 0in 0in 0.0001pt; font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Shatanawi, Sawsan Muhammad <sawsan.shatanawi@wsu.edu><br>
<b>Sent:</b> Tuesday, January 16, 2024 10:43 AM<br>
<b>To:</b> Junchao Zhang <junchao.zhang@gmail.com><br>
<b>Cc:</b> Barry Smith <bsmith@petsc.dev>; Matthew Knepley <knepley@gmail.com>; Mark Adams <mfadams@lbl.gov>; petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div class="x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Hello all,<br>
<br>
Thank you for your valuable help. I will do your recommendations and hope it will run without any issues.</div>
<div class="x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Bests,<br>
Sawsan</div>
<div class="x_elementToProof" style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
</div>
<div id="x_appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Junchao Zhang <junchao.zhang@gmail.com><br>
<b>Sent:</b> Friday, January 12, 2024 8:46 AM<br>
<b>To:</b> Shatanawi, Sawsan Muhammad <sawsan.shatanawi@wsu.edu><br>
<b>Cc:</b> Barry Smith <bsmith@petsc.dev>; Matthew Knepley <knepley@gmail.com>; Mark Adams <mfadams@lbl.gov>; petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<div>
<p style="border:1px; border-style:solid; border-color:#FFCACA; background-color:#FFEB9C; padding:1px">
<span style="font-size:10.0pt; color:black"><b>[EXTERNAL EMAIL]</span></b></p>
<div>
<div dir="ltr">Hi, Sawsan,
<div> First in test_main.F90, you need to call VecGetArrayF90(temp_solution, H_vector, ierr) and VecRestoreArrayF90 (temp_solution, H_vector, ierr) as Barry mentioned.<br>
<div> Secondly, in the loop of test_main.F90, it calls GW_solver(). Within it, it calls PetscInitialize()/PetscFinalize(). But without MPI being initialized, PetscInitialize()/PetscFinalize()<b> can only be called once.
</b></div>
<div><font face="monospace"> do timestep =2 , NTSP<br>
call GW_boundary_conditions(timestep-1)<br>
!print *,HNEW(1,1,1)<br>
call GW_elevation()<br>
! print *, GWTOP(2,2,2)<br>
call GW_conductance()<br>
! print *, CC(2,2,2)<br>
call GW_recharge()<br>
! print *, B_Rech(5,4)<br>
call GW_pumping(timestep-1)<br>
! print *, B_pump(2,2,2)<br>
call GW_SW(timestep-1)<br>
print *,B_RIVER (2,2,2)<br>
call GW_solver(timestep-1,N)<br>
call GW_deallocate_loop()<br>
end do</font><br>
</div>
<div><br>
</div>
A solution is to delete PetscInitialize()/PetscFinalize() in GW_solver_try.F90 and move it to test_main.F90, outside the do loop.
<div><br>
</div>
<div><font face="monospace">diff --git a/test_main.F90 b/test_main.F90<br>
index b5997c55..107bd3ee 100644<br>
--- a/test_main.F90<br>
+++ b/test_main.F90<br>
@@ -1,5 +1,6 @@<br>
program test_GW<br>
<br>
+#include <petsc/finclude/petsc.h><br>
use petsc<br>
use GW_constants<br>
use GW_param_by_user<br>
@@ -8,6 +9,9 @@ program test_GW<br>
implicit none<br>
integer :: N<br>
integer :: timestep<br>
+ PetscErrorCode ierr<br>
+<br>
+ call PetscInitialize(ierr)<br>
call GW_domain(N)<br>
!print *, "N=",N<br>
!print *, DELTAT<br>
@@ -37,4 +41,5 @@ program test_GW<br>
end do<br>
print *, HNEW(NCOL,3,2)<br>
call GW_deallocate ()<br>
+ call PetscFinalize(ierr)<br>
end program test_GW</font><br>
</div>
<div><br>
</div>
<div>With that, the MPI error will be fixed. The code could run to gw_deallocate () before abort. There are other memory errors. You can install/use valgrind to fix them. Run it with valgrind ./GW.exe and look through the output</div>
<div><br>
</div>
<div><br>
</div>
<div>Thanks.</div>
<div>--Junchao Zhang<br>
</div>
<div><br>
</div>
</div>
</div>
<br>
<div class="x_x_gmail_quote">
<div dir="ltr" class="x_x_gmail_attr">On Thu, Jan 11, 2024 at 10:49 PM Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu">sawsan.shatanawi@wsu.edu</a>> wrote:<br>
</div>
<blockquote class="x_x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div class="x_x_msg6313238004422759188">
<div dir="ltr">
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Hello,<br>
<br>
Thank you all for your help.<br>
<br>
I have changed VecGetArray to VecGetArrayF90, and the location of destory call. but I want to make sure that VecGet ArrayF90 is to make a new array( vector) that I can use in the rest of my Fortran code?<br>
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
when I run it and debugged it, I got </div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"> 5.2000000E-03</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"> 50.00000</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"> 10.00000</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"> 0.0000000E+00</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">PETSC: Attaching gdb to /weka/data/lab/richey/sawsan/GW_CODE/code2024/SS_GWM/./GW.exe of pid 33065 on display :0.0 on
machine sn16</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Unable to start debugger in xterm: No such file or directory</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"> 0.0000000E+00</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,0)">Attempting to use an MPI routine after finalizing MPICH</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">srun: error: sn16: task 0: Exited with exit code 1</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">[sawsan.shatanawi@login-p2n02 SS_GWM]$ gdb ./GW/exe</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">GNU gdb (GDB) Red Hat Enterprise Linux 7.6.1-100.el7</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Copyright (C) 2013 Free Software Foundation, Inc.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">License GPLv3+: GNU GPL version 3 or later <<a href="https://urldefense.com/v3/__http://gnu.org/licenses/gpl.html__;!!JmPEgBY0HMszNaDT!o4qS1zaFLg2L8PlawWJVyYsJpnwYHuL6SIIZsmzTO98RCzbP26HwTi_0-ipCS_D2SBG0X4gtEnM13-nbzFEFKvtraKaM$" originalsrc="https://urldefense.com/v3/__http://gnu.org/licenses/gpl.html__;!!JmPEgBY0HMszNaDT!o4qS1zaFLg2L8PlawWJVyYsJpnwYHuL6SIIZsmzTO98RCzbP26HwTi_0-ipCS_D2SBG0X4gtEnM13-nbzFEFKvtraKaM$" shash="Yw2mN7hWQC3K1XPN8YH8YZb/IB3/69wogiYM/Q2Oeph0cFQlR51C3pnzjEPHROff2FTZCjqqtwQpGY12XJr5vO3hEUCfZWZBRUYhXnzSLAwMl4gF09n+Fk3SUDSllJFN7TfNpf9fQSDm5OSeo9HMqCPZ2oLbnuIrw2RyK86IEO8=" target="_blank">http://gnu.org/licenses/gpl.html</a>></span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">This is free software: you are free to change and redistribute it.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">There is NO WARRANTY, to the extent permitted by law. Type "show copying"</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">and "show warranty" for details.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">This GDB was configured as "x86_64-redhat-linux-gnu".</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">For bug reporting instructions, please see:</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)"><<a href="https://urldefense.com/v3/__http://www.gnu.org/software/gdb/bugs/__;!!JmPEgBY0HMszNaDT!o4qS1zaFLg2L8PlawWJVyYsJpnwYHuL6SIIZsmzTO98RCzbP26HwTi_0-ipCS_D2SBG0X4gtEnM13-nbzFEFKjO8jNfs$" originalsrc="https://urldefense.com/v3/__http://www.gnu.org/software/gdb/bugs/__;!!JmPEgBY0HMszNaDT!o4qS1zaFLg2L8PlawWJVyYsJpnwYHuL6SIIZsmzTO98RCzbP26HwTi_0-ipCS_D2SBG0X4gtEnM13-nbzFEFKjO8jNfs$" shash="eB5DijySrhtlvnFYSQBciMY5NHYo7ew/4cvNNnc5CbzloZpMQ1xAt95GpbAk3w3i40tTCiTbSqPR9QanBpZjU8AaBluFOUXK2jfr/kLK59hx4ggWNpYEzs88OyNndPkzh4yObyoED/m3Pi/fQJQ6XDv9HUZZu0MHlhFTnB6yMOI=" target="_blank">http://www.gnu.org/software/gdb/bugs/</a>>...</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">./GW/exe: No such file or directory.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">(gdb) run</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Starting program:</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">No executable file specified.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Use the "file" or "exec-file" command.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">(gdb) bt</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">No stack.</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">(gdb)<br>
<br>
If the highlighted line is the error, I don't know why when I write gdb , it does not show me the location of error</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">The code :
<a href="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!o4qS1zaFLg2L8PlawWJVyYsJpnwYHuL6SIIZsmzTO98RCzbP26HwTi_0-ipCS_D2SBG0X4gtEnM13-nbzFEFKp8i33ur$" originalsrc="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!o4qS1zaFLg2L8PlawWJVyYsJpnwYHuL6SIIZsmzTO98RCzbP26HwTi_0-ipCS_D2SBG0X4gtEnM13-nbzFEFKp8i33ur$" shash="rfeyfrK1lPD4FCFzH+Njnu35ia+RENa/KPfyS0ZTMara8F344H2ia3D6XRdf6aJw9RPj8VQtMqKKLplfCUGMn58GGtOGyy/esthW6uvBPhZWeX38gRatJOtUNtfZEBZJoRvG9U9TOdXBM/T4JgrIEZX94b5BZ55STMC0SxpUSRE=" id="x_x_m_6313238004422759188OWAf2a88c8d-1aa9-f2d9-963c-456c02baba7d" target="_blank">
sshatanawi/SS_GWM (github.com)</a> <br>
<br>
I really appreciate your helps<br>
<br>
</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Sawsan</span></div>
<div id="x_x_m_6313238004422759188appendonsend"></div>
<hr style="display:inline-block; width:98%">
<div id="x_x_m_6313238004422759188divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>><br>
<b>Sent:</b> Wednesday, January 10, 2024 5:35 PM<br>
<b>To:</b> Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>><br>
<b>Cc:</b> Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>>; Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>>;
<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<div>
<p style="border:1px solid rgb(255,202,202); background-color:rgb(255,235,156); padding:1px">
<span style="font-size:10pt; color:black"><b>[EXTERNAL EMAIL]</b></span></p>
<div><br id="x_x_m_6313238004422759188x_lineBreakAtBeginningOfMessage">
<div><br>
<blockquote type="cite">
<div>On Jan 10, 2024, at 6:49 PM, Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>> wrote:</div>
<br>
<div>
<div dir="ltr">Hi, Sawsan,
<div> I could build your code and I also could gdb it.</div>
<div><br>
</div>
<div>$ gdb ./GW.exe</div>
<div>...</div>
<div>$ Thread 1 "GW.exe" received signal SIGSEGV, Segmentation fault.</div>
0x00007ffff1e6d44f in vecgetarray_ (x=0x7fffffffa718, fa=0x0, ia=0x7fffffffa75c, ierr=0x0) at /scratch/jczhang/petsc/src/vec/vec/interface/ftn-custom/zvectorf.c:257<br>
257 *ierr = VecGetArray(*x, &lx);<br>
(gdb) bt<br>
#0 0x00007ffff1e6d44f in vecgetarray_ (x=0x7fffffffa718, fa=0x0, ia=0x7fffffffa75c, ierr=0x0) at /scratch/jczhang/petsc/src/vec/vec/interface/ftn-custom/zvectorf.c:257<br>
#1 0x000000000040b6e3 in gw_solver (t_s=1.40129846e-45, n=300) at GW_solver_try.F90:169<br>
#2 0x000000000040c6a8 in test_gw () at test_main.F90:35
<div> </div>
<div><font color="#ff0000">ierr=0x0 </font><font>caused the segfault. See </font><a href="https://urldefense.com/v3/__https://petsc.org/release/manualpages/Vec/VecGetArray/*vecgetarray__;Iw!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ9Km12jA$" originalsrc="https://urldefense.com/v3/__https://petsc.org/release/manualpages/Vec/VecGetArray/*vecgetarray__;Iw!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ9Km12jA$" shash="eTMT/2IvtfZNSLkFumxiDqpRSfF/dgzmXotN+byvDZRxq/oemvkHqAoRDJ6VBQDDgt8NDoy4uaJUa/2nLvLg63RpYKZTYAE/FXqA5XMXVcQghCk2cdxNU3znqka7CkTfwKg7j36WcSIvZae3GsncanykV1g7dy693xt9SULqdzU=" target="_blank">https://petsc.org/release/manualpages/Vec/VecGetArray/#vecgetarray</a>,
you should use VecGetArrayF90 instead.<br>
</div>
<div><br>
</div>
<div>BTW, Barry, the code <a href="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM/blob/main/GW_solver_try.F90*L169__;Iw!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZh2eAi4o$" originalsrc="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM/blob/main/GW_solver_try.F90*L169__;Iw!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZh2eAi4o$" shash="ke6goHMtptjx+CF9yL/R/sT9aCkk5ug5hGph6iJ3ifOSAGLKfOEA6sk77fVkMpX/xyZMIUB5Gy9b03uElW8dhgqfO5/VzeFcnvak8tbHE0vTSkjXrE32O5UoQRLeAuU/Yv/w+5MShST2VC6pH0V0rrJ0g8hbBidUrHGFbsOmCrU=" target="_blank">
https://github.com/sshatanawi/SS_GWM/blob/main/GW_solver_try.F90#L169</a> has "call VecGetArray(temp_solution, H_vector, ierr)". I don't find petsc Fortran examples doing VecGetArray. Do we still support it?</div>
</div>
</div>
</blockquote>
<div><br>
</div>
This is not the correct calling sequence for VecGetArray() from Fortran. </div>
<div><br>
</div>
<div> Regardless, definitely should not be writing any new code that uses VecGetArray() from Fortran. Should use VecGetArrayF90().</div>
<div><br>
<blockquote type="cite">
<div>
<div dir="ltr">
<div><br>
</div>
<div>
<div>
<div dir="ltr">
<div dir="ltr">--Junchao Zhang</div>
</div>
</div>
<br>
</div>
</div>
<br>
<div>
<div dir="ltr">On Wed, Jan 10, 2024 at 2:38 PM Shatanawi, Sawsan Muhammad via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">Hello all,</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"><br>
</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">I hope you are doing well.</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"><br>
</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">Generally, I use gdb <the name of my exe.file> to debug the
code.</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"> I got the attached error message.</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"><br>
</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">I have tried to add the flag -start_in_debugger in the make
file, but it didn't work, so it seems I was doing it in the wrong way</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"><br>
</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">This is the link for the whole code: <a href="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ1Veab3M$" originalsrc="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ1Veab3M$" shash="FLwlmhpwFDsmWuDZID6fKEI46YwJ6CEVHJBvcIEUDIW3AOTXlxCYofuA8lQTJ/tebHAtGpNGIKiUMXt+GL3gvn0RL+lKsrkJYWP4XO9HSc5oCGbRUsyUaWGo69Ae9SpaoDYvzBFPkxPatWaBsvBBPI8tAHTouaaxR68Mm8Dqfsg=" id="x_x_m_6313238004422759188x_m_5694742774502394814OWA37645568-445c-55da-2f91-abc541b1f4d6" target="_blank" style="margin:0px">sshatanawi/SS_GWM
(github.com)</a></span></div>
<div style="font-style:normal; font-weight:400; font-size:medium; margin:0px; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; word-spacing:0px">
<div id="x_x_m_6313238004422759188x_m_5694742774502394814LPBorder_GTaHR0cHM6Ly9naXRodWIuY29tL3NzaGF0YW5hd2kvU1NfR1dN" style="margin:16px 0px; max-width:800px; min-width:424px">
<table id="x_x_m_6313238004422759188x_m_5694742774502394814LPContainer867037" style="padding:12px 36px 12px 12px; width:600px; border-width:1px; border-style:solid; border-color:rgb(200,200,200); border-radius:2px">
<tbody>
<tr valign="top" style="border-spacing:0px">
<td>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814LPImageContainer867037" style="margin:0px 12px 0px 0px; height:120px; overflow:hidden">
<a id="x_x_m_6313238004422759188x_m_5694742774502394814LPImageAnchor867037" href="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ1Veab3M$" originalsrc="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ1Veab3M$" shash="isDxSGgS3nK2npqU3rzu1ZTz5Q/oQ2lSy6glpeM8s50e1l5wVug1hDei/E/gAENeGg3ZD6oRzEBrS1uuEwpL+IcrTwxAFhlu1rl6HVxwSgrIX6XJwCeEeGGLSfBeZt/nKSIAR13CM4hNhthq9lSmjaDHb6fQLKcg3pDMu+szwlM=" target="_blank" style="margin:0px"><img id="x_x_m_6313238004422759188x_m_5694742774502394814LPThumbnailImageId867037" alt="" height="120" width="240" style="margin:0px; display:block" src="https://opengraph.githubassets.com/9eb6cd14baf12f04848ed209b6f502415eb531bdd7b3a5f9696af68663b870c0/sshatanawi/SS_GWM"></a></div>
</td>
<td style="width:289.6px">
<div id="x_x_m_6313238004422759188x_m_5694742774502394814LPTitle867037" style="font-weight:300; font-size:21px; font-family:wf_segoe-ui_light,"Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; margin:0px 8px 12px 0px">
<a id="x_x_m_6313238004422759188x_m_5694742774502394814LPUrlAnchor867037" href="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ1Veab3M$" originalsrc="https://urldefense.com/v3/__https://github.com/sshatanawi/SS_GWM__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ1Veab3M$" shash="R7iq0vLmqA5uzT4Yhkx/uS6vZ5JeBm9WqUdYDe25+4uB228z+MCFgXsSQ2j195jfqKTS5V/9rLKo+/ySIhA6VCOqkxNyU58WZgg2LWLIsnZKfqQzL1+mLfc9Bfrs83c6EG/nHas1NxcUOM87IunNjf3O/z8RZs0nS998tLOrRH4=" target="_blank" style="margin:0px; text-decoration:none">GitHub
- sshatanawi/SS_GWM</a></div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814LPDescription867037" style="font-size:14px; font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; margin:0px 8px 12px 0px; max-height:100px; overflow:hidden; color:rgb(102,102,102)">
Contribute to sshatanawi/SS_GWM development by creating an account on GitHub.</div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814LPMetadata867037" style="font-size:14px; font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; margin:0px; color:rgb(166,166,166)">
<a href="https://urldefense.com/v3/__http://github.com/__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ8rcrPiA$" originalsrc="https://urldefense.com/v3/__http://github.com/__;!!JmPEgBY0HMszNaDT!tqBApprMfYxwNz4Zvnk8coNE5AeWjA9wSdAM7QJcIIVP1z0VDsVIalo4Sew2b0fW3bZtTAbPh-h0MUsZ8rcrPiA$" shash="mWGnWyvZz+gmGJRm3ZSa1fIl1Oe1uv/FraARdnPINna3UyQH+8OfZPzylj8Xkpr9/vmfmclaZI88i1EyEK5JpeUUI0zun2K1pn2hTZCTavsyGdnAY10ZH9vJ9HRFwHIBqgulfIWog50KH7ng1HFoHMDET6Hwblr8wg46NkF9ORw=" target="_blank">github.com</a></div>
</td>
</tr>
</tbody>
</table>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814LPCloseButtonContainer867037" title="Remove link preview" role="button">
<i aria-hidden="true" id="x_x_m_6313238004422759188x_m_5694742774502394814LPCloseButton867037"></i></div>
</div>
</div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"><br>
</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">You can read the description of the code in " </span><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400; color:rgb(12,100,192)">Model
Desprciption.pdf</span><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">"</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">the compiling file is </span><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400; color:rgb(12,100,192)">makefile_f90</span><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"> where
you can find the linked code files</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400"><br>
</span></div>
<div style="text-align:left; margin:0px"><span style="letter-spacing:normal; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt; font-weight:400">I really appreciate your help<br>
<br>
Bests,<br>
Sawsan</span></div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814appendonsend"></div>
<hr style="display:inline-block; width:98%">
<div id="x_x_m_6313238004422759188x_m_5694742774502394814divRplyFwdMsg" dir="ltr">
<font face="Calibri, sans-serif" style="font-size:11pt"><b>From:</b> Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>><br>
<b>Sent:</b> Friday, January 5, 2024 4:53 AM<br>
<b>To:</b> Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>><br>
<b>Cc:</b> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>>;
<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<div>
<p style="border:1px solid rgb(255,202,202); background-color:rgb(255,235,156); padding:1px">
<span style="font-size:10pt"><b>[EXTERNAL EMAIL]</b></span></p>
<div>
<div dir="ltr">This is a segv. As Matt said, you need to use a debugger for this or add print statements to narrow down the place where this happens.
<div><br>
</div>
<div>You will need to learn how to use debuggers to do your project so you might as well start now.</div>
<div><br>
</div>
<div>If you have a machine with a GUI debugger that is easier but command line debuggers are good to learn anyway.</div>
<div><br>
</div>
<div>I tend to run debuggers directly (eg, lldb ./a.out -- program-args ...) and use a GUI debugger (eg, Totalview or DDT) if available.</div>
<div><br>
</div>
<div>Mark</div>
<div><br>
</div>
</div>
<br>
<div>
<div dir="ltr">On Wed, Dec 20, 2023 at 10:02 PM Shatanawi, Sawsan Muhammad via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
Hello Matthew,<br>
<br>
Thank you for your help. I am sorry that I keep coming back with my error messages, but I reached a point that I don't know how to fix them, and I don't understand them easily.<br>
The list of errors is getting shorter, now I am getting the attached error messages </div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
Thank you again,<br>
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
Sawsan</div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399appendonsend">
</div>
<hr style="display:inline-block; width:98%">
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399divRplyFwdMsg" dir="ltr">
<font face="Calibri, sans-serif" style="font-size:11pt"><b>From:</b> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
<b>Sent:</b> Wednesday, December 20, 2023 6:54 PM<br>
<b>To:</b> Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>><br>
<b>Cc:</b> Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>>;
<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<div>
<p style="border:1px solid rgb(255,202,202); background-color:rgb(255,235,156); padding:1px">
<span style="font-size:10pt"><b>[EXTERNAL EMAIL]</b></span></p>
<div>
<div dir="ltr">
<div dir="ltr">On Wed, Dec 20, 2023 at 9:49 PM Shatanawi, Sawsan Muhammad via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
</div>
<div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif">Hello Barry,<br>
<br>
Thank you a lot for your help, Now I am getting the attached error message.</span></div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Do not destroy the PC from KSPGetPC()</div>
<div><br>
</div>
<div> THanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif">Bests,<br>
Sawsan</span></div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543appendonsend">
</div>
<hr style="display:inline-block; width:98%">
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543divRplyFwdMsg" dir="ltr">
<font face="Calibri, sans-serif" style="font-size:11pt"><b>From:</b> Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>><br>
<b>Sent:</b> Wednesday, December 20, 2023 6:32 PM<br>
<b>To:</b> Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>><br>
<b>Cc:</b> Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>>;
<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<div>
<p style="border:1px solid rgb(255,202,202); background-color:rgb(255,235,156); padding:1px">
<span style="font-size:10pt"><b>[EXTERNAL EMAIL]</b></span></p>
<div>
<div><br>
</div>
Instead of
<div><br>
</div>
<div>
<div> call PCCreate(PETSC_COMM_WORLD, pc, ierr)</div>
<div> call PCSetType(pc, PCILU,ierr) ! Choose a preconditioner type (ILU)</div>
<div> call KSPSetPC(ksp, pc,ierr) ! Associate the preconditioner with the KSP solver</div>
<div><br>
</div>
<div>do</div>
<div><br>
</div>
<div> call KSPGetPC(ksp,pc,ierr)</div>
<div> call PCSetType(pc, PCILU,ierr)</div>
<div><br>
</div>
<div>Do not call KSPSetUp(). It will be taken care of automatically during the solve</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div><br>
<blockquote type="cite">
<div>On Dec 20, 2023, at 8:52 PM, Shatanawi, Sawsan Muhammad via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:</div>
<br>
<div>
<div style="font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
Hello, <br>
I don't think that I set preallocation values when I created the matrix, would you please have look at my code. It is just the petsc related part from my code.<br>
I was able to fix some of the error messages. Now I have a new set of error messages related to the KSP solver (attached)<br>
<br>
I appreciate your help</div>
<div style="font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
<br>
</div>
<div style="font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
Sawsan</div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_appendonsend" style="font-family:Helvetica; font-size:18px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
</div>
<hr style="font-family:Helvetica; font-size:18px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; display:inline-block; width:934.906px">
<span style="font-family:Helvetica; font-size:18px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; float:none; display:inline"></span>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_divRplyFwdMsg" dir="ltr" style="font-family:Helvetica; font-size:18px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<font face="Calibri, sans-serif" style="font-size:11pt"><b>From:</b><span> </span>Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>><br>
<b>Sent:</b><span> </span>Wednesday, December 20, 2023 6:44 AM<br>
<b>To:</b><span> </span>Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>><br>
<b>Cc:</b><span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b><span> </span>Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code</font>
<div> </div>
</div>
<div style="font-family:Helvetica; font-size:18px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div style="margin-top:0px; margin-bottom:0px; border:1px solid rgb(255,202,202); background-color:rgb(255,235,156); padding:1px">
<span style="font-size:10pt"><b>[EXTERNAL EMAIL]</b></span></div>
<div>
<div dir="ltr">Did you set preallocation values when you created the matrix?
<div>Don't do that.</div>
</div>
<br>
<div>
<div dir="ltr">On Wed, Dec 20, 2023 at 9:36 AM Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div dir="ltr">Hello, </div>
<div dir="ltr"><br>
</div>
<div dir="ltr">I am trying to create a sparse matrix( which is as I believe a zero matrix) then adding some nonzero elements to it over a loop, then assembling it </div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_x_m_-677654408871201876ms-outlook-mobile-signature">
<div><br>
</div>
Get<span> </span><a href="https://urldefense.com/v3/__https://aka.ms/o0ukef__;!!JmPEgBY0HMszNaDT!uUJ_jeYf45gcXDGR_PeMjhU7hbd_fKcXJPn0pM9eb-YQihKNYuXMYM9x-hglsbXsCFIwNBWgHXdetHODupsOloE$" originalsrc="https://urldefense.com/v3/__https://aka.ms/o0ukef__;!!JmPEgBY0HMszNaDT!uUJ_jeYf45gcXDGR_PeMjhU7hbd_fKcXJPn0pM9eb-YQihKNYuXMYM9x-hglsbXsCFIwNBWgHXdetHODupsOloE$" shash="zqhJ/KQrH0HxFZ5lBgQjyvfWQbG9RQK4AaDjsXWVp9mxi3ZyFa2CzjZDijjfC66UU6a1ZipCkxazY/kPXaJzYRaCNUZEKSnAZkQtelKepLSwHHQQProEnWRFu4uwYvMWNHOzjF0jR2qv6qWSUTCDTB2rLpoU5ql6F2nnc1ohbw8=" target="_blank">Outlook
for iOS</a></div>
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_x_m_-677654408871201876mail-editor-reference-message-container">
<hr style="display:inline-block; width:917.344px">
<div id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_x_m_-677654408871201876divRplyFwdMsg" dir="ltr">
<font face="Calibri, sans-serif"><b>From:</b><span> </span>Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>><br>
<b>Sent:</b><span> </span>Wednesday, December 20, 2023 2:48 AM<br>
<b>To:</b><span> </span>Shatanawi, Sawsan Muhammad <<a href="mailto:sawsan.shatanawi@wsu.edu" target="_blank">sawsan.shatanawi@wsu.edu</a>><br>
<b>Cc:</b><span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b><span> </span>Re: [petsc-users] Help with Integrating PETSc into Fortran Groundwater Flow Simulation Code
<div> </div>
</font></div>
<div style="margin-top:0px; margin-bottom:0px; border:1px solid rgb(255,202,202); padding:1px; background-color:rgb(255,235,156)">
<span style="font-size:10pt"><b>[EXTERNAL EMAIL]</b></span></div>
<div>
<div dir="ltr">I am guessing that you are creating a matrix, adding to it, finalizing it ("assembly"), and then adding to it again, which is fine, but you are adding new non-zeros to the sparsity pattern.
<div>If this is what you want then you can tell the matrix to let you do that.</div>
<div>Otherwise you have a bug.</div>
<div><br>
</div>
<div>Mark</div>
</div>
<br>
<div>
<div dir="ltr">On Tue, Dec 19, 2023 at 9:50 PM Shatanawi, Sawsan Muhammad via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
</div>
Hello everyone,
<div><br>
</div>
<div>I hope this email finds you well.</div>
<div><span><br>
</span></div>
<div><span> My Name is Sawsan Shatanawi,</span><span><span style="display:inline; color:rgb(55,65,81)"> and I am currently working on developing a Fortran code for simulating groundwater flow in a 3D system. The code involves solving a nonlinear system, and
I have created the matrix to be solved using the PCG solver and Picard iteration. However, </span></span><span>when I tried to assign it as a PETSc matrix I started getting a lot of error messages.</span></div>
<div><br>
</div>
<div>I am kindly asking if someone can help me,<span style="display:inline; color:rgb(55,65,81)"> I would be happy to share my code with him/her.</span></div>
<div><span><br>
</span></div>
<div><span>Please find the attached file contains a list of errors I have gotten</span></div>
<div><span style="color:rgb(55,65,81)"><br>
</span></div>
<div><span style="color:rgb(55,65,81)">Thank you in advance for your time and assistance.</span><br>
</div>
<div><span>
<p style="margin:1.25em 0px; border:0px solid rgb(217,217,227); box-sizing:border-box; color:rgb(55,65,81)">
Best regards,</p>
<p style="margin:1.25em 0px; border:0px solid rgb(217,217,227); box-sizing:border-box; color:rgb(55,65,81)">
Sawsan</p>
<br>
</span></div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif; font-size:12pt">
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
<span id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_cid:66357E60-6F0B-4E28-9883-5421A4F75577"><Matrix_RHS.F90></span><span id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_cid:FB95EDCC-3170-4E85-BBFB-5804756A9022"><out.txt></span><span id="x_x_m_6313238004422759188x_m_5694742774502394814x_m_1291163558337743399x_m_6596943672973964543x_cid:8AF6D377-8CA5-4FFF-81AC-12EFB1837766"><solver.F90></span></div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
<span>-- </span><br>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener</div>
<div><br>
</div>
<div><a href="https://urldefense.com/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!JmPEgBY0HMszNaDT!uskvAyF0pMMWDbMIexr9g4qN46V7Rea17GQdNIVG2vH_HMaX7mXgie4ZYgusmPpss_DS7H1_8vn8arGQNSkC$" originalsrc="https://urldefense.com/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!JmPEgBY0HMszNaDT!uskvAyF0pMMWDbMIexr9g4qN46V7Rea17GQdNIVG2vH_HMaX7mXgie4ZYgusmPpss_DS7H1_8vn8arGQNSkC$" shash="VGO9igiXq9J+bviDFWJOUVLsYQFFKhC1Alfb9jdMmC2HqnMvrePDWYz7H3sFjdiNAcfFwDBloP4Uutc3w2cyjcRjyjpza6wWjTsAS4F59aE3au/jXXdje8BFTjjsQgL6F/8dO/rISfyxkeVuhZ2uWiHRpDttXyOKulUc2punTlU=" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</body>
</html>