1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84
|
function testh
%TESTH test Householder reflections
%
% Example:
% testh
% See also: testall
% CXSparse, Copyright (c) 2006-2022, Timothy A. Davis. All Rights Reserved.
% SPDX-License-Identifier: LGPL-2.1+
format long e
fprintf ('-------------------------------------------------\n') ;
x = [-3 4 5]' ;
disp (x) ;
[v, beta, s] = house (x) ;
disp ('v = ') ;
disp (v) ;
disp ('beta = ') ;
disp (beta) ;
disp ('s = ') ;
disp (s) ;
x = x - v*(beta*(v'*x)) ;
disp (x) ;
fprintf ('-------------------------------------------------\n') ;
x = [3 4 5]' ;
disp (x) ;
[v, beta, s] = house (x) ;
disp ('v = ') ;
disp (v) ;
disp ('beta = ') ;
disp (beta) ;
disp ('s = ') ;
disp (s) ;
x = x - v*(beta*(v'*x)) ;
disp (x) ;
fprintf ('-------------------------------------------------\n') ;
x = [1 eps]' ;
disp (x) ;
[v, beta, s] = house (x) ;
disp ('v = ') ;
disp (v) ;
disp ('beta = ') ;
disp (beta) ;
disp ('s = ') ;
disp (s) ;
x = x - v*(beta*(v'*x)) ;
disp (x) ;
fprintf ('-------------------------------------------------\n') ;
x = pi ;
disp (x) ;
[v, beta, s] = house (x) ;
disp ('v = ') ;
disp (v) ;
disp ('beta = ') ;
disp (beta) ;
disp ('s = ') ;
disp (s) ;
x = x - v*(beta*(v'*x)) ;
disp (x) ;
fprintf ('-------------------------------------------------\n') ;
x = -pi ;
disp (x) ;
[v, beta, s] = house (x) ;
disp ('v = ') ;
disp (v) ;
disp ('beta = ') ;
disp (beta) ;
disp ('s = ') ;
disp (s) ;
x = x - v*(beta*(v'*x)) ;
disp (x) ;
fprintf ('-------------------------------------------------\n') ;
x = [1 0 0]' ;
disp (x) ;
[v, beta, s] = house (x) ; %#ok
x = x - v*(beta*(v'*x)) ;
disp (x) ;
fprintf ('-------------------------------------------------\n') ;
|