Ефективність елементарних неунітарних перетворень в поєднанні з віртуальними матрицями в організації обчислень матричної алгебри

Abstract

Сучасні прикладні задачі матричної алгебри потребують алгоритмів ефективної та надійної роботи з великою кількістю рівнянь. Ціль таких алгоритмів полягає в перетворенні цих рівнянь з широким діапазоном чисельних значень їх коефіцієнтів. Перетворення спрямовані на спрощення форми матриць з метою надання їм однієї з канонічних. Розроблено багато чисельних методів, схем, алгоритмів для досягнення цих цілей. На якість результатів перетворень безумовно впливає не лише сам метод а також і порядок виконання окремих операцій. В даній роботі запропоновано прийнятний ( на наш погляд ) варіант роботи з відомими методами та схему їх застосування коли пріоритетні для обчислень питання точності та стійкості, економії оперативної пам’яті та раціональної кількості операцій узгоджуються між собою. Базовою в цих перетвореннях прийнята елементарна неунітарна матриця. Для обмеження безконтрольного росту коефіцієнтів призначено певний порядок операцій з використанням Virtual матриць, в основі яких управління перестановками рядків / стовбців. Ланцюги перетворень з використанням елементарних матриць формують циклічні процеси та можуть бути універсальними. В роботі наведені приклади використання запропонованих схем для двох важливих типових матричних перетворень.

Modern applied problems of matrix algebra require algorithms for efficient and reliable operation with a large number of equations. The purpose of such algorithms is to transform these equations with a wide range of numerical values of their coefficients. The transformations are aimed at simplifying the form of matrices in order to make them one of the canonical ones. Many numerical methods, schemes, algorithms have been developed to achieve these goals. The quality of the transformation results is certainly influenced not only by the method itself, but also by the order in which individual operations are performed. In this paper, an acceptable (in our opinion) version of working with known methods and a scheme for their use is proposed when the issues of accuracy and stability, saving RAM and a rational number of operations are consistent with each other. An elementary non-unitary matrix is accepted as the basic one in these transformations. To limit the uncontrolled growth of the coefficients, a certain order of operations is assigned using Virtual matrices, which are based on the management of row / column permutations. Chains of transformations using elementary matrices form cyclical processes and can be universal. The paper gives examples of using the proposed schemes for two important typical matrix transformations.