[petsc-users] GCREATEF.C error

Sanjay Govindjee s_g at berkeley.edu
Sun Mar 23 22:10:08 CDT 2025


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 5264
s_g at berkeley.edu
https://urldefense.us/v3/__http://faculty.ce.berkeley.edu/sanjay__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLJDGn-rKw$ 
-------------------------------------------------------------------

Books:

Introduction to Mechanics of Solid Materials
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 Solids
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 Solids
https://urldefense.us/v3/__https://www.amazon.com/dp/1083047361/__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLIyAaK4Eg$ 

Engineering Mechanics of Deformable Solids
https://urldefense.us/v3/__https://www.amazon.com/dp/0199651647__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLLf5ewSKQ$ 

Engineering Mechanics 3 (Dynamics) 2nd Edition
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!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLJV9YZdRA$ 

-------------------------------------------------------------------
NSF NHERI SimCenter
https://urldefense.us/v3/__https://simcenter.designsafe-ci.org/__;!!G_uCfscf7eWS!d3a3dMZHNHbpo2pUZUkrJ_sVG2VKWJenpA8SfSld478I1qG5iDDI_A5HMLojdmimArbDAt-E50BvvLJMVyG6GQ$ 
-------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20250323/7f00c0c7/attachment.html>


More information about the petsc-users mailing list