Otimizando a solução inversa para um sistema linear triangular inferior esparso
Eu tenho a representação da coluna esparsa compactada (csc) da matriz triangular inferior nxn A com zeros na diagonal principal e gostaria de resolver b em (A + I)' * x = b Esta é a rotina que tenho para calcular isso: void backsolve(const int*__restrict__ Lp, const int*__restrict__ Li, const...