On Wed, Mar 25, 2009 at 1:05 AM, Khan, Irfan <span dir="ltr"><<a href="mailto:irfan.khan@gatech.edu">irfan.khan@gatech.edu</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi<br>
Can the petsc parallel vectors be used with two different communicators? For instance, I have created two different communicators called FEA_Comm and FSI_Comm. The total number of processes are x+y. FSI_Comm works on x+y but FEA_Comm works only on x.<br>
<br>
Now I am trying to create parallel vectors a1 and a2 such that a1 has entries from x+y processes but a2 has entries from only y processes.<br>
<br>
After splitting the communicators I assign PETSC_COMM_WORLD to FEA_Comm which works on only x processes. Subsequently petsc is initialized (PetscInitialize()). But when the parallel vectors are created, the processes hang.</blockquote>
<div><br>PETSC_COMM_WORLD should encompass all processes you wish to use in PETSc, so that means x+y. You can create Vec<br>objects on subcommunicators, like x.<br><br> Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Any suggestions will be helpful<br>
<br>
Thankyou<br>
Irfan<br>
Graduate Research Assistant<br>
Woodruff school of Mechanical Engineering<br>
Atlanta, GA (30307)<br>
</blockquote></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<br>