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