14
edits
User:Bumi: Difference between revisions
Jump to navigation
Jump to search
m
→Computing left eigenvectors as a third output
| (4 intermediate revisions by the same user not shown) | |||
| Line 33: | Line 33: | ||
* Please state the commits and patches you already contributed to Octave. | * Please state the commits and patches you already contributed to Octave. | ||
** [https://savannah.gnu.org/bugs/?36437 bug #36437] | ** [https://savannah.gnu.org/bugs/?36437 bug #36437] | ||
** '''Currently I am working on [http://wiki.octave.org/User:Bumi#Preliminary_balancing | ** '''Currently I am working on [http://wiki.octave.org/User:Bumi#Preliminary_balancing Preliminary balancing]''' | ||
*** '''[https://savannah.gnu.org/patch/?8960 patch #8960]''' | |||
*** [http://pastebin.com/CBBUYYDW EIG.cc] | |||
*** [http://pastebin.com/9met8n18 EIG.h] | |||
*** [http://pastebin.com/9PZLTd48 fEIG.cc] | |||
*** [http://pastebin.com/Z6RifG2i fEIG.h] | |||
*** [http://pastebin.com/RJ6Fi04f eig.cc] | |||
*** And I also created a diff file | |||
**** [http://pastebin.com/C7PhbzVS diff] | |||
== F: Feeling fine == | == F: Feeling fine == | ||
| Line 112: | Line 120: | ||
=== Timeline === | === Timeline === | ||
* '''Community Bonding period''' (Until May 22) | * '''Community Bonding period''' (Until May 22) | ||
** Get acquainted with the code and LAPACK | ** Get acquainted with the code and LAPACK, <u>finish most of preliminary balancing</u> | ||
* '''Week 1-2''' (May 23 - Jun 5) | * '''Week 1-2''' (May 23 - Jun 5) | ||
** Finals, non-coding time | ** Finals, non-coding time | ||
| Line 120: | Line 128: | ||
* '''Midterm evaluations''' (Jun 20 - Jun 27) | * '''Midterm evaluations''' (Jun 20 - Jun 27) | ||
* '''Week 6''' (Jun 27 - Jul 3) | * '''Week 6''' (Jun 27 - Jul 3) | ||
** Implementing preliminary balancing, testing | ** <s>Implementing preliminary balancing, testing</s> | ||
* '''Week 7-10''' (Jul 4 - Jul 31) | * '''Week 7-10''' (Jul 4 - Jul 31) | ||
** algorithm choosing for eigenvalue calculation (chol or qz) | ** algorithm choosing for eigenvalue calculation (chol or qz) | ||
| Line 158: | Line 166: | ||
* The Matlab documentation does not mention whether there is balancing in the generalised case, but if needed the *ggevx could be used same. | * The Matlab documentation does not mention whether there is balancing in the generalised case, but if needed the *ggevx could be used same. | ||
** '''EDIT:''' In Matlab there is no balance option in the generalised case: | |||
* IN: | |||
{{Code|Matlab's eig|<syntaxhighlight lang="matlab"> | |||
A = [1, 2 ; 3, 4] | |||
B = [5, 6 ; 7, 8] | |||
e = eig(A, B, 'nobalance') | |||
</syntaxhighlight>}} | |||
* OUT: | |||
Error using eig | |||
For generalized eigenproblem EIG(A,B), flag must be 'vector', 'matrix', 'qz', or 'chol'. | |||
==== Computing left eigenvectors as a third output ==== | ==== Computing left eigenvectors as a third output ==== | ||
* *geevx and *ggevx could also be used as these can compute not just right but left | * *geevx and *ggevx could also be used as these can compute not just right but left eigenvectors also. | ||
==== Choosing among generalized eigenvalue algorithms ==== | ==== Choosing among generalized eigenvalue algorithms ==== | ||