Add license header (overloading module)
[scilab.git] / scilab / modules / overloading / macros / %sp_inv.sci
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) INRIA
3 // 
4 // This file must be used under the terms of the CeCILL.
5 // This source file is licensed as described in the file COPYING, which
6 // you should have received as part of this distribution.  The terms
7 // are also available at    
8 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
9
10 function x=%sp_inv(a)
11
12 [ma,na]=size(a)
13 if ma<>na then error(20,1),end
14 [hand,rk]=lufact(a)
15 if rk<na then ludel(hand);error(19),end
16 x=[]
17 for k=1:ma
18   b=0*ones(ma,1);b(k)=1;
19   x=[x,sparse(lusolve(hand,b))]
20 end
21 ludel(hand);
22 endfunction