Hi all, Is there an equivalent of MatCreateMPIAIJWithSplitArrays function for symmetric matrices? Or does this function work with symmetric matrices?
I work with symmetric matrices that are extremely expensive to compute and I don't know their structure in advance. So to do preallocation, it is not practical to compute a matrix twice - firstly only to get d_nnz and o_nnz arrays and secondly to set values. Now, I solve this problem by storing a matrix to a file system, which is faster but still far from ideal. The best solution I see is to compute a matrix only once, store it in memory as split CSR arrays and then just call something like MatCreateMPISBAIJWithSplitArrays. Thanks, Daniel
