Приложение 1_2.
Описания глобальных параметров подпрограмм
при обработке ленточных и трехдиагональных матриц.

В отличие от плотных матриц ленточные матрицы распределяются по одномерной решетке процессов (1 * P) блочно - столбцовым способом (см. во второй части методического пособия п.5). В подпрограммы обработки ленточных матриц также можно передать не всю исходную распределенную глобальную матрицу, а ее подматрицу sub (A), которая должна быть квадратной (N * N). При этом предполагается, что первый элемент sub (A) имеет в исходной матрице координаты (1, JA). Таким образом, параметры, передаваемые в подпрограммы обработки ленточных матриц, следующие.

Для подматрицы sub(A):

• число строк и столбцов подматрицы N;
• локальный массив A, содержащий локальную часть исходной распределенной матрицы;
• номер первого столбца подматрицы JA;
• дескриптор DESCA исходной глобальной распределенной матрицы;
• число наддиагоналей BWU и число поддиагоналей BWL (для несимметричных матриц) или ширина ленты BW (для симметричных матриц).

Трехдиагональные матрицы, как и ленточные, распределяются по одномерной решетке процессов (1 * P) (см. во второй части методического пособия п.6). Локальные части таких матриц представляются не в виде единого двумерного массива, а в виде трех (для несимметричных матриц) или двух (для симметричных матриц) векторов, которые являются локальными частями трех (или двух) глобальных векторов, содержащих соответственно элементы трех (или двух) диагоналей исходной трехдиагональной матрицы. Таким образом, параметры, передаваемые в подпрограммы обработки трехдиагональных матриц, следующие.

Для подматрицы sub(A):

• число строк и столбцов подматрицы N;
• локальный массив D, содержащий локальную часть исходной главной диагонали матрицы A;
• локальный массив DL, содержащий локальную часть исходного вектора с элементами поддигонали матрицы A;
• локальный массив DU, содержащий локальную часть исходного вектора с элементами наддиагонали матрицы A;
• номер первого столбца подматрицы JA;
• дескриптор DESCA исходной глобальной распределенной матрицы.

Правила описания параметров ленточных и трехдиагональных матриц аналогичны приведенным в п.5.3 правилам описания параметров плотных матриц.