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
|
// Copyright Bruno Pinçon, ESIAL-IECN, Inria CORIDA project
// <bruno.pincon@iecn.u-nancy.fr>
//
// This set of scilab 's macros provide a few sparse utilities.
//
// This software is governed by the CeCILL license under French law and
// abiding by the rules of distribution of free software. You can use,
// modify and/ or redistribute the software under the terms of the CeCILL
// license as circulated by CEA, CNRS and INRIA at the following URL
// "http://www.cecill.info".
//
// As a counterpart to the access to the source code and rights to copy,
// modify and redistribute granted by the license, users are provided only
// with a limited warranty and the software's author, the holder of the
// economic rights, and the successive licensors have only limited
// liability.
//
// In this respect, the user's attention is drawn to the risks associated
// with loading, using, modifying and/or developing or reproducing the
// software by the user in light of its specific status of free software,
// that may mean that it is complicated to manipulate, and that also
// therefore means that it is reserved for developers and experienced
// professionals having in-depth computer knowledge. Users are therefore
// encouraged to load and test the software's suitability as regards their
// requirements in conditions enabling the security of their systems and/or
// data to be ensured and, more generally, to use and operate it in the
// same conditions as regards security.
//
// The fact that you are presently reading this means that you have had
// knowledge of the CeCILL license and that you accept its terms.
function [X] = test_on_columns(X,Xold)
//
// X and Xold are (n,t) matrix , Xold being an optional argument
// verify and force all columns of the new X to be non paralleles
// and also to be non paralleles to the columns of Xold if Xold
// is given
// 2 columns j1 and j2 are // if | Cj1'*Cj2 | = n because
// all elements are in {+1,-1}
//
[lhs,rhs] = argn()
[n,t] = size(X)
if rhs == 1 then
for j=2:t
while %t
res_test = abs(X(:,j)'*X(:,1:j-1))
if or(res_test == n) then
X(:,j) = fsign(rand(n,1)-0.5)
else
break
end
end
end
else // rhs = 2
for j=1:t
while %t
res_test = abs([X(:,j)'*Xold X(:,j)'*X(:,1:j-1)])
if or(res_test == n) then
X(:,j) = fsign(rand(n,1)-0.5)
else
break
end
end
end
end
endfunction
|
