Changeset 9ee2ce3 for thirdparty/SZ/sz/src/szd_float_pwr.c

Timestamp:

09/28/18 16:32:55 (6 years ago)

Author:

Hal Finkel <hfinkel@…>

Branches:

master, pympi

Children:

e6aa0eb

Parents:

abca157

git-author:

Hal Finkel <hfinkel@…> (09/28/18 16:32:55)

git-committer:

Hal Finkel <hfinkel@…> (09/28/18 16:32:55)

Message:

importing new SZ files

File:

• thirdparty/SZ/sz/src/szd_float_pwr.c (modified) (2 diffs)

thirdparty/SZ/sz/src/szd_float_pwr.c

@@ -17 +17 @@
 #include "Huffman.h"
 #include "sz_float_pwr.h"
+#include "utility.h"
 //#include "rw.h"
 //
@@ -1350 +1351 @@
         free(groupID);
 }
+
+void decompressDataSeries_float_1D_pwr_pre_log(float** data, size_t dataSeriesLength, TightDataPointStorageF* tdps) {
+
+        decompressDataSeries_float_1D(data, dataSeriesLength, tdps);
+        float threshold = tdps->minLogValue;
+        if(tdps->pwrErrBoundBytes_size > 0){
+                unsigned char * signs;
+                sz_lossless_decompress(confparams_dec->losslessCompressor, tdps->pwrErrBoundBytes, tdps->pwrErrBoundBytes_size, &signs, dataSeriesLength);
+                for(size_t i=0; i<dataSeriesLength; i++){
+                        if((*data)[i] < threshold) (*data)[i] = 0;
+                        else (*data)[i] = exp2((*data)[i]);
+                        if(signs[i]) (*data)[i] = -((*data)[i]);
+                }
+                free(signs);
+        }
+        else{
+                for(size_t i=0; i<dataSeriesLength; i++){
+                        if((*data)[i] < threshold) (*data)[i] = 0;
+                        else (*data)[i] = exp2((*data)[i]);
+                }
+        }
+
+}
+
+void decompressDataSeries_float_2D_pwr_pre_log(float** data, size_t r1, size_t r2, TightDataPointStorageF* tdps) {
+
+        size_t dataSeriesLength = r1 * r2;
+        decompressDataSeries_float_2D(data, r1, r2, tdps);
+        float threshold = tdps->minLogValue;
+        if(tdps->pwrErrBoundBytes_size > 0){
+                unsigned char * signs;
+                sz_lossless_decompress(confparams_dec->losslessCompressor, tdps->pwrErrBoundBytes, tdps->pwrErrBoundBytes_size, &signs, dataSeriesLength);
+                for(size_t i=0; i<dataSeriesLength; i++){
+                        if((*data)[i] < threshold) (*data)[i] = 0;
+                        else (*data)[i] = exp2((*data)[i]);
+                        if(signs[i]) (*data)[i] = -((*data)[i]);
+                }
+                free(signs);
+        }
+        else{
+                for(size_t i=0; i<dataSeriesLength; i++){
+                        if((*data)[i] < threshold) (*data)[i] = 0;
+                        else (*data)[i] = exp2((*data)[i]);
+                }
+        }
+
+}
+
+void decompressDataSeries_float_3D_pwr_pre_log(float** data, size_t r1, size_t r2, size_t r3, TightDataPointStorageF* tdps) {
+
+        size_t dataSeriesLength = r1 * r2 * r3;
+        decompressDataSeries_float_3D(data, r1, r2, r3, tdps);
+        float threshold = tdps->minLogValue;
+        if(tdps->pwrErrBoundBytes_size > 0){
+                unsigned char * signs;
+                sz_lossless_decompress(confparams_dec->losslessCompressor, tdps->pwrErrBoundBytes, tdps->pwrErrBoundBytes_size, &signs, dataSeriesLength);
+                for(size_t i=0; i<dataSeriesLength; i++){
+                        if((*data)[i] < threshold) (*data)[i] = 0;
+                        else (*data)[i] = exp2((*data)[i]);
+                        if(signs[i]) (*data)[i] = -((*data)[i]);
+                }
+                free(signs);
+        }
+        else{
+                for(size_t i=0; i<dataSeriesLength; i++){
+                        if((*data)[i] < threshold) (*data)[i] = 0;
+                        else (*data)[i] = exp2((*data)[i]);
+                }
+        }
+}
 #pragma GCC diagnostic pop