* Bug #13810 fixed - householder(v, k*v) returned column of %nan. Input parameters...
[scilab.git] / scilab / modules / linear_algebra / tests / nonreg_tests / bug_13810.tst
1 // ============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2015 - Samuel GOUGEON
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // ============================================================================
7 // <-- Non-regression test for bug 13810 -->
8 //
9 // <-- Bugzilla URL -->
10 //     http://bugzilla.scilab.org/13810
11 //
12 // <-- Short Description -->
13 // householder(v, k*v) returned [%nan %nan %nan]'
14 assert_checkequal(householder([%e ; 0 ; 0]), [0 ; 1 ; 0]);
15 v = [ 2 3 1]';
16 ref = [ 0.8320503  -0.5547002  0 ]';
17 assert_checkalmostequal(householder(v,v), ref, 1e-6);
18 assert_checkalmostequal(householder(v,2*v), ref, 1e-6);
19 ref = [  0.5345225  0.8017837  0.2672612 ]';
20 assert_checkalmostequal(householder(v,-2*v), ref, 1e-6);
21 v = [ 2 -%i 1+%i].';
22 ref = [0.7559289 -0.3779645*%i  0.3779645 + 0.3779645*%i].';
23 assert_checkalmostequal(householder(v,v), ref, 1e-6);