.. include:: defs.rst .. _hl-api: High-Level C API ================ The |libzfp| C API provides functionality for sequentially compressing and decompressing whole integer and floating-point arrays or single blocks. It is broken down into a :ref:`high-level API ` and a :ref:`low-level API `. The high-level API handles compression of entire arrays and supports a variety of back-ends (e.g., serial, OpenMP). The low-level API exists for processing individual, possibly partial blocks as well as reduced-precision integer data less than 32 bits wide. Both C APIs are declared in :file:`zfp.h`. The following sections are available: * :ref:`hl-macros` * :ref:`hl-types` * :ref:`hl-data` * :ref:`hl-functions` * :ref:`hl-func-bitstream` * :ref:`hl-func-stream` * :ref:`hl-func-exec` * :ref:`hl-func-config` * :ref:`hl-func-field` * :ref:`hl-func-codec` .. _hl-macros: Macros ------ .. _version-id: .. c:macro:: ZFP_VERSION_MAJOR .. c:macro:: ZFP_VERSION_MINOR .. c:macro:: ZFP_VERSION_PATCH .. c:macro:: ZFP_VERSION_TWEAK Macros identifying the |zfp| library version (*major*.\ *minor*.\ *patch*.\ *tweak*). :c:macro:`ZFP_VERSION_TWEAK` is new as of |zfp| |verrelease| and is used to mark intermediate develop versions (unofficial releases). ---- .. c:macro:: ZFP_VERSION_DEVELOP Macro signifying that the current version is an intermediate version that differs from the last official release. This macro is undefined for official releases; when defined, its value equals 1. Note that this macro may be defined even if the four :ref:`version identifiers ` have not changed. Available as of |zfp| |verrelease|. ---- .. c:macro:: ZFP_VERSION A single integer constructed from the four :ref:`version identifiers `. This integer can be generated by :c:macro:`ZFP_MAKE_VERSION` or :c:macro:`ZFP_MAKE_FULLVERSION`. Its value equals the global constant :c:data:`zfp_library_version`. .. note:: Although :c:macro:`ZFP_VERSION` increases monotonically with release date and with the four :ref:`version identifiers ` it depends on, the mapping to :c:macro:`ZFP_VERSION` changed with the introduction of :c:macro:`ZFP_VERSION_TWEAK` in |zfp| |verrelease|. Going forward, we recommend using :c:macro:`ZFP_MAKE_VERSION` or :c:macro:`ZFP_MAKE_FULLVERSION` in conditional code that depends on :c:macro:`ZFP_VERSION`, e.g., :code:`#if ZFP_VERSION >= ZFP_MAKE_VERSION(1, 0, 0)`. Note that such constructions should not be used with older versions of |zfp|, e.g., :code:`if (zfp_library_version == ZFP_MAKE_VERSION(0, 5, 5))` will not give the expected result with binary versions of |libzfp| before version |verrelease|. ---- .. c:macro:: ZFP_VERSION_STRING :c:macro:`ZFP_VERSION_STRING` is a string literal composed of the four :ref:`version identifiers `. It is a component of :c:data:`zfp_version_string`. ---- .. c:macro:: ZFP_MAKE_VERSION(major, minor, patch) .. c:macro:: ZFP_MAKE_VERSION_STRING(major, minor, patch) Utility macros for constructing :c:macro:`ZFP_VERSION` and :c:macro:`ZFP_VERSION_STRING`, respectively. Available as of |zfp| |verrelease|, these macros may be used by applications to test for a certain |zfp| version number, e.g., :code:`#if ZFP_VERSION >= ZFP_MAKE_VERSION(1, 0, 0)`. ---- .. c:macro:: ZFP_MAKE_FULLVERSION(major, minor, patch, tweak) .. c:macro:: ZFP_MAKE_FULLVERSION_STRING(major, minor, patch, tweak) Utility macros for constructing :c:macro:`ZFP_VERSION` and :c:macro:`ZFP_VERSION_STRING`, respectively. Includes tweak version used by intermediate develop versions. Available as of |zfp| |verrelease|, these macros may be used by applications to test for a certain |zfp| version number, e.g., :code:`#if ZFP_VERSION >= ZFP_MAKE_FULLVERSION(1, 0, 0, 2)`. ---- .. c:macro:: ZFP_CODEC Macro identifying the version of the compression CODEC. See also :c:data:`zfp_codec_version`. ---- .. c:macro:: ZFP_MIN_BITS .. c:macro:: ZFP_MAX_BITS .. c:macro:: ZFP_MAX_PREC .. c:macro:: ZFP_MIN_EXP Default compression parameter settings that impose no constraints. The largest possible compressed block size, corresponding to 4D blocks of doubles, is given by :c:macro:`ZFP_MAX_BITS`. See also :c:type:`zfp_stream`. ---- .. c:macro:: ZFP_META_NULL Null representation of the 52-bit encoding of field metadata. This value is returned by :c:func:`zfp_field_metadata` when the field metadata cannot be encoded in 64 bits, such as when the array dimensions are too large (see :ref:`limitations`). In addition to signaling error, this value is guaranteed not to represent valid metadata. ---- .. _header-macros: The :code:`ZFP_HEADER` bit mask specifies which portions of a header to output (if any). The constants below should be bitwise ORed together. Use :c:macro:`ZFP_HEADER_FULL` to output all header information available. The compressor and decompressor must agree on which parts of the header to read/write. See :c:func:`zfp_read_header` and :c:func:`zfp_write_header` for how to read and write header information. .. c:macro:: ZFP_HEADER_MAGIC Magic constant that identifies the data as a |zfp| stream compressed using a particular CODEC version. .. c:macro:: ZFP_HEADER_META Array size and scalar type information stored in the :c:type:`zfp_field` struct. .. c:macro:: ZFP_HEADER_MODE Compression mode and parameters stored in the :c:type:`zfp_stream` struct. .. c:macro:: ZFP_HEADER_FULL Full header information (bitwise OR of all :code:`ZFP_HEADER` constants). ---- .. c:macro:: ZFP_MAGIC_BITS .. c:macro:: ZFP_META_BITS .. c:macro:: ZFP_MODE_SHORT_BITS .. c:macro:: ZFP_MODE_LONG_BITS .. c:macro:: ZFP_HEADER_MAX_BITS .. c:macro:: ZFP_MODE_SHORT_MAX Number of bits used by each portion of the header. These macros are primarily informational and should not be accessed by the user through the high-level API. For most common compression parameter settings, only :c:macro:`ZFP_MODE_SHORT_BITS` bits of header information are stored to encode the mode (see :c:func:`zfp_stream_mode`). ---- .. _data-macros: The :code:`ZFP_DATA` bit mask specifies which portions of array data structures to compute total storage size for. These constants should be bitwise ORed together. Use :c:macro:`ZFP_DATA_ALL` to count all storage used. .. c:macro:: ZFP_DATA_UNUSED Allocated but unused data. .. c:macro:: ZFP_DATA_PADDING Padding for alignment purposes. .. c:macro:: ZFP_DATA_META Class members and other fixed-size storage. .. c:macro:: ZFP_DATA_MISC Miscellaneous uncategorized storage. .. c:macro:: ZFP_DATA_PAYLOAD Compressed data encoding array elements. .. c:macro:: ZFP_DATA_INDEX Block :ref:`index ` information. .. c:macro:: ZFP_DATA_CACHE Uncompressed :ref:`cached ` data. .. C:macro:: ZFP_DATA_HEADER :ref:`header` information. .. c:macro:: ZFP_DATA_ALL All storage (bitwise OR of all :code:`ZFP_DATA` constants). ---- .. c:macro:: ZFP_ROUND_FIRST .. c:macro:: ZFP_ROUND_NEVER .. c:macro:: ZFP_ROUND_LAST Available rounding modes for :c:macro:`ZFP_ROUNDING_MODE`, which specifies at build time how |zfp| performs rounding in lossy compression mode. .. _hl-types: Types ----- .. c:type:: zfp_stream The :c:type:`zfp_stream` struct encapsulates all information about the compressed stream for a single block or a collection of blocks that represent an array. See the section on :ref:`compression modes ` for a description of the members of this struct. :: typedef struct { uint minbits; // minimum number of bits to store per block uint maxbits; // maximum number of bits to store per block uint maxprec; // maximum number of bit planes to store int minexp; // minimum floating point bit plane number to store bitstream* stream; // compressed bit stream zfp_execution exec; // execution policy and parameters } zfp_stream; ---- .. c:type:: zfp_execution The :c:type:`zfp_stream` also stores information about how to execute compression, e.g., sequentially or in parallel. The execution is determined by the policy and any policy-specific parameters such as number of threads. :: typedef struct { zfp_exec_policy policy; // execution policy (serial, omp, cuda, ...) void* params; // execution parameters } zfp_execution; .. warning:: As of |zfp| |verrelease| :c:type:`zfp_execution` replaces the former :code:`zfp_exec_params` with a :code:`void*` to the associated :code:`zfp_exec_params` type (e.g., :c:type:`zfp_exec_params_omp`) to limit ABI-breaking changes due to future extensions to |zfp| execution policies. ---- .. c:type:: zfp_exec_policy Currently three execution policies are available: serial, OpenMP parallel, and CUDA parallel. :: typedef enum { zfp_exec_serial = 0, // serial execution (default) zfp_exec_omp = 1, // OpenMP multi-threaded execution zfp_exec_cuda = 2 // CUDA parallel execution } zfp_exec_policy; ---- .. c:type:: zfp_exec_params_omp Execution parameters for OpenMP parallel compression. These are initialized to default values. When nonzero, they indicate the number of threads to request for parallel compression and the number of consecutive blocks to assign to each thread. :: typedef struct { uint threads; // number of requested threads uint chunk_size; // number of blocks per chunk } zfp_exec_params_omp; ---- .. _mode_struct: .. c:type:: zfp_mode Enumerates the compression modes. :: typedef enum { zfp_mode_null = 0, // an invalid configuration of the 4 params zfp_mode_expert = 1, // expert mode (4 params set manually) zfp_mode_fixed_rate = 2, // fixed rate mode zfp_mode_fixed_precision = 3, // fixed precision mode zfp_mode_fixed_accuracy = 4, // fixed accuracy mode zfp_mode_reversible = 5 // reversible (lossless) mode } zfp_mode; ---- .. _config_struct: .. c:type:: zfp_config Encapsulates compression mode and parameters (if any). :: typedef struct { zfp_mode mode; // compression mode */ union { double rate; // compressed bits/value (negative for word alignment) uint precision; // uncompressed bits/value double tolerance; // absolute error tolerance struct { uint minbits; // min number of compressed bits/block uint maxbits; // max number of compressed bits/block uint maxprec; // max number of uncompressed bits/value int minexp; // min floating point bit plane number to store } expert; // expert mode arguments } arg; // arguments corresponding to compression mode } zfp_config; ---- .. c:type:: zfp_type Enumerates the scalar types supported by the compressor and describes the uncompressed array. The compressor and decompressor must use the same :c:type:`zfp_type`, e.g., one cannot compress doubles and decompress to floats or integers. :: typedef enum { zfp_type_none = 0, // unspecified type zfp_type_int32 = 1, // 32-bit signed integer zfp_type_int64 = 2, // 64-bit signed integer zfp_type_float = 3, // single precision floating point zfp_type_double = 4 // double precision floating point } zfp_type; ---- .. _field: .. index:: single: Strided Arrays .. c:type:: zfp_field The uncompressed array is described by the :c:type:`zfp_field` struct, which encodes the array's scalar type, dimensions, and memory layout. :: typedef struct { zfp_type type; // scalar type (e.g., int32, double) size_t nx, ny, nz, nw; // sizes (zero for unused dimensions) ptrdiff_t sx, sy, sz, sw; // strides (zero for contiguous array a[nw][nz][ny][nx]) void* data; // pointer to array data } zfp_field; For example, a static multidimensional C array declared as :: double array[n1][n2][n3][n4]; would be described by a :c:type:`zfp_field` with members :: type = zfp_type_double; nx = n4; ny = n3; nz = n2; nw = n1; sx = 1; sy = n4; sz = n3 * n4; sw = n2 * n3 * n4; data = &array[0][0][0][0]; The strides, when nonzero, specify how the array is laid out in memory. Strides can be used in case multiple fields are stored interleaved via "array of struct" (AoS) rather than "struct of array" (SoA) storage, or if the dimensions should be transposed during (de)compression. Strides may even be negative, allowing one or more dimensions to be traversed in reverse order. Given 4D array indices (*x*, *y*, *z*, *w*), the corresponding array element is stored at :: data[x * sx + y * sy + z * sz + w * sw] where :code:`data` is a pointer to the first array element. .. _new-field: .. warning:: The :c:type:`zfp_field` struct was modified in |zfp| |64bitrelease| to use :code:`size_t` and :code:`ptrdiff_t` for array dimensions and strides, respectively, to support 64-bit addressing of very large arrays (previously, :code:`uint` and :code:`int` were used). This ABI incompatible change may require rebuilding applications that use |zfp| and may in some cases also require code changes to handle pointers to :code:`size_t` instead of pointers to :code:`uint` (see :c:func:`zfp_field_size`, for instance). .. _indexing: .. index:: single: C order single: Fortran order .. warning:: It is paramount that the field dimensions, *nx*, *ny*, *nz*, and *nw*, and strides, *sx*, *sy*, *sz*, and *sw*, be correctly mapped to how the uncompressed array is laid out in memory. Although compression will still succeed if array dimensions are accidentally transposed, compression ratio and/or accuracy may suffer greatly. Since the leftmost index, *x*, is assumed to vary fastest, |zfp| can be thought of as assuming Fortran ordering. For C ordered arrays, the user should transpose the dimensions or specify strides to properly describe the memory layout. See :ref:`this FAQ ` for further details. .. c:type:: zfp_bool :c:type:`zfp_bool` is new as of |zfp| |boolrelease|. Although merely an alias for :code:`int`, this type serves to document that a return value or function parameter should be treated as Boolean. Two enumerated constants are available:: enum { zfp_false = 0, zfp_true = !zfp_false }; The reason why :c:type:`zfp_bool` is not an enumerated type itself is that in C++ this would require an explicit cast between the :code:`bool` type resulting from logical expressions, e.g., :code:`zfp_bool done = static_cast(queue.empty() && work == 0)`. Such casts from :code:`bool` to a non-enumerated :code:`int` are not necessary. The |zfp| |boolrelease| API has changed to use :c:type:`zfp_bool` in place of :code:`int` where appropriate; this change should not affect existing code. .. _hl-data: Constants --------- .. c:var:: const uint zfp_codec_version The version of the compression CODEC implemented by this version of the |zfp| library. The library can decompress files generated by the same CODEC only. To ensure that the :file:`zfp.h` header matches the binary library linked to, :c:data:`zfp_codec_version` should match :c:macro:`ZFP_CODEC`. ---- .. c:var:: const uint zfp_library_version The library version. The binary library and headers are compatible if :c:data:`zfp_library_version` matches :c:macro:`ZFP_VERSION`. ---- .. c:var:: const char* const zfp_version_string A constant string representing the |zfp| library version and release date. One can search for this string in executables and libraries that link to |libzfp| when built as a static library. .. _hl-functions: Functions --------- .. c:function:: size_t zfp_type_size(zfp_type type) Return byte size of the given scalar type, e.g., :code:`zfp_type_size(zfp_type_float) = 4`. .. _hl-func-bitstream: Compressed Stream ^^^^^^^^^^^^^^^^^ .. c:function:: zfp_stream* zfp_stream_open(bitstream* stream) Allocate compressed stream and associate it with bit stream for reading and writing bits to/from memory. *stream* may be :c:macro:`NULL` and attached later via :c:func:`zfp_stream_set_bit_stream`. ---- .. c:function:: void zfp_stream_close(zfp_stream* stream) Close and deallocate compressed stream. This does not affect the attached bit stream. ---- .. c:function:: void zfp_stream_rewind(zfp_stream* stream) Rewind bit stream to beginning for compression or decompression. ---- .. c:function:: bitstream* zfp_stream_bit_stream(const zfp_stream* stream) Return bit stream associated with compressed stream (see :c:func:`zfp_stream_set_bit_stream`). ---- .. c:function:: void zfp_stream_set_bit_stream(zfp_stream* stream, bitstream* bs) Associate bit stream with compressed stream. ---- .. c:function:: size_t zfp_stream_compressed_size(const zfp_stream* stream) Number of bytes of compressed storage. This function returns the current byte offset within the bit stream from the beginning of the bit stream memory buffer. To ensure all buffered compressed data has been output call :c:func:`zfp_stream_flush` first. ---- .. c:function:: size_t zfp_stream_maximum_size(const zfp_stream* stream, const zfp_field* field) Conservative estimate of the compressed byte size for the compression parameters stored in *stream* and the array whose scalar type and dimensions are given by *field*. This function may be used to determine how large a memory buffer to allocate to safely hold the entire compressed array. The buffer may then be resized (using :code:`realloc()`) after the actual number of bytes is known, as returned by :c:func:`zfp_compress`. .. _hl-func-stream: Compression Parameters ^^^^^^^^^^^^^^^^^^^^^^ .. c:function:: zfp_mode zfp_stream_compression_mode(const zfp_stream* stream) Return compression mode associated with compression parameters. Return :code:`zfp_mode_null` when compression parameters are invalid. ---- .. c:function:: void zfp_stream_set_reversible(zfp_stream* stream) Enable :ref:`reversible ` (lossless) compression. ---- .. c:function:: double zfp_stream_rate(const zfp_stream* stream, uint dims) Return rate in compressed bits per value if *stream* is in :ref:`fixed-rate mode ` (see :c:func:`zfp_stream_set_rate`), else zero. *dims* is the dimensionality of the compressed data. ---- .. c:function:: double zfp_stream_set_rate(zfp_stream* stream, double rate, zfp_type type, uint dims, zfp_bool align) Set *rate* for :ref:`fixed-rate mode ` in compressed bits per value. The target scalar *type* and array *dimensionality* are needed to correctly translate the rate to the number of bits per block. The Boolean *align* should be :code:`zfp_true` if :ref:`word alignment ` is needed, e.g., to support random access writes of blocks for |zfp|'s :ref:`compressed arrays `. Such alignment may further constrain the rate. The closest supported rate is returned, which may differ from the requested rate. ---- .. c:function:: uint zfp_stream_precision(const zfp_stream* stream) Return precision in uncompressed bits per value if *stream* is in :ref:`fixed-precision mode ` (see :c:func:`zfp_stream_set_precision`), else zero. ---- .. c:function:: uint zfp_stream_set_precision(zfp_stream* stream, uint precision) Set *precision* for :ref:`fixed-precision mode `. The precision specifies how many uncompressed bits per value to store, and indirectly governs the relative error. The actual precision is returned, e.g., in case the desired precision is out of range. To preserve a certain floating-point mantissa or integer precision in the decompressed data, see FAQ :ref:`#21 `. ---- .. c:function:: double zfp_stream_accuracy(const zfp_stream* stream) Return accuracy as an absolute error tolerance if *stream* is in :ref:`fixed-accuracy mode ` (see :c:func:`zfp_stream_set_accuracy`), else zero. ---- .. c:function:: double zfp_stream_set_accuracy(zfp_stream* stream, double tolerance) Set absolute error *tolerance* for :ref:`fixed-accuracy mode `. The tolerance ensures that values in the decompressed array differ from the input array by no more than this tolerance (in all but exceptional circumstances; see FAQ :ref:`#17 `). This compression mode should be used only with floating-point (not integer) data. ---- .. c:function:: uint64 zfp_stream_mode(const zfp_stream* stream) Return compact encoding of compression parameters. If the return value is no larger than :c:macro:`ZFP_MODE_SHORT_MAX`, then the least significant :c:macro:`ZFP_MODE_SHORT_BITS` (12 in the current version) suffice to encode the parameters. Otherwise all 64 bits are needed, and the low :c:macro:`ZFP_MODE_SHORT_BITS` bits will be all ones. Thus, this variable-length encoding can be used to economically encode and decode the compression parameters, which is especially important if the parameters are to vary spatially over small regions. Such spatially adaptive coding would have to be implemented via the :ref:`low-level API `. ---- .. c:function:: zfp_mode zfp_stream_set_mode(zfp_stream* stream, uint64 mode) Set all compression parameters from compact integer representation. See :c:func:`zfp_stream_mode` for how to encode the parameters. Return the mode associated with the newly-set compression parameters. If the decoded compression parameters are invalid, they are not set and the function returns :code:`zfp_mode_null`. ---- .. c:function:: void zfp_stream_params(const zfp_stream* stream, uint* minbits, uint* maxbits, uint* maxprec, int* minexp) Query :ref:`compression parameters `. For any parameter not needed, pass :c:macro:`NULL` for the corresponding pointer. ---- .. c:function:: zfp_bool zfp_stream_set_params(zfp_stream* stream, uint minbits, uint maxbits, uint maxprec, int minexp) Set all compression parameters directly. See the section on :ref:`expert mode ` for a discussion of the parameters. The return value is :code:`zfp_true` upon success. .. _hl-func-exec: Execution Policy ^^^^^^^^^^^^^^^^ .. c:function:: zfp_exec_policy zfp_stream_execution(const zfp_stream* stream) Return current :ref:`execution policy `. ---- .. c:function:: uint zfp_stream_omp_threads(const zfp_stream* stream) Return number of OpenMP threads to request for compression. See :c:func:`zfp_stream_set_omp_threads`. ---- .. c:function:: uint zfp_stream_omp_chunk_size(const zfp_stream* stream) Return number of blocks to compress together per OpenMP thread. See :c:func:`zfp_stream_set_omp_chunk_size`. ---- .. c:function:: zfp_bool zfp_stream_set_execution(zfp_stream* stream, zfp_exec_policy policy) Set :ref:`execution policy `. If different from the previous policy, initialize the execution parameters to their default values. :code:`zfp_true` is returned if the execution policy is supported. ---- .. c:function:: zfp_bool zfp_stream_set_omp_threads(zfp_stream* stream, uint threads) Set the number of OpenMP threads to use during compression. If *threads* is zero, then the number of threads is given by the value of the OpenMP *nthreads-var* internal control variable when :c:func:`zfp_compress` is called (usually the maximum number available). This function also sets the execution policy to OpenMP. Upon success, :code:`zfp_true` is returned. ---- .. c:function:: zfp_bool zfp_stream_set_omp_chunk_size(zfp_stream* stream, uint chunk_size) Set the number of consecutive blocks to compress together per OpenMP thread. If zero, use one chunk per thread. This function also sets the execution policy to OpenMP. Upon success, :code:`zfp_true` is returned. .. _hl-func-config: Compression Configuration ^^^^^^^^^^^^^^^^^^^^^^^^^ These functions encode a desired compression mode and associated parameters (if any) in a single struct, e.g., for configuring |zfp|'s :ref:`read-only array classes `. .. c:function:: zfp_config zfp_config_none() Unspecified configuration. ---- .. c:function:: zfp_config zfp_config_rate(double rate, zfp_bool align) :ref:`Fixed-rate mode ` using *rate* compressed bits per value. When *align* is true, word alignment is enforced to further constrain the rate (see :c:func:`zfp_stream_set_rate`). ---- .. c:function:: zfp_config zfp_config_precision(uint precision) :ref:`Fixed-precision mode ` using *precision* uncompressed bits per value (see also :c:func:`zfp_stream_set_precision`). ---- .. c:function:: zfp_config zfp_config_accuracy(double tolerance) :ref:`Fixed-accuracy mode ` with absolute error no larger than *tolerance* (see also :c:func:`zfp_stream_set_accuracy`). ---- .. c:function:: zfp_config zfp_config_reversible() :ref:`Reversible (lossless) mode ` (see also :c:func:`zfp_stream_set_reversible`). ---- .. c:function:: zfp_config zfp_config_expert(uint minbits, uint maxbits, uint maxprec, int minexp) :ref:`Expert mode ` with given parameters (see also :c:func:`zfp_stream_set_params`). .. _hl-func-field: Array Metadata ^^^^^^^^^^^^^^ .. c:function:: zfp_field* zfp_field_alloc() Allocates and returns a default initialized :c:type:`zfp_field` struct. The caller must free this struct using :c:func:`zfp_field_free`. ---- .. c:function:: zfp_field* zfp_field_1d(void* pointer, zfp_type type, size_t nx) Allocate and return a field struct that describes an existing 1D array, :code:`a[nx]`, of *nx* uncompressed scalars of given *type* stored at *pointer*, which may be :c:macro:`NULL` and specified later. ---- .. c:function:: zfp_field* zfp_field_2d(void* pointer, zfp_type type, size_t nx, size_t ny) Allocate and return a field struct that describes an existing 2D array, :code:`a[ny][nx]`, of *nx* |times| *ny* uncompressed scalars of given *type* stored at *pointer*, which may be :c:macro:`NULL` and specified later. ---- .. c:function:: zfp_field* zfp_field_3d(void* pointer, zfp_type type, size_t nx, size_t ny, size_t nz) Allocate and return a field struct that describes an existing 3D array, :code:`a[nz][ny][nx]`, of *nx* |times| *ny* |times| *nz* uncompressed scalars of given *type* stored at *pointer*, which may be :c:macro:`NULL` and specified later. ---- .. c:function:: zfp_field* zfp_field_4d(void* pointer, zfp_type type, size_t nx, size_t ny, size_t nz, size_t nw) Allocate and return a field struct that describes an existing 4D array, :code:`a[nw][nz][ny][nx]`, of *nx* |times| *ny* |times| *nz* |times| *nw* uncompressed scalars of given *type* stored at *pointer*, which may be :c:macro:`NULL` and specified later. ---- .. c:function:: void zfp_field_free(zfp_field* field) Free :c:type:`zfp_field` struct previously allocated by one of the functions above. ---- .. c:function:: void* zfp_field_pointer(const zfp_field* field) Return pointer to the first scalar in the field with index *x* = *y* = *z* = *w* = 0. ---- .. c:function:: void* zfp_field_begin(const zfp_field* field) Return pointer to the lowest memory address occupied by the field. Equals :c:func:`zfp_field_pointer` if all strides are positive. Available since |zfp| |fieldrelease|. ---- .. c:function:: zfp_type zfp_field_type(const zfp_field* field) Return array scalar type. ---- .. c:function:: uint zfp_field_precision(const zfp_field* field) Return scalar precision in number of bits, e.g., 32 for :code:`zfp_type_float`. ---- .. c:function:: uint zfp_field_dimensionality(const zfp_field* field) Return array dimensionality (1, 2, 3, or 4). ---- .. c:function:: size_t zfp_field_size(const zfp_field* field, size_t* size) Return total number of scalars stored in the array, e.g., *nx* |times| *ny* |times| *nz* for a 3D array. If *size* is not :c:macro:`NULL`, then store the number of scalars for each dimension, e.g., :code:`size[0] = nx; size[1] = ny; size[2] = nz` for a 3D array. ---- .. c:function:: size_t zfp_field_size_bytes(const zfp_field* field) Return number of bytes spanned by the field payload data. This includes gaps in memory in case the field layout, as given by the strides, is not contiguous (see :c:func:`zfp_field_is_contiguous`). Available since |zfp| |fieldrelease|. ---- .. c:function:: size_t zfp_field_blocks(const zfp_field* field) Return total number of *d*-dimensional blocks (whether partial or whole) spanning the array. Each whole block consists of |4powd| scalars. Available since |zfp| |fieldrelease|. ---- .. c:function:: zfp_bool zfp_field_stride(const zfp_field* field, ptrdiff_t* stride) Return :code:`zfp_false` if the array is stored contiguously as :code:`a[nx]`, :code:`a[ny][nx]`, :code:`a[nz][ny][nx]`, or :code:`a[nw][nz][ny][nx]` depending on dimensionality. Return :code:`zfp_true` if the array is strided and laid out differently in memory. If *stride* is not :c:macro:`NULL`, then store the stride for each dimension, e.g., :code:`stride[0] = sx; stride[1] = sy; stride[2] = sz;` for a 3D array. See :c:type:`zfp_field` for more information on strides. Return false if the array is stored contiguously (the default) as :code:`a[nx]`, :code:`a[ny][nx]`, :code:`a[nz][ny][nx]`, or :code:`a[nw][nz][ny][nx]` depending on dimensionality. Return true if nonzero strides have been specified. ---- .. c:function:: zfp_bool zfp_field_is_contiguous(const zfp_field* field) Return true if the field occupies a contiguous portion of memory. Note that the field layout may be contiguous even if a raster order traversal does not visit memory in a monotonically increasing or decreasing order, e.g., if the layout is simply a permutation of the default layout. Available since |zfp| |fieldrelease|. ---- .. c:function:: uint64 zfp_field_metadata(const zfp_field* field) Return 52-bit compact encoding of the scalar type and array dimensions. This function returns :c:macro:`ZFP_META_NULL` on failure, e.g., if the array dimensions are :ref:`too large ` to be encoded in 52 bits. ---- .. _zfp_field_set: .. c:function:: void zfp_field_set_pointer(zfp_field* field, void* pointer) Set pointer to first scalar in the array. ---- .. c:function:: zfp_type zfp_field_set_type(zfp_field* field, zfp_type type) Set array scalar type. ---- .. c:function:: void zfp_field_set_size_1d(zfp_field* field, size_t nx) Specify dimensions of 1D array :code:`a[nx]`. ---- .. c:function:: void zfp_field_set_size_2d(zfp_field* field, size_t nx, size_t ny) Specify dimensions of 2D array :code:`a[ny][nx]`. ---- .. c:function:: void zfp_field_set_size_3d(zfp_field* field, size_t nx, size_t ny, size_t nz) Specify dimensions of 3D array :code:`a[nz][ny][nx]`. ---- .. c:function:: void zfp_field_set_size_4d(zfp_field* field, size_t nx, size_t ny, size_t nz, size_t nw) Specify dimensions of 4D array :code:`a[nw][nz][ny][nx]`. ---- .. c:function:: void zfp_field_set_stride_1d(zfp_field* field, ptrdiff_t sx) Specify stride for 1D array: :code:`sx = &a[1] - &a[0]`. ---- .. c:function:: void zfp_field_set_stride_2d(zfp_field* field, ptrdiff_t sx, ptrdiff_t sy) Specify strides for 2D array: :code:`sx = &a[0][1] - &a[0][0]; sy = &a[1][0] - &a[0][0]`. ---- .. c:function:: void zfp_field_set_stride_3d(zfp_field* field, ptrdiff_t sx, ptrdiff_t sy, ptrdiff_t sz) Specify strides for 3D array: :code:`sx = &a[0][0][1] - &a[0][0][0]; sy = &a[0][1][0] - &a[0][0][0]; sz = &a[1][0][0] - &a[0][0][0]`. ---- .. c:function:: void zfp_field_set_stride_4d(zfp_field* field, ptrdiff_t sx, ptrdiff_t sy, ptrdiff_t sz, ptrdiff_t sw) Specify strides for 4D array: :code:`sx = &a[0][0][0][1] - &a[0][0][0][0]; sy = &a[0][0][1][0] - &a[0][0][0][0]; sz = &a[0][1][0][0] - &a[0][0][0][0]; sw = &a[1][0][0][0] - &a[0][0][0][0]`. ---- .. c:function:: zfp_bool zfp_field_set_metadata(zfp_field* field, uint64 meta) Specify array scalar type and dimensions from compact 52-bit representation. Return :code:`zfp_true` upon success. See :c:func:`zfp_field_metadata` for how to encode *meta*. .. _hl-func-codec: Compression and Decompression ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ .. c:function:: size_t zfp_compress(zfp_stream* stream, const zfp_field* field) Compress the whole array described by *field* using parameters given by *stream*. Then flush the stream to emit any buffered bits and align the stream on a word boundary. The resulting byte offset within the bit stream is returned, which equals the total number of bytes of compressed storage if the stream was rewound before the :c:func:`zfp_compress` call. Zero is returned if compression failed. ---- .. c:function:: size_t zfp_decompress(zfp_stream* stream, zfp_field* field) Decompress from *stream* to array described by *field* and align the stream on the next word boundary. Upon success, the nonzero return value is the same as would be returned by a corresponding :c:func:`zfp_compress` call, i.e., the current byte offset or the number of compressed bytes consumed. Zero is returned if decompression failed. ---- .. _zfp-header: .. c:function:: size_t zfp_write_header(zfp_stream* stream, const zfp_field* field, uint mask) Write an optional variable-length header to the stream that encodes compression parameters, array metadata, etc. The header information written is determined by the bit *mask* (see :c:macro:`macros `). Unlike in :c:func:`zfp_compress`, no word alignment is enforced. See the :ref:`limitations ` section for limits on the maximum array size supported by the header. The return value is the number of bits written, or zero upon failure. ---- .. c:function:: size_t zfp_read_header(zfp_stream* stream, zfp_field* field, uint mask) Read header if one was previously written using :c:func:`zfp_write_header`. The *stream* and *field* data structures are populated with the information stored in the header, as specified by the bit *mask* (see :c:macro:`macros `). The caller must ensure that *mask* agrees between header read and write calls. The return value is the number of bits read, or zero upon failure.