Probably someone can help me with a relatively simple question. I need to convert a matrix to a full symmetric matrix again, after using for example pdsyrk. So in general copying the upper triangular part of the matrix to the lower (or vice versa). Is there a simpler way then using pdgeadd and manipulating the diagonal again? If not, can I use pdgeadd with C==A to do this operation in place, so that I do not need the matrix twice in memory?
If someone has a suggestion it would be great.
Thanks a lot,