How to code block lanczos algorithm?
7 views (last 30 days)
Show older comments
Thanks ....
Answers (1)
Steven Lord
on 29 Jun 2021
function [T]=bsymLanczos2(A,v,m)
m=13;
n=100;
r = -5 + (5+5).*rand(n);
A = triu(r) + triu(r,1)';
v=rand(n,2); % 2*2 block vector
You ask the user to pass A, v, and m into your function then promptly throw the user-specified values in the garbage and use other values. If you're creating these to test that your function works correctly move them out of the bsymLanczos2 function and pass the variables created by running these lines into bsymLanczos2 as input arguments.
I'm also assuming there are some lines of code beyond what you've posted, as the code as written will only return an all zero matrix. All the code that calculates the alpha, delta, and V arrays has no impact on what bsymLanczos2 returns.
Since I'm guessing this is homework, to check that your code works as you expect see if your textbook has any examples worked out step by step. If it does run your code on the data from one or more of those examples using the debugging tools in MATLAB to step through the code line by line. Check that at each step you get the same answer as in the textbook example and if you don't that you understand why not.
See Also
Categories
Find more on Sources in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!