mat_vec_multiply.m

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

function mat_vec_multiply(filename_1,vec)
[n,nnz,va,ja,ia]=read_2(filename_1);

%-----------------------------------------------------
%======================================
k=1;t=1;i=1;j=1;valM=0; va_out=[]; ja_out=[]; root_out=1;
for i=1:n 
  vec(i)=i;
endfor

for i=1:n

    [vl,vic]=linii_2(i,n,nnz,va,ja,ia); % vl(i)=elem nenule de pe linia i,respectiv vic indice elementelor(coloana)
  
    k=1; t=1;
      for k=1:length(vl)   %parcurge linie cu linie
	for t=1:length(vec) %parcurge vectorul vec
	  if(vic(k) == t) %daca exista elem pe linie corespunzator cu elem din vector vec => inmulteste
	    valM=valM+vl(k)*vec(t); %elementul i din vect rezultat 
	  endif
	endfor
      endfor
	  
	  va_out=[va_out,valM]; %vector ce contine elem nenule din matricea rezultata
	
	
	  valM=0;
  
 endfor
va_out'


%================================================================

endfunction