krylov.m File Reference

Functions
	if (isa(A,"single")||isa(V,"single")) defeps
endif	if (nargin< 3||nargin > 5) print_usage()
	elseif (nargin< 5) pflg=0
endif	if (nargin< 4) eps1
endif	if (isempty(eps1)) eps1
	if (!na) error("A(%d x %d) must be square"
	rows (A)
	columns (A))
	if (m!=na) error("A(%d x %d)
endif	if (columns(V) > 0) Q
	Q (idx,:)
	H (i, nu-columns(V)+(1:columns(V)))
	U (:, i)
	U (idx(1), i)
	U (idx, i:j1)
Variables
	function [Uret, H, nu] = krylov (A, V, k, eps1, pflg)
endif	na = issquare (A)
	endif [m, kb] = size (V)
V(%d x%d) submit a bug	report
for	kk
endfor for	ii
	hv = U(idx,ii)
	av = alpha(ii)
endfor endif	V = A*Q
	nu = length (alpha)
for	i
endfor endwhile	j1 = columns (U)
	idx = pivot_vec(i:na)
	Uret = U

---

Function Documentation

columns

(

A

)

elseif

(

)

[pure virtual]

H	(	i	,
		nu-	columnsV)+(1:columns(V)
	)

endif if	(	columns(V)	,
		0
	)

if

(

m!

= na

)

if

(

!

na

)

endif if

(

isempty(eps1)

)

endif if

(

)

endif if

(

nargin< 3||nargin >

5

)

if

(

isa(A,"single")||isa(V,"single")

)

Q	(	idx	,
		:
	)

rows

(

A

)

U	(	idx	,
		i:j1
	)

U	(	idx(1)	,
		i
	)

U	(	:	,
		i
	)

---

Variable Documentation

av = alpha(ii)

endif[m, kb] = size (V)

function[Uret, H, nu] = krylov (A, V, k, eps1, pflg)

hv = U(idx,ii)

for i

Initial value:

 1:nu
      hv = U(:,i)

idx = pivot_vec(i:na)

endfor for ii

Initial value:

 nu:-1:1
        idx = pivot_vec(ii:na)

endfor endwhile j1 = columns (U)

for kk

Initial value:

 1:columns (Q)
        Q(pivot_vec(nu-columns(Q)+kk),kk) = 1

return endif na = issquare (A)

endfor nu = length (alpha)

V (%d x %d) submit a bug report

Uret = U

V = A*Q