Hello i have a conceptual idea of a sparse matvec implementation where i have multiple matrices, how would i go about implementing something such as<div><br></div><div>for i = ...</div><div>   for k = ...</div><div>      w_{ik} = K_i * u_k</div>
<div>   end</div><div>end</div><div><br></div><div>Where each of the K_i are sparse matrices... the K_i are various stiffness matrices whose size can range (although they are all the same size).  The u_k are reused</div><div>
<br></div><div>Thanks!</div><div><br></div><div>Mark</div>