Changeset 16560 for trunk

Timestamp:

06/01/13 19:52:12 (11 years ago)

Author:

tbretz

Message:

Performance improved and heavily tested version of the calibration algorithm.

File:

• trunk/Mars/mcore/factfits.h (modified) (1 diff)

trunk/Mars/mcore/factfits.h

@@ -58 +58 @@
          int16_t *data      = reinterpret_cast<int16_t*>(fBufferRow.data() + offset + fOffsetData);
 
+         /*
          for (uint32_t i=0; i<1440*1024; i+=1024, startCell++)
          {
+             if (*startCell < 0)
+                 continue;
+
              for (uint32_t j=0; j<fNumRoi; j++, data++)
+                 *data += fOffsetCalibration[i + (*startCell+j)%1024];
+         }
+         */
+
+         // This version is faster because the compilers optimization
+         // is not biased by the evaluation of %1024
+         for (int ch=0; ch<1440; ch++)
+         {
+             if (startCell[ch]<0)
              {
-                 if ((*startCell) < 0) continue;
-                 *data += fOffsetCalibration[i + (*startCell+j)%1024];
+                 data += fNumRoi;
+                 continue;
              }
+
+             const int16_t *off = fOffsetCalibration.data() + ch*1024;
+
+             const int16_t *cal = off+startCell[ch];
+             const int16_t *end = data+fNumRoi;
+
+             if (startCell[ch]+fNumRoi>1024)
+             {
+                 while (cal<off+1024)
+                     *data++ += *cal++;
+
+                 cal = off;
+             }
+             while (data<end)
+                 *data++ += *cal++;
          }
     }