Я хочу сделать факторизацию Холецкого в распределенной среде. Для этого я использую pdpotrf (). Тем не менее, я изо всех сил пытаюсь понять параметры, необходимые для функции, и они не предоставляют пример C на том, как его использовать (и пример был бы действительно здорово иметь).
Предположим, у меня есть матрица NxX, которую я хочу разложить. Тогда какие значения должны иметь параметры? uplo
, a
а также info
хорошо определены в моем уме. Как насчет отдыха?
n
должен быть равен N, я бы сказал. Тем не мение, desca
, ia
а также ja
это те, которые смущают меня. Более того, desca
является глобальным и локальным, то, что я не могу понять. Может кто-нибудь объяснить? Может быть, с подходом-примером (не обязательно кодом, поскольку я хочу понять, а не просто копировать, вставить).
я делаю не знаю любой фортран.
Из этого ответ, у нас есть:
IA
а такжеJA
: эти параметры предназначены для обеспечения начальной строки и столбца вашей глобальной матрицы внутри большей матрицы. Они актуальны только в том случае, если у вас большая матрица, и вы хотите только декомплекс Холецкого для подматрицы. В твоем случае,IA
а такжеJA
оба должны быть1
!
Опять же, из этого ответа существует ссылка на сайт, это описывает desca
: