blob: 37c1e9f77a0b02ec3a261513046a0b889e1b65e6 [file] [log] [blame]
Yann Collet4856a002015-01-24 01:58:16 +01001/*
Yann Collet99909862016-04-09 16:17:18 +02002 bench.c - open-source compression benchmark module
3 Copyright (C) Yann Collet 2012-2016
Yann Collet4856a002015-01-24 01:58:16 +01004
5 GPL v2 License
6
7 This program is free software; you can redistribute it and/or modify
8 it under the terms of the GNU General Public License as published by
9 the Free Software Foundation; either version 2 of the License, or
10 (at your option) any later version.
11
12 This program is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU General Public License for more details.
16
17 You should have received a copy of the GNU General Public License along
18 with this program; if not, write to the Free Software Foundation, Inc.,
19 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
20
21 You can contact the author at :
Yann Collet99909862016-04-09 16:17:18 +020022 - zstd homepage : http://www.zstd.net
Yann Collet4856a002015-01-24 01:58:16 +010023 - zstd source repository : https://github.com/Cyan4973/zstd
Yann Collet4856a002015-01-24 01:58:16 +010024*/
25
Yann Colletf3eca252015-10-22 15:31:46 +010026/* **************************************
Yann Collet4856a002015-01-24 01:58:16 +010027* Compiler Options
Yann Colletf3eca252015-10-22 15:31:46 +010028****************************************/
Yann Collet4856a002015-01-24 01:58:16 +010029/* Disable some Visual warning messages */
Yann Collet6c8b9252015-12-16 02:44:56 +010030#ifdef _MSC_VER
31# define _CRT_SECURE_NO_WARNINGS /* fopen */
32# pragma warning(disable : 4127) /* disable: C4127: conditional expression is constant */
33#endif
Yann Collet4856a002015-01-24 01:58:16 +010034
Yann Colletf3eca252015-10-22 15:31:46 +010035/* Unix Large Files support (>4GB) */
Yann Collet4856a002015-01-24 01:58:16 +010036#define _FILE_OFFSET_BITS 64
Yann Colletf3eca252015-10-22 15:31:46 +010037#if (defined(__sun__) && (!defined(__LP64__))) /* Sun Solaris 32-bits requires specific definitions */
Yann Collet4856a002015-01-24 01:58:16 +010038# define _LARGEFILE_SOURCE
Yann Colletf3eca252015-10-22 15:31:46 +010039#elif ! defined(__LP64__) /* No point defining Large file for 64 bit */
Yann Collet4856a002015-01-24 01:58:16 +010040# define _LARGEFILE64_SOURCE
41#endif
42
Yann Collet4856a002015-01-24 01:58:16 +010043
Yann Colletf3eca252015-10-22 15:31:46 +010044/* *************************************
Yann Collet4856a002015-01-24 01:58:16 +010045* Includes
Yann Colletf3eca252015-10-22 15:31:46 +010046***************************************/
Yann Collet99909862016-04-09 16:17:18 +020047#define _POSIX_C_SOURCE 199309L /* before <time.h> - needed for nanosleep() */
Yann Collet4856a002015-01-24 01:58:16 +010048#include <stdlib.h> /* malloc, free */
49#include <string.h> /* memset */
Yann Colleteeb8ba12015-10-22 16:55:40 +010050#include <stdio.h> /* fprintf, fopen, ftello64 */
51#include <sys/types.h> /* stat64 */
52#include <sys/stat.h> /* stat64 */
Yann Collet99909862016-04-09 16:17:18 +020053#include <time.h> /* clock_t, nanosleep, clock, CLOCKS_PER_SEC */
Yann Collet4856a002015-01-24 01:58:16 +010054
Yann Collet27d3dad2016-03-11 13:41:20 +010055/* sleep : posix - windows - others */
56#if !defined(_WIN32) && (defined(__unix__) || defined(__unix) || (defined(__APPLE__) && defined(__MACH__)))
57# include <unistd.h>
inikep472638c2016-03-23 12:28:28 +010058# include <sys/resource.h> /* setpriority */
Yann Collet27d3dad2016-03-11 13:41:20 +010059# define BMK_sleep(s) sleep(s)
inikep4c12f232016-03-29 14:52:13 +020060# define mili_sleep(mili) { struct timespec t; t.tv_sec=0; t.tv_nsec=mili*1000000ULL; nanosleep(&t, NULL); }
inikep19bd48f2016-04-04 12:10:00 +020061# define SET_HIGH_PRIORITY setpriority(PRIO_PROCESS, 0, -20)
Yann Collet27d3dad2016-03-11 13:41:20 +010062#elif defined(_WIN32)
63# include <windows.h>
64# define BMK_sleep(s) Sleep(1000*s)
inikep472638c2016-03-23 12:28:28 +010065# define mili_sleep(mili) Sleep(mili)
inikep19bd48f2016-04-04 12:10:00 +020066# define SET_HIGH_PRIORITY SetPriorityClass(GetCurrentProcess(), REALTIME_PRIORITY_CLASS)
Yann Collet4856a002015-01-24 01:58:16 +010067#else
Yann Collet27d3dad2016-03-11 13:41:20 +010068# define BMK_sleep(s) /* disabled */
inikep472638c2016-03-23 12:28:28 +010069# define mili_sleep(mili) /* disabled */
inikep19bd48f2016-04-04 12:10:00 +020070# define SET_HIGH_PRIORITY /* disabled */
Yann Collet4856a002015-01-24 01:58:16 +010071#endif
72
inikep4c12f232016-03-29 14:52:13 +020073#if !defined(_WIN32) && (defined(__unix__) || defined(__unix) || (defined(__APPLE__) && defined(__MACH__)))
inikep4c12f232016-03-29 14:52:13 +020074 typedef clock_t BMK_time_t;
inikep1c556a32016-03-30 10:59:48 +020075# define BMK_initTimer(ticksPerSecond) ticksPerSecond=0
inikep4c12f232016-03-29 14:52:13 +020076# define BMK_getTime(x) x = clock()
77# define BMK_getSpanTimeMicro(ticksPerSecond, clockStart, clockEnd) (1000000ULL * (clockEnd - clockStart) / CLOCKS_PER_SEC)
78# define BMK_getSpanTimeNano(ticksPerSecond, clockStart, clockEnd) (1000000000ULL * (clockEnd - clockStart) / CLOCKS_PER_SEC)
inikep4c12f232016-03-29 14:52:13 +020079#elif defined(_WIN32)
80 typedef LARGE_INTEGER BMK_time_t;
inikep4c12f232016-03-29 14:52:13 +020081# define BMK_initTimer(x) if (!QueryPerformanceFrequency(&x)) { fprintf(stderr, "ERROR: QueryPerformance not present\n"); }
82# define BMK_getTime(x) QueryPerformanceCounter(&x)
83# define BMK_getSpanTimeMicro(ticksPerSecond, clockStart, clockEnd) (1000000ULL*(clockEnd.QuadPart - clockStart.QuadPart)/ticksPerSecond.QuadPart)
84# define BMK_getSpanTimeNano(ticksPerSecond, clockStart, clockEnd) (1000000000ULL*(clockEnd.QuadPart - clockStart.QuadPart)/ticksPerSecond.QuadPart)
85#else
86 typedef int BMK_time_t;
inikep1c556a32016-03-30 10:59:48 +020087# define BMK_initTimer(ticksPerSecond) ticksPerSecond=0
inikep4c12f232016-03-29 14:52:13 +020088# define BMK_getTimeMicro(clockStart) clockStart=1
89# define BMK_getSpanTimeMicro(ticksPerSecond, clockStart, clockEnd) (TIMELOOP_S*1000000ULL+clockEnd-clockStart)
90# define BMK_getSpanTimeNano(ticksPerSecond, clockStart, clockEnd) (TIMELOOP_S*1000000000ULL+clockEnd-clockStart)
91#endif
92
Yann Colletf3eca252015-10-22 15:31:46 +010093#include "mem.h"
Yann Collet31683c02015-12-18 01:26:48 +010094#include "zstd_static.h"
inikepeaba91a2016-03-23 20:30:26 +010095#include "zstd_internal.h" /* ZSTD_compressBegin_targetSrcSize */
96#include "datagen.h" /* RDG_genBuffer */
Yann Collet4856a002015-01-24 01:58:16 +010097#include "xxhash.h"
98
99
Yann Colletf3eca252015-10-22 15:31:46 +0100100/* *************************************
Yann Collet4856a002015-01-24 01:58:16 +0100101* Compiler specifics
Yann Colletf3eca252015-10-22 15:31:46 +0100102***************************************/
Yann Collet4856a002015-01-24 01:58:16 +0100103#if !defined(S_ISREG)
104# define S_ISREG(x) (((x) & S_IFMT) == S_IFREG)
105#endif
106
Yann Collet99909862016-04-09 16:17:18 +0200107#if defined(_MSC_VER)
108# define snprintf sprintf_s
109#elif defined (__cplusplus) || (defined (__STDC_VERSION__) && (__STDC_VERSION__ >= 199901L))
110 /* part of <stdio.h> */
111#else
112 extern int snprintf (char* s, size_t maxlen, const char* format, ...); /* not declared in <stdio.h> when C version < c99 */
Yann Colleted699e62015-12-16 02:37:24 +0100113#endif
114
Yann Collet4856a002015-01-24 01:58:16 +0100115
Yann Colletf3eca252015-10-22 15:31:46 +0100116/* *************************************
Yann Collet4856a002015-01-24 01:58:16 +0100117* Constants
Yann Colletf3eca252015-10-22 15:31:46 +0100118***************************************/
inikep44af12d2016-03-14 15:59:04 +0100119#ifndef ZSTD_VERSION
120# define ZSTD_VERSION ""
121#endif
122
Yann Collet27d3dad2016-03-11 13:41:20 +0100123#define NBLOOPS 3
Yann Collet699b14d2016-03-17 19:37:33 +0100124#define TIMELOOP_S 1
Yann Collet27d3dad2016-03-11 13:41:20 +0100125#define ACTIVEPERIOD_S 70
126#define COOLPERIOD_S 10
Yann Collet4856a002015-01-24 01:58:16 +0100127
128#define KB *(1 <<10)
129#define MB *(1 <<20)
130#define GB *(1U<<30)
131
Yann Colletd062f132015-12-01 01:31:17 +0100132static const size_t maxMemory = (sizeof(size_t)==4) ? (2 GB - 64 MB) : (size_t)(1ULL << ((sizeof(size_t)*8)-31));
Yann Collet4856a002015-01-24 01:58:16 +0100133
134static U32 g_compressibilityDefault = 50;
Yann Collet4856a002015-01-24 01:58:16 +0100135
136
Yann Colletf3eca252015-10-22 15:31:46 +0100137/* *************************************
Yann Colleted699e62015-12-16 02:37:24 +0100138* console display
Yann Colletf3eca252015-10-22 15:31:46 +0100139***************************************/
Yann Colleted699e62015-12-16 02:37:24 +0100140#define DISPLAY(...) fprintf(stderr, __VA_ARGS__)
141#define DISPLAYLEVEL(l, ...) if (g_displayLevel>=l) { DISPLAY(__VA_ARGS__); }
142static U32 g_displayLevel = 2; /* 0 : no display; 1: errors; 2 : + result + interaction + warnings; 3 : + progression; 4 : + information */
143
144
145/* *************************************
146* Exceptions
147***************************************/
148#ifndef DEBUG
149# define DEBUG 0
150#endif
151#define DEBUGOUTPUT(...) if (DEBUG) DISPLAY(__VA_ARGS__);
152#define EXM_THROW(error, ...) \
153{ \
154 DEBUGOUTPUT("Error defined at %s, line %i : \n", __FILE__, __LINE__); \
155 DISPLAYLEVEL(1, "Error %i : ", error); \
156 DISPLAYLEVEL(1, __VA_ARGS__); \
157 DISPLAYLEVEL(1, "\n"); \
158 exit(error); \
159}
Yann Collet4856a002015-01-24 01:58:16 +0100160
161
Yann Colletf3eca252015-10-22 15:31:46 +0100162/* *************************************
Yann Collet4856a002015-01-24 01:58:16 +0100163* Benchmark Parameters
Yann Colletf3eca252015-10-22 15:31:46 +0100164***************************************/
Yann Collet1d1ae402016-03-17 19:51:02 +0100165static U32 g_nbIterations = NBLOOPS;
Yann Collet1c00dc32015-10-21 08:22:25 +0100166static size_t g_blockSize = 0;
inikepeaba91a2016-03-23 20:30:26 +0100167int g_additionalParam = 0;
Yann Collet4856a002015-01-24 01:58:16 +0100168
inikep4e26bb62016-03-14 12:48:51 +0100169void BMK_setNotificationLevel(unsigned level) { g_displayLevel=level; }
170
inikepeaba91a2016-03-23 20:30:26 +0100171void BMK_setAdditionalParam(int additionalParam) { g_additionalParam=additionalParam; }
172
Yann Collet1d1ae402016-03-17 19:51:02 +0100173void BMK_SetNbIterations(unsigned nbLoops)
Yann Collet4856a002015-01-24 01:58:16 +0100174{
Yann Collet1d1ae402016-03-17 19:51:02 +0100175 g_nbIterations = nbLoops;
inikep2872b6f2016-03-22 14:38:34 +0100176 DISPLAYLEVEL(2, "- %i iterations -\n", g_nbIterations);
Yann Collet4856a002015-01-24 01:58:16 +0100177}
178
Yann Collet1c00dc32015-10-21 08:22:25 +0100179void BMK_SetBlockSize(size_t blockSize)
180{
181 g_blockSize = blockSize;
inikep4e26bb62016-03-14 12:48:51 +0100182 DISPLAYLEVEL(2, "using blocks of size %u KB \n", (U32)(blockSize>>10));
Yann Collet1c00dc32015-10-21 08:22:25 +0100183}
184
Yann Collet4856a002015-01-24 01:58:16 +0100185
Yann Colletf3eca252015-10-22 15:31:46 +0100186/* ********************************************************
Yann Collet4856a002015-01-24 01:58:16 +0100187* Private functions
Yann Colletf3eca252015-10-22 15:31:46 +0100188**********************************************************/
inikep4c12f232016-03-29 14:52:13 +0200189/* returns time span in microseconds */
190static U64 BMK_clockSpan( BMK_time_t clockStart, BMK_time_t ticksPerSecond )
Yann Collet4856a002015-01-24 01:58:16 +0100191{
inikep4c12f232016-03-29 14:52:13 +0200192 BMK_time_t clockEnd;
Yann Collet99909862016-04-09 16:17:18 +0200193
inikep4c12f232016-03-29 14:52:13 +0200194 (void)ticksPerSecond;
195 BMK_getTime(clockEnd);
196 return BMK_getSpanTimeMicro(ticksPerSecond, clockStart, clockEnd);
inikep06f793a2016-03-29 11:17:58 +0200197}
198
Yann Colleted699e62015-12-16 02:37:24 +0100199static U64 BMK_getFileSize(const char* infilename)
200{
201 int r;
202#if defined(_MSC_VER)
203 struct _stat64 statbuf;
204 r = _stat64(infilename, &statbuf);
205#else
206 struct stat statbuf;
207 r = stat(infilename, &statbuf);
208#endif
209 if (r || !S_ISREG(statbuf.st_mode)) return 0; /* No good... */
210 return (U64)statbuf.st_size;
211}
212
Yann Collet4856a002015-01-24 01:58:16 +0100213
Yann Colletf3eca252015-10-22 15:31:46 +0100214/* ********************************************************
Yann Collet4856a002015-01-24 01:58:16 +0100215* Bench functions
Yann Colletf3eca252015-10-22 15:31:46 +0100216**********************************************************/
Yann Collet1c00dc32015-10-21 08:22:25 +0100217typedef struct
Yann Collet4856a002015-01-24 01:58:16 +0100218{
Yann Colleted699e62015-12-16 02:37:24 +0100219 const char* srcPtr;
Yann Collet1c00dc32015-10-21 08:22:25 +0100220 size_t srcSize;
221 char* cPtr;
222 size_t cRoom;
223 size_t cSize;
224 char* resPtr;
225 size_t resSize;
226} blockParam_t;
227
inikep4e26bb62016-03-14 12:48:51 +0100228typedef struct
229{
inikepc034b732016-03-14 13:13:42 +0100230 double ratio;
inikep4e26bb62016-03-14 12:48:51 +0100231 size_t cSize;
inikepc034b732016-03-14 13:13:42 +0100232 double cSpeed;
233 double dSpeed;
inikep4e26bb62016-03-14 12:48:51 +0100234} benchResult_t;
235
Yann Collet99909862016-04-09 16:17:18 +0200236
Yann Colletbe2010e2015-10-31 12:57:14 +0100237#define MIN(a,b) ((a)<(b) ? (a) : (b))
Yann Collet2ce49232016-02-02 14:36:49 +0100238#define MAX(a,b) ((a)>(b) ? (a) : (b))
Yann Collet1c00dc32015-10-21 08:22:25 +0100239
Yann Colleted699e62015-12-16 02:37:24 +0100240static int BMK_benchMem(const void* srcBuffer, size_t srcSize,
inikep472638c2016-03-23 12:28:28 +0100241 const char* displayName, int cLevel,
Yann Collet31683c02015-12-18 01:26:48 +0100242 const size_t* fileSizes, U32 nbFiles,
inikep4e26bb62016-03-14 12:48:51 +0100243 const void* dictBuffer, size_t dictBufferSize, benchResult_t *result)
Yann Collet1c00dc32015-10-21 08:22:25 +0100244{
Yann Colletd64f4352016-03-21 00:07:42 +0100245 size_t const blockSize = (g_blockSize ? g_blockSize : srcSize) + (!srcSize); /* avoid div by 0 */
246 U32 const maxNbBlocks = (U32) ((srcSize + (blockSize-1)) / blockSize) + nbFiles;
Yann Colleted699e62015-12-16 02:37:24 +0100247 blockParam_t* const blockTable = (blockParam_t*) malloc(maxNbBlocks * sizeof(blockParam_t));
Yann Colletb9151402016-03-26 17:18:11 +0100248 size_t const maxCompressedSize = ZSTD_compressBound(srcSize) + (maxNbBlocks * 1024); /* add some room for safety */
Yann Collet1c00dc32015-10-21 08:22:25 +0100249 void* const compressedBuffer = malloc(maxCompressedSize);
250 void* const resultBuffer = malloc(srcSize);
Yann Collet2630a5e2016-01-14 19:13:22 +0100251 ZSTD_CCtx* refCtx = ZSTD_createCCtx();
Yann Collet31683c02015-12-18 01:26:48 +0100252 ZSTD_CCtx* ctx = ZSTD_createCCtx();
Yann Collet7b51a292016-01-26 15:58:49 +0100253 ZSTD_DCtx* refDCtx = ZSTD_createDCtx();
Yann Collet31683c02015-12-18 01:26:48 +0100254 ZSTD_DCtx* dctx = ZSTD_createDCtx();
Yann Colletde406ee2016-03-20 15:46:10 +0100255 U32 nbBlocks;
inikep4c12f232016-03-29 14:52:13 +0200256 BMK_time_t ticksPerSecond;
Yann Colletde406ee2016-03-20 15:46:10 +0100257
258 /* checks */
Yann Collet7b51a292016-01-26 15:58:49 +0100259 if (!compressedBuffer || !resultBuffer || !blockTable || !refCtx || !ctx || !refDCtx || !dctx)
Yann Colleted699e62015-12-16 02:37:24 +0100260 EXM_THROW(31, "not enough memory");
Yann Collet4856a002015-01-24 01:58:16 +0100261
Yann Collet4856a002015-01-24 01:58:16 +0100262 /* init */
263 if (strlen(displayName)>17) displayName += strlen(displayName)-17; /* can only display 17 characters */
inikep4c12f232016-03-29 14:52:13 +0200264 BMK_initTimer(ticksPerSecond);
265
Yann Collet1c00dc32015-10-21 08:22:25 +0100266 /* Init blockTable data */
Yann Colletde406ee2016-03-20 15:46:10 +0100267 { const char* srcPtr = (const char*)srcBuffer;
Yann Collet1c00dc32015-10-21 08:22:25 +0100268 char* cPtr = (char*)compressedBuffer;
269 char* resPtr = (char*)resultBuffer;
Yann Collet699b14d2016-03-17 19:37:33 +0100270 U32 fileNb;
Yann Colletde406ee2016-03-20 15:46:10 +0100271 for (nbBlocks=0, fileNb=0; fileNb<nbFiles; fileNb++) {
Yann Collet70611352015-12-16 03:01:03 +0100272 size_t remaining = fileSizes[fileNb];
Yann Collet699b14d2016-03-17 19:37:33 +0100273 U32 const nbBlocksforThisFile = (U32)((remaining + (blockSize-1)) / blockSize);
274 U32 const blockEnd = nbBlocks + nbBlocksforThisFile;
Yann Colletfd416f12016-01-30 03:14:15 +0100275 for ( ; nbBlocks<blockEnd; nbBlocks++) {
Yann Colletde406ee2016-03-20 15:46:10 +0100276 size_t const thisBlockSize = MIN(remaining, blockSize);
Yann Colleted699e62015-12-16 02:37:24 +0100277 blockTable[nbBlocks].srcPtr = srcPtr;
278 blockTable[nbBlocks].cPtr = cPtr;
279 blockTable[nbBlocks].resPtr = resPtr;
280 blockTable[nbBlocks].srcSize = thisBlockSize;
281 blockTable[nbBlocks].cRoom = ZSTD_compressBound(thisBlockSize);
282 srcPtr += thisBlockSize;
283 cPtr += blockTable[nbBlocks].cRoom;
284 resPtr += thisBlockSize;
285 remaining -= thisBlockSize;
Yann Colletfd416f12016-01-30 03:14:15 +0100286 } } }
Yann Collet1c00dc32015-10-21 08:22:25 +0100287
Yann Collet4856a002015-01-24 01:58:16 +0100288 /* warmimg up memory */
Yann Colletd062f132015-12-01 01:31:17 +0100289 RDG_genBuffer(compressedBuffer, maxCompressedSize, 0.10, 0.50, 1);
Yann Collet4856a002015-01-24 01:58:16 +0100290
291 /* Bench */
Yann Colleta5b66e32016-03-26 01:48:27 +0100292 { double fastestC = 100000000., fastestD = 100000000.;
Yann Colletb9151402016-03-26 17:18:11 +0100293 U64 const crcOrig = XXH64(srcBuffer, srcSize, 0);
Yann Collet4856a002015-01-24 01:58:16 +0100294 U64 crcCheck = 0;
inikep4c12f232016-03-29 14:52:13 +0200295 BMK_time_t coolTime;
Yann Colletde406ee2016-03-20 15:46:10 +0100296 U32 testNb;
inikep19bd48f2016-04-04 12:10:00 +0200297 size_t cSize = 0;
298 double ratio = 0.;
Yann Collet4856a002015-01-24 01:58:16 +0100299
inikep4c12f232016-03-29 14:52:13 +0200300 BMK_getTime(coolTime);
inikep4e26bb62016-03-14 12:48:51 +0100301 DISPLAYLEVEL(2, "\r%79s\r", "");
Yann Colletde406ee2016-03-20 15:46:10 +0100302 for (testNb = 1; testNb <= (g_nbIterations + !g_nbIterations); testNb++) {
inikep4c12f232016-03-29 14:52:13 +0200303 BMK_time_t clockStart, clockEnd;
inikep19bd48f2016-04-04 12:10:00 +0200304 U64 clockLoop = g_nbIterations ? TIMELOOP_S*1000000ULL : 10;
Yann Collet4856a002015-01-24 01:58:16 +0100305
Yann Collet27d3dad2016-03-11 13:41:20 +0100306 /* overheat protection */
inikep4c12f232016-03-29 14:52:13 +0200307 if (BMK_clockSpan(coolTime, ticksPerSecond) > ACTIVEPERIOD_S*1000000ULL) {
Yann Collet27d3dad2016-03-11 13:41:20 +0100308 DISPLAY("\rcooling down ... \r");
309 BMK_sleep(COOLPERIOD_S);
inikep4c12f232016-03-29 14:52:13 +0200310 BMK_getTime(coolTime);
Yann Collet27d3dad2016-03-11 13:41:20 +0100311 }
Yann Collet4856a002015-01-24 01:58:16 +0100312
313 /* Compression */
inikep2872b6f2016-03-22 14:38:34 +0100314 DISPLAYLEVEL(2, "%2i-%-17.17s :%10u ->\r", testNb, displayName, (U32)srcSize);
Yann Collet699b14d2016-03-17 19:37:33 +0100315 memset(compressedBuffer, 0xE5, maxCompressedSize); /* warm up and erase result buffer */
Yann Collet4856a002015-01-24 01:58:16 +0100316
inikep472638c2016-03-23 12:28:28 +0100317 mili_sleep(1); /* give processor time to other processes */
inikep4c12f232016-03-29 14:52:13 +0200318 BMK_getTime(clockStart);
319 do { BMK_getTime(clockEnd); }
320 while (BMK_getSpanTimeNano(ticksPerSecond, clockStart, clockEnd) == 0);
321 BMK_getTime(clockStart);
Yann Colletde406ee2016-03-20 15:46:10 +0100322
Yann Colleta5b66e32016-03-26 01:48:27 +0100323 { U32 nbLoops;
inikep19bd48f2016-04-04 12:10:00 +0200324 for (nbLoops = 0 ; BMK_clockSpan(clockStart, ticksPerSecond) < clockLoop ; nbLoops++) {
Yann Colleta5b66e32016-03-26 01:48:27 +0100325 U32 blockNb;
inikep19bd48f2016-04-04 12:10:00 +0200326 ZSTD_compressBegin_targetSrcSize(refCtx, dictBuffer, dictBufferSize, blockSize, cLevel);
327 // ZSTD_compressBegin_usingDict(refCtx, dictBuffer, dictBufferSize, cLevel);
Yann Colleta5b66e32016-03-26 01:48:27 +0100328 for (blockNb=0; blockNb<nbBlocks; blockNb++) {
329 size_t const rSize = ZSTD_compress_usingPreparedCCtx(ctx, refCtx,
330 blockTable[blockNb].cPtr, blockTable[blockNb].cRoom,
331 blockTable[blockNb].srcPtr,blockTable[blockNb].srcSize);
332 if (ZSTD_isError(rSize)) EXM_THROW(1, "ZSTD_compress_usingPreparedCCtx() failed : %s", ZSTD_getErrorName(rSize));
333 blockTable[blockNb].cSize = rSize;
334 } }
inikep19bd48f2016-04-04 12:10:00 +0200335 { U64 const clockSpan = BMK_clockSpan(clockStart, ticksPerSecond);
Yann Colleta5b66e32016-03-26 01:48:27 +0100336 if ((double)clockSpan < fastestC*nbLoops) fastestC = (double)clockSpan / nbLoops;
Yann Colletde406ee2016-03-20 15:46:10 +0100337 } }
Yann Collet4856a002015-01-24 01:58:16 +0100338
Yann Collet1c00dc32015-10-21 08:22:25 +0100339 cSize = 0;
Yann Colletde406ee2016-03-20 15:46:10 +0100340 { U32 blockNb; for (blockNb=0; blockNb<nbBlocks; blockNb++) cSize += blockTable[blockNb].cSize; }
Yann Collet2acb5d32015-10-29 16:49:43 +0100341 ratio = (double)srcSize / (double)cSize;
inikep2872b6f2016-03-22 14:38:34 +0100342 DISPLAYLEVEL(2, "%2i-%-17.17s :%10u ->%10u (%5.3f),%6.1f MB/s\r",
Yann Colletde406ee2016-03-20 15:46:10 +0100343 testNb, displayName, (U32)srcSize, (U32)cSize, ratio,
inikep06f793a2016-03-29 11:17:58 +0200344 (double)srcSize / fastestC );
Yann Collet4856a002015-01-24 01:58:16 +0100345
Yann Colleta5b66e32016-03-26 01:48:27 +0100346 (void)fastestD; (void)crcOrig; /* unused when decompression disabled */
Yann Collete93d6ce2016-01-31 00:58:06 +0100347#if 1
Yann Collet4856a002015-01-24 01:58:16 +0100348 /* Decompression */
Yann Collet7b51a292016-01-26 15:58:49 +0100349 memset(resultBuffer, 0xD6, srcSize); /* warm result buffer */
Yann Collet4856a002015-01-24 01:58:16 +0100350
inikep472638c2016-03-23 12:28:28 +0100351 mili_sleep(1); /* give processor time to other processes */
inikep4c12f232016-03-29 14:52:13 +0200352 BMK_getTime(clockStart);
353 do { BMK_getTime(clockEnd); }
354 while (BMK_getSpanTimeNano(ticksPerSecond, clockStart, clockEnd) == 0);
355 BMK_getTime(clockStart);
Yann Collet7b51a292016-01-26 15:58:49 +0100356
Yann Colleta5b66e32016-03-26 01:48:27 +0100357 { U32 nbLoops;
inikep19bd48f2016-04-04 12:10:00 +0200358 for (nbLoops = 0 ; BMK_clockSpan(clockStart, ticksPerSecond) < clockLoop ; nbLoops++) {
Yann Colleta5b66e32016-03-26 01:48:27 +0100359 U32 blockNb;
360 ZSTD_decompressBegin_usingDict(refDCtx, dictBuffer, dictBufferSize);
361 for (blockNb=0; blockNb<nbBlocks; blockNb++) {
362 size_t const regenSize = ZSTD_decompress_usingPreparedDCtx(dctx, refDCtx,
363 blockTable[blockNb].resPtr, blockTable[blockNb].srcSize,
364 blockTable[blockNb].cPtr, blockTable[blockNb].cSize);
365 if (ZSTD_isError(regenSize)) {
366 DISPLAY("ZSTD_decompress_usingPreparedDCtx() failed on block %u : %s \n",
367 blockNb, ZSTD_getErrorName(regenSize));
inikep19bd48f2016-04-04 12:10:00 +0200368 clockLoop = 0; /* force immediate test end */
Yann Colleta5b66e32016-03-26 01:48:27 +0100369 break;
370 }
371 blockTable[blockNb].resSize = regenSize;
372 } }
inikep19bd48f2016-04-04 12:10:00 +0200373 { U64 const clockSpan = BMK_clockSpan(clockStart, ticksPerSecond);
Yann Colleta5b66e32016-03-26 01:48:27 +0100374 if ((double)clockSpan < fastestD*nbLoops) fastestD = (double)clockSpan / nbLoops;
Yann Collet7b51a292016-01-26 15:58:49 +0100375 } }
376
inikep2872b6f2016-03-22 14:38:34 +0100377 DISPLAYLEVEL(2, "%2i-%-17.17s :%10u ->%10u (%5.3f),%6.1f MB/s ,%6.1f MB/s\r",
Yann Colletde406ee2016-03-20 15:46:10 +0100378 testNb, displayName, (U32)srcSize, (U32)cSize, ratio,
inikep06f793a2016-03-29 11:17:58 +0200379 (double)srcSize / fastestC,
380 (double)srcSize / fastestD );
Yann Collet4856a002015-01-24 01:58:16 +0100381
382 /* CRC Checking */
inikep19bd48f2016-04-04 12:10:00 +0200383 { crcCheck = XXH64(resultBuffer, srcSize, 0);
Yann Colleta5b66e32016-03-26 01:48:27 +0100384 if (crcOrig!=crcCheck) {
385 size_t u;
386 DISPLAY("!!! WARNING !!! %14s : Invalid Checksum : %x != %x \n", displayName, (unsigned)crcOrig, (unsigned)crcCheck);
387 for (u=0; u<srcSize; u++) {
388 if (((const BYTE*)srcBuffer)[u] != ((const BYTE*)resultBuffer)[u]) {
389 U32 segNb, bNb, pos;
390 size_t bacc = 0;
391 DISPLAY("Decoding error at pos %u ", (U32)u);
392 for (segNb = 0; segNb < nbBlocks; segNb++) {
393 if (bacc + blockTable[segNb].srcSize > u) break;
394 bacc += blockTable[segNb].srcSize;
395 }
396 pos = (U32)(u - bacc);
397 bNb = pos / (128 KB);
398 DISPLAY("(block %u, sub %u, pos %u) \n", segNb, bNb, pos);
399 break;
Yann Collet03a6dab2016-01-21 02:21:17 +0100400 }
Yann Colleta5b66e32016-03-26 01:48:27 +0100401 if (u==srcSize-1) { /* should never happen */
402 DISPLAY("no difference detected\n");
403 } }
404 break;
405 } } /* CRC Checking */
Yann Collete8c6bb12015-07-26 00:23:57 +0100406#endif
Yann Colletde406ee2016-03-20 15:46:10 +0100407 } /* for (testNb = 1; testNb <= (g_nbIterations + !g_nbIterations); testNb++) */
Yann Collet4856a002015-01-24 01:58:16 +0100408
inikep2872b6f2016-03-22 14:38:34 +0100409 if (crcOrig == crcCheck) {
inikep4e26bb62016-03-14 12:48:51 +0100410 result->ratio = ratio;
411 result->cSize = cSize;
Yann Collet99909862016-04-09 16:17:18 +0200412 result->cSpeed = (double)srcSize / fastestC;
inikep06f793a2016-03-29 11:17:58 +0200413 result->dSpeed = (double)srcSize / fastestD;
inikep4e26bb62016-03-14 12:48:51 +0100414 }
inikep2872b6f2016-03-22 14:38:34 +0100415 DISPLAYLEVEL(2, "%2i#\n", cLevel);
Yann Colletde406ee2016-03-20 15:46:10 +0100416 } /* Bench */
Yann Collet4856a002015-01-24 01:58:16 +0100417
Yann Colleted699e62015-12-16 02:37:24 +0100418 /* clean up */
Yann Collet4856a002015-01-24 01:58:16 +0100419 free(compressedBuffer);
420 free(resultBuffer);
Yann Collet2630a5e2016-01-14 19:13:22 +0100421 ZSTD_freeCCtx(refCtx);
Yann Collet31683c02015-12-18 01:26:48 +0100422 ZSTD_freeCCtx(ctx);
Yann Collet7b51a292016-01-26 15:58:49 +0100423 ZSTD_freeDCtx(refDCtx);
Yann Collet31683c02015-12-18 01:26:48 +0100424 ZSTD_freeDCtx(dctx);
Yann Collet4856a002015-01-24 01:58:16 +0100425 return 0;
426}
427
428
Yann Collet4856a002015-01-24 01:58:16 +0100429static size_t BMK_findMaxMem(U64 requiredMem)
430{
Yann Colletde406ee2016-03-20 15:46:10 +0100431 size_t const step = 64 MB;
Yann Collet4856a002015-01-24 01:58:16 +0100432 BYTE* testmem = NULL;
433
434 requiredMem = (((requiredMem >> 26) + 1) << 26);
Yann Colletde406ee2016-03-20 15:46:10 +0100435 requiredMem += step;
Yann Collet050efba2015-11-03 09:49:30 +0100436 if (requiredMem > maxMemory) requiredMem = maxMemory;
Yann Collet4856a002015-01-24 01:58:16 +0100437
Yann Colletde406ee2016-03-20 15:46:10 +0100438 do {
Yann Collet4856a002015-01-24 01:58:16 +0100439 testmem = (BYTE*)malloc((size_t)requiredMem);
Yann Colletde406ee2016-03-20 15:46:10 +0100440 requiredMem -= step;
441 } while (!testmem);
442
Yann Collet4856a002015-01-24 01:58:16 +0100443 free(testmem);
Yann Colletde406ee2016-03-20 15:46:10 +0100444 return (size_t)(requiredMem);
Yann Collet4856a002015-01-24 01:58:16 +0100445}
446
Yann Colleted699e62015-12-16 02:37:24 +0100447static void BMK_benchCLevel(void* srcBuffer, size_t benchedSize,
inikepeaba91a2016-03-23 20:30:26 +0100448 const char* displayName, int cLevel, int cLevelLast,
Yann Collet31683c02015-12-18 01:26:48 +0100449 const size_t* fileSizes, unsigned nbFiles,
450 const void* dictBuffer, size_t dictBufferSize)
Yann Collet4856a002015-01-24 01:58:16 +0100451{
inikep4e26bb62016-03-14 12:48:51 +0100452 benchResult_t result, total;
inikepe9554b72016-03-14 18:10:30 +0100453 int l;
inikep4e26bb62016-03-14 12:48:51 +0100454
inikep19bd48f2016-04-04 12:10:00 +0200455 SET_HIGH_PRIORITY;
inikep472638c2016-03-23 12:28:28 +0100456
inikepe9554b72016-03-14 18:10:30 +0100457 const char* pch = strrchr(displayName, '\\'); /* Windows */
458 if (!pch) pch = strrchr(displayName, '/'); /* Linux */
459 if (pch) displayName = pch+1;
inikep4e26bb62016-03-14 12:48:51 +0100460
inikepe9554b72016-03-14 18:10:30 +0100461 memset(&result, 0, sizeof(result));
462 memset(&total, 0, sizeof(total));
inikep5fdd0b42016-03-14 19:51:11 +0100463
inikepeaba91a2016-03-23 20:30:26 +0100464 if (g_displayLevel == 1 && !g_additionalParam)
inikep2872b6f2016-03-22 14:38:34 +0100465 DISPLAY("bench %s: input %u bytes, %i iterations, %u KB blocks\n", ZSTD_VERSION, (U32)benchedSize, g_nbIterations, (U32)(g_blockSize>>10));
inikepe9554b72016-03-14 18:10:30 +0100466
467 if (cLevelLast < cLevel) cLevelLast = cLevel;
468
Yann Collet99909862016-04-09 16:17:18 +0200469 for (l=cLevel; l <= cLevelLast; l++) {
inikepe9554b72016-03-14 18:10:30 +0100470 BMK_benchMem(srcBuffer, benchedSize,
inikep472638c2016-03-23 12:28:28 +0100471 displayName, l,
inikepe9554b72016-03-14 18:10:30 +0100472 fileSizes, nbFiles,
473 dictBuffer, dictBufferSize, &result);
474 if (g_displayLevel == 1) {
inikepeaba91a2016-03-23 20:30:26 +0100475 if (g_additionalParam)
476 DISPLAY("%-3i%11i (%5.3f) %6.1f MB/s %6.1f MB/s %s (param=%d)\n", -l, (int)result.cSize, result.ratio, result.cSpeed, result.dSpeed, displayName, g_additionalParam);
inikepd700a1a2016-03-15 12:18:44 +0100477 else
478 DISPLAY("%-3i%11i (%5.3f) %6.1f MB/s %6.1f MB/s %s\n", -l, (int)result.cSize, result.ratio, result.cSpeed, result.dSpeed, displayName);
inikepe9554b72016-03-14 18:10:30 +0100479 total.cSize += result.cSize;
480 total.cSpeed += result.cSpeed;
481 total.dSpeed += result.dSpeed;
482 total.ratio += result.ratio;
Yann Collet99909862016-04-09 16:17:18 +0200483 } }
484 if (g_displayLevel == 1 && cLevelLast > cLevel) {
inikepe9554b72016-03-14 18:10:30 +0100485 total.cSize /= 1+cLevelLast-cLevel;
486 total.cSpeed /= 1+cLevelLast-cLevel;
487 total.dSpeed /= 1+cLevelLast-cLevel;
488 total.ratio /= 1+cLevelLast-cLevel;
Yann Collet99909862016-04-09 16:17:18 +0200489 DISPLAY("avg%11i (%5.3f) %6.1f MB/s %6.1f MB/s %s\n", (int)total.cSize, total.ratio, total.cSpeed, total.dSpeed, displayName);
inikepe9554b72016-03-14 18:10:30 +0100490 }
Yann Colleted699e62015-12-16 02:37:24 +0100491}
492
493static U64 BMK_getTotalFileSize(const char** fileNamesTable, unsigned nbFiles)
494{
495 U64 total = 0;
496 unsigned n;
497 for (n=0; n<nbFiles; n++)
498 total += BMK_getFileSize(fileNamesTable[n]);
499 return total;
500}
501
Yann Colleta5b66e32016-03-26 01:48:27 +0100502/*! BMK_loadFiles() :
503 Loads `buffer` with content of files listed within `fileNamesTable`.
504 At most, fills `buffer` entirely */
Yann Collet70611352015-12-16 03:01:03 +0100505static void BMK_loadFiles(void* buffer, size_t bufferSize,
506 size_t* fileSizes,
Yann Colleta5b66e32016-03-26 01:48:27 +0100507 const char** fileNamesTable, unsigned nbFiles)
Yann Colleted699e62015-12-16 02:37:24 +0100508{
Yann Colleted699e62015-12-16 02:37:24 +0100509 size_t pos = 0;
Yann Colleted699e62015-12-16 02:37:24 +0100510
Yann Collet699b14d2016-03-17 19:37:33 +0100511 unsigned n;
Yann Colletfd416f12016-01-30 03:14:15 +0100512 for (n=0; n<nbFiles; n++) {
Yann Colleted699e62015-12-16 02:37:24 +0100513 U64 fileSize = BMK_getFileSize(fileNamesTable[n]);
Yann Colleta5b66e32016-03-26 01:48:27 +0100514 FILE* const f = fopen(fileNamesTable[n], "rb");
Yann Colleted699e62015-12-16 02:37:24 +0100515 if (f==NULL) EXM_THROW(10, "impossible to open file %s", fileNamesTable[n]);
516 DISPLAYLEVEL(2, "Loading %s... \r", fileNamesTable[n]);
Yann Colleta5b66e32016-03-26 01:48:27 +0100517 if (fileSize > bufferSize-pos) fileSize = bufferSize-pos, nbFiles=n; /* buffer too small - stop after this file */
518 { size_t const readSize = fread(((char*)buffer)+pos, 1, (size_t)fileSize, f);
519 if (readSize != (size_t)fileSize) EXM_THROW(11, "could not read %s", fileNamesTable[n]);
520 pos += readSize; }
Yann Colleta52c98d2015-12-16 03:12:31 +0100521 fileSizes[n] = (size_t)fileSize;
Yann Colleted699e62015-12-16 02:37:24 +0100522 fclose(f);
523 }
524}
525
Yann Collet31683c02015-12-18 01:26:48 +0100526static void BMK_benchFileTable(const char** fileNamesTable, unsigned nbFiles,
inikepeaba91a2016-03-23 20:30:26 +0100527 const char* dictFileName, int cLevel, int cLevelLast)
Yann Colleted699e62015-12-16 02:37:24 +0100528{
529 void* srcBuffer;
530 size_t benchedSize;
Yann Collet31683c02015-12-18 01:26:48 +0100531 void* dictBuffer = NULL;
532 size_t dictBufferSize = 0;
533 size_t* fileSizes = (size_t*)malloc(nbFiles * sizeof(size_t));
Yann Colleted699e62015-12-16 02:37:24 +0100534 U64 totalSizeToLoad = BMK_getTotalFileSize(fileNamesTable, nbFiles);
535 char mfName[20] = {0};
536 const char* displayName = NULL;
537
Yann Collet31683c02015-12-18 01:26:48 +0100538 if (!fileSizes) EXM_THROW(12, "not enough memory for fileSizes");
539
540 /* Load dictionary */
Yann Colletfd416f12016-01-30 03:14:15 +0100541 if (dictFileName != NULL) {
Yann Collet31683c02015-12-18 01:26:48 +0100542 U64 dictFileSize = BMK_getFileSize(dictFileName);
543 if (dictFileSize > 64 MB) EXM_THROW(10, "dictionary file %s too large", dictFileName);
544 dictBufferSize = (size_t)dictFileSize;
545 dictBuffer = malloc(dictBufferSize);
546 if (dictBuffer==NULL) EXM_THROW(11, "not enough memory for dictionary (%u bytes)", (U32)dictBufferSize);
547 BMK_loadFiles(dictBuffer, dictBufferSize, fileSizes, &dictFileName, 1);
548 }
549
Yann Colleted699e62015-12-16 02:37:24 +0100550 /* Memory allocation & restrictions */
551 benchedSize = BMK_findMaxMem(totalSizeToLoad * 3) / 3;
552 if ((U64)benchedSize > totalSizeToLoad) benchedSize = (size_t)totalSizeToLoad;
553 if (benchedSize < totalSizeToLoad)
554 DISPLAY("Not enough memory; testing %u MB only...\n", (U32)(benchedSize >> 20));
555 srcBuffer = malloc(benchedSize);
556 if (!srcBuffer) EXM_THROW(12, "not enough memory");
557
558 /* Load input buffer */
Yann Collet70611352015-12-16 03:01:03 +0100559 BMK_loadFiles(srcBuffer, benchedSize, fileSizes, fileNamesTable, nbFiles);
Yann Colleted699e62015-12-16 02:37:24 +0100560
561 /* Bench */
562 snprintf (mfName, sizeof(mfName), " %u files", nbFiles);
563 if (nbFiles > 1) displayName = mfName;
564 else displayName = fileNamesTable[0];
565
Yann Collet31683c02015-12-18 01:26:48 +0100566 BMK_benchCLevel(srcBuffer, benchedSize,
inikepeaba91a2016-03-23 20:30:26 +0100567 displayName, cLevel, cLevelLast,
Yann Collet31683c02015-12-18 01:26:48 +0100568 fileSizes, nbFiles,
569 dictBuffer, dictBufferSize);
Yann Collet4856a002015-01-24 01:58:16 +0100570
Yann Colleteeb8ba12015-10-22 16:55:40 +0100571 /* clean up */
Yann Collet4856a002015-01-24 01:58:16 +0100572 free(srcBuffer);
Yann Collet31683c02015-12-18 01:26:48 +0100573 free(dictBuffer);
Yann Collet70611352015-12-16 03:01:03 +0100574 free(fileSizes);
Yann Collet4856a002015-01-24 01:58:16 +0100575}
576
577
inikepeaba91a2016-03-23 20:30:26 +0100578static void BMK_syntheticTest(int cLevel, int cLevelLast, double compressibility)
Yann Collet4856a002015-01-24 01:58:16 +0100579{
Yann Colleted699e62015-12-16 02:37:24 +0100580 char name[20] = {0};
Yann Collet4856a002015-01-24 01:58:16 +0100581 size_t benchedSize = 10000000;
582 void* srcBuffer = malloc(benchedSize);
Yann Collet4856a002015-01-24 01:58:16 +0100583
Yann Collet4856a002015-01-24 01:58:16 +0100584 /* Memory allocation */
Yann Colleted699e62015-12-16 02:37:24 +0100585 if (!srcBuffer) EXM_THROW(21, "not enough memory");
Yann Collet4856a002015-01-24 01:58:16 +0100586
587 /* Fill input buffer */
Yann Colletd062f132015-12-01 01:31:17 +0100588 RDG_genBuffer(srcBuffer, benchedSize, compressibility, 0.0, 0);
Yann Collet4856a002015-01-24 01:58:16 +0100589
590 /* Bench */
Yann Colleted699e62015-12-16 02:37:24 +0100591 snprintf (name, sizeof(name), "Synthetic %2u%%", (unsigned)(compressibility*100));
inikepeaba91a2016-03-23 20:30:26 +0100592 BMK_benchCLevel(srcBuffer, benchedSize, name, cLevel, cLevelLast, &benchedSize, 1, NULL, 0);
Yann Collet4856a002015-01-24 01:58:16 +0100593
Yann Colleted699e62015-12-16 02:37:24 +0100594 /* clean up */
Yann Collet4856a002015-01-24 01:58:16 +0100595 free(srcBuffer);
Yann Collet4856a002015-01-24 01:58:16 +0100596}
597
598
Yann Collet31683c02015-12-18 01:26:48 +0100599int BMK_benchFiles(const char** fileNamesTable, unsigned nbFiles,
inikepeaba91a2016-03-23 20:30:26 +0100600 const char* dictFileName, int cLevel, int cLevelLast)
Yann Collet4856a002015-01-24 01:58:16 +0100601{
Yann Collet699b14d2016-03-17 19:37:33 +0100602 double const compressibility = (double)g_compressibilityDefault / 100;
Yann Collet4856a002015-01-24 01:58:16 +0100603
604 if (nbFiles == 0)
inikepeaba91a2016-03-23 20:30:26 +0100605 BMK_syntheticTest(cLevel, cLevelLast, compressibility);
Yann Collet4856a002015-01-24 01:58:16 +0100606 else
inikepeaba91a2016-03-23 20:30:26 +0100607 BMK_benchFileTable(fileNamesTable, nbFiles, dictFileName, cLevel, cLevelLast);
Yann Collet4856a002015-01-24 01:58:16 +0100608 return 0;
609}
610