Floating point exception

[pierre-24]

Dear developers,

I would like to report an error that I get when I try to use HDF5 files in XDS using neggia.

The stacktrace looks like this:

...

 UNTRUSTED_RECTANGLE=   1030.  1041.       0.  3270.   (X1,X2, Y1,Y2)
 UNTRUSTED_RECTANGLE=   2070.  2081.       0.  3270.   (X1,X2, Y1,Y2)
 UNTRUSTED_RECTANGLE=      0.  3111.     514.   552.   (X1,X2, Y1,Y2)
 UNTRUSTED_RECTANGLE=      0.  3111.    1065.  1103.   (X1,X2, Y1,Y2)
 UNTRUSTED_RECTANGLE=      0.  3111.    1616.  1654.   (X1,X2, Y1,Y2)
 UNTRUSTED_RECTANGLE=      0.  3111.    2167.  2205.   (X1,X2, Y1,Y2)
 UNTRUSTED_RECTANGLE=      0.  3111.    2718.  2756.   (X1,X2, Y1,Y2)
 PROCESS #  1 COVERS THE IMAGE RANGE        1 ...      50
 NUMBER OF TRUSTED DETECTOR PIXELS  9486379

 DETERMINATION OF THE DARK CURRENT LOOK-UP TABLE "BLANK.cbf"
 -----------------------------------------------------------

 DARK CURRENT LOOK-UP TABLE IS SET CONSTANT TO        0

 DETERMINATION OF THE GAIN LOOK-UP TABLE "GAIN.cbf"
 --------------------------------------------------
 GAIN LOOK-UP TABLE IS DETERMINED FROM DATA IMAGES       1      50
 AND SAMPLED AT   445 X   467 EQUIDISTENT POINTS ALONG X- AND Y.

 *** DEFINITION OF SYMBOLS ***
  GAIN     = MEAN OF THE GAIN-VALUES ON DATA IMAGE
             "GAIN" RELATES POISSON-DISTRIBUTED X-RAY COUNTS
             TO DATA IMAGE PIXEL VALUES SUCH THAT
             GAIN * X-RAY COUNTS = PIXEL VALUE - DARK CURRENT
 MGAIN     = NUMBER OF SAMPLING POINTS ON DATA IMAGE
             USED FOR DETERMINATION OF GAIN
 I/O-FLAG  = ERROR CODE AFTER ACCESSING DATA IMAGE 
              0: NO ERROR
             -1: CANNOT OPEN OR READ IMAGE FILE
             -3: WRONG DATA FORMAT

 FRAME #    GAIN     MGAIN

 [generic_data_plugin] - INFO - generic_open
       + library          = <
 /opt/cecisw/arch/easybuild/2016b/software/Neggia/1.0.1-foss-2016b/lib/dectris-n
 eggia.so>
       + template_name    = <CPS2305_M6_1_data_000001.h5>
 [generic_data_plugin] - ERROR - 'handle' not null
 generic_getfrm: could not open 
 /opt/cecisw/arch/easybuild/2016b/software/Neggia/1.0.1-foss-2016b/lib/dectris-n
 eggia.so ier=          -1
 [generic_data_plugin] - INFO - generic_get_header
NEGGIA ERROR: CANNOT READ PIXEL MASK FROM CPS2305_M6_1_data_000001.h5
 nx,ny,nbyte,qx,qy,number_of_frames,info_array(1:5)=           0           0
           0  0.0000000E+00  0.0000000E+00           0           1           0
           5           3          -1
INFO(1:5)=vendor/major version/minor version/patch/timestamp=   1   0   5   3          -1
 generic_getfrm: data are from Dectris
forrtl: severe (71): integer divide by zero
Image              PC                Routine            Line        Source             
xds                000000000060CF09  Unknown               Unknown  Unknown
xds                000000000060B880  Unknown               Unknown  Unknown
xds                00000000005BFE52  Unknown               Unknown  Unknown
xds                000000000056CB53  Unknown               Unknown  Unknown
xds                0000000000573951  Unknown               Unknown  Unknown
libpthread.so.0    00007F5A0FE0E370  Unknown               Unknown  Unknown
dectris-neggia.so  00007F5A104F1149  Unknown               Unknown  Unknown
dectris-neggia.so  00007F5A104F1DD8  Unknown               Unknown  Unknown
xds                00000000005502FE  generic_getfrm_         24930  MAIN_XDS.f90
xds                000000000054D38E  getfrm_                 24993  MAIN_XDS.f90
xds                000000000056043B  mkgtbl_                 10733  MAIN_XDS.f90
xds                00000000005681BE  init_                    8337  MAIN_XDS.f90
xds                0000000000408271  xds_                    20918  MAIN_XDS.f90
xds                0000000000406256  MAIN__                      1  MAIN_XDS.f90
xds                00000000004036B6  Unknown               Unknown  Unknown
libc.so.6          00007F5A0F85BB35  Unknown               Unknown  Unknown
xds                00000000004035A9  Unknown               Unknown  Unknown

The version corresponding to the tag "v1.0.1" and with the standard CMAKE procedure.

Is it possible that this is coming from your side ?

Thanks in advance !

[pilipp]

I've found the error message "NEGGIA ERROR: CANNOT READ PIXEL MASK FROM CPS2305_M6_1_data_000001.h5" in your output. The pixel mask is usually found in the master file (e.g. CPS2305_M6_1_master.h5). Did you process from a data file?

[pierre-24]

I use the script XDS_to_H5.py, as reported here. The generated XDS.INP is to be found on this gist, but maybe I need to change it.

[dirkboye]

closing this issue due to inactivity. if problem should still remain please re-open with further information.