array/Array-util.cc:563:            idx = idx * dvx(i) + idxa(i)(0);
array/Array-util.cc:575:              start = start * dvx(i) + xstart;
array/Array-util.cc:576:              step = step * dvx(i) + xstep;
array/CNDArray.cc:897:    a = a * s;
array/MArray.cc:223:    a = a + s;
array/MArray.cc:234:    a = a - s;
array/MArray.cc:245:    a = a * s;
array/MArray.cc:256:    a = a / s;
array/MArray.cc:269:    a = a + b;
array/MArray.cc:280:    a = a - b;
array/dim-vector.h:461:      k = k * rep[i] + idx[i];
array/dim-vector.h:472:      k = k * rep[i] + idx[i];
array/fCNDArray.cc:908:    a = a * s;
array/fCNDArray.cc:917:    a = a / s;
numeric/CollocWt.cc:178:              z = z * z;
numeric/CollocWt.cc:180:              y = y * (ap + y);
numeric/CollocWt.cc:221:                zc = zc - z / (x - root[j-1]);
numeric/CollocWt.cc:224:          z = z / zc;
numeric/CollocWt.cc:225:          x = x - z;
numeric/CollocWt.cc:238:      x = x + sqrt (std::numeric_limits<double>::epsilon ());
numeric/CollocWt.cc:358:            ax = ax / x / x;
numeric/CollocWt.cc:361:            ax = ax / (1.0 - x) / (1.0 - x);
numeric/CollocWt.cc:365:          y = y + vect[j];
numeric/eigs-base.cc:460:              AminusSigmaB = AminusSigmaB - sigma * tmp *
numeric/eigs-base.cc:464:            AminusSigmaB = AminusSigmaB - sigma *
numeric/eigs-base.cc:468:        AminusSigmaB = AminusSigmaB - sigma * b;
numeric/eigs-base.cc:483:      AminusSigmaB = AminusSigmaB - sigmat;
numeric/eigs-base.cc:555:            AminusSigmaB = AminusSigmaB - tmp;
numeric/eigs-base.cc:558:        AminusSigmaB = AminusSigmaB - sigma * b;
numeric/eigs-base.cc:625:              AminusSigmaB = AminusSigmaB - tmp * b.hermitian () * b *
numeric/eigs-base.cc:629:            AminusSigmaB = AminusSigmaB - sigma * b.hermitian () * b;
numeric/eigs-base.cc:632:        AminusSigmaB = AminusSigmaB - sigma * b;
numeric/eigs-base.cc:647:      AminusSigmaB = AminusSigmaB - sigmat;
numeric/eigs-base.cc:719:            AminusSigmaB = AminusSigmaB - tmp;
numeric/eigs-base.cc:722:        AminusSigmaB = AminusSigmaB - sigma * b;
numeric/lo-specfun.cc:3253:      term = term * temp * rx / (pp + ai);
numeric/lo-specfun.cc:3254:      value = value + term;
numeric/lo-specfun.cc:3259:          value = value * exp (pp * gnulib::log (xx)
numeric/lo-specfun.cc:3269:      ai = ai + 1.0;
numeric/lo-specfun.cc:3270:      ns = ns - 1;
numeric/lo-specfun.cc:3283:          psq = psq + 1.0;
numeric/lo-specfun.cc:3455:              g = g / 3.0;
numeric/lo-specfun.cc:3481:          g = g / 3.0;
numeric/lo-specfun.cc:3907:        ans = ans / (pow (-1.0, n + 1) / xgamma (double (n+1)));
numeric/oct-norm.cc:441:  x = x / vector_norm (x, p);
numeric/oct-norm.cc:454:      z = z * m;
numeric/oct-rand.cc:137:  i = i > 0 ? i : -i;
numeric/oct-rand.cc:142:    i = i % hi;
numeric/randgamma.c:76:  r = r / sum(r)
numeric/randpoisson.c:328:      P = P*lambda/(double)tableidx;
numeric/randpoisson.c:359:          P = P*lambda/(double)tableidx;
numeric/randpoisson.c:390:      P = P*lambda/(double)tableidx;
numeric/randpoisson.c:408:          P = P*lambda/(double)tableidx;
util/oct-inttypes.cc:703:            retval = retval * a_val;
util/oct-inttypes.cc:705:          b_val = b_val >> 1;
util/oct-inttypes.cc:708:            a_val = a_val * a_val;