Sorting a vector

[tog]

Thanks Matthew,

Well I am interested by the parallel case unfortunately ;)
Any reference paper that could help me to implement this ?

Guillaume

On Tue, Jun 16, 2009 at 8:56 PM, Matthew Knepley <knepley at gmail.com> wrote:

> On Tue, Jun 16, 2009 at 5:38 AM, tog <guillaume.alleon at gmail.com> wrote:
>
>> Hi there
>>
>> This is a kind of newbie question :)
>> 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 ?
>
>
> You can get p=1 with VecMax().
>
>
>>
>> What if I want to sort my vector to get the indices in ascending order of
>> their values ?
>
>
> The p=2+ and permutation cases are sophisticated in parallel and not
> present. In serial, you can use
> VecGetArray() and PetscSortReal()/SortRealWithPermutation().
>
>   Matt
>
>
>>
>> Thanks
>> Guillaume
>>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>



-- 
PGP KeyID: 1024D/69B00854  subkeys.pgp.net

http://cheztog.blogspot.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090616/9e6bfb95/attachment.htm>