BeBOP Optimized Sparse Kernel Interface Library
1.0.1h
|
This routine duplicates the functionality of oski_CountZeroRowsCSR(), and could be eliminated.
This routine duplicates the functionality of oski_CountZeroRowsCSR(), and could be eliminated.
Check symmetric and Hermitian full storage cases.
Add symmetric/Hermitian pattern check
The current implementation does not try to re-tune if already tuned.
Check that the new data structure really is faster than the old.
The output properties data structure actually defines a more general property about the diagonal, namely, that it is all ones. However, the available input matrix properties only allow the user to specify whether or not there is an implicit unit diagonal. Thus, it is possible that the user could create an input matrix with an explicit unit diagonal, but this condition is not checked when wrapping the data structure. It might be desirable to do this to make optimized triangular solve for the unit diagonal case more efficient.
Similarly, the oski_matCSR_t data structure has "is_upper" and "is_lower" flags, which could be set even if the user asserts that the matrix has a "general" pattern.