On Tue, Jun 16, 2009 at 8:30 AM, tog <span dir="ltr"><<a href="mailto:guillaume.alleon@gmail.com">guillaume.alleon@gmail.com</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;">
Thanks Matthew,<br><br>Well I am interested by the parallel case unfortunately ;)<br>Any reference paper that could help me to implement this ?</blockquote><div><br>Maybe you can use STAPL (<a href="http://parasol.tamu.edu/stapl/">http://parasol.tamu.edu/stapl/</a>)<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>Guillaume<div><div></div><div class="h5"><br><br><div class="gmail_quote">
On Tue, Jun 16, 2009 at 8:56 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div>On Tue, Jun 16, 2009 at 5:38 AM, tog <span dir="ltr"><<a href="mailto:guillaume.alleon@gmail.com" target="_blank">guillaume.alleon@gmail.com</a>></span> wrote:<br>
</div><div class="gmail_quote"><div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi there<br><br>This is a kind of newbie question :)<br>I have a vector which is the solution of my simulation, I want then to extract the p indices corresponding to the p biggest values in my vector. Is there a way to do this with PETSc ?</blockquote>
</div><div><br>You can get p=1 with VecMax().<br> </div><div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
What if I want to sort my vector to get the indices in ascending order of their values ?</blockquote></div><div><br>The p=2+ and permutation cases are sophisticated in parallel and not present. In serial, you can use<br>
VecGetArray() and PetscSortReal()/SortRealWithPermutation().<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>Thanks<br>Guillaume<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>
</blockquote></div><br><br clear="all"><br></div></div><font color="#888888">-- <br>PGP KeyID: 1024D/69B00854 <a href="http://subkeys.pgp.net" target="_blank">subkeys.pgp.net</a><br><br><a href="http://cheztog.blogspot.com" target="_blank">http://cheztog.blogspot.com</a><br>
</font></blockquote></div><br><br clear="all"><br>-- <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>