Shouldn't the NR::sprsm function prototype look like:

void NR::sprspm(... Vec_O_INT &ijc) and not void NR::sprspm(... Vec_I_INT &ijc)?

This is how NR:sprstm is prototyped.

I noticed this first in the book and checked the distribution files also.

Pls advise

regards

Jeff

Hi Jeff,

No, the code is correct. Look at the main comment for the routine in the book: only those components of the matrix product that are pre-specified in the input array ijc are returned.

Best wishes,
Saul Teukolsky