<!DOCTYPE html>
<!-- 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;">
Another issue arising from the tighter type checking. I have fixed the issue in the branch, again I rebased the branch so you cannot just git pull it (sorry). Also you will now need to change PETSC_NULL_OBJECT with PETSC_NULL_MAT for the second
</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="pfptBannerscrfq3j" 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="pfptBannerscrfq3j" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBannerscrfq3j" 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="pfptBannerscrfq3j" 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>
#pfptBannerscrfq3j { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButtonscrfq3j:hover, .pfptPrimaryButtonscrfq3j:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButtonscrfq3j:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
<meta charset="utf-8"></head><body><pre style="font-family: sans-serif; font-size: 100%; white-space: pre-wrap; word-wrap: break-word">
Another issue arising from the tighter type checking. I have fixed the issue in the branch, again I rebased the branch so you cannot just git pull it (sorry).
Also you will now need to change PETSC_NULL_OBJECT with PETSC_NULL_MAT for the second matrix argument for SNESGetJacobian().
Barry
> On Jun 26, 2024, at 12:49 AM, Adrian Croucher <a.croucher@auckland.ac.nz> wrote:
>
>
> On 24/06/24 1:39 pm, Adrian Croucher wrote:
>>
>> If this kind of usage is no longer going to be supported I may be able to change it to something more like what ex1f.F90 does. My SNES_Jacobian() is basically just a wrapper for SNESComputeJacobianDefaultColor() now - I think it may have had some other things in it once but it doesn't now. However, this issue might cause trouble for other people as well.
>
> I've just been experimenting with altering my code so that the SNESSetJacobian() works the same as in snes/tests/ex1f.F90, i.e. using something like:
>
> call SNESSetJacobian(solver, J, J, SNESComputeJacobianDefaultColor, fd_coloring, ierr)
>
> instead of using a custom Jacobian routine that calls SNESComputeJacobianDefaultColor(). That seems to build ok using your PETSc branch and I don't get a runtime error calling SNESComputeJacobianDefaultColor().
>
> However I get one when another part of the code subsequently calls SNESGetJacobian(). I call it like this:
>
> call SNESGetJacobian(solver, J, PETSC_NULL_MAT, PETSC_NULL_FUNCTION, &
> PETSC_NULL_INTEGER, ierr)
>
> and the error is below. This has worked up until now- is there something wrong with it?
>
> I tried putting a similar line into snes/tests/ex1f.F90 and it gets the same error.
>
> - Adrian
>
> [0]PETSC ERROR: Corrupt argument: <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/*valgrind__;Iw!!G_uCfscf7eWS!b1Gx5q-NMarPE0nExBMOr0s4mG_Z0OJaWf6V7NyRuywDwA4IJ75v2lgJRkc-mbZbytwXA1WY28Bf_8usFwBOWg$">https://urldefense.us/v3/__https://petsc.org/release/faq/*valgrind__;Iw!!G_uCfscf7eWS!b1Gx5q-NMarPE0nExBMOr0s4mG_Z0OJaWf6V7NyRuywDwA4IJ75v2lgJRkc-mbZbytwXA1WY28Bf_8usFwBOWg$</a>
> [0]PETSC ERROR: Fortran callback Id invalid
> [0]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!b1Gx5q-NMarPE0nExBMOr0s4mG_Z0OJaWf6V7NyRuywDwA4IJ75v2lgJRkc-mbZbytwXA1WY28Bf_8uqReGPEw$">https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!b1Gx5q-NMarPE0nExBMOr0s4mG_Z0OJaWf6V7NyRuywDwA4IJ75v2lgJRkc-mbZbytwXA1WY28Bf_8uqReGPEw$</a> for trouble shooting.
> [0]PETSC ERROR: Petsc Development GIT revision: v3.21.2-274-g55d8e0067c3 GIT Date: 2024-06-20 10:02:08 -0400
> [0]PETSC ERROR: /home/acro018/software/geothermal/waiwera/code/testing-debug-PETSc-main-barry/timestepper_test with PETSC_ARCH main-debug-barry on EN438880 by acro018 Wed Jun 26 16:30:51 2024
> [0]PETSC ERROR: Configure options: --with-x --download-hdf5 --download-zlib --download-netcdf --download-pnetcdf --download-exodusii --download-triangle --download-ptscotch --download-chaco --download-hypre
> [0]PETSC ERROR: #1 PetscObjectGetFortranCallback() at /home/acro018/software/PETSc/code/src/sys/objects/inherit.c:316
> [0]PETSC ERROR: #2 ../src/timestepper.F90:610
>
> --
> Dr Adrian Croucher
> Senior Research Fellow
> Department of Engineering Science
> Waipapa Taumata Rau / University of Auckland, New Zealand
> email: a.croucher@auckland.ac.nz
> tel: +64 (0)9 923 4611
>
</pre></body></html>