For the iterative methods, look at iml.ps.
For the basic matrix subroutines, look at sp1_5c (which includes spblas and mv)
For direct sparse solves, use SuperLU.