I have also written the following test:
   pele::Array<double> v1(6,1);
   pele::Array<double> v2(6,2);
   pele::Array<double> v4(6,4);
   
   pele::Array<double> v = v4/(v1*v2);
    v -= 2;

and everything works as it should.