ZMat

% Try different compression methods
data = rand(100, 100);

% Zlib (default, balanced)
[zlib_out, info] = zmat(data, 1, 'zlib');

% Gzip (similar to zlib)
[gzip_out, info] = zmat(data, 1, 'gzip');

% LZ4 (fastest)
[lz4_out, info] = zmat(data, 1, 'lz4');

% LZMA (highest compression)
[lzma_out, info] = zmat(data, 1, 'lzma');

% Zstd (modern, flexible)
[zstd_out, info] = zmat(data, 1, 'zstd');

% Blosc2 meta-compressor optimized for N-D arrays
data = uint8(rand(100, 100, 50) * 255);

% Blosc2 with different internal codecs
[b_lz4, info] = zmat(data, 1, 'blosc2lz4');
[b_zlib, info] = zmat(data, 1, 'blosc2zlib');
[b_zstd, info] = zmat(data, 1, 'blosc2zstd');

% With multi-threading and byte-shuffling
[compressed, info] = zmat(data, -6, 'blosc2lz4', ...
    'nthread', 4, 'shuffle', 1, 'typesize', 1);

% Decompress
original = zmat(compressed, info);

/* test_zmat.c - Compress data using libzmat */
#include "zmatlib.h"
#include <stdio.h>
#include <string.h>
#include <stdlib.h>

int main() {
    unsigned char input[] = "ZMat compression test data";
    size_t inputsize = strlen((char*)input);
    unsigned char *output = NULL;
    size_t outputsize = 0;
    int status = 0;
    
    /* Compress: zipid=0 (zlib), clevel=-6 */
    zmat_run(inputsize, input, &outputsize, &output, 
             0, &status, -6);
    
    printf("Compressed: %zu -> %zu bytes\n", inputsize, outputsize);
    
    /* Decompress: clevel=0 for decompression */
    unsigned char *decoded = NULL;
    size_t decodedsize = 0;
    zmat_run(outputsize, output, &decodedsize, &decoded,
             0, &status, 0);
    
    printf("Decompressed: %s\n", decoded);
    
    free(output);
    free(decoded);
    return 0;
}

! test_zmat.f90 - Use ZMat in Fortran90
program zmat_demo
    use iso_c_binding
    use zmat_module
    implicit none
    
    character(len=100) :: input_str = "ZMat Fortran compression test"
    integer(c_size_t) :: input_len, output_len, decoded_len
    integer(c_int) :: zipid, status, clevel
    type(c_ptr) :: output_ptr, decoded_ptr
    character(kind=c_char), pointer :: output(:), decoded(:)
    integer :: i
    
    ! Set compression parameters
    input_len = len_trim(input_str)
    zipid = 0          ! 0=zlib, 1=gzip, 3=lzma, 5=lz4, 7=zstd
    clevel = -6        ! compression level (-1 to -9)
    
    ! Compress data
    call zmat_run(input_len, input_str, output_len, &
                  output_ptr, zipid, status, clevel)
    
    call c_f_pointer(output_ptr, output, [output_len])
    print *, "Compressed:", input_len, "->", output_len, "bytes"
    
    ! Decompress (clevel=0 for decompression)
    call zmat_run(output_len, output, decoded_len, &
                  decoded_ptr, zipid, status, 0)
    
    call c_f_pointer(decoded_ptr, decoded, [decoded_len])
    print *, "Decompressed:", (decoded(i), i=1,decoded_len)
    
end program zmat_demo

% Base64 encoding and decoding
str = 'ZMat compression test';

% Encode to base64 (iscompress=2 removes trailing newline)
encoded = char(zmat(str, 2, 'base64'));

% Decode from base64
decoded = char(zmat(encoded, 0, 'base64'));

% Keep all newlines with iscompress=3
encoded_full = char(zmat(str, 3, 'base64'));

% Combine compression + base64 encoding
data = rand(50, 50);
[compressed, info] = zmat(data, 1, 'zlib');
base64_compressed = char(zmat(compressed, 2, 'base64'));