4 Feb
2010
4 Feb
'10
6:09 p.m.
on 04.02.2010 at 20:46 joel falcou wrote :
forgot to answer: it throws an exception
DE wrote: triggering a static assert makes more sense as the type of matrix is known at CT
if you meant that we try to assign to a symmetric matrix a matrix that is not symmetric as in the example m = m1*m2; //m is symmetric then the code does not compile to force the assignment we can write m = reshape<symmetric>(m1*m2); or m = reshape<symmetric, lower_tr>(m1*m2); to get mirrored lower triangular part -- Pavel