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 // Copyright (C) 2018 - Samuel GOUGEON
   5 //
   6 // Copyright (C) 2012 - 2016 - Scilab Enterprises
   7 //
   8 // This file is hereby licensed under the terms of the GNU GPL v2.0,
   9 // pursuant to article 5.3.4 of the CeCILL v.2.1.
  10 // This file was originally licensed under the terms of the CeCILL v2.1,
  11 // and continues to be available under such terms.
  12 // For more information, see the COPYING file which you should have received
  13 // along with this program.
  14
  15 function x = %sp_inv(a)
  16     [ma, na] = size(a)
  17     if ma<>na then
  18         msg = _("%s: Argument #%d: Square matrix expected.\n")
  19         error(msprintf(msg, "%sp_inv", 1))
  20     end
  21
  22     oldw = warning("query");
  23     warning("off")
  24     Lup = umf_lufact(a);
  25     warning(oldw)
  26
  27     x = sparse(umf_lusolve(Lup, eye(ma, na)));
  28     umf_ludel(Lup)
  29 endfunction