Dear Matt and Hong:<br><br>Based what you said, it looks like a little difficult to evalute the matrix in PETSc, especailly regarding a big dimension. However, when I select iterative methods, how to select a suitable one based on some evaluation? Could you give me some advice? thanks a lot.
<br><br>Regards,<br>Yujie<br><br><div><span class="gmail_quote">On 1/14/08, <b class="gmail_sendername">Hong Zhang</b> &lt;<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>&gt; wrote:</span><blockquote class="gmail_quote" style="margin-top: 0; margin-right: 0; margin-bottom: 0; margin-left: 0; margin-left: 0.80ex; border-left-color: #cccccc; border-left-width: 1px; border-left-style: solid; padding-left: 1ex">
<br>If you want few selected eigen solutions of sparse matrix,<br>you should use sparse eigen solver. Take a look at&#39;<br>slepc (<a href="http://www.grycap.upv.es/slepc/">http://www.grycap.upv.es/slepc/</a>)<br>or use slepc interface with arpack.
<br><br>Hong<br><br><br>On Mon, 14 Jan 2008, Yujie wrote:<br><br>&gt; Thank you for your advice.<br>&gt; I have used -ksp_compute_eigenvalues_explicitly to get the eigen values.<br>&gt; However, it is very very<br>&gt; slow because the dimension of the matrix is about ten thousand.
<br>&gt;<br>&gt; Yujie<br>&gt;<br>&gt; On 1/14/08, Matthew Knepley &lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt; wrote:<br>&gt;&gt;<br>&gt;&gt; You can use<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; <a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/KSP/KSPComputeEigenvaluesExplicitly.html">
http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/KSP/KSPComputeEigenvaluesExplicitly.html</a><br>&gt;&gt;<br>&gt;&gt; with and without a preconditioner. We have not coded the SVD<br>&gt;&gt; counterpart, but you can use
<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; <a href="http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/KSP/KSPComputeExplicitOperator.html">http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/KSP/KSPComputeExplicitOperator.html
</a><br>&gt;&gt;<br>&gt;&gt; and then call the LAPACK yourself.<br>&gt;&gt;<br>&gt;&gt;&nbsp;&nbsp; Matt<br>&gt;&gt;<br>&gt;&gt; On Jan 13, 2008 11:23 PM, Yujie &lt;<a href="mailto:recrusader@gmail.com">recrusader@gmail.com</a>&gt; wrote:
<br>&gt;&gt;&gt; Hi, everyone<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; I want to select iterative methods by observing the singular values<br>&gt;&gt;&gt; decompostion of the matrix. However, I don&#39;t know how to get all the<br>
&gt;&gt;&gt; singular values of the matrix in PETSc. I know the command<br>&gt;&gt;&gt; &quot;-ksp_monitor_singular_value&quot; may get the max and min singular values at<br>&gt;&gt;&gt; each iteration. How to get the singular values of the matrix I want to
<br>&gt;&gt;&gt; solve? In addition, when I use the preconditioned iterative method, how<br>&gt;&gt; to<br>&gt;&gt;&gt; get the singular values of the preconditioned iterative operator?<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; thanks a lot.
<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Regards,<br>&gt;&gt;&gt; Yujie<br>&gt;&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; --<br>&gt;&gt; What most experimenters take for granted before they begin their<br>&gt;&gt; experiments is infinitely more interesting than any results to which
<br>&gt;&gt; their experiments lead.<br>&gt;&gt; -- Norbert Wiener<br>&gt;&gt;<br>&gt;&gt;<br>&gt;<br><br></blockquote></div><br>