housh.m and krylov.m documentations

John W. Eaton jwe at bevo.che.wisc.edu
Tue Dec 4 11:17:32 CST 2007


On 28-Nov-2007, Marco Caliari wrote:

| On Tue, 27 Nov 2007, John W. Eaton wrote:
| 
| > On 27-Nov-2007, Marco Caliari wrote:
| > 
| > | In krylov.m, there is a
| > | 
| > | If B is a vector...
| > | 
| > | but B is not an input argument.
| > 
| > I'm not sure what the correct fix is for this.
| 
| Not sure how to write:
| 
| If V is a vector, then H contains the Hessenberg matrix such that
| A*U == U*H+RK*EK', where 
| 
| RK = A*U(:,K)-U*H(:,K) and EK' is the vector of length K [0 0 ... 0 1].
| Otherwise, H is meaningless.
| 
| ...
| 
| If V is a vector and K is greater than length(A)-1, then H contains the
| Hessenberg matrix such that A*U == U*H.

I made these changes.

Thanks,

jwe


More information about the Bug-octave mailing list