v1 ublas_elementwise_div(const v1& a, const v2& b)
 {
   typedef typename v1::value_type D;
   v1 c(a.size());
   std::transform(a.begin(),a.end(),b.begin(),c.begin(),std::divides<D>());
   return c;
 }