1// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab 2// Copyright (C) INRIA 3// 4// Copyright (C) 2012 - 2016 - Scilab Enterprises 5// 6// This file is hereby licensed under the terms of the GNU GPL v2.0, 7// pursuant to article 5.3.4 of the CeCILL v.2.1. 8// This file was originally licensed under the terms of the CeCILL v2.1, 9// and continues to be available under such terms. 10// For more information, see the COPYING file which you should have received 11// along with this program. 12 13function r=%sp_k_sp(a,b) 14 // a.*.b with a and b sparse 15 16 [ija,va,mna]=spget(a) 17 [ijb,vb,mnb]=spget(b) 18 ia=ija(:,1);ja=ija(:,2) 19 ib=ijb(:,1);jb=ijb(:,2) 20 21 ij=[((ia-ones(ia))*mnb(1)).*.ones(ib)+ones(ia).*.ib,.. 22 ((ja-ones(ja))*mnb(2)).*.ones(jb)+ones(ia).*.jb] 23 r=sparse(ij,va.*.vb,mna.*mnb) 24endfunction 25