commit 147763bfa562a00c043b1af31207e278c11f6636
parent e3a783fa575560abdd9259107db5fe946c4e923f
Author: Vincent Forest <vincent.forest@meso-star.com>
Date: Tue, 1 Oct 2019 13:45:18 +0200
Test the accel struct configuration
Diffstat:
2 files changed, 253 insertions(+), 0 deletions(-)
diff --git a/cmake/CMakeLists.txt b/cmake/CMakeLists.txt
@@ -130,6 +130,7 @@ if(NOT NO_TEST)
register_test(${_name} ${_name})
endfunction()
+ new_test(test_s3d_accel_struct_conf)
new_test(test_s3d_device)
new_test(test_s3d_sampler)
new_test(test_s3d_sample_sphere)
diff --git a/src/test_s3d_accel_struct_conf.c b/src/test_s3d_accel_struct_conf.c
@@ -0,0 +1,252 @@
+/* Copyright (C) 2015-2019 |Meso|Star> (contact@meso-star.com)
+ *
+ * This software is a computer program whose purpose is to describe a
+ * virtual 3D environment that can be ray-traced and sampled both robustly
+ * and efficiently.
+ *
+ * This software is governed by the CeCILL license under French law and
+ * abiding by the rules of distribution of free software. You can use,
+ * modify and/or redistribute the software under the terms of the CeCILL
+ * license as circulated by CEA, CNRS and INRIA at the following URL
+ * "http://www.cecill.info".
+ *
+ * As a counterpart to the access to the source code and rights to copy,
+ * modify and redistribute granted by the license, users are provided only
+ * with a limited warranty and the software's author, the holder of the
+ * economic rights, and the successive licensors have only limited
+ * liability.
+ *
+ * In this respect, the user's attention is drawn to the risks associated
+ * with loading, using, modifying and/or developing or reproducing the
+ * software by the user in light of its specific status of free software,
+ * that may mean that it is complicated to manipulate, and that also
+ * therefore means that it is reserved for developers and experienced
+ * professionals having in-depth computer knowledge. Users are therefore
+ * encouraged to load and test the software's suitability as regards their
+ * requirements in conditions enabling the security of their systems and/or
+ * data to be ensured and, more generally, to use and operate it in the
+ * same conditions as regards security.
+ *
+ * The fact that you are presently reading this means that you have had
+ * knowledge of the CeCILL license and that you accept its terms. */
+
+#include "s3d.h"
+#include "test_s3d_utils.h"
+
+#include <rsys/clock_time.h>
+#include <rsys/math.h>
+#include <string.h>
+
+/*******************************************************************************
+ * Mesh functions and data structure
+ ******************************************************************************/
+struct mesh {
+ double* pos;
+ size_t* ids;
+ size_t nverts;
+ size_t ntris;
+ struct mem_allocator* allocator;
+};
+
+static void
+mesh_init_sphere
+ (struct mesh* sphere,
+ struct mem_allocator* allocator,
+ const size_t nthetas)
+{
+ const size_t nphis = (size_t)(((double)nthetas + 0.5) * 0.5);
+ const double step_theta = 2*PI / (double)nthetas;
+ const double step_phi = PI / (double)nphis;
+ size_t itheta, iphi;
+ size_t i;
+
+ CHK(sphere && allocator && nthetas);
+ memset(sphere, 0, sizeof(*sphere));
+
+ sphere->nverts = nthetas * (nphis-1)/*#contour verts*/ + 2 /*polar verts*/;
+ sphere->ntris = 2*nthetas * (nphis-2)/*#contour tris*/ + 2*nthetas/*#polar tris*/;
+ sphere->allocator = allocator;
+
+ CHK(sphere->pos = MEM_CALLOC(allocator, sphere->nverts, sizeof(double[3])));
+ CHK(sphere->ids = MEM_CALLOC(allocator, sphere->ntris, sizeof(size_t[3])));
+
+ /* Build the contour vertices */
+ i = 0;
+ FOR_EACH(itheta, 0, nthetas) {
+ const double theta = -PI + (double)itheta * step_theta;
+ const double cos_theta = cos(theta);
+ const double sin_theta = sin(theta);
+ FOR_EACH(iphi, 0, nphis-1) {
+ const double phi = -PI*0.5 + (double)(iphi + 1) * step_phi;
+ const double cos_phi = cos(phi);
+ const double sin_phi = sin(phi);
+ sphere->pos[i++] = cos_phi * cos_theta;
+ sphere->pos[i++] = cos_phi * sin_theta;
+ sphere->pos[i++] = sin_phi;
+ }
+ }
+ /* polar vertices */
+ sphere->pos[i++] = 0.0; sphere->pos[i++] = 0.0; sphere->pos[i++] =-1.0;
+ sphere->pos[i++] = 0.0; sphere->pos[i++] = 0.0; sphere->pos[i++] = 1.0;
+ CHK(i == sphere->nverts*3);
+
+ /* Define the indices of the contour primitives */
+ i = 0;
+ FOR_EACH(itheta, 0, nthetas) {
+ const size_t itheta0 = itheta * (nphis - 1);
+ const size_t itheta1 = ((itheta + 1) % nthetas) * (nphis - 1);
+ FOR_EACH(iphi, 0, nphis-2) {
+ const size_t iphi0 = iphi + 0;
+ const size_t iphi1 = iphi + 1;
+ sphere->ids[i++] = itheta0 + iphi0; /* First triangle */
+ sphere->ids[i++] = itheta0 + iphi1;
+ sphere->ids[i++] = itheta1 + iphi0;
+ sphere->ids[i++] = itheta1 + iphi0; /* Second triangle */
+ sphere->ids[i++] = itheta0 + iphi1;
+ sphere->ids[i++] = itheta1 + iphi1;
+ }
+ }
+ /* Define the indices of the polar primitives */
+ FOR_EACH(itheta, 0, nthetas) {
+ const size_t itheta0 = itheta * (nphis - 1);
+ const size_t itheta1 = ((itheta + 1) % nthetas) * (nphis - 1);
+ sphere->ids[i++] = nthetas * (nphis - 1);
+ sphere->ids[i++] = itheta0;
+ sphere->ids[i++] = itheta1;
+ sphere->ids[i++] = nthetas * (nphis - 1) + 1;
+ sphere->ids[i++] = itheta1 + (nphis - 2);
+ sphere->ids[i++] = itheta0 + (nphis - 2);
+ }
+ CHK(i == sphere->ntris*3);
+}
+
+static void
+mesh_release(struct mesh* mesh)
+{
+ CHK(mesh);
+ MEM_RM(mesh->allocator, mesh->pos);
+ MEM_RM(mesh->allocator, mesh->ids);
+}
+
+static INLINE void
+mesh_dump(const struct mesh* mesh, FILE* stream)
+{
+ size_t i;
+ CHK(mesh && stream);
+ FOR_EACH(i, 0, mesh->nverts) {
+ fprintf(stream, "v %g %g %g\n",
+ mesh->pos[i*3+0],
+ mesh->pos[i*3+1],
+ mesh->pos[i*3+2]);
+ }
+ FOR_EACH(i, 0, mesh->ntris) {
+ fprintf(stream, "f %lu %lu %lu\n",
+ (unsigned long)mesh->ids[i*3+0]+1,
+ (unsigned long)mesh->ids[i*3+1]+1,
+ (unsigned long)mesh->ids[i*3+2]+1);
+ }
+}
+
+static void
+mesh_get_pos(const unsigned ivert, float pos[3], void* ctx)
+{
+ const struct mesh* mesh = ctx;
+ CHK(pos && ctx && ivert < mesh->nverts);
+ pos[0] = (float)mesh->pos[ivert*3+0];
+ pos[1] = (float)mesh->pos[ivert*3+1];
+ pos[2] = (float)mesh->pos[ivert*3+2];
+}
+
+static void
+mesh_get_tri(const unsigned itri, unsigned ids[3], void* ctx)
+{
+ const struct mesh* mesh = ctx;
+ CHK(ids && ctx && itri < mesh->ntris);
+ ids[0] = (unsigned)mesh->ids[itri*3+0];
+ ids[1] = (unsigned)mesh->ids[itri*3+1];
+ ids[2] = (unsigned)mesh->ids[itri*3+2];
+}
+
+/*******************************************************************************
+ * Helper functions
+ ******************************************************************************/
+static void
+time_scene_view_creation
+ (struct s3d_scene* scn,
+ const struct s3d_accel_struct_conf* cfg,
+ const char* string)
+{
+ char dump[128];
+ struct time t0, t1;
+ struct s3d_scene_view* view;
+ CHK(scn);
+
+ time_current(&t0);
+ CHK(s3d_scene_view_create2(scn, S3D_TRACE, cfg, &view) == RES_OK);
+ time_sub(&t0, time_current(&t1), &t0);
+ time_dump(&t0, TIME_ALL, NULL, dump, sizeof(dump));
+ printf("%s: %s\n", string, dump);
+ CHK(s3d_scene_view_ref_put(view) == RES_OK);
+}
+
+/*******************************************************************************
+ * Main test function
+ ******************************************************************************/
+int
+main(int argc, char** argv)
+{
+ struct mem_allocator allocator;
+ struct mesh sphere;
+ struct s3d_device* dev;
+ struct s3d_shape* shape;
+ struct s3d_scene* scn;
+ struct s3d_scene_view* view;
+ struct s3d_vertex_data vdata = S3D_VERTEX_DATA_NULL;
+ struct s3d_accel_struct_conf cfg = S3D_ACCEL_STRUCT_CONF_DEFAULT;
+ (void)argc, (void)argv;
+
+ mem_init_proxy_allocator(&allocator, &mem_default_allocator);
+
+ mesh_init_sphere(&sphere, &allocator, 256);
+ /*mesh_dump(&sphere, stdout);*/
+
+ CHK(s3d_device_create(NULL, &allocator, 1, &dev) == RES_OK);
+ CHK(s3d_scene_create(dev, &scn) == RES_OK);
+ CHK(s3d_shape_create_mesh(dev, &shape) == RES_OK);
+ CHK(s3d_scene_attach_shape(scn, shape) == RES_OK);
+
+ vdata.usage = S3D_POSITION;
+ vdata.type = S3D_FLOAT3;
+ vdata.get = mesh_get_pos;
+ CHK(s3d_mesh_setup_indexed_vertices(shape, (unsigned)sphere.ntris, mesh_get_tri,
+ (unsigned)sphere.nverts, &vdata, 1, &sphere) == RES_OK);
+
+ CHK(s3d_scene_view_create2(NULL, S3D_TRACE, NULL, &view) == RES_BAD_ARG);
+ CHK(s3d_scene_view_create2(scn, S3D_TRACE, NULL, NULL) == RES_BAD_ARG);
+
+ time_scene_view_creation(scn, NULL, "All default");
+
+ cfg.quality = S3D_ACCEL_STRUCT_QUALITY_LOW;
+ cfg.mask = S3D_ACCEL_STRUCT_FLAG_ROBUST | S3D_ACCEL_STRUCT_FLAG_DYNAMIC;
+ time_scene_view_creation(scn, &cfg, "Low quality, robust & dynamic");
+
+ cfg.quality = S3D_ACCEL_STRUCT_QUALITY_MEDIUM;
+ cfg.mask = S3D_ACCEL_STRUCT_FLAG_COMPACT;
+ time_scene_view_creation(scn, &cfg, "Medium quality, compact");
+
+ cfg.quality = S3D_ACCEL_STRUCT_QUALITY_HIGH;
+ cfg.mask = S3D_ACCEL_STRUCT_FLAG_ROBUST | S3D_ACCEL_STRUCT_FLAG_COMPACT;
+ time_scene_view_creation(scn, &cfg, "High quality, compact & robust");
+
+ CHK(s3d_shape_ref_put(shape) == RES_OK);
+ CHK(s3d_scene_ref_put(scn) == RES_OK);
+ CHK(s3d_device_ref_put(dev) == RES_OK);
+
+ mesh_release(&sphere);
+
+ check_memory_allocator(&allocator);
+ mem_shutdown_proxy_allocator(&allocator);
+ CHK(mem_allocated_size() == 0);
+ return 0;
+}
+
