c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab c Copyright (C) INRIA c c This file must be used under the terms of the CeCILL. c This source file is licensed as described in the file COPYING, which c you should have received as part of this distribution. The terms c are also available at c http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt subroutine wspmsp(p,q,r,ar,ai,nela,inda,br,bi,nelb,indb, $ cr,ci,nelc,indc, $ ib,ic,xr,xi,xb,ita,itb,ierr) c multiply complex sparse matrices by the method of gustafson,acm t.o.m.s. c vol 4 (1978) p250. c*** input c p number of rows in a, c q number of columns in a c r number of columns in b c a a one-dimensional array containing the non-zero elements c of the first matrix,arranged row-wise, but not c necessarily in order within rows. c nela number of non-zero elements in a c inda(i) 1<=i<=p number of non-zero elements in row i of a. c inda(p+i) 1<=i