C Output from Public domain Ratfor, version 1.0
      subroutine spwy(n,m,k,nnza,a,ja,ia,ao,jao,iao,nnzdmax,d,jd,id, dsu
     *b,jdsub,nnzemax,e,je,ie,nsubmax,lindx,xlindx, nnzlmax,lnz,xlnz,iw,
     *iwmax,iwork,xsuper,tmpmax, tmpvec, wwm,wwn,cachsz,level,x,s,u,c,y,
     *b,small, ierr,maxit,time)
      integer nnza,m,n,nnzdmax,nnzemax,iwmax, nnzlmax,nsubmax,cachsz,lev
     *el,tmpmax,ierr,maxit,time(7)
      ja(nnza),jao(nnza),jdsub(nnzemax+1),jd(nnzdmax), ia(n+1),iao(m+1),
     *id(m+1),lindx(nsubmax),xlindx(m+1), iw(m,5),xlnz(m+1),iwork(iwmax)
     *,xsuper(m+1),je(nnzemax),ie(m+1)
      double precision small, a(nnza),ao(nnza),dsub(nnzemax+1),d(nnzdmax
     *), lnz(nnzlmax),c(n),y(m),wwm(m,3),tmpvec(tmpmax), wwn(n,14),x(n),
     *s(n),u(n),e(nnzemax),b(m,k)
      do23000 i=1,k
      call dcopy(m,b(1,i),1,a(nnza - m),1)
      call slpfn(n,m,nnza,a,ja,ia,ao,jao,iao,nnzdmax,d,jd,id, dsub,jdsub
     *,nsubmax,lindx,xlindx,nnzlmax,lnz, xlnz,iw(1,1),iw(1,2),iwmax,iwor
     *k,iw(1,3),iw(1,4), xsuper,iw(1,5),tmpmax,tmpvec,wwm(1,2),cachsz, l
     *evel,x,s,u,c,y,b,wwn(1,1),wwn(1,2),wwn(1,3), wwn(1,4),nnzemax,e,je
     *,ie,wwm(1,3),wwn(1,5),wwn(1,6), wwn(1,7),wwn(1,8),wwn(1,9),wwn(1,1
     *0),wwn(1,11), wwn(1,12),wwn(1,13),wwn(1,14),wwm(1,1),small,ierr,ma
     *xit, time)
      call dcopy(m,y,1,b(1,i),1)
23000 continue
23001 continue
      return
      end