* Bug #12121 fixed - inv function did not work for complex sparse matrices.
[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 // Copyright (C) 2013 - Scilab Enterprises - Adeline CARNIS
4 //
5 // This file must be used under the terms of the CeCILL.
6 // This source file is licensed as described in the file COPYING, which
7 // you should have received as part of this distribution.  The terms
8 // are also available at
9 // http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
10
11 function x = %sp_inv(a)
12     
13     [ma, na] = size(a)
14     if ma<>na then 
15         error(20,1)
16     end
17     warning("off")
18     Lup = umf_lufact(a);
19     warning("on")
20     x = sparse(umf_lusolve(Lup, eye(ma, na)));
21     umf_ludel(Lup)
22
23 endfunction