Подпрограмма:  PAR_SYGVZ

Назначение

Выбор пaраметров распараллеливания для целевых программ комплекса PARALG, вычисления собственных значений и собственных векторов в обобщенной проблеме собственных значений для вещественных симметричных матриц

Описание

Эта служебная подпрограмма выбирает и вычисляет за пользователя параметры распараллеливания, которые конкретизируют распределение вычислений между параллельными процессами и значения которых необходимо передавать целевым программам вычисления собственных значений и собственных векторов в обобщенной проблеме собственных значений для вещественных симметричных матриц

Литература:

http://num_anal.srcc.msu.ru/par_prog/
http://software.intel.com/en-us/articles/intel-math-kernel-library-documentation/
http://www.netlib.org/scalapack/slug/index.html

Использование 

         CALL  PAR_SYGVZ(NAME, N, NPROCS, MMIN, NB, NPROW, NPCOL)

   Параметры
NAME - имя целевой подпрограммы, с помощью которой пользователь собирается решать свою задачу (одно из PDSYGV2, PDSYGV4, PDSYGV6); (входной параметр, тип символьный);
N - порядок исходных матриц (входной параметр, тип целый);
NPROCS - число процессов, которое может выделить пользователь для решения задачи (входной параметр, тип целый);
MMIN - минимальное количество локальной памяти, которое необходимо выделить на каждом из параллельных процессов для решения задачи (выходной параметр, тип целый);
NB - величина блока, на которые делятся исходные матрицы; NB полагается равным 0, если пользователь доверяет выбор его значения этой подпрограмме (входной и выходной параметр, тип целый);
NPROW - число строк в решетке процессов (выходной параметр, тип целый);
NPCOL - число столбцов в решетке процессов (выходной параметр, тип целый); 
   Вызываемые подпрограммы: нет 
   Замечания по использованию
  1.  Подпрограмма является последовательной фортранной программой и не использует распределенную память (требуется только один процессор)
  2.  Используются подпрограммы ICEIL, NUMROC ( из библиотеки ScaLAPACK_TOOLS)

Пример использования

Необходимо выбрать размер блока NB, число строк и столбцов решетки процессов (NPROW,NPCOL) и вычислить размер локальной памяти, необходимой для вычисления собственных значений и собственных векторов в обобщенной проблеме собственных значений для вещественных симметричных матриц с помощью подпрограммы PDSYGV2.
Матрицы - квадратные симметричные порядка 5.
Выделяемое число процессов NPROCS = 4. NB = 0.

Фрагмент фортранного текста вызывающей программы
(полный текст теста можно получить в tpr_sygvz.zip) 

      PROGRAM TPR_SYGVZ
*
      INTEGER                 N, NPROCS, MMIN, NB, NPROW, NPCOL
      CHARACTER*(*)   NAME
      PARAMETER          ( NAME = 'PDSYGV2')
      EXTERNAL             PAR_SYGVZ
*
      N = 5
      NPROCS = 4
      NB=0
*
      CALL PAR_SYGVZ(NAME, N, NPROCS, MMIN, NB, NPROW, NPCOL)
      STOP
      END

Результаты:

      NB = 4
      NPROW = 2
      NPCOL = 2
      MMIN = 5613