Trisurf Monte Carlo simulator
Samo Penic
2016-02-16 f0bcead40eea062fb556e142f6f1f7a7240bc144
Added base64 encoding to zipped data
2 files modified
163 ■■■■■ changed files
src/snapshot.c 156 ●●●● patch | view | raw | blame | history
src/snapshot.h 7 ●●●●● patch | view | raw | blame | history
src/snapshot.c
@@ -4,6 +4,8 @@
#include<stdlib.h>
#include<stdarg.h>
#include <zlib.h>
#include<inttypes.h>
ts_uint ts_sprintf(ts_string *str, char *fmt, ...){
    va_list ap;
    va_start(ap,fmt);
@@ -26,22 +28,10 @@
    xml_trisurf_vtx_neigh(data,vesicle->vlist);    
    xml_trisurf_vtx_tristar(data,vesicle->vlist);
#ifdef COMPRESSION
    z_stream defstream;
    defstream.zalloc = Z_NULL;
    defstream.zfree = Z_NULL;
    defstream.opaque = Z_NULL;
    defstream.avail_in = data->beg+1;
    defstream.next_in = (unsigned char *)data->string;
    char *compr=(char *)malloc((data->beg+1)*sizeof(char *));
    defstream.avail_out = data->beg+1;
    defstream.next_out = (unsigned char *)compr;
    deflateInit(&defstream, Z_BEST_COMPRESSION);
        deflate(&defstream, Z_FINISH);
        deflateEnd(&defstream);
    fwrite(compr, sizeof(unsigned char), defstream.total_out, fh);
//    fprintf(fh,"%s",compr);
//printf("Uncompressed size is: %lu\n", defstream.total_out);
    free(compr);
    char *compressed;
    ts_uint nbytes=ts_compress_string64(data->string, data->beg, &compressed);
    fwrite(compressed, sizeof(unsigned char), nbytes, fh);
    free (compressed);
#else
    fprintf(fh,"%s", data->string);
#endif
@@ -52,6 +42,12 @@
}
ts_bool xml_trisurf_header(FILE *fh, ts_vesicle *vesicle){
    fprintf(fh, "<trisurfversion></trisurfversion>\n");
    fprintf(fh, "<trisurfdate></trisurfdate>\n");
    fprintf(fh, "<dumpdate></dumpdate>\n");
    fprintf(fh, "<tape>\n");
    fprintf(fh, "</tape>\n");
    fprintf(fh, "<trisurf nvtx=\"%u\" npoly=\"%u\" nfono=\"%u\">\n", vesicle->vlist->n, vesicle->poly_list->n, vesicle->poly_list->poly[0]->vlist->n);
    return TS_SUCCESS;
}
@@ -105,3 +101,131 @@
    return TS_SUCCESS;
}
/* UTILITIES */
/* zlib compression base64 encoded */
/* compressed must not be pre-malloced */
ts_uint ts_compress_string64(char *data, ts_uint data_len, char **compressed){
    z_stream defstream;
    defstream.zalloc = Z_NULL;
    defstream.zfree = Z_NULL;
    defstream.opaque = Z_NULL;
    defstream.avail_in = data_len+1;
    defstream.next_in = (unsigned char *)data;
    char *compr=(char *)malloc(data_len*sizeof(char *));
    defstream.avail_out = data_len+1;
    defstream.next_out = (unsigned char *)compr;
    deflateInit(&defstream, Z_BEST_COMPRESSION);
        deflate(&defstream, Z_FINISH);
        deflateEnd(&defstream);
    /*base64 encode*/
    size_t nbase;
    *compressed=base64_encode((unsigned char *)compr,(size_t)defstream.total_out,&nbase);
    //fwrite(base64, sizeof(unsigned char), nbase, fh);
    free(compr);
    return nbase;
}
/* base64 encoding, taken from http://stackoverflow.com/questions/342409/how-do-i-base64-encode-decode-in-c */
static char encoding_table[] = {'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H',
                                'I', 'J', 'K', 'L', 'M', 'N', 'O', 'P',
                                'Q', 'R', 'S', 'T', 'U', 'V', 'W', 'X',
                                'Y', 'Z', 'a', 'b', 'c', 'd', 'e', 'f',
                                'g', 'h', 'i', 'j', 'k', 'l', 'm', 'n',
                                'o', 'p', 'q', 'r', 's', 't', 'u', 'v',
                                'w', 'x', 'y', 'z', '0', '1', '2', '3',
                                '4', '5', '6', '7', '8', '9', '+', '/'};
