R: explicitly call Matrix::solve in bonpow and alpha.centrality, otherwise the dense solvers are called. (Why???)