Методы решения систем с разреженными матрицами. Способы хранения и представления разреженных матриц, операции над ними. Блатов И.А - 18 стр.

                                        - 18 -

  k =1; С (1) =0 ;  I A(2) =I A(1)
  i =1,...,3 : i =1 :IX ( JA(1)) =IX (1) =1 ≠0
                    C (1) =0 +1 ⋅ BN (1) =1
             i =2 : IX ( JA(2)) =IX (4) =3 ≠0
                    C (1) =1 +3 ⋅ BN (3) =1 +3 ⋅ 4 =13
             i =3 : IX ( JA(3)) =IX (6) =0
  k =2 ; С (2) =0 ; I A(3) ≠I A(2)
  i =4,...,8 :
             i =4 : IX ( JA(4)) =IX (3) =0
             i =5 : IX ( JA(5)) =IX (2) =2
                    C (2) =0 +9 ⋅ BN (2) =9 ⋅ 2 =18
  . . . . . . . . . . . . . . . . . . . . . . . . . . .
      3.6.   Умножение РМ на заполненный вектор (результат –
             заполненный вектор)
      Этот алгоритм незначительно отличается от предыдущего. Разница в том,
что РМЦУ заполняется по k -й строке матрицы A , а не по вектору B .
Здесь, по-видимому, можно обойтись без РМЦУ следующим образом.
                                Алгоритм
1. k =1.
2. Просматриваем содержимое        k –й    строки в JA и умножаем
   соответствующие компоненты AN на компоненты, накапливая результат в
   соответствующей ячейке C (k ) вектора C .
                         Схема       решения
                                    k =1
                                 С (k ) =0
                              IA(k +1) =IA(k )
                        нет                        да
          i =I A(k ), IA(k +1) −1                k =k +1
      C (k ) =C (k ) +AN (i ) ⋅ B ( JA(i ))
        k =k +1
  Задача 23. Умножить РМ A из задачи 22 на заполненный вектор
  B =(2 4 6 8 10 9 )T .
  k =1; С (1) =0 ; IA(2) ≠IA(1)
  i =1,...,3 :
               i =1 : C (1) =0 +1×2 =2