GammaLib

Several cfitsio functions are used widely in the FITS module, hence they could be define in a single utility function to reduce the number of lines in the code. An example is the “move to HDU” code

    // Move to HDU
    int status = 0;
    status     = __ffmahd(FPTR(vptr), (FPTR(vptr)->HDUposition)+1, NULL,
                          &status);
    if (status != 0) {
        std::string msg = "Unable to move FITS file pointer to extension" 
                          " number "+
                          gammalib::str((FPTR(vptr)->HDUposition))+".";
        throw GException::fits_error(G_READ_STR, status, msg);
    }

Here, the NULL pointer could be replaced by an integer that holds the HDU type, and this type could be returned by the function when necessary.