Hello,

I want to solve many symmetric linear systems one after another in parallel 
using boomerAMG + KSPCG  and need to make the matrix transfer more efficient. 
Matrices are symmetric in structure and values. boomerAMG + KSPCG work fine.

So far I have been loading the entire matrices but I read in a paper, that it's 
sufficient to load the upper part only and tell petsc that the matrix is 
symmetric using MatSetOption(A,MAT_SYMMETRIC,PETSC_TRUE); Unfortunately all 
computations fail if I load only the upper values and use 
MatSetOption(A,MAT_SYMMETRIC,PETSC_TRUE);
The idea is:
    if (matrix_.symmetric())
    {
        MatSetOption(A,MAT_SYMMETRIC,PETSC_TRUE);
       //load only upper part of the matrix MatSetValues(...)

   }else //asymmetric matrix
           {               //load the entire matrix MatSetValues(...)
           }

Is it possible at all? 

Klaus

Reply via email to