Matrix-Vektor-Produkt bei gegebenem CRS-Format

Bachelor Algorithmik kontinuierlicher Systeme
1

Disclaimer: Dieser Thread wurde aus dem alten Forum importiert. Daher werden eventuell nicht alle Formatierungen richtig angezeigt. Der ursprüngliche Thread beginnt im zweiten Post dieses Threads.

2

Matrix-Vektor-Produkt bei gegebenem CRS-Format
Hallo,
weiß jemand wie man das Matrix-Vektor-Produkt bei gegebenem CRS-Format berechnet, ohne die Matrix zu rekonstruieren?

Grüße,
Leo

3

Matrix-Vektor-Produkt bei gegebenem CRS-Format


def f(val,row,col,v)
  x = zeros(v.shape)
    for i in 0..len(row)-1:
       for j in row[i]...row[i+1]:
         x[i] += val[j] * v[col[j]]

  return x