citibeth
commented
5 years ago I had no idea that tinyvector and array ever even mixed like this. Can we get rid of tinyvector in favor of std::array? I think that would be an improvement
On Fri, Mar 29, 2019 at 07:33 AnderOne notifications@github.com wrote:
There is a problem with to unwrap of expression templates when we use multicomponent Array on the left side and a vector expression on the right.
include <blitz/array.h>
include
int main() {
blitz::Array<blitz::TinyVector<double, 2>, 1> DAT(5);
blitz::TinyVector<double, 2> TMP(3, 7);
if false
//Here we construct TinyVector from expression and fill Array. It's right!
DAT = blitz::TinyVector<double, 2> (TMP * 2 - 1);
else
//This code leads to incorrect behavior!
DAT = TMP * 2 - 1;
endif
std::cout << DAT << std::endl;
return 0;
}
If N-dimensional array is used then it leads to compilation error:
include <blitz/array.h>
include
int main() {
blitz::Array<blitz::TinyVector<double, 2>, 2> DAT(5, 4);
blitz::TinyVector<double, 2> TMP(3, 7);
if false
//Here we construct TinyVector from expression and fill Array. It's right!
DAT = blitz::TinyVector<double, 2> (TMP * 2 - 1);
else
//This code leads to compilation error!
/./blitz/blitz/array/expr.h:194:81: error: no match for call to ‘(const T_expr {aka const blitz::FastTV2Iterator<double, 2>}) (const blitz::TinyVector<int, 2>&)’ T_result operator()(const TinyVector<int, Nrank>& i) const { return iter(i); } /
DAT = TMP * 2 - 1;
endif
std::cout << DAT << std::endl;
return 0;
}
— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/blitzpp/blitz/issues/108, or mute the thread https://github.com/notifications/unsubscribe-auth/AB1cdydM8D_tsZUuEoR6_dPvAhsrXdz9ks5vbhYrgaJpZM4cSZjn .
AnderOne
lutorm
There is a problem with to unwrap of expression templates when we use multicomponent Array on the left side and a vector expression on the right.
If N-dimensional array is used then it leads to compilation error: