... Thanks again for pointing this out so clear. What I now wonder about is : are'nt column_major and column_major_tag redundant ? Both are types and denote...
... Yes, and just for this reason it is good to have both. Construction of an empty object is less expensive. So, if you have to overload a template function,...
Hi Joerg, Thanks for the reply. I have done a little more work on the noalias syntax support to address the issues you raised. ... I think the dangling...
... I must admit that I didn't look into the definition of row_major. But even if you are right (and you are ;o), can you be sure that it will be empty in the...
Hi Toon, hi Kresimir, you both made good points. My main rationale not to change things is that we have three storage categories currently (row major, column...
Hello, I need to compute the condition number of a matrix X with respect to inversion. Matlab is using NORM(X) * NORM(INV(X)) with NORM as max(svd(X)) which...
Hi Alex, ... Ouch! Another bug. I've added typedef matrix_reference<self_type> closure_type; to identity_matrix<> (and zero_matrix<> and something similar to ...
Hello, The recently posted lu works fine (thanks Joerg) but the underlying outer_product does not keep the sparsity of the matrix, it fills in lots of zeroes...
Gunter Winkler
guwi17@...
May 10, 2003 10:03 am
675
Hi all, Geert Ceuppens privately reported a 1.30.0 regression for packed adaptors ... #include <boost/numeric/ublas/matrix.hpp> #include...
Hi all, It may be worthwhile to report a problem with gcc 3.3 to the gcc team. This problem affects uBLAS, and can be seen by trying to compile glucat-0.1.2...
Hi Ian, ... John Maddock reported and explained this problem to me some time ago (thanks!) and I've (hopefully) implemented a *ugly* workaround for it. The ...
Hi Zhang, ... unsigned ... I'd try to add specializations of the form template<> struct promote_traits<std::size_t, unsigned char> { typedef std::size_t...
Hi all, [snip] ... what i ... Just a few days ago I was trying to do the same thing and got the same error. Pretty much by accident, I found that by defining ...
Hello all, today I encountered a strange problem. Please regard the code below (file included as well). It produces a warning ublas.cpp(9) : warning C4172:...
Hello, Any idea if there is out there a library that uses uBLAS' interface, and provides algorithms for eigenvectors calculations? matrix decomposition?...
Hello, after Joergs latest changes the LU of a banded matrix has (finally) a computational complexity of O(n) which is optimal. [You can count the number of...
Gunter Winkler
guwi17@...
May 14, 2003 10:24 am
683
... interface, ... Well - there's some excellent work being done by Kresimir Fresl and Toon Knapen; and by now Joerg Walter is an expert as well ( ;-) to Joerg...
Hi Gunter (hi all), ... I hadn't even the time to announce them ;-) ... You privately pointed out that performance problems could arise for certain computed...
Hi all, ... certain ... assignment ... it2; ... it2; ... I applied a similar patch for sparse proxy evaluators. Gunter reported another problem in some...
Hi Patrick, ... included as well). It produces a warning ... temporary ... 'const double ... Did I overlook ... No. This is an idiosyncrasy of current ublas....
Hi Andy, ... Interesting. ... The former: maybe. The latter: probably. Indexed iterators are an artefact of the time when I believed ublas wouldn't need sparse...
Dear Joerg, Dear developers! The following program (Sorry, I had to include the part from lu.hpp, since the new axpy_lu_factorize produced some kind of...
Hi Paul, ... This ... with ... gcc ... regression ... Yep, I've noticed this one. ... does ... Ouch. Isn't that a regression w.r.t. GCC 3.2 then (or is the...
Dear Joerg, Dear Developers! My last example was not perfect, since it used a matrix, but the same is of course true for vectors: <snip> int main(void) { ...
Hi Benedikt, ... and ... Changed to std::rand() in my local copy (seems like I've been programming too much C in the past ;-) ... implicitely. And yes,...