[petsc-dev] New error in SNESComputeJacobianDefaultColor()

[Barry Smith]

   Adrian,

    Sorry for the delay. I think it is now fixed in the branch 

    Since I rebased the branch it is simplest if you do in the PETSc directory

    git checkout main
    git fetch
    git branch -D barry/2024-06-14/fix-fortransnessetjacobianmatfdcoloring
  git checkout barry/2024-06-14/fix-fortransnessetjacobianmatfdcoloring

to get the changes

Barry


> On Jun 16, 2024, at 7:30 PM, Adrian Croucher <a.croucher at auckland.ac.nz> wrote:
> 
> hi Barry,
> 
> On 15/06/24 7:56 am, Barry Smith wrote:
>> 
>>    Thanks for reporting the problem. It was an oversight on my part when working on the Fortran interfaces.
>> 
>>     It should work now in https://urldefense.us/v3/__https://gitlab.com/petsc/petsc/-/merge_requests/7631__;!!G_uCfscf7eWS!YpdXxzYPyXZ9AndH2YvFeoREgCWxzvg8SJuEPYxwVPskWkDOUVW1pZ32XZta3qnJcxaPOA_RDMxbaoWadUZvZw$ 
>> 
> Thanks, I think that's fixed that problem, but it looks like there might be another one now. When I build my code I'm getting an error when I call SNESSetJacobian():
> 
> 1565 |          SNES_Jacobian, self%context, ierr)
>       |                                           1
> Error: There is no specific subroutine for the generic ‘snessetjacobian’ at (1)
> 
> The self%context object I'm passing in for the ctx argument is a derived type, which used to work ok.
> 
> - Adrian
> 
>>    Barry
>> 
>> 
>>> On Jun 14, 2024, at 1:29 AM, Adrian Croucher <a.croucher at auckland.ac.nz> <mailto:a.croucher at auckland.ac.nz> wrote:
>>> 
>>> This Message Is From an External Sender
>>> This message came from outside your organization.
>>> hi,
>>> 
>>> I just tried updating to the latest PETSc 'main' branch. Since I'm using 
>>> Fortran I had to modify a few bits of my code for the new type checking 
>>> etc. recently introduced here:
>>> 
>>> https://urldefense.us/v3/__https://gitlab.com/petsc/petsc/-/commit/5d83a8b16d06840f96948f1a43aa9c83c769a60a__;!!G_uCfscf7eWS!agn9OOuCAtlXUe5B5aW7NZLJZ-8jp-mPz7HxVYUCk8_iplZ6EstbNYXU3HX2ndPpAI2v6OW7qKKmtSXaUlaN_K3CI_9LIOA$
>>> 
>>> That was ok and my code now builds again but I am getting a new run-time 
>>> error (see below). It seems to be complaining about the type of the 
>>> argument 'ctx' in SNESComputeJacobianDefaultColor(). I am passing a 
>>> MatFDColoring object in for that. Up to now that has been ok.
>>> 
>>> Has something changed with the way that function has to be called?
>>> 
>>> - Adrian
>>> 
>>> --
>>> 
>>> [0]PETSC ERROR: Invalid argument
>>> [0]PETSC ERROR: Wrong type of object: Parameter # 5
>>> [0]PETSC ERROR: See https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!agn9OOuCAtlXUe5B5aW7NZLJZ-8jp-mPz7HxVYUCk8_iplZ6EstbNYXU3HX2ndPpAI2v6OW7qKKmtSXaUlaN_K3CbpQVJDw$ for trouble shooting.
>>> [0]PETSC ERROR: Petsc Development GIT revision: 
>>> v3.21.2-195-g3719fecf4fb  GIT Date: 2024-06-03 17:57:55 +0000
>>> [0]PETSC ERROR: ./timestepper_test with PETSC_ARCH main-debug on 
>>> EN438880 by acro018 Fri Jun 14 16:32:19 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 SNESComputeJacobianDefaultColor() at 
>>> /home/acro018/software/PETSc/code/src/snes/interface/snesj2.c:72
>>> [0]PETSC ERROR: #2 ../src/timestepper.F90:639
>>> [0]PETSC ERROR: #3 oursnesjacobian() at 
>>> /home/acro018/software/PETSc/code/src/snes/interface/ftn-custom/zsnesf.c:176
>>> [0]PETSC ERROR: #4 SNESComputeJacobian() at 
>>> /home/acro018/software/PETSc/code/src/snes/interface/snes.c:2888
>>> [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() at 
>>> /home/acro018/software/PETSc/code/src/snes/impls/ls/ls.c:217
>>> [0]PETSC ERROR: #6 SNESSolve() at 
>>> /home/acro018/software/PETSc/code/src/snes/interface/snes.c:4755
>>> [0]PETSC ERROR: #7 ../src/timestepper.F90:2295
>>> 
>>> -- 
>>> Dr Adrian Croucher
>>> Senior Research Fellow
>>> Department of Engineering Science
>>> Waipapa Taumata Rau / University of Auckland, New Zealand
>>> email: a.croucher at auckland.ac.nz <mailto:a.croucher at auckland.ac.nz>
>>> tel: +64 (0)9 923 4611
>>> 
>>> 
>> 
> -- 
> Dr Adrian Croucher
> Senior Research Fellow
> Department of Engineering Science
> Waipapa Taumata Rau / University of Auckland, New Zealand
> email: a.croucher at auckland.ac.nz <mailto:a.croucher at auckland.ac.nz>
> tel: +64 (0)9 923 4611

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20240619/41e50ffa/attachment-0001.html>