Запутался с аргументами pdfotrf

Я хочу сделать факторизацию Холецкого в распределенной среде. Для этой цели я использую pdpotrf(). Тем не менее, я изо всех сил пытаюсь понять параметры, необходимые для функции, и они не предоставляют пример C о том, как ее использовать (и пример был бы действительно замечательным).

Предположим, у меня есть матрица NxX, которую я хочу разложить на множители. Тогда какие значения должны быть у параметров? uplo, a и info хорошо определены в моей голове. Как насчет остальных?

n должно быть равно N, я бы сказал. Однако меня смущают desca, ia и ja. Более того, desca бывает глобальным и локальным, чего я не могу понять. Может кто-нибудь объяснить? Может быть с подходом-примером (не обязательно кодом, так как я хочу понять, а не просто копипастить).

Я не знаю Фортран.


person gsamaras    schedule 19.06.2015    source источник
comment
Вам сначала нужно изучить, как этот пакет распределяет матрицы, после этого вы будете знать, что означают аргументы. Вам не нужен Фортран. Вам необходимо прочитать руководство ScaLAPACK.   -  person Vladimir F    schedule 20.06.2015
comment
Если бы я мог найти это, было бы здорово! Нельзя ли по параметрам @VladimirF сказать, как происходит распределение?   -  person gsamaras    schedule 20.06.2015
comment
Вы должны rtfm netlib.org/scalapack/slug/node35.html   -  person Vladimir F    schedule 20.06.2015
comment
@VladimirF Верьте или нет, вчера вечером (когда меня не было дома) я подумал, что мне следует погуглить матрицу распространения в ScaLAPACK. И да, должен РТФМ! Изучу ссылку, которую вы дали. Вы также можете ответить на это.   -  person gsamaras    schedule 21.06.2015
comment
@VladimirF Я прочитал это, но все еще не могу понять. Я продолжу поиски и отпишусь, если нужно.   -  person gsamaras    schedule 22.06.2015


Ответы (1)


Из этого ответа мы имеем:

IA и JA: эти параметры предназначены для предоставления начальной строки и столбца вашей глобальной матрицы внутри большей матрицы. Они актуальны только в том случае, если у вас есть большая матрица и вам нужно только разложение Холецкого подматрицы. В вашем случае IA и JA оба должны быть 1!

Опять же, из этого ответа существует файл ссылка, описывающая desca:

введите здесь описание изображения

person gsamaras    schedule 30.06.2015