... This should map to a matrix_expression<const_self_type> which makes sense. ... But what should row(S,k) = x; mean? Does it assign x to row k and column k...
Gunter Winkler
guwi17@...
Jul 1, 2004 4:39 pm
1631
I think Dave's evaluation of uBlas should be closely looked at and we should address the issues mentioned. So I would like to comment so that we can further...
... I just reposted my own reply as well. I think we have about 100% agreement which is nice to see! Michael -- ___________________________________ Michael...
Salut, (note that I only use the dense matrices of ublas) ... I prefer size1 and size2, since it fits to the sheme of having a tensor structure with a size3. ...
Salut, ... Actually, there"s more that zou can out into matrices, eg. a algabra based on small matrices. This imples, that the elements of a matrix need not be...
Hi Toon, ... I would have like to, but was swamped with other work when uBLAS went through review (I was probably writing the BGL book). ... No, there are...
... Mixing up abstractions is not a good idea. Matrices and tensors are different creatures. True, sometimes it is useful to think of a tensor as a matrix, and...
Hi Jeremy, Thanks for chipping in on this. I don't think we have conversed before, apart from maybe a few patches I submitted for MTL. Although it would be...
Hi Michael, ... Thanks! ... Although automagically inserting temporaries in non-trivial (especially since it must all be done with template meta-programming) I...
Hi Aleksey, You are confusing the matrix abstraction with the multi-dimensional array abstraction. A *linear algebra* algorithm will never need to be...
Hi all, Reply to a sub-issue below. Best regards ... Can anyone explain to me why operator* implies the commutativity of elements? The GluCat library (...
Hi Paul, I agree with you. The use of operator* does not imply commutativity. Cheers, Jeremy ... _______________________________________________ Jeremy Siek...
Salut, ... I meant it the other way round. With operator* it's often assumed, that it is commutative. A while ago we discussed that on ublas, and Jörg said, ...
Hello all, I am not sure whether this has already been discussed. Resizing a matrix<T,F,A> with A=unbounded_array<T> uses reallocation. I am now facing a...
Karl Meerbergen
Karl.Meerbergen@...
Jul 6, 2004 8:51 am
1646
Salut Karl, ... This would not be nice for small vectors, which don't need resize. ... So, I'd prefer this option. But then, why not use std::vector<>. Is...
Salut, ... I took a look at storage.hpp, and it seems to me, that also for std::vector<>, resize() will copy the elements, if resized to smaller arrays. Have I...
... I tried a small test problem with std::vector . First I created a 5x5 matrix , then resized it to 5x4. Then I resized to 5x6. The capacity() and begin()...
Karl Meerbergen
Karl.Meerbergen@...
Jul 6, 2004 9:18 am
1649
Salut, ... it's starting in line 1073 in storage.hpp, at least in the boost_1_31_0 version: template<class T> BOOST_UBLAS_INLINE void resize ( std::vector<T>...
Hi Jeremy, ... OK, this 'eager' evaluation logic seems sound to me. I assume it operates in a bottom up fashion on the expressions replacing operators with...
On Tuesday 06 July 2004 01:30, "Aleksey Gurtovoy" <agurtovoy@...> ... Aleksey, Thanks for point this out. I had forgotten about Vesa contribution. It...
Salut, ... why not, the fusion aka tensor product is a linear operation, B: a NxN matrix C = A ° B , C_{i*N+j,k*N + m} = A[ i , k] * B[ j, m] A ° ( q B +...
... I already made a syntax adaption from uBLAS to ITL. So we should have a look at the ITL defininitions as well. mfg Gunter...
Gunter Winkler
guwi17@...
Jul 6, 2004 11:29 am
1656
Hi all, ... I do not think this is odd. I'd liked very often to have a matrix of matrices, but I purged it from my wishlist due to practical reasons: The...
An example use of a block diagonal matrix is when you apply decomposition to model decentralized planning. Each of the independent subsystems can be...
Hi Deane, I will answer just to ublas-dev as this is getting a bit detailed. I not sure how to express this in a very clear fashion. A block diagonal matrix is...