[petsc-users] GCREATEF.C error
Sanjay Govindjee
s_g at berkeley.edu
Mon Mar 24 13:40:47 CDT 2025
Thanks! I look forward to seeing the Change Note :)
Time now to perform some search and destroy through the code!
-sanjay
-------------------------------------------------------------------
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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssWC_usPow$
-------------------------------------------------------------------
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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssXfRsAazA$
Continuum Mechanics of Solids
https://urldefense.us/v3/__https://global.oup.com/academic/product/continuum-mechanics-of-solids-9780198864721__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssXoPG4Icg$
Example Problems for Continuum Mechanics of Solids
https://urldefense.us/v3/__https://www.amazon.com/dp/1083047361/__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssV55pN4oA$
Engineering Mechanics of Deformable Solids
https://urldefense.us/v3/__https://www.amazon.com/dp/0199651647__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssUkQt68JQ$
Engineering Mechanics 3 (Dynamics) 2nd Edition
https://urldefense.us/v3/__http://www.amazon.com/dp/3642537111__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssVLab2hTw$
Engineering Mechanics 3, Supplementary Problems: Dynamics
https://urldefense.us/v3/__http://www.amzn.com/B00SOXN8JU__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssWTSu_3Cw$
-------------------------------------------------------------------
NSF NHERI SimCenter
https://urldefense.us/v3/__https://simcenter.designsafe-ci.org/__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssVRW5b4BA$
-------------------------------------------------------------------
On 3/24/25 11:33 AM, Barry Smith wrote:
> Ahh, you cannot use this type of construct
>
> if (x .eq. PETSC_NULL_XXX) then
>
> anymore.
>
> Nor can you ever assign x = PETSC_NULL_XXXX
>
> Instead use
>
> if (PetscObjectIsNull(x)) then
>
> or
>
> if (.not. PetscObjectIsNull(x)) then
>
> Also all PETSc objects are now null when they are declared so just
> write
>
> Mat Kmat
> if (PetscObjectIsNull(Kmat)) then
> call MatCreate(PETSC_COMM_WORLD,Kmat,ierr)
>> etc...
>
> Barry
>
>
>
>> On Mar 24, 2025, at 11:14 AM, Sanjay Govindjee <s_g at berkeley.edu> wrote:
>>
>> 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 <bsmith at 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 at 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 5264
>>> s_g at berkeley.edu
>>> https://urldefense.us/v3/__http://faculty.ce.berkeley.edu/sanjay__;!!G_uCfscf7eWS!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssWC_usPow$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssXfRsAazA$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssXoPG4Icg$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssV55pN4oA$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssUkQt68JQ$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssVLab2hTw$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssWTSu_3Cw$ <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!eHYMl71AKju95WDkDP_equtIeYH7RhHSzQqTS35mJGW2jkTzRvu62AtY2nbXhngzrBkB78YbsJ62ssVRW5b4BA$ <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/20250324/8d6715d8/attachment-0001.html>
More information about the petsc-users
mailing list