Thanks Thomas, I needed that. &nbsp;<div><br></div><div>I had used trevc as you had, but I hadn&#39;t done the multiplication with the original matrix and the eigenvalue and see how the results compare. &nbsp;I was just comparing my eigenvectors with eigenvectors calculated with Mathematica and another application. &nbsp;I too get a bunch of zeros indicating the eigenvectors are good, but my eigenvectors don&#39;t match with Mathematica. &nbsp;I have tried different methods of normalizing the eigenvectors to see if that is the difference, but have had no such luck in getting them to match up.</div>
<div><br></div><div>Jeremy<br><br><div class="gmail_quote">On Sat, Jun 21, 2008 at 6:44 PM, Thomas Klimpel &lt;<a href="mailto:Thomas.Klimpel@synopsys.com">Thomas.Klimpel@synopsys.com</a>&gt; wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hello Jeremy,<br>
<br>
I modified your file hseqr.cpp slightly (just update libs/numeric/bindings), to create an example of trecv in use. Hope it helps.<br>
<br>
Regards,<br>
Thomas<br>
<br>
Index: bindings/lapack/test/hseqr.cpp<br>
===================================================================<br>
--- bindings/lapack/test/hseqr.cpp &nbsp; &nbsp; &nbsp;(Revision 46588)<br>
+++ bindings/lapack/test/hseqr.cpp &nbsp; &nbsp; &nbsp;(Arbeitskopie)<br>
@@ -90,6 +90,30 @@<br>
 &nbsp; &nbsp; cout &lt;&lt; &quot;&#39;New&#39; original matrix:\n&quot; &lt;&lt; origG &lt;&lt; endl;<br>
 &nbsp; &nbsp; cout &lt;&lt; &quot;==================================&quot; &lt;&lt; endl;<br>
<br>
+ &nbsp; &nbsp;ublas::matrix&lt;complex&lt;double&gt;, ublas::column_major&gt; cVL(cZ);<br>
+ &nbsp; &nbsp;ublas::matrix&lt;complex&lt;double&gt;, ublas::column_major&gt; cVR(cZ);<br>
+ &nbsp; &nbsp;boost::numeric::bindings::traits::detail::array&lt;complex&lt;double&gt; &gt; work(3*n);<br>
+ &nbsp; &nbsp;lapack::trevc(&#39;B&#39;,&#39;B&#39;,G,cVL,cVR,work);<br>
+<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;\n==================================&quot; &lt;&lt; endl;<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;Testing left &amp; right eigenvectors...&quot; &lt;&lt; endl;<br>
+ &nbsp; &nbsp;Hessenberg(G);<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;Many &#39;zeros&#39;:&quot; &lt;&lt; endl;<br>
+ &nbsp; &nbsp;for ( int i=0; i&lt;n; ++i){<br>
+ &nbsp; &nbsp; &nbsp; &nbsp;cout &lt;&lt; (ublas::prod( G, column(cVR, i) ) - values(i) * column(cVR, i)) &lt;&lt; endl;<br>
+ &nbsp; &nbsp; &nbsp; &nbsp;cout &lt;&lt; (ublas::prod( ublas::herm(G), column(cVL, i) ) - conj( values(i) ) * column(cVL, i)) &lt;&lt; endl;<br>
+ &nbsp; &nbsp;}<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;==================================&quot; &lt;&lt; endl;<br>
+<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;\n==================================&quot; &lt;&lt; endl;<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;Verifying diagonal matrix...&quot; &lt;&lt; endl;<br>
+ &nbsp; &nbsp;ublas::matrix&lt;complex&lt;double&gt;, ublas::column_major&gt; cG(n,n,0);<br>
+ &nbsp; &nbsp;Hessenberg(G);<br>
+ &nbsp; &nbsp;G = ublas::prod( G, cVR );<br>
+ &nbsp; &nbsp;G = ublas::prod( ublas::herm(cVL), G );<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;&#39;diagonal&#39; matrix:\n&quot; &lt;&lt; G &lt;&lt; endl;<br>
+ &nbsp; &nbsp;cout &lt;&lt; &quot;==================================&quot; &lt;&lt; endl;<br>
+<br>
&nbsp;}<br>
<br>
<br>
<br>
<br>
-----Ursprüngliche Nachricht-----<br>
Von: <a href="mailto:ublas-bounces@lists.boost.org">ublas-bounces@lists.boost.org</a> im Auftrag von Jeremy Conlin<br>
Gesendet: Do 19.06.2008 21:41<br>
An: ublas mailing list<br>
Betreff: [ublas] Example of using trevc.hpp?<br>
<div><div></div><div class="Wj3C7c"><br>
I am trying to use the trevc.hpp header with my new &nbsp;hseqr.hpp header. &nbsp;I&#39;m<br>
not getting the results I expected and I&#39;m trying to track down the error.<br>
&nbsp;Everything compiles and runs, but the result is not what I expected. &nbsp;I can<br>
confirm the results of hseqr by a call to Mathematica, but I can&#39;t figure<br>
out an equivalent command for trevc.<br>
Does anyone have an example of trevc in use?<br>
<br>
Thanks,<br>
Jeremy<br>
<br>
</div></div><br>_______________________________________________<br>
ublas mailing list<br>
<a href="mailto:ublas@lists.boost.org">ublas@lists.boost.org</a><br>
<a href="http://lists.boost.org/mailman/listinfo.cgi/ublas" target="_blank">http://lists.boost.org/mailman/listinfo.cgi/ublas</a><br>
<br></blockquote></div><br></div>