Re: [petsc-users] Problem with composite DM index sets

2025-06-23 Thread Randall Mackie
Hi Barry, That was a good question, and I didn’t know the answer but I did a little digging, and in fact, in the 2018 Fortran standard, an extended form of implicit none was introduced that disable also implicit procedures: implicit none (type, external) Indeed, when adding that to the little

Re: [petsc-users] Problem with composite DM index sets

2025-06-22 Thread Barry Smith
implicit none prevents using undeclared variables. Is there a way to avoid calling any functions/suboutines that don't have an interface declared in a module? That would have found this problem. Barry > On Jun 22, 2025, at 8:06 PM, Randall Mackie wrote: > > Thanks Barry! > > > >> On

Re: [petsc-users] Problem with composite DM index sets

2025-06-22 Thread Randall Mackie
Thanks Barry! > On Jun 22, 2025, at 4:50 PM, Barry Smith wrote: > > >I'm sorry for not getting back to you sooner. I have attached a working > version of your code. Since you were missing > > use petscdmcomposite > > the compiler could not generate the correct call to DMCompositeGet

Re: [petsc-users] Problem with composite DM index sets

2025-06-22 Thread Barry Smith
   I'm sorry for not getting back to you sooner. I have attached a working version of your code. Since you were missing   use petscdmcompositethe compiler could not generate the correct call to DMCompositeGetGlobalISs()  Barry test.F90 Description: Binary data On Jun 17, 2025, at 6:39 PM, Randall

Re: [petsc-users] Problem with composite DM index sets

2025-06-20 Thread Randall Mackie
Dear PETSc, On Tuesday I sent a small reproducer program that shows a problem with the Fortran version of DMCompositeGetGlobalISs, with the valgrind stack trace below. Not having yet received a reply, I’m just wondering if perhaps my email slipped through the cracks. Would appreciate some guid

Re: [petsc-users] Problem with composite DM index sets

2025-06-18 Thread Randall Mackie
Follow up: running with valgrind shows the following issues….is this a bug in PETSc? ==5216== Use of uninitialised value of size 8 ==5216==at 0x49ED69C: f90array1dcreatefortranaddr_ (f90_fwrap.F90:52) ==5216==by 0x4D7EA94: F90Array1dCreate (f90_cwrap.c:140) ==5216==by 0x6B7295A: dmcom

[petsc-users] Problem with composite DM index sets

2025-06-17 Thread Randall Mackie
Dear Petsc users -I am trying to upgrade my code to petsc-3.23 (from 3.19), and I seem to have run into a problem with DMCompositeGetGlobalISs.The example program listed on the man page for DMCompositeGetGlobalISs, https://petsc.org/release/src/snes/tutorials/ex73f90t.F90.html, seems to indicate th