“…For example, the algebraic multigrid (AMG) solver [13], an iterative algorithm widely used in both laser fusion and climate modeling, reports above 90% SpMV operation of its overall iterations. Since 1970s, plenty of researches have been dedicated to optimizing SpMV performance for its fundamental importance, which are generally separated into two paths, one for developing new application-specific storage formats [7,21,[29][30][31][32]36], and the other for tuning on emerging processor architectures [10,11,16,22,27,28,34,35]. This separation leads to low performance and low productivity in SpMV solvers and applications:…”