Hi, Is it correct that vector_temporary_type is not defined for expression types? Some old code that used matrix_range of an expression does no longer compile....
2318
Karl Meerbergen
Karl.Meerbergen@...
Nov 2, 2004 2:50 pm
Hello again, I found some more differences with previous uBlas versions. In vector_expression.hpp, vector_binary::expression1() and ...
2319
Michael Stevens
mmssau
Nov 2, 2004 3:52 pm
Hi Karl, ... Intersting. The vector/matrix_temporary_types are designed as a property of the container concept. Proxies are designed to propagate this via...
2320
Michael Stevens
mmssau
Nov 2, 2004 3:58 pm
... Looks like a silly mistake. ... Hmmm sounds bad. Certain expression2() should not disappear! Thanks for the info, Michael...
2321
Karl Meerbergen
Karl.Meerbergen@...
Nov 2, 2004 4:03 pm
... It is probably possible with a traits class too, that should then be specialized for all expression types!! :-( In any case, I agree that adding too much...
2322
Karl Meerbergen
Karl.Meerbergen@...
Nov 2, 2004 4:07 pm
... If I had access to boost, I could do the changes myself. I am currently using ublas_pure to get our stuff compiled. I suppose some other surprises will...
2323
Gunter Winkler
guwi17@...
Nov 2, 2004 4:26 pm
... Ok. It works, but resize is broken. I changed it to respect the preserve parameter correctly. New timings are on my page. I expected that...
2324
Michael Stevens
mmssau
Nov 2, 2004 6:09 pm
... Sounds good. If it can't wait Toon can commit to CVS Michael...
2325
Karl Meerbergen
Karl.Meerbergen@...
Nov 2, 2004 6:12 pm
... Well, Toon is on holiday. But it can wait till he's back. I can easily correct the errors related to expression1() etc, but I am not familiar with the...
2326
Michael Stevens
mmssau
Nov 2, 2004 7:01 pm
... It tell algorithms what type they should create when the need a temporary vector container before assigning to an expression. Because proxies are not ...
2327
Karl Meerbergen
Karl.Meerbergen@...
Nov 2, 2004 7:06 pm
... I see. I see some difficulties when a proxy is used on a binary_expression. But I think the storage_category could be used. Is that right? The unknown...
2328
Michael Stevens
mmssau
Nov 2, 2004 7:37 pm
... Ahh. I should have made it clear that the above temporaries only occur in assignments. Therefore there is always a LHS which is either a Container or an...
2329
Karl Meerbergen
Karl.Meerbergen@...
Nov 2, 2004 8:36 pm
... So, instead of having a member temporary_vector_type in the proxies, would it not be better to define a traits class to determine temporary_vector_type? ...
2330
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 10:57 am
Hi all, I think we need a discussion about which members in expressions should be available to public. Currently, the expressions do not allow access of their ...
2331
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 11:20 am
... Hi Michael, I added to my version of the code a temporary_traits that computes vector_temporary_type and matrix_temporary_type for containers and proxies. ...
2332
Gunter Winkler
guwi17@...
Nov 3, 2004 11:37 am
... I think, its better to have less public functions. I use constructions like this in my code: // (v * t) [i] = v [i] * t template<class E1, class T2> ...
2333
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 11:47 am
... For the BLAS bindings, access to the data is inevitable, because you need the containers, constants etc. Karl...
2334
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 11:53 am
Hello I have another problem with a piece of code that does no longer compile. It does not seem to be possible to assign a range into a range. Is there a ...
2335
Patrick Kowalzick
patrick_kowa...
Nov 3, 2004 11:59 am
Hi Karl, Do you include "matrix_proxy" and "vector_proxy" ? For vectors it works fine for me. Patrick...
2336
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 12:01 pm
... In order to limit the public functions, we can create an accessor object for each expression that allows to access the data. Karl...
2337
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 12:06 pm
... Patrick I mean, ublas::range a( 0, 1 ) ; ublas::range b ; b = a; Karl...
2338
Patrick Kowalzick
patrick_kowa...
Nov 3, 2004 12:14 pm
... This works fine for vectors: rp(ublas::range(0,3)) = p_l(ublas::range(3,6)); rp(ublas::range(3,6)) = p_l(ublas::range(0,3)); I will try yours.... Patrick...
2339
Patrick Kowalzick
patrick_kowa...
Nov 3, 2004 12:22 pm
... Not working here as well. But copy ctor works. Patrick...
2340
Gunter Winkler
guwi17@...
Nov 3, 2004 12:46 pm
... a range cannot be assigned: from storage2.hpp: class basic_range: public storage_array<basic_range<Z, D> > ... and template<class E> class storage_array: ...
2341
Gunter Winkler
guwi17@...
Nov 3, 2004 12:47 pm
... I can do things like matrix<double>; A(10,10); A = scalar_matrix<double>(10,10,1.0); project(A, range(1,4), range(1,4)) = 2 * project(A, range(3,6),...
2342
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 12:59 pm
... So it appears, but it was possible in the past! Karl...
2343
Karl Meerbergen
Karl.Meerbergen@...
Nov 3, 2004 2:18 pm
... Hello all, I think we have to make a decision to what we allow to do with expressions. 1) Do we want the user to interpret expressions and so have access...
2344
Karl Meerbergen
Karl.Meerbergen@...
Nov 4, 2004 8:14 am
Hello all, I used uBLAS_pure on IBM. After the modifications I did yesterday on my local version to get old stuff working, everything seems to work well on ...
2345
Karl Meerbergen
Karl.Meerbergen@...
Nov 4, 2004 1:35 pm
Hello again, Sorry for continuing reporting problems. I think, we have to do a little exercise about what is the size_type in the sparse matrix types. For...
2346
Ian McCulloch
ianmcc@...
Nov 4, 2004 6:51 pm
(from the boost-user group) ... is ... I don't know why any of the ublas devs have commented on this at boost-user (or maybe they have but I missed it?). There...
