Changeset 9ee2ce3 for thirdparty/SZ/sz/src/sz_double_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/sz_double_pwr.c (modified) (2 diffs)

thirdparty/SZ/sz/src/sz_double_pwr.c

@@ -24 +24 @@
 #include "zlib.h"
 #include "rw.h"
+#include "utility.h"
 
 void compute_segment_precisions_double_1D(double *oriData, size_t dataLength, double* pwrErrBound, unsigned char* pwrErrBoundBytes, double globalPrecision)
@@ -1773 +1774 @@
         free_TightDataPointStorageD(tdps);
 }
+
+#include <stdbool.h>
+
+void SZ_compress_args_double_NoCkRngeNoGzip_1D_pwr_pre_log(unsigned char** newByteData, double *oriData, double pwrErrRatio, size_t dataLength, size_t *outSize, double min, double max){
+
+        double * log_data = (double *) malloc(dataLength * sizeof(double));
+
+        unsigned char * signs = (unsigned char *) malloc(dataLength);
+        memset(signs, 0, dataLength);
+        // preprocess
+        double max_abs_log_data;
+    if(min == 0) max_abs_log_data = fabs(log2(fabs(max)));
+    else if(max == 0) max_abs_log_data = fabs(log2(fabs(min)));
+    else max_abs_log_data = fabs(log2(fabs(min))) > fabs(log2(fabs(max))) ? fabs(log2(fabs(min))) : fabs(log2(fabs(max)));
+    double min_log_data = max_abs_log_data;
+        bool positive = true;
+        for(size_t i=0; i<dataLength; i++){
+                if(oriData[i] < 0){
+                        signs[i] = 1;
+                        log_data[i] = -oriData[i];
+                        positive = false;
+                }
+                else
+                        log_data[i] = oriData[i];
+                if(log_data[i] > 0){
+                        log_data[i] = log2(log_data[i]);
+                        if(log_data[i] > max_abs_log_data) max_abs_log_data = log_data[i];
+                        if(log_data[i] < min_log_data) min_log_data = log_data[i];
+                }
+        }
+
+        double valueRangeSize, medianValue_f;
+        computeRangeSize_double(log_data, dataLength, &valueRangeSize, &medianValue_f); 
+        if(fabs(min_log_data) > max_abs_log_data) max_abs_log_data = fabs(min_log_data);
+        double realPrecision = log2(1.0 + pwrErrRatio) - max_abs_log_data * 2.23e-16;
+        for(size_t i=0; i<dataLength; i++){
+                if(oriData[i] == 0){
+                        log_data[i] = min_log_data - 2.0001*realPrecision;
+                }
+        }
+    TightDataPointStorageD* tdps = SZ_compress_double_1D_MDQ(log_data, dataLength, realPrecision, valueRangeSize, medianValue_f);
+    tdps->minLogValue = min_log_data - 1.0001*realPrecision;
+    free(log_data);
+    if(!positive){
+            unsigned char * comp_signs;
+                // compress signs
+                unsigned long signSize = sz_lossless_compress(confparams_cpr->losslessCompressor, confparams_cpr->gzipMode, signs, dataLength, &comp_signs);
+                tdps->pwrErrBoundBytes = comp_signs;
+                tdps->pwrErrBoundBytes_size = signSize;
+        }
+        else{
+                tdps->pwrErrBoundBytes = NULL;
+                tdps->pwrErrBoundBytes_size = 0;
+        }
+        free(signs);
+
+    convertTDPStoFlatBytes_double(tdps, newByteData, outSize);
+    if(*outSize>dataLength*sizeof(double))
+            SZ_compress_args_double_StoreOriData(oriData, dataLength+2, tdps, newByteData, outSize);
+
+    free_TightDataPointStorageD(tdps);
+}
+
+void SZ_compress_args_double_NoCkRngeNoGzip_2D_pwr_pre_log(unsigned char** newByteData, double *oriData, double pwrErrRatio, size_t r1, size_t r2, size_t *outSize, double min, double max){
+
+        size_t dataLength = r1 * r2;
+        double * log_data = (double *) malloc(dataLength * sizeof(double));
+
+        unsigned char * signs = (unsigned char *) malloc(dataLength);
+        memset(signs, 0, dataLength);
+        // preprocess
+        double max_abs_log_data;
+    if(min == 0) max_abs_log_data = fabs(log2(fabs(max)));
+    else if(max == 0) max_abs_log_data = fabs(log2(fabs(min)));
+    else max_abs_log_data = fabs(log2(fabs(min))) > fabs(log2(fabs(max))) ? fabs(log2(fabs(min))) : fabs(log2(fabs(max)));
+    double min_log_data = max_abs_log_data;
+        bool positive = true;
+        for(size_t i=0; i<dataLength; i++){
+                if(oriData[i] < 0){
+                        signs[i] = 1;
+                        log_data[i] = -oriData[i];
+                        positive = false;
+                }
+                else
+                        log_data[i] = oriData[i];
+                if(log_data[i] > 0){
+                        log_data[i] = log2(log_data[i]);
+                        if(log_data[i] > max_abs_log_data) max_abs_log_data = log_data[i];
+                        if(log_data[i] < min_log_data) min_log_data = log_data[i];
+                }
+        }
+
+        double valueRangeSize, medianValue_f;
+        computeRangeSize_double(log_data, dataLength, &valueRangeSize, &medianValue_f); 
+        if(fabs(min_log_data) > max_abs_log_data) max_abs_log_data = fabs(min_log_data);
+        double realPrecision = log2(1.0 + pwrErrRatio) - max_abs_log_data * 2.23e-16;
+        for(size_t i=0; i<dataLength; i++){
+                if(oriData[i] == 0){
+                        log_data[i] = min_log_data - 2.0001*realPrecision;
+                }
+        }
+    TightDataPointStorageD* tdps = SZ_compress_double_2D_MDQ(log_data, r1, r2, realPrecision, valueRangeSize, medianValue_f);
+    tdps->minLogValue = min_log_data - 1.0001*realPrecision;
+    free(log_data);
+
+    if(!positive){
+            unsigned char * comp_signs;
+                // compress signs
+                unsigned long signSize = sz_lossless_compress(confparams_cpr->losslessCompressor, confparams_cpr->gzipMode, signs, dataLength, &comp_signs);
+                tdps->pwrErrBoundBytes = comp_signs;
+                tdps->pwrErrBoundBytes_size = signSize;
+        }
+        else{
+                tdps->pwrErrBoundBytes = NULL;
+                tdps->pwrErrBoundBytes_size = 0;
+        }
+        free(signs);
+
+    convertTDPStoFlatBytes_double(tdps, newByteData, outSize);
+    if(*outSize>dataLength*sizeof(double))
+            SZ_compress_args_double_StoreOriData(oriData, dataLength+2, tdps, newByteData, outSize);
+
+    free_TightDataPointStorageD(tdps);
+}
+
+void SZ_compress_args_double_NoCkRngeNoGzip_3D_pwr_pre_log(unsigned char** newByteData, double *oriData, double pwrErrRatio, size_t r1, size_t r2, size_t r3, size_t *outSize, double min, double max){
+
+        size_t dataLength = r1 * r2 * r3;
+        double * log_data = (double *) malloc(dataLength * sizeof(double));
+
+        unsigned char * signs = (unsigned char *) malloc(dataLength);
+        memset(signs, 0, dataLength);
+        // preprocess
+        double max_abs_log_data;
+    if(min == 0) max_abs_log_data = fabs(log2(fabs(max)));
+    else if(max == 0) max_abs_log_data = fabs(log2(fabs(min)));
+    else max_abs_log_data = fabs(log2(fabs(min))) > fabs(log2(fabs(max))) ? fabs(log2(fabs(min))) : fabs(log2(fabs(max)));
+    double min_log_data = max_abs_log_data;
+        bool positive = true;
+        for(size_t i=0; i<dataLength; i++){
+                if(oriData[i] < 0){
+                        signs[i] = 1;
+                        log_data[i] = -oriData[i];
+                        positive = false;
+                }
+                else
+                        log_data[i] = oriData[i];
+                if(log_data[i] > 0){
+                        log_data[i] = log2(log_data[i]);
+                        if(log_data[i] > max_abs_log_data) max_abs_log_data = log_data[i];
+                        if(log_data[i] < min_log_data) min_log_data = log_data[i];
+                }
+        }
+
+        double valueRangeSize, medianValue_f;
+        computeRangeSize_double(log_data, dataLength, &valueRangeSize, &medianValue_f); 
+        if(fabs(min_log_data) > max_abs_log_data) max_abs_log_data = fabs(min_log_data);
+        double realPrecision = log2(1.0 + pwrErrRatio) - max_abs_log_data * 2.23e-16;
+        for(size_t i=0; i<dataLength; i++){
+                if(oriData[i] == 0){
+                        log_data[i] = min_log_data - 2.0001*realPrecision;
+                }
+        }
+    TightDataPointStorageD* tdps = SZ_compress_double_3D_MDQ(log_data, r1, r2, r3, realPrecision, valueRangeSize, medianValue_f);
+    tdps->minLogValue = min_log_data - 1.0001*realPrecision;
+    free(log_data);
+    if(!positive){
+            unsigned char * comp_signs;
+                // compress signs
+                unsigned long signSize = sz_lossless_compress(confparams_cpr->losslessCompressor, confparams_cpr->gzipMode, signs, dataLength, &comp_signs);
+                tdps->pwrErrBoundBytes = comp_signs;
+                tdps->pwrErrBoundBytes_size = signSize;
+        }
+        else{
+                tdps->pwrErrBoundBytes = NULL;
+                tdps->pwrErrBoundBytes_size = 0;
+        }
+        free(signs);
+
+    convertTDPStoFlatBytes_double(tdps, newByteData, outSize);
+    if(*outSize>dataLength*sizeof(double))
+            SZ_compress_args_double_StoreOriData(oriData, dataLength+2, tdps, newByteData, outSize);
+
+    free_TightDataPointStorageD(tdps);
+}