#1015 Matrix<T> operator*(Matrix<T>,Matrix<T>) docu inaccurate
Metadata
| Property | Value |
|---|---|
| Reported by | Uli Sack (usack@math.fu-berlin.de) |
| Reported at | Jan 9, 2012 11:08 |
| Type | Feature Request |
| Version | Git (pre2.4) [autotools] |
| Operating System | Unspecified / All |
Description
For a quick and dirty test I wanted to use Matrix operator*(Matrix,Matrix) but had to realize that it won't work for any T with block_size other than one due to the row_type of Matrix and the operator* in block_vector_unmanaged. In general this came as a surprise, as I would expect a method for a block matrix to be able handle block matrices of arbitrary block sizes.
Therefore my question is, is there any reasonable way to use that method? Why is it there at all? Up to now I thought that ISTL-philosophy was "do not provide methods that are very costly (here:copying around matrices) or [cannot be | are not] implemented generically". Should it not be removed maybe?
In any case I think that this shortcoming should be documented, now documentation says "Generic matrix multiplication." I don't have to argue that 'generic' doesn't quite hit the mark. Anyway, there is a patch attached that slightly enhances documentation in that respect.