I want to do a Cholesky factorization in a distributed environment. For that purpose, I use pdpotrf(). However, I am struggling understanding the parameters needed by the function and they provide no C example on how to use it (and an example would be really great to have).
Assume I have a NxX matrix I want to factorize. Then, what values should the parameters have? uplo
, a
and info
are well defined in my mind. How about the rest?
n
should be equal to N I would say. However, desca
, ia
and ja
are the ones that confuse me. Moreover, desca
is global and local, something that I can't understand. Can someone explain? Maybe with an approach-example (not necessarily code, since I want to understand, not just copy paste).
I do not know any Fortran.
From this answer, we have:
Again, from this answer, there exists a link, that describes
desca
: