Registered Member
|
Hello,
I have a problem concerning sparse matrix multiplication and storage order. I have tried the following code:
To unknown reasons the storage order is changed from column-major to row-major when evaluating the expression U*P*U.T. I know the combination of auto keyword and template expressions is difficult, but as far as I know the above expression should be safe and valid. When I replace auto by sparse_type, the storage order is unchanged. However in the end I would like to do something like:
in order to remove the creation of temporary objects, which is not possible due to mismatching storage orders. Probably, I missed something in the docs concerning triple matrix products. Thank you for your help. Matthias // Edit: added case D, which also changes the storage order. |
Registered users: Bing [Bot], Google [Bot], Sogou [Bot]