Thanks to everyone! Sorry if I did not answer before but on the we I did not check the group. I think I understand now the mechanism, what I did not understand was why it was upcasting the 2 to Float64 if the matrix was Float32. But probably I'll give numarray a try. thanks again. curzio