[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: how to get all the singular values of the matrix?
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: how to get all the singular values of the matrix?
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Mon, 14 Jan 2008 00:06:19 -0600
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=uk6MC1zsV59EWlw4BYykyLGfB62FLG/DDFQVtg3+/SY=; b=hHm2D09upgr7KQ2k7GiSI9ByOiUU90j+gPAdwNoSLtnmYvu3y0TC1UnSzCrEBQaoMiJ/qu7ETQe+id7HXSTm7rxU5ODjyRAFA4gIGqUaGWrKsVDBZxZgZqSNeVScgT9vMF8JKLSQqXiEl6z2SSBQPHruXN2EKaEwRNdUd/i6Nys=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=CbBWbkrD/EkgEig+8f8vuIAbOahM8AT54KzzAne+FGRgqniKh/6XMZh91RJLP7Gy0nTxLNEE2kf91xX1qeOFfefipOMPn4Gm8nFPqb5oecaCNJAAX+JRw8ZNc13o6t4+3NbdP/gZxLCA06DVSTZYJnbN1Td/TmAAJTR/J0nmO3A=
- In-reply-to: <7ff0ee010801132123n5cc5194cld73b2d958be957b9@mail.gmail.com>
- References: <7ff0ee010801132123n5cc5194cld73b2d958be957b9@mail.gmail.com>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
You can use
http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/KSP/KSPComputeEigenvaluesExplicitly.html
with and without a preconditioner. We have not coded the SVD
counterpart, but you can use
http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/KSP/KSPComputeExplicitOperator.html
and then call the LAPACK yourself.
Matt
On Jan 13, 2008 11:23 PM, Yujie <recrusader@xxxxxxxxx> wrote:
> Hi, everyone
>
> I want to select iterative methods by observing the singular values
> decompostion of the matrix. However, I don't know how to get all the
> singular values of the matrix in PETSc. I know the command
> "-ksp_monitor_singular_value" may get the max and min singular values at
> each iteration. How to get the singular values of the matrix I want to
> solve? In addition, when I use the preconditioned iterative method, how to
> get the singular values of the preconditioned iterative operator?
>
> thanks a lot.
>
> Regards,
> Yujie
>
--
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