linii_1.m

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

function [vl,vic]=linii_1(x,n,nalloc,val,col,link)
  k=0;
  j=1;
  for i=1:nalloc;
    if ((link(i) == -1) &&( link(i-1) != -1));
      k(j)=link(i-1);
      j++;
    endif;
  endfor;
  k(1:j-1); %ultimele elemente nenule de pe fiecare linie
%--------------------------------------------------
  j=1; i=1;
  while (i<=n)
    if(i == 1 && x==1)
	vl=val(i:k(i)); %elemem nenule pt linia 1
	vic=col(i:k(i)); %indicii elem

    endif

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

    endif
    i++;
  endwhile;


endfunction