Re: [eigen] Using Triangular Views

[Christoph Hertzberg <chtz@xxxxxxxxxxxxxxxxxxxxxxxx>]

On 2016-08-30 15:03, Gael Guennebaud wrote:
> The situation here is different because tri*tri never compiled (on
> purpose), [...]

Ok, fair enough. At least not worth blocking 3.3
But I'm afraid I don't recall why tri*tri is forbidden on purpose. Isn't 
this one of the most obvious things one wants to do with triangular 
matrices?

I actually just found the following in FullPivLU::reconstructedMatrix():
>   // FIXME the .toDenseMatrix() should not be needed...
>   res = m_lu.leftCols(smalldim)
>             .template triangularView<UnitLower>().toDenseMatrix()
>       * m_lu.topRows(smalldim)
>             .template triangularView<Upper>().toDe

I guess there are other cases for that.

Christoph






--
 Dipl. Inf., Dipl. Math. Christoph Hertzberg

 Universität Bremen
 FB 3 - Mathematik und Informatik
 AG Robotik
 Robert-Hooke-Straße 1
 28359 Bremen, Germany

 Zentrale: +49 421 178 45-6611

 Besuchsadresse der Nebengeschäftsstelle:
 Robert-Hooke-Straße 5
 28359 Bremen, Germany

 Tel.:    +49 421 178 45-4021
 Empfang: +49 421 178 45-6600
 Fax:     +49 421 178 45-4150
 E-Mail:  chtz@xxxxxxxxxxxxxxxxxxxxxxxx

 Weitere Informationen: http://www.informatik.uni-bremen.de/robotik