Hi Barry,
The call sequence happens across several routines. I believe it is as
follows:
call PetscInitialize(PETSC_NULL_CHARACTER, ierr)
call MPI_Comm_rank(PETSC_COMM_WORLD, rank, ierr)
call MPI_Comm_size(PETSC_COMM_WORLD, ntasks, ierr)
Kmat = PETSC_NULL_MAT
if(Kmat.eq.PETSC_NULL_MAT) then
call MatCreate(PETSC_COMM_WORLD,Kmat,ierr)
etc...
Kmat itself is declared in a module
module pfeapc
# include <petsc/finclude/petscksp.h>
use petscksp
implicit none
Vec :: rhs, sol, xvec
Vec :: yvec, zvec
Vec :: Mdiag, Msqrt
Mat :: Kmat, Mmat, Pmat
KSP :: kspsol
end module pfeapc
-
On Mon, Mar 24, 2025 at 7:15 AM Barry Smith <bsm...@petsc.dev> wrote:
>
> How do you declare and initialize the matrix in usolve.F before
> calling MatCreate()? You should not initialize it with any value before the
> call.
>
> Barry
>
>
> On Mar 23, 2025, at 11:10 PM, Sanjay Govindjee via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
> Barry,
>
> I now have a compiled version of my code using the main branch. When I
> run however I am getting an error in matcreate_( ) when I try to solve
> (actually just set up the matrix). The console window reports
>
> [0]PETSC ERROR: matcreate_() at
> /Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c:14 Cannot create
> PETSC_NULL_XXX object
> [3]PETSC ERROR: matcreate_() at
> /Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c:14 Cannot create
> PETSC_NULL_XXX object
> [2]PETSC ERROR: matcreate_() at
> /Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c:14 Cannot create
> PETSC_NULL_XXX object
> [1]PETSC ERROR: matcreate_() at
> /Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c:14 Cannot create
> PETSC_NULL_XXX object
>
> The debugger windows all report (modulo the pid):
>
> (lldb) process attach --pid 90952
> Process 90952 stopped
> * thread #1, queue = 'com.apple.main-thread', stop reason = signal SIGSTOP
> frame #0: 0x00007fff69d92746 libsystem_kernel.dylib`__semwait_signal +
> 10
> libsystem_kernel.dylib`__semwait_signal:
> -> 0x7fff69d92746 <+10>: jae 0x7fff69d92750 ; <+20>
> 0x7fff69d92748 <+12>: movq %rax, %rdi
> 0x7fff69d9274b <+15>: jmp 0x7fff69d9121d ; cerror
> 0x7fff69d92750 <+20>: retq
> Target 0: (feap) stopped.
> frame
>
> The debugger reports for the stack:
>
> (lldb) thread backtrace
> * thread #1, queue = 'com.apple.main-thread', stop reason = signal SIGSTOP
> * frame #0: 0x00007fff69d92746 libsystem_kernel.dylib`__semwait_signal +
> 10
> frame #1: 0x00007fff69d15eea libsystem_c.dylib`nanosleep + 196
> frame #2: 0x00007fff69d15d52 libsystem_c.dylib`sleep + 41
> frame #3: 0x0000000111acb04c libpetsc.3.022.dylib`PetscSleep(s=10) at
> psleep.c:48:5
> frame #4: 0x0000000111722961 libpetsc.3.022.dylib`PetscAttachDebugger
> at adebug.c:458:5
> frame #5: 0x00000001162ec7c8
> libpetsc.3.022.dylib`PetscAttachDebuggerErrorHandler(comm=0x000000011788fde8,
> line=14, fun="matcreate_",
> file="/Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c", num=62,
> p=PETSC_ERROR_INITIAL, mess="Cannot create PETSC_NULL_XXX object",
> ctx=0x0000000000000000) at adebug.c:522:9
> frame #6: 0x00000001162ecdb0
> libpetsc.3.022.dylib`PetscError(comm=0x000000011788fde8, line=14,
> func="matcreate_",
> file="/Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c", n=62,
> p=PETSC_ERROR_INITIAL, mess="Cannot create PETSC_NULL_XXX object") at
> err.c:409:15
> frame #7: 0x000000011233ecad
> libpetsc.3.022.dylib`matcreate_(a=0x00000001165fdc74, b=0x000000010fcde8c8,
> ierr=0x00007ffee0460308) at gcreatef.c:14:3
> frame #8: 0x000000010f7cf072 feap`usolve_ at usolve.F:138:72
> frame #9: 0x000000010f942de2 feap`presol_ at presol.f:181:72
> frame #10: 0x000000010f8cb8d8 feap`pmacr1_ at pmacr1.f:555:72
> frame #11: 0x000000010f8c60ad feap`pmacr_ at pmacr.f:614:72
> frame #12: 0x000000010f86ae4f feap`pcontr_ at pcontr.f:1375:72
> frame #13: 0x000000010fc1215e feap`main at feap87.f:173:72
> frame #14: 0x00007fff69c4ecc9 libdyld.dylib`start + 1
> frame #15: 0x00007fff69c4ecc9 libdyld.dylib`start + 1
>
> Here is a peek at the frame stack:
>
> frame #3: 0x0000000105c7104c libpetsc.3.022.dylib`PetscSleep(s=10) at
> psleep.c:48:5
> 45
> 46 #if defined(PETSC_HAVE_SLEEP)
> 47 else
> -> 48 sleep((int)s);
> 49 #elif defined(PETSC_HAVE__SLEEP) &&
> defined(PETSC_HAVE__SLEEP_MILISEC)
> 50 else _sleep((int)(s * 1000));
> 51 #elif defined(PETSC_HAVE__SLEEP)
> (lldb) up
> frame #4: 0x00000001058c8961 libpetsc.3.022.dylib`PetscAttachDebugger at
> adebug.c:458:5
> 455 while (left > 0) left = PetscSleep(left) - 1;
> 456 }
> 457 #else
> -> 458 PetscCall(PetscSleep(sleeptime));
> 459 #endif
> 460 }
> 461 #endif
> (lldb) up
> frame #5: 0x000000010a4927c8
> libpetsc.3.022.dylib`PetscAttachDebuggerErrorHandler(comm=0x000000010c7f1de8,
> line=14, fun="matcreate_",
> file="/Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c", num=62,
> p=PETSC_ERROR_INITIAL, mess="Cannot create PETSC_NULL_XXX object",
> ctx=0x0000000000000000) at adebug.c:522:9
> 519 if (fun) (void)(*PetscErrorPrintf)("%s() at %s:%d %s\n", fun,
> file, line, mess);
> 520 else (void)(*PetscErrorPrintf)("%s:%d %s\n", file, line,
> mess);
> 521
> -> 522 (void)PetscAttachDebugger();
> 523 abort(); /* call abort because don't want to kill other MPI
> ranks that may successfully attach to debugger */
> 524 PetscFunctionReturn(PETSC_SUCCESS);
> 525 }
> (lldb) up
> frame #6: 0x000000010a492db0
> libpetsc.3.022.dylib`PetscError(comm=0x000000010c7f1de8, line=14,
> func="matcreate_",
> file="/Users/sg/petsc-3.22.4main/gnu/ftn/mat/utils/gcreatef.c", n=62,
> p=PETSC_ERROR_INITIAL, mess="Cannot create PETSC_NULL_XXX object") at
> err.c:409:15
> 406 if (p == PETSC_ERROR_INITIAL && n != PETSC_ERR_MEMC)
> (void)PetscMallocValidate(__LINE__, PETSC_FUNCTION_NAME, __FILE__);
> 407
> 408 if (!eh) ierr = PetscTraceBackErrorHandler(comm, line, func,
> file, n, p, lbuf, NULL);
> -> 409 else ierr = (*eh->handler)(comm, line, func, file, n, p,
> lbuf, eh->ctx);
> 410 PetscStackClearTop;
> 411
> 412 /*
> (lldb) up
> frame #7: 0x00000001064e4cad
> libpetsc.3.022.dylib`matcreate_(a=0x000000010a7a3c74, b=0x0000000104e3e8c8,
> ierr=0x00007ffeeb300308) at gcreatef.c:14:3
> 11 PETSC_EXTERN void matcreate_(MPI_Fint *a, Mat *b,
> PetscErrorCode *ierr)
> 12 {
> 13 PetscBool null_b = !*(void**) b ? PETSC_TRUE : PETSC_FALSE;
> -> 14 PETSC_FORTRAN_OBJECT_CREATE(b);
> 15 CHKFORTRANNULLOBJECT(b);
> 16 *ierr = MatCreate(MPI_Comm_f2c(*(a)), b);
> 17 if (*ierr) return;
> (lldb) up
> frame #8: 0x000000010492f072 feap`usolve_ at usolve.F:138:72
> 135 onnz => mr(np(246):np(246)+ilist(2,246)-1)
> 136 dnnz => mr(np(247):np(247)+ilist(2,247)-1)
> 137
> -> 138 call MatCreate(PETSC_COMM_WORLD,Kmat,ierr)
> 139 call MatSetSizes(Kmat,numpeq,numpeq,PETSC_DETERMINE,
> 140 & PETSC_DETERMINE,ierr)
> 141 if(pfeap_bcin) call MatSetBlockSize(Kmat,nsbk,ierr)
>
>
> --
> -------------------------------------------------------------------
> Sanjay Govindjee, PhD, PE
> Horace, Dorothy, and Katherine Johnson Professor in Engineering
> Distinguished Professor of Civil and Environmental Engineering
>
> 779 Davis Hall
> University of California
> Berkeley, CA 94720-1710
>
> Voice: +1 510 642 6060
> FAX: +1 510 643 5264s_g@berkeley.eduhttp://faculty.ce.berkeley.edu/sanjay
> <https://urldefense.us/v3/__http://faculty.ce.berkeley.edu/sanjay__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLJDGn-rKw$>
> -------------------------------------------------------------------
>
> Books:
>
> Introduction to Mechanics of Solid
> Materialshttps://global.oup.com/academic/product/introduction-to-mechanics-of-solid-materials-9780192866080
>
> <https://urldefense.us/v3/__https://global.oup.com/academic/product/introduction-to-mechanics-of-solid-materials-9780192866080__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLK6eJ2xMA$>
>
> Continuum Mechanics of
> Solidshttps://global.oup.com/academic/product/continuum-mechanics-of-solids-9780198864721
>
> <https://urldefense.us/v3/__https://global.oup.com/academic/product/continuum-mechanics-of-solids-9780198864721__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLL3Y4t2YA$>
>
> Example Problems for Continuum Mechanics of
> Solidshttps://www.amazon.com/dp/1083047361/
> <https://urldefense.us/v3/__https://www.amazon.com/dp/1083047361/__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLIyAaK4Eg$>
>
> Engineering Mechanics of Deformable
> Solidshttps://www.amazon.com/dp/0199651647
> <https://urldefense.us/v3/__https://www.amazon.com/dp/0199651647__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLLf5ewSKQ$>
>
> Engineering Mechanics 3 (Dynamics) 2nd
> Editionhttp://www.amazon.com/dp/3642537111
> <https://urldefense.us/v3/__http://www.amazon.com/dp/3642537111__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLIJH2xWfA$>
>
> Engineering Mechanics 3, Supplementary Problems: Dynamics
> https://urldefense.us/v3/__http://www.amzn.com/B00SOXN8JU__;!!G_uCfscf7eWS!Ze1e1pn5iR1dyxVghSG5nRjhvjR5ODwW_5cnNeDlkxHfiBm7LC3cDe1n-nOt8RCETQzWjV0DIc1Il5Sei2SjZQ$
>
> <https://urldefense.us/v3/__http://www.amzn.com/B00SOXN8JU__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLJV9YZdRA$>
>
> -------------------------------------------------------------------
> NSF NHERI SimCenterhttps://simcenter.designsafe-ci.org/
> <https://urldefense.us/v3/__https://simcenter.designsafe-ci.org/__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLJMVyG6GQ$>
> -------------------------------------------------------------------
>
>
>
>
