Re: [boost] linear algebra library... again

Subject: Re: [boost] linear algebra library... again
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2010-02-03 18:06:59

DE wrote:
> as far as i can tell now the lib i propose provide a broader variety
> of concepts such as shape of the matrix (rectangular. symmetric etc.).
> this implies that a product of symmetric matrix is symmetric and so on.
> 
>   matrix<double, symmetric> m1, m2, m;
>   //...
>   m = m1*m2;

What does this last statement do in case m1 and m2 don't have the same size? Does it throw an exception, or does it simply use the smaller of the two sizes? Same question for the case that m1*m2 is not symmetric, as in the following example:

m1 =
0 1
1 0

m2 =
1 0
0 2

m1*m2 =
0 2
1 0

Does it throw an exception, or does it silently set
m=
0 2
2 0
(or m=
0 1
1 0) ?

Regards,
Thomas