Added Python doc strings

Author: Martin D. Weinberg

expui/BiorthBasis.H

@@ -120,6 +120,7 @@ namespace BasisClasses
     unsigned H5compress  = 5;
     unsigned H5chunk     = 1024*1024; // (1 MB)
     bool     H5shuffle   = true;
+    bool     H5szip      = false;
     //@}
 
     //! Round time key to emulated fixed-point arithmetic
@@ -290,11 +291,12 @@ namespace BasisClasses
     void setCovarFloatType(bool flt) { floatType = flt; }
 
     //! HDF5 compression settings
-    void setCovarH5Compress(unsigned level, unsigned chunksize, bool shuffle)
+    void setCovarH5Compress(unsigned level, unsigned chunksize, bool shuffle, bool szip=false)
     {
       H5compress = level;
       H5chunk    = chunksize;
       H5shuffle  = shuffle;
+      H5szip     = szip;
     }
 
     //! Make covariance after accumulation

expui/BiorthBasis.cc

@@ -4890,14 +4890,22 @@ namespace BasisClasses
       ("rankSize", HighFive::DataSpace::From(nmax)).write(nmax);
 
     if (H5compress) {
+      // Szip parameters
+      const int options_mask = H5_SZIP_NN_OPTION_MASK;
+      const int pixels_per_block = 8;
+
       // Properties for coefficients
       //
       unsigned int csz2 = nmax * ltot * sampleSize;
       HighFive::Chunking data_dims2{std::min<unsigned>(csz2, H5chunk), 1};
 
       dcpl2.add(data_dims2);
       if (H5shuffle) dcpl2.add(HighFive::Shuffle());
-      dcpl2.add(HighFive::Deflate(H5compress));
+      if (H5szip) {
+	dcpl2.add(HighFive::Szip(options_mask, pixels_per_block));
+      } else {
+	dcpl2.add(HighFive::Deflate(H5compress));
+      }
 
       // Properties for  covariance
       //
@@ -4906,7 +4914,11 @@ namespace BasisClasses
 
       dcpl3.add(data_dims3);
       if (H5shuffle) dcpl3.add(HighFive::Shuffle());
-      dcpl3.add(HighFive::Deflate(H5compress));
+      if (H5szip) {
+	dcpl3.add(HighFive::Szip(options_mask, pixels_per_block));
+      } else {
+	dcpl3.add(HighFive::Deflate(H5compress));
+      }
     }
 
     // Pack the coefficient data

pyEXP/BasisWrappers.cc

@@ -1434,21 +1434,26 @@ void BasisFactoryClasses(py::module &m)
          )")
     .def("setCovarH5Compress", &BasisClasses::BiorthBasis::setCovarH5Compress,
          R"(
-	 Set the HDF5 compression level for covariance storage in HDF5.
+	 Set the HDF5 compression level for covariance storage in HDF5.  The Szip
+     	 compression algorithm may also be enabled but seems to not have better
+         performance for these data than the standard Gzip compression.  Gzip with
+      	 level 5 and shuffling is enabled by default.
 
          Parameters
          ----------
          compress : int
-		 HDF5 compression level 0-9 (default: 5)
+		 HDF5 Gzip compression level 0-9 (default: 5)
          chunkSize : int
 		 HDF5 chunk size for dataset storage (default: 1024*1024)
          shuffle : bool
 		 Use shuffle filter if true (default: true)
+         szip : bool
+		 Use Szip compression algorithm (default: false)
 
          Returns
          -------
          None
-         )", py::arg("compress")=5, py::arg("chunkSize")=1024*1024, py::arg("shuffle")=true)
+         )", py::arg("compress")=5, py::arg("chunkSize")=1024*1024, py::arg("shuffle")=true, py::arg("azip")=false)
     .def("getFields", &BasisClasses::BiorthBasis::getFields,
 	 R"(
          Return the field evaluations for a given cartesian position. The