<div dir="ltr">On Sat, Feb 16, 2013 at 5:46 AM, Hui Zhang <span dir="ltr"><<a href="mailto:mike.hui.zhang@hotmail.com" target="_blank">mike.hui.zhang@hotmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I updated petsc-dev just now with the last changeset:   26190:c5c1560f6126<br>
<br>
I reconfigure and rebuild petsc.  I modified src/sys/threadcomm/examples/tutorials/ex4.c to include a VecDuplicate like the following ediff file (see behind, ! indicating the modifications).<br>
<br>
Then I got error messages (using mpiexec -np 1)<br>
<br>
[0]PETSC ERROR: PetscCommGetThreadComm() line 111 in /Users/huizhang/Software/petsc-dev/src/sys/threadcomm/interface/threadcomm.c<br>
[0]PETSC ERROR: PetscThreadCommRunKernel0() line 792 in /Users/huizhang/Software/petsc-dev/src/sys/threadcomm/interface/threadcomm.c<br>
[0]PETSC ERROR: PetscThreadCommStackDestroy() line 197 in /Users/huizhang/Software/petsc-dev/src/sys/threadcomm/interface/threadcomm.c<br>
[0]PETSC ERROR: PetscThreadCommDestroy() line 242 in /Users/huizhang/Software/petsc-dev/src/sys/threadcomm/interface/threadcomm.c<br>
[0]PETSC ERROR: Petsc_DelThreadComm() line 71 in /Users/huizhang/Software/petsc-dev/src/sys/threadcomm/interface/dlregisthreadcomm.c<br>
[0]PETSC ERROR: PetscFinalize() line 1269 in /Users/huizhang/Software/petsc-dev/src/sys/objects/pinit.c<br></blockquote><div><br></div><div style>You failed to VecDestroy the vector.</div><div style><br></div><div style>
   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Using more processes gives such error messages for each process (or thread?).  I use scalar-tye complex number.<br>
<br>
I also got the same error messages from my application programs (still with scalar-type complex), in which I do not include petscthreadcomm.h.<br>
<br>
<br>
*** 12,18 ****<br>
  {<br>
    PetscErrorCode ierr;<br>
    PetscScalar    dot=0.0,v;<br>
!   Vec            x,y;<br>
    PetscInt       N  =8;<br>
    PetscScalar    one=1.0,two=2.0,alpha=2.0;<br>
<br>
--- 12,18 ----<br>
  {<br>
    PetscErrorCode ierr;<br>
    PetscScalar    dot=0.0,v;<br>
!   Vec            x,y,z;<br>
    PetscInt       N  =8;<br>
    PetscScalar    one=1.0,two=2.0,alpha=2.0;<br>
<br>
*** 26,32 ****<br>
    ierr = VecSetFromOptions(x);CHKERRQ(ierr);<br>
    ierr = VecSet(x,one);CHKERRQ(ierr);<br>
    ierr = PetscPrintf(PETSC_COMM_WORLD,"x = %lf\n",PetscRealPart(one));CHKERRQ(ierr);<br>
!<br>
    ierr = VecCreate(PETSC_COMM_WORLD,&y);CHKERRQ(ierr);<br>
    ierr = VecSetSizes(y,PETSC_DECIDE,N);CHKERRQ(ierr);<br>
    ierr = VecSetFromOptions(y);CHKERRQ(ierr);<br>
--- 26,32 ----<br>
    ierr = VecSetFromOptions(x);CHKERRQ(ierr);<br>
    ierr = VecSet(x,one);CHKERRQ(ierr);<br>
    ierr = PetscPrintf(PETSC_COMM_WORLD,"x = %lf\n",PetscRealPart(one));CHKERRQ(ierr);<br>
!   ierr = VecDuplicate(x,&z); CHKERRQ(z);<br>
    ierr = VecCreate(PETSC_COMM_WORLD,&y);CHKERRQ(ierr);<br>
    ierr = VecSetSizes(y,PETSC_DECIDE,N);CHKERRQ(ierr);<br>
    ierr = VecSetFromOptions(y);CHKERRQ(ierr);<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>