add fits_inconsistent_tdim exception

knodlseder
1 parent 92099a5d
Showing 2 changed files with 53 additions and 0 deletions
include/GException.hpp
src/fits/GException_fits.cpp
@@ -29,6 +29,7 @@
 /* __ Includes ___________________________________________________________ */
 #include <string>                             // string
+#include <vector>                             // string
 #include <sstream>                            // ostringstream
 #include <stdexcept>                          // exception
 #include "GTime.hpp"
@@ -310,6 +311,13 @@ public:
                            int         max_rows,
                            std::string message = "");
     };
+    class fits_inconsistent_tdim : public GExceptionHandler {
+    public:
+        fits_inconsistent_tdim(std::string      origin,
+                               std::vector<int> tdim,
+                               int              number,
+                               std::string      message = "");
+    };
     // GSkymap exceptions
@@ -565,3 +565,48 @@ GException::fits_invalid_nrows::fits_invalid_nrows(std::string origin,
 }
+/***********************************************************************//**
+ * @brief FITS error: Inconsistent TDIM information
+ *
+ * @param[in] origin Method that throws the error.
+ * @param[in] tdim Number of rows.
+ * @param[in] number Expected total number of elements.
+ * @param[in] message Error message.
+ ***************************************************************************/
+GException::fits_inconsistent_tdim::fits_inconsistent_tdim(std::string      origin,
+                                                           std::vector<int> tdim,
+                                                           int              number,
+                                                           std::string      message)
+{
+    // Set origin
+    m_origin  = origin;
+
+    // Set message
+    if (tdim.size() < 1) {
+        m_message  = "Empty TDIM keyword encountered.";
+    }
+    else {
+        // Compute expectation
+        std::string sdim = "("+str(tdim[0]);
+        int         num  = tdim[0];
+        for (int k = 1; k < tdim.size(); ++k) {
+            sdim += ","+str(tdim[k]);
+            num  *= tdim[k];
+        }
+        sdim += ")";
+        
+        // Set message
+        m_message  = "TDIM keyword "+sdim+" predicts "+str(num)+" column"
+                     " elements, while there are "+str(number)+" elements"
+                     " in the column.";
+    }
+
+    // Set optional message
+    if (message.length() > 0)
+        m_message += " "+message;
+
+    // Return
+    return;
+}
+
+