From 782283a3d067dd5a856de4fda91030e3394cb598 Mon Sep 17 00:00:00 2001 From: Sam Hocevar Date: Mon, 29 Aug 2011 23:26:58 +0000 Subject: [PATCH] test: add matrix code to the benchmark program. --- test/lol-bench.cpp | 178 ++++++++++++++++++++++++++++++++++----------- 1 file changed, 135 insertions(+), 43 deletions(-) diff --git a/test/lol-bench.cpp b/test/lol-bench.cpp index d3dc480e..5f95fd19 100644 --- a/test/lol-bench.cpp +++ b/test/lol-bench.cpp @@ -20,75 +20,154 @@ using namespace std; using namespace lol; +static size_t const MATRIX_TABLE_SIZE = 64 * 1024; +static size_t const MATRIX_RUNS = 100; + static size_t const HALF_TABLE_SIZE = 1024 * 1024; -static size_t const HALF_RUNS = 100; +static size_t const HALF_RUNS = 50; +static void bench_matrix(int mode); static void bench_half(int mode); int main(int argc, char **argv) { - Log::Info("-----------------------------------\n"); - Log::Info("Half precision floats (random bits)\n"); - Log::Info("-----------------------------------\n"); + Log::Info("----------------------------\n"); + Log::Info(" Float matrices [-2.0, 2.0]\n"); + Log::Info("----------------------------\n"); + bench_matrix(1); + + Log::Info("-------------------------------------\n"); + Log::Info(" Half precision floats (random bits)\n"); + Log::Info("-------------------------------------\n"); bench_half(1); - Log::Info("---------------------------------\n"); - Log::Info("Half precision floats [-2.0, 2.0]\n"); - Log::Info("---------------------------------\n"); + Log::Info("-----------------------------------\n"); + Log::Info(" Half precision floats [-2.0, 2.0]\n"); + Log::Info("-----------------------------------\n"); bench_half(2); return EXIT_SUCCESS; } -static void bench_half(int mode) +static void bench_matrix(int mode) { - float result[8] = { 0.0f }; + float result[5] = { 0.0f }; Timer timer; /* Set up tables */ - float *pf = new float[HALF_TABLE_SIZE]; - half *ph = new half[HALF_TABLE_SIZE]; + mat4 *pm = new mat4[MATRIX_TABLE_SIZE + 1]; + float *pf = new float[MATRIX_TABLE_SIZE]; - switch (mode) + for (size_t run = 0; run < MATRIX_RUNS; run++) { - case 1: - for (size_t i = 0; i < HALF_TABLE_SIZE; i++) - ph[i] = half::makebits(rand()); - break; - case 2: - for (size_t i = 0; i < HALF_TABLE_SIZE; i++) - ph[i] = RandF(-2.0f, 2.0f); - break; + switch (mode) + { + case 1: + for (size_t i = 0; i < MATRIX_TABLE_SIZE; i++) + for (int j = 0; j < 4; j++) + for (int k = 0; k < 4; k++) + pm[i][j][k] = RandF(-2.0f, 2.0f); + break; + } + + /* Copy matrices */ + timer.GetMs(); + for (size_t i = 0; i < MATRIX_TABLE_SIZE; i++) + pm[i] = pm[i + 1]; + result[0] += timer.GetMs(); + + /* Determinant */ + timer.GetMs(); + for (size_t i = 0; i < MATRIX_TABLE_SIZE; i++) + pf[i] = pm[i].det(); + result[1] += timer.GetMs(); + + /* Multiply matrices */ + timer.GetMs(); + for (size_t i = 0; i < MATRIX_TABLE_SIZE; i++) + pm[i] *= pm[i + 1]; + result[2] += timer.GetMs(); + + /* Add matrices */ + timer.GetMs(); + for (size_t i = 0; i < MATRIX_TABLE_SIZE; i++) + pm[i] += pm[i + 1]; + result[3] += timer.GetMs(); + + /* Invert matrix */ + timer.GetMs(); + for (size_t i = 0; i < MATRIX_TABLE_SIZE; i++) + pm[i] = pm[i].invert(); + result[4] += timer.GetMs(); } + delete[] pm; + delete[] pf; + + for (size_t i = 0; i < sizeof(result) / sizeof(*result); i++) + result[i] *= 1000000.0f / (MATRIX_TABLE_SIZE * MATRIX_RUNS); + + Log::Info(" ns/elem\n"); + Log::Info("mat4 = mat4 %7.3f\n", result[0]); + Log::Info("float = mat4.det() %7.3f\n", result[1]); + Log::Info("mat4 *= mat4 %7.3f\n", result[2]); + Log::Info("mat4 += mat4 %7.3f\n", result[3]); + Log::Info("mat4 = mat4.invert() %7.3f\n", result[4]); +} + +static void bench_half(int mode) +{ + float result[10] = { 0.0f }; + Timer timer; + + /* Set up tables */ + float *pf = new float[HALF_TABLE_SIZE + 1]; + half *ph = new half[HALF_TABLE_SIZE + 1]; + for (size_t run = 0; run < HALF_RUNS; run++) { + switch (mode) + { + case 1: + for (size_t i = 0; i < HALF_TABLE_SIZE + 1; i++) + ph[i] = half::makebits(rand()); + break; + case 2: + for (size_t i = 0; i < HALF_TABLE_SIZE + 1; i++) + ph[i] = RandF(-2.0f, 2.0f); + break; + } +#if defined __EXCEPTIONS +#error lol +#endif + + /* Copy float */ + timer.GetMs(); + for (size_t i = 0; i < HALF_TABLE_SIZE; i++) + pf[i] = pf[i + 1]; + result[0] += timer.GetMs(); + /* Convert half to float (array) */ timer.GetMs(); half::convert(pf, ph, HALF_TABLE_SIZE); - result[0] += timer.GetMs(); + result[1] += timer.GetMs(); /* Convert half to float (fast) */ timer.GetMs(); for (size_t i = 0; i < HALF_TABLE_SIZE; i++) pf[i] = (float)ph[i]; - result[1] += timer.GetMs(); - - /* Convert float to half (array) */ - timer.GetMs(); - half::convert(ph, pf, HALF_TABLE_SIZE); result[2] += timer.GetMs(); - /* Convert float to half (fast) */ + /* Add a half to every float */ timer.GetMs(); for (size_t i = 0; i < HALF_TABLE_SIZE; i++) - ph[i] = (half)pf[i]; + pf[i] += ph[i]; result[3] += timer.GetMs(); - /* Convert float to half (accurate) */ + /* Copy half */ timer.GetMs(); for (size_t i = 0; i < HALF_TABLE_SIZE; i++) - ph[i] = half::makeaccurate(pf[i]); + ph[i] = ph[i + 1]; result[4] += timer.GetMs(); /* Change sign of every half */ @@ -97,33 +176,46 @@ static void bench_half(int mode) ph[i] = -ph[i]; result[5] += timer.GetMs(); - /* Add a half to every float */ + /* Convert float to half (array) */ timer.GetMs(); - for (size_t i = 0; i < HALF_TABLE_SIZE; i++) - pf[i] += ph[i]; + half::convert(ph, pf, HALF_TABLE_SIZE); result[6] += timer.GetMs(); + /* Convert float to half (fast) */ + timer.GetMs(); + for (size_t i = 0; i < HALF_TABLE_SIZE; i++) + ph[i] = (half)pf[i]; + result[7] += timer.GetMs(); + + /* Convert float to half (accurate) */ + timer.GetMs(); + for (size_t i = 0; i < HALF_TABLE_SIZE; i++) + ph[i] = half::makeaccurate(pf[i]); + result[8] += timer.GetMs(); + /* Add a float to every half */ timer.GetMs(); for (size_t i = 0; i < HALF_TABLE_SIZE; i++) ph[i] += pf[i]; - result[7] += timer.GetMs(); + result[9] += timer.GetMs(); } - delete[]pf; - delete[]ph; + delete[] pf; + delete[] ph; for (size_t i = 0; i < sizeof(result) / sizeof(*result); i++) result[i] *= 1000000.0f / (HALF_TABLE_SIZE * HALF_RUNS); Log::Info(" ns/elem\n"); - Log::Info("float[] = half[] %7.3f\n", result[0]); - Log::Info("float = half %7.3f\n", result[1]); - Log::Info("half[] = float[] %7.3f\n", result[2]); - Log::Info("half = float (fast) %7.3f\n", result[3]); - Log::Info("half = float (accurate) %7.3f\n", result[4]); + Log::Info("float = float %7.3f\n", result[0]); + Log::Info("float = half (array) %7.3f\n", result[1]); + Log::Info("float = half (fast) %7.3f\n", result[2]); + Log::Info("float += half %7.3f\n", result[3]); + Log::Info("half = half %7.3f\n", result[4]); Log::Info("half = -half %7.3f\n", result[5]); - Log::Info("float += half %7.3f\n", result[6]); - Log::Info("half += float %7.3f\n", result[7]); + Log::Info("half = float (array) %7.3f\n", result[6]); + Log::Info("half = float (fast) %7.3f\n", result[7]); + Log::Info("half = float (accurate) %7.3f\n", result[8]); + Log::Info("half += float %7.3f\n", result[9]); }