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

UptoLike

ВУЗ: 

ВГУ | Воронеж

Составители: 

Рубрика: 

- 21 -
a) определяем с помощью
IB
участок массива
JB
, описывающий
k
-
ю строку массива
B
;
b) просматриваем выделенный участок
JB
и для каждого просмат -
риваемого элемента
j
, стоящего в
n
-й позиции массива
JB
, полагаем
(
)
(
)
(
)
(
)
nBNmANjXjX
⋅
.
0
6
. Выбираем из РВН значения
i
-й строки матрицы
С
и присваиваем
их соответствующим компонентам
CN
.
0
7
. Если просмотр
i
-й строки завершен, то полагаем
1
i
i
и
возвращаемся на
0
2
.
Алгоритм
1
i
)
1
(
)
(
i
IC
i
IC
нет да
1
i
i
для 1)1(),( −+= iICiCIj :
0
))
(
(
j
JC
X
)
1
(
)
(
i
IA
i
IA
нет да
1
i
i
для
1)1(),( −+= iIAiAIm
:
)
(
m
JA
k
)
(
m
AN
A
для
1)1(),( −+= kIBkBIn
:
)
(
n
JB
j
)
(
)
(
)
(
n
BN
A
j
X
j
X
⋅
для
1)1(),( −+= iICiCIj
:
))
(
(
)
(
j
JC
X
j
CN
1
i
i
Задача 24. Найти произведение двух РМ
A
и
B
.
,
7531:
514253:
543211:
IA
JA
AN
                                      - 21 -

    a) определяем с помощью IB участок массива JB , описывающий          k-
       ю строку массива B ;
    b) просматриваем выделенный участок JB и для каждого просмат-
       риваемого элемента j , стоящего в n -й позиции массива JB , полагаем
       X ( j ) =X ( j ) +AN (m )⋅ BN (n ).
  6 0 . Выбираем из РВН значения i -й строки матрицы С и присваиваем
  их соответствующим компонентам CN .
  70 . Если просмотр i -й строки завершен, то полагаем i =i +1 и
                     0
  возвращаемся на 2 .



                                     Алгоритм
                                       i =1
                                IC (i ) =IC (i +1)
                                       нет              да
                                               i =i +1
                   для j =I C (i ), IC (i +1) −1: X ( JC ( j )) =0
                       IA(i ) =IA(i +1)
                              нет              да
                                             i =i +1
           для m =I A(i ), IA(i +1) −1 : k =JA(m)
                                       A =AN (m)
                       для n =I B ( k ), IB (k +1) −1: j =JB(n)
                                           X ( j ) = X ( j ) +A ⋅ BN (n)
                    для     j =I C (i ), IC (i +1) −1: CN ( j ) =X ( JC ( j ))
                                                         i =i +1
Задача 24. Найти произведение двух РМ A и B .
� AN : 1 1 2 3 4 5
 �
   � JA : 3 5 2 4 1 5 ,
� IA : 1 3 5 7
 �

Страницы