| fa046f9f | 12-Sep-2021 |
Junchao Zhang <jczhang@mcs.anl.gov> |
MatProduct: need to check symmetry between symbolic and numeric
Suppose we have code:
MatTransposeMatMult(A,B,MAT_INITIAL_MATRIX,fill,C);
MatTransposeMatMult(E,B,MAT_REUSE_MATRIX,fill,C);
If A
MatProduct: need to check symmetry between symbolic and numeric
Suppose we have code:
MatTransposeMatMult(A,B,MAT_INITIAL_MATRIX,fill,C);
MatTransposeMatMult(E,B,MAT_REUSE_MATRIX,fill,C);
If A is symmetric and E is not, but E has the same nonzero pattern as A, then above code is legitimate.
If we change the MatProductType in MAT_INITIAL_MATRIX from MATPRODUCT_AtB to MATPRODUCT_AB, then
in MAT_REUSE_MATRIX, we need to redo symbolic otherwise C=E^B would be wronlgy computed as C=EB.
show more ...
|