<!DOCTYPE html>
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body>
My odessy to update my code to the main branch continues...I am now
encountering an error with VecSetValue:<br>
<br>
<blockquote>[0]PETSC ERROR: VecSetValues() at
/Users/sg/petsc-3.22.4main/src/vec/vec/interface/rvector.c:926
Null Pointer: Parameter # 4<br>
</blockquote>
This is on a single processor run -- I got tired of closing all
those debug windows. My call looks like<br>
<blockquote>call VecSetValue(rhs, k-1, b(j), ADD_VALUES, ierr)<br>
</blockquote>
and the first two data types are Vec and integer, and for the third<br>
<blockquote>real (kind=8) :: b(:) ! comes through the argument list
for the subroutine<br>
</blockquote>
vecsetvalue_ however seems to be receiving va=0x0000000000000000. I
see that formally the value is supposed to be PetscScalar; is that
the issue? The thread backtrace is below.<br>
<blockquote>(lldb) process attach --pid 51640<br>
Process 51640 stopped<br>
* thread #1, queue = 'com.apple.main-thread', stop reason = signal
SIGSTOP<br>
frame #0: 0x00007fff69d92746
libsystem_kernel.dylib`__semwait_signal + 10<br>
libsystem_kernel.dylib`__semwait_signal:<br>
-> 0x7fff69d92746 <+10>: jae
0x7fff69d92750 ; <+20><br>
0x7fff69d92748 <+12>: movq %rax, %rdi<br>
0x7fff69d9274b <+15>: jmp 0x7fff69d9121d ;
cerror<br>
0x7fff69d92750 <+20>: retq <br>
Target 0: (feap) stopped.<br>
<br>
Executable module set to "/Users/sg/Feap/ver87/parfeap/feap".<br>
Architecture set to: x86_64h-apple-macosx-.<br>
(lldb) thread backtrace<br>
* thread #1, queue = 'com.apple.main-thread', stop reason = signal
SIGSTOP<br>
* frame #0: 0x00007fff69d92746
libsystem_kernel.dylib`__semwait_signal + 10<br>
frame #1: 0x00007fff69d15eea libsystem_c.dylib`nanosleep + 196<br>
frame #2: 0x00007fff69d15d52 libsystem_c.dylib`sleep + 41<br>
frame #3: 0x000000010a1fb04c
libpetsc.3.022.dylib`PetscSleep(s=10) at psleep.c:48:5<br>
frame #4: 0x0000000109e52961
libpetsc.3.022.dylib`PetscAttachDebugger at adebug.c:458:5<br>
frame #5: 0x000000010ea1c7c8
libpetsc.3.022.dylib`PetscAttachDebuggerErrorHandler(comm=0x0000000110e78de8,
line=926, fun="VecSetValues",
file="/Users/sg/petsc-3.22.4main/src/vec/vec/interface/rvector.c",
num=85, p=PETSC_ERROR_INITIAL, mess="Null Pointer: Parameter # 4",
ctx=0x0000000000000000) at adebug.c:522:9<br>
frame #6: 0x000000010ea1cdb0
libpetsc.3.022.dylib`PetscError(comm=0x0000000110e78de8, line=926,
func="VecSetValues",
file="/Users/sg/petsc-3.22.4main/src/vec/vec/interface/rvector.c",
n=85, p=PETSC_ERROR_INITIAL, mess="Null Pointer: Parameter # %d")
at err.c:409:15<br>
frame #7: 0x000000010a866bf2
libpetsc.3.022.dylib`VecSetValues(x=0x00007fdb53068450, ni=1,
ix=0x00007ffee6d74098, y=0x0000000000000000, iora=ADD_VALUES) at
rvector.c:926:3<br>
frame #8: 0x000000010a82c1dc
libpetsc.3.022.dylib`vecsetvalue_(v=0x00000001093ca8f8,
i=0x00007ffee6d74098, va=0x0000000000000000,
mode=0x00007ffee6d74094, ierr=0x00007ffee6d74238) at
zvectorf.c:20:11<br>
frame #9: 0x0000000108ebbb48 feap`usolve_ at usolve.F:256:72<br>
frame #10: 0x000000010904cd9c feap`psolve_ at psolve.f:232:72<br>
frame #11: 0x0000000108fb7edc feap`pmacr1_ at pmacr1.f:667:72<br>
frame #12: 0x0000000108fb20dd feap`pmacr_ at pmacr.f:614:72<br>
frame #13: 0x0000000108f56e7f feap`pcontr_ at pcontr.f:1375:72<br>
frame #14: 0x00000001092fe18e feap`main at feap87.f:173:72<br>
frame #15: 0x00007fff69c4ecc9 libdyld.dylib`start + 1<br>
(lldb) up<br>
frame #1: 0x00007fff69d15eea libsystem_c.dylib`nanosleep + 196<br>
libsystem_c.dylib`nanosleep:<br>
-> 0x7fff69d15eea <+196>: testl %eax, %eax<br>
0x7fff69d15eec <+198>: jns 0x7fff69d15eb5
; <+143><br>
0x7fff69d15eee <+200>: callq 0x7fff69d1e008
; symbol stub for: __error<br>
0x7fff69d15ef3 <+205>: cmpl $0x3c, (%rax)<br>
(lldb) up<br>
frame #2: 0x00007fff69d15d52 libsystem_c.dylib`sleep + 41<br>
libsystem_c.dylib`sleep:<br>
-> 0x7fff69d15d52 <+41>: movl %eax, %ecx<br>
0x7fff69d15d54 <+43>: xorl %eax, %eax<br>
0x7fff69d15d56 <+45>: cmpl $-0x1, %ecx<br>
0x7fff69d15d59 <+48>: jne 0x7fff69d15d85 ;
<+92><br>
(lldb) up<br>
frame #3: 0x000000010a1fb04c libpetsc.3.022.dylib`PetscSleep(s=10)
at psleep.c:48:5<br>
45 <br>
46 #if defined(PETSC_HAVE_SLEEP)<br>
47 else<br>
-> 48 sleep((int)s);<br>
49 #elif defined(PETSC_HAVE__SLEEP) &&
defined(PETSC_HAVE__SLEEP_MILISEC)<br>
50 else _sleep((int)(s * 1000));<br>
51 #elif defined(PETSC_HAVE__SLEEP)<br>
(lldb) up<br>
frame #4: 0x0000000109e52961
libpetsc.3.022.dylib`PetscAttachDebugger at adebug.c:458:5<br>
455 while (left > 0) left = PetscSleep(left) - 1;<br>
456 }<br>
457 #else<br>
-> 458 PetscCall(PetscSleep(sleeptime));<br>
459 #endif<br>
460 }<br>
461 #endif<br>
(lldb) up<br>
frame #5: 0x000000010ea1c7c8
libpetsc.3.022.dylib`PetscAttachDebuggerErrorHandler(comm=0x0000000110e78de8,
line=926, fun="VecSetValues",
file="/Users/sg/petsc-3.22.4main/src/vec/vec/interface/rvector.c",
num=85, p=PETSC_ERROR_INITIAL, mess="Null Pointer: Parameter # 4",
ctx=0x0000000000000000) at adebug.c:522:9<br>
519 if (fun) (void)(*PetscErrorPrintf)("%s() at %s:%d
%s\n", fun, file, line, mess);<br>
520 else (void)(*PetscErrorPrintf)("%s:%d %s\n", file,
line, mess);<br>
521 <br>
-> 522 (void)PetscAttachDebugger();<br>
523 abort(); /* call abort because don't want to kill
other MPI ranks that may successfully attach to debugger */<br>
524 PetscFunctionReturn(PETSC_SUCCESS);<br>
525 }<br>
(lldb) up<br>
frame #6: 0x000000010ea1cdb0
libpetsc.3.022.dylib`PetscError(comm=0x0000000110e78de8, line=926,
func="VecSetValues",
file="/Users/sg/petsc-3.22.4main/src/vec/vec/interface/rvector.c",
n=85, p=PETSC_ERROR_INITIAL, mess="Null Pointer: Parameter # %d")
at err.c:409:15<br>
406 if (p == PETSC_ERROR_INITIAL && n !=
PETSC_ERR_MEMC) (void)PetscMallocValidate(__LINE__,
PETSC_FUNCTION_NAME, __FILE__);<br>
407 <br>
408 if (!eh) ierr = PetscTraceBackErrorHandler(comm,
line, func, file, n, p, lbuf, NULL);<br>
-> 409 else ierr = (*eh->handler)(comm, line, func,
file, n, p, lbuf, eh->ctx);<br>
410 PetscStackClearTop;<br>
411 <br>
412 /*<br>
(lldb) up<br>
frame #7: 0x000000010a866bf2
libpetsc.3.022.dylib`VecSetValues(x=0x00007fdb53068450, ni=1,
ix=0x00007ffee6d74098, y=0x0000000000000000, iora=ADD_VALUES) at
rvector.c:926:3<br>
923 PetscValidHeaderSpecific(x, VEC_CLASSID, 1);<br>
924 if (!ni) PetscFunctionReturn(PETSC_SUCCESS);<br>
925 PetscAssertPointer(ix, 3);<br>
-> 926 PetscAssertPointer(y, 4);<br>
927 PetscValidType(x, 1);<br>
928 <br>
929 PetscCall(PetscLogEventBegin(VEC_SetValues, x, 0, 0,
0));<br>
(lldb) up<br>
frame #8: 0x000000010a82c1dc
libpetsc.3.022.dylib`vecsetvalue_(v=0x00000001093ca8f8,
i=0x00007ffee6d74098, va=0x0000000000000000,
mode=0x00007ffee6d74094, ierr=0x00007ffee6d74238) at
zvectorf.c:20:11<br>
17 PETSC_EXTERN void vecsetvalue_(Vec *v, PetscInt *i,
PetscScalar *va, InsertMode *mode, PetscErrorCode *ierr)<br>
18 {<br>
19 /* cannot use VecSetValue() here since that uses
PetscCall() which has a return in it */<br>
-> 20 *ierr = VecSetValues(*v, 1, i, va, *mode);<br>
21 }<br>
22 <br>
23 PETSC_EXTERN void vecsetvaluelocal_(Vec *v, PetscInt
*i, PetscScalar *va, InsertMode *mode, PetscErrorCode *ierr)<br>
(lldb) up<br>
frame #9: 0x0000000108ebbb48 feap`usolve_ at usolve.F:256:72<br>
253 if( k .gt. 0 ) then<br>
254 j = j + 1<br>
255 ! write(*,*) 'rank: ',rank,' k,j,b(j)
',k,j,b(j)<br>
-> 256 call VecSetValue(rhs, k-1, b(j),
ADD_VALUES, ierr)<br>
257 endif<br>
258 end do<br>
</blockquote>
<pre class="moz-signature" cols="72">--
-------------------------------------------------------------------
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
<a class="moz-txt-link-abbreviated" href="mailto:s_g@berkeley.edu">s_g@berkeley.edu</a>
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__http://faculty.ce.berkeley.edu/sanjay__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w9azVgSHQ$">http://faculty.ce.berkeley.edu/sanjay</a>
-------------------------------------------------------------------
Books:
Introduction to Mechanics of Solid Materials
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__https://global.oup.com/academic/product/introduction-to-mechanics-of-solid-materials-9780192866080__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w_VG5HHMA$">https://global.oup.com/academic/product/introduction-to-mechanics-of-solid-materials-9780192866080</a>
Continuum Mechanics of Solids
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__https://global.oup.com/academic/product/continuum-mechanics-of-solids-9780198864721__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w9T8MhzbQ$">https://global.oup.com/academic/product/continuum-mechanics-of-solids-9780198864721</a>
Example Problems for Continuum Mechanics of Solids
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__https://www.amazon.com/dp/1083047361/__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w_GqMlonw$">https://www.amazon.com/dp/1083047361/</a>
Engineering Mechanics of Deformable Solids
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__https://www.amazon.com/dp/0199651647__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w9troczOg$">https://www.amazon.com/dp/0199651647</a>
Engineering Mechanics 3 (Dynamics) 2nd Edition
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__http://www.amazon.com/dp/3642537111__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w_g7GDa1w$">http://www.amazon.com/dp/3642537111</a>
Engineering Mechanics 3, Supplementary Problems: Dynamics
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__http://www.amzn.com/B00SOXN8JU__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w97Vbz6Dw$">http://www.amzn.com/B00SOXN8JU</a>
-------------------------------------------------------------------
NSF NHERI SimCenter
<a class="moz-txt-link-freetext" href="https://urldefense.us/v3/__https://simcenter.designsafe-ci.org/__;!!G_uCfscf7eWS!eYHZJmdbbAvMW_dwF5Gl5kuFLWS9Tllmp4OCyVBwo0hVI25qi3ZiV-maVgXQ7wmZ8dKjEGu1hb5u7w_SCyLZOg$">https://simcenter.designsafe-ci.org/</a>
-------------------------------------------------------------------
</pre>
</body>
</html>