static char *decoding_table = NULL;
static int mod_table[] = {0, 2, 1};
char *base64_encode(const unsigned char *data,
                    size_t input_length,
                    size_t *output_length) {
    *output_length = 4 * ((input_length + 2) / 3);
    int i,j;
    char *encoded_data = malloc(*output_length);
    if (encoded_data == NULL) return NULL;
    for (i = 0, j = 0; i < input_length;) {
            uint32_t octet_a = i < input_length ? (unsigned char)data[i++] : 0;
            uint32_t octet_b = i < input_length ? (unsigned char)data[i++] : 0;
            uint32_t octet_c = i < input_length ? (unsigned char)data[i++] : 0;
            uint32_t triple = (octet_a << 0x10) + (octet_b << 0x08) + octet_c;
            encoded_data[j++] = encoding_table[(triple >> 3 * 6) & 0x3F];
            encoded_data[j++] = encoding_table[(triple >> 2 * 6) & 0x3F];
            encoded_data[j++] = encoding_table[(triple >> 1 * 6) & 0x3F];
            encoded_data[j++] = encoding_table[(triple >> 0 * 6) & 0x3F];
    }
    for (i = 0; i < mod_table[input_length % 3]; i++)
            encoded_data[*output_length - 1 - i] = '=';
    return encoded_data;
}
unsigned char *base64_decode(const char *data,
                             size_t input_length,
                             size_t *output_length) {
    int i,j;
    if (decoding_table == NULL) build_decoding_table();
    if (input_length % 4 != 0) return NULL;
    *output_length = input_length / 4 * 3;
    if (data[input_length - 1] == '=') (*output_length)--;
    if (data[input_length - 2] == '=') (*output_length)--;
    unsigned char *decoded_data = malloc(*output_length);
    if (decoded_data == NULL) return NULL;
    for (i = 0, j = 0; i < input_length;) {
            uint32_t sextet_a = data[i] == '=' ? 0 & i++ : decoding_table[(int)data[i++]];
            uint32_t sextet_b = data[i] == '=' ? 0 & i++ : decoding_table[(int)data[i++]];
               uint32_t sextet_c = data[i] == '=' ? 0 & i++ : decoding_table[(int)data[i++]];
            uint32_t sextet_d = data[i] == '=' ? 0 & i++ : decoding_table[(int)data[i++]];
            uint32_t triple = (sextet_a << 3 * 6)
            + (sextet_b << 2 * 6)
            + (sextet_c << 1 * 6)
            + (sextet_d << 0 * 6);
            if (j < *output_length) decoded_data[j++] = (triple >> 2 * 8) & 0xFF;
            if (j < *output_length) decoded_data[j++] = (triple >> 1 * 8) & 0xFF;
            if (j < *output_length) decoded_data[j++] = (triple >> 0 * 8) & 0xFF;
    }
    return decoded_data;
}
void build_decoding_table() {
    decoding_table = malloc(256);
    int i;
    for (i = 0; i < 64; i++)
        decoding_table[(unsigned char) encoding_table[i]] = i;
}
void base64_cleanup() {
    free(decoding_table);
}
src/snapshot.h
@@ -17,4 +17,11 @@
ts_bool xml_trisurf_tria_neigh(ts_string *data, ts_triangle_list *tlist);
ts_bool xml_trisurf_vtx_neigh(ts_string *data, ts_vertex_list *vlist);
ts_bool xml_trisurf_vtx_tristar(ts_string *data, ts_vertex_list *vlist);
/* UTILITIES */
char *base64_encode(const unsigned char *data, size_t input_length, size_t *output_length);
unsigned char *base64_decode(const char *data, size_t input_length, size_t *output_length);
void build_decoding_table();
void base64_cleanup();
ts_uint ts_compress_string64(char *data, ts_uint data_len, char **compressed);
#endif