Thanks Matthew,<br><br>Well I am interested by the parallel case unfortunately ;)<br>Any reference paper that could help me to implement this ?<br><br>Guillaume<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">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 class="im">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 class="im"><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 class="im"><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>-- <br>PGP KeyID: 1024D/69B00854 <a href="http://subkeys.pgp.net">subkeys.pgp.net</a><br><br><a href="http://cheztog.blogspot.com">http://cheztog.blogspot.com</a><br>