1. ´Ł - ¨ 14-02-00326 Ł4ª47 ß Æ ŁŁ ae ae ª Œ ı 11.38.241.2015. ´ß ŁaeºŁ º ßØ ˇÆˆ .
ŁŁ ae º
Ø ŒaeŁ ae Œ ı ae Ł º ß ŁºŁ
ae Œ Ł
Ł º ae Ł ª º ƺ Œ ae ŁØ Ø Œ ae Ł,
Ł ae Œ aeŁae
Ł Ø. ´ Œ Ł ª Ł Ø Œ Ł ŁŒº Ø ŁŒŁ Ł ß ŁaeºŁ º ºŁ Ø ßı ºª Æ Ł aeŒŁı
ae ae Øł ª Ł Ø Ł º ª Ł aeº ƺ ŒŁ
Ł Ł ª º Ø, . ´ ¸ Æ ß Æº ŒŁ Ł Ł
. ºª Ł º ª ŒŁ [1]. ˜ ƺ Ł . ˝ ºº ºŁ ŁŁ ae Œ Ø Ø Œ Ł ßØ Ł ƺ
ƺ ae Œ
ßı ae ae ºº ºŁ ŁŁ º Ł Ł Æ ø Ł ƺ Œ . ł Ł ı Ł ª º ßı ¸
Œ . ˚ ß ˇ.. `º, ¯.—. ˝ª, .¸. Œº Œ-ˇÆªaeŒŁØ ˆaeaeßØ ŁaeŁ ˇ ae º
ºº º aeŒŁı
Ł Ø ŁØ Ł ŁŁ. ˇ Œ º ºº º ßı aeae ß ß ß ŁaeºŁ º aeŒ Ł . ´ß ŁaeºŁ º ß ß ŁaeºŁ º ª aeŒ Ł .
ae º Ł Ø Œ ł Ł ƺ - ı Ł ª ŁØ ( ¸ ). ´ß ŁaeºŁ º ª ŒŁ ºŁ Ł aeŒŁ Ł º ßı Ł aeº , aeŒ ºº º ŒaeŁ º . ˜º Ø aeŁae Œae Ł
ª ß, Ł Œ ß Ł ŁŁ Æ ae º Ø aeŁae aeŒ Ł Ł º ßı ae Łª
Łaeı ßı ae Łª Łaeº Ł øŁØ Œ Ł ß ºŁ Ø ßı ºª Æ Łł Ł Œ º ßı ae Ø Œ ŁaeŒŁ ae ß Ł Ł ƺ
ł Ł ßØ Ł
º ŁŁ Ł Ł Łae º
ae ae Πae
´ Ø ¸ . Œ Ł º Œ Ł ae ºÆ º Ł Ø ı Ł ª º Ø ae ı Ł aeŒ Ø Æ ae Ł Łaeı
ª Œ ŁŁ [2], ı [4]. ˜ ı Ł ª º º ßı ¸
Œ Ł ß
Œ ´ ª [3], Œ ß ß ÆߺŁ Æ Æø ß Ł Ł Ł [5]. ˇ Ł Æߺ Æ Æ º Ł ŁŁ ƺ ae Ł. ˇ Ł
ŁaeŁ ß Æº ŒŁ Ł ae ßı Æ ı, Ł [7],
ScaLAPACK [8].
ß Æ ÆøŁºŁ Œ Ł ŁŁ Ł ae Œ Æ Ł Łaeı Ø Æº ae ae º Ł Ø Ø Ø
Ø aeŁae ß. ´ º , Ł , Œº ø Ø Œ ß ae -ƺ ŒŁ Ł º Ł ae º º Æ ß ºŁae ae Œ ae Ł łŁ Œ Ł ßØ ae ł Ł ß º Łae º Ł ae º
Łaeº ßı , º Œ Ł ßØ ºª Ł , ae Œ ¸ , Œº øŁØ ƺ Œ [6]. ı ßı Ł Ø Æº ı Ł ª º Ł ß aeŁae Ł ª
ae ß º ªŁ aeŒŁ ºŁ, aeŒ ı Ł ª Ø ¸ ŒŁ ƺ Ł øŁ
º Œ
ßı Ł ae ˇÆˆ Œ ´ Æ
C :
A
S 1WRh ˙ ae S ae ae Ł Ł º ß Æ ß ŁaeŁ ßı ae -ƺ Œ Łae Œ 1. — ł Ł aeŁae Ł ª ß (2) ºŁSk, k = 1; : : : ; M . ˛ae º ae º Ł ºº º Ł ª ae
ae ø Ł ae
ß ŁaeºŁ º º ae Ł ae ª ae º Ł (3)
aeae ƺ Œ Sk º ae ß Łaeº ŁØ Æ ł ª . ´ Ł ŁØ Œ ª Œ h = (h1; : : : ; hM 1)T ,
ºº º º
ß ae ae -ƺ Œ Sk Ł º Ø 4ª48ŒŁ, ı º Łı Ł ŁŁ º . —Łae. 1. ˇ Ł Łaeı ¸ ¸ Ø ¸ ae ae Æ º Ł Ł Ø ae Ł Ø ae Ł ¸ ƺ Œ , Œ º Ł Ø º Ł
Ł ´ ae ´ß ŁaeºŁ º aeŒ Ł ae º Ł º Ł Ø ª ŒŁ ł Ł ¸ ª ŒŁ Œ ae ae
ß Łaeº Ł ae ae ae ae Łaeº . ß ae ae Ł ºŁ Œ ºŁ ae aeº º ßı º Ł Ø Œ Ł ŁŁ Ł, Œ Œ aeº ae Ł , Ł ß ŁŁ Ł ß ºŁae , ae Ł , Ł
ª Ł. ´ ae ae ŁŁ ae [11], ß ae Ł , Œ Œ º n n Æ ae n3 º Ł ºŁŒ Ł ßı Œ Ø Œ ºŁ ae ŁØ Æø Ł Ø
LAPACK. ª Œ ºŁ ae º Ł ºŁŒ Ł ºª Ł Ł Ø ª ŒŁ Ø Œ Ł ŁŁ Ł Œ Œ ł Ł ø Ł aeº º º Ł ºŁŒ Ł ŁºŁ ß ŁaeºŁ º Æ Ł º ßı ßı ł 7Nkn3 + 5Nkn2l + (3M 5)(n3 + n2l); 3. ´ßŁaeºŁº aeŒŁ
OAD = ((N
M + 1) =M ) (7n3 + 5n2l) + (3M 5)(n3 + n2l) ª ß ŁaeºŁ º aeŒ Ł , Œ Œ ł Ł S = OMS=OAD, º ªŒ ae : ŒŁ Ł ª Æ ºŁ Ł º ª Œ ºŁ ae , Ł Łae º ßÆŁ . ˛ Ł º
ƺ Œ N 4. ˇŒŁaeŒŁ ºß
ŁŁ Ł Œ Œ ºŁ Ł ª
ae ª ae ºŁ ¸ º Ł
, .Œ. ºº º ae 5. ˙ŒºŁ Œ ºŁº
, ßØ Ł . ºŁ Ł aeŒŁ º Ł Ø ª ŒŁ. — ºŁ ae Ł ß ŁaeºŁ º ß Œae Œ ŁŒ . ¨ae º Ł ae
ł Ł ¸ aeŒ º Œ ae Œ ß Łaeº Ł ŁaeŒ ßı ºŁ Ł ŁŒº
Ł º Ł aeŒŁØ . , ˝ŁŒ º ¯. .
ß ł Ł ae ßı ŁØ, .: ˝ 2. Hockney R.W. The potential calculation and some applications // Meth. Comput. Phys.
1970. Vol. 9. P. 136 211. 3. Wang H.H. A parallel method for tridiagonal equations // ACM Trans. Math. Software.
1981. Vol. 7. P. 170 183. 4. ´ º ı .´., ˙ º
ł Ł aeŁae ß ´ß ŁaeºŁ º ß 5. Johnsson S.L. Solving narrow banded systems on ensemble architectures // ACM Trans.
Math. Software. 1985. Vol. 11. P. 271 288. 6. Ortega J.M. Introduction to parallel and vector solution of linear systems,
New York: Springer, 1988. 7. ŒŁ ª 8. Cleary A., Dongarra J.J. Implementation in ScaLAPACK of divide-and-conquer algorithms for banded and tridiagonal linear systems // Report No. UT-CS-97-358 (LAWN No. 125), University of Tennessee, 1997.
aeŒŁØ . ., ˆ ºŁ .´. Łaeº ß ß, .: ˝ 12. Dongarra J.J., et al. LINPACK Users’ Guide, Philadelphia:SIAM, 1979. 13. Belov P.A., Yakovlev S.L. Asymptotic method for determining the amplitude for three-particle breakup: Neutron-deuteron scattering // Phys. Atom. Nucl. 2013. Vol. 76. P. 126 138. 452 Arrowhead decomposition for a block-tridiagonal system of linear equations Pavel Belov, Eduard Nugumanov and Sergey Yakovlev Keywords: arrowhead decomposition method, block-tridiagonal matrix, system of linear equations, parallel solution, computational speedup, matrix sweeping algorithm The arrowhead decomposition method which allows efficient parallel solution of a blocktridiagonal system of linear equations is presented. The computational speedup with respect to the matrix sweeping algorithm is analytically estimated by taking into account the number of elementary operations of multiplication for the parallel and serial parts of the decomposition method. It is shown that the maximal speedup is achieved for the finite number of parallel processors. For a given size of the initial system of linear equations, the parameters of the computational system which give the maximal speedup are obtained. Computational experiments confirm the analytical estimations of the computational speedup.