linii_2.m

来自「Format Linked List for sparse Matrix」· M 代码 · 共 29 行

function [vl,vic]=linii_2(x,n,nnz,va,ja,ia)
  k=0;
  j=1;
  for i=2:n+1;
    %if ((link(i) == -1) &&( link(i-1) != -1));
      
      k(j)=ia(i)-1;
      j++;
  endfor;
  k(1:n); %indicii ultimelor elemente nenule de pe fiecare linie
%--------------------------------------------------
  j=1; i=1;
  while (i<=n)
    if(i == 1 && x==1)
	vl=va(i:k(i)); %elemem nenule pt linia 1
	vic=ja(i:k(i)); %indicii elem

    endif

    if(i == x && x != 1);
	vl=va(k(i-1)+1:k(i)); %elemen nenule pt linia x
	vic=ja(k(i-1)+1:k(i)); %indicii elemen nenule

    endif
    i++;
  endwhile;


endfunction