[petsc-users] MatSetValues for nonzeros of a sparse matrix

[Umut Tabak]

Jed Brown wrote:
>   MatSetValues(A, 1, &idxm[r], 1, &idxn[r], &valuesMat[r*4+r], INSERT_VALUES);
>
> or
>
>   MatSetValue(A, idxm[r], idxn[r], valuesMat[r*4+r], INSERT_VALUES);
>
>   
Thanks for the quick replies, the above, I guess, still means to keep 
the zeros but below explanation is more practical from an efficiency 
point of view.
>> Also I keep the row and column indices of the non-zero values and as
>> well the non-zeros values in arrays of proper type.
>>     
>
> you should loop over rows and insert all entries of that row at once.
>   
This was where I was confused :), I should keep the non-zeros of the row 
and set them with the loop.

Thanks Jed.

Umut