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