- double x, y; // Rectang coords of center of pixel
- int ix, iy;
- for (x = xMin + xInc / 2, ix = 0; ix < nx; x += xInc, ix++)
- for (y = yMin + yInc / 2, iy = 0; iy < ny; y += yInc, iy++) {
+ double xstart = xMin + xInc / 2;
+
+#if HAVE_OPENMP
+ #pragma omp parallel for
+#endif
+ for (int ix = 0; ix < nx; ix++) {
+ double x = xstart + (ix * xInc);
+ double y = yMin + yInc / 2;
+ for (int iy = 0; iy < ny; iy++, y += yInc) {