star-sf

commit 8ce539a4586708eceb22dcd07b03436ea8778aa3
parent c773354e904af2919f9077c712cdd12c1d8af395
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Mon, 19 Sep 2016 12:55:08 +0200

Test the BSDF view data structure

Diffstat:
M
cmake/CMakeLists.txt
 | 
1
+
A
src/test_ssf_bsdf_view.c
 | 
94
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

2 files changed, 95 insertions(+), 0 deletions(-)
diff --git a/cmake/CMakeLists.txt b/cmake/CMakeLists.txt
@@ -94,6 +94,7 @@ if(NOT NO_TEST)
   endfunction()
 
   new_test(test_ssf_bsdf)
+  new_test(test_ssf_bsdf_view)
   new_test(test_ssf_bxdf)
   new_test(test_ssf_fresnel)
   new_test(test_ssf_fresnel_dielectric_conductor)
diff --git a/src/test_ssf_bsdf_view.c b/src/test_ssf_bsdf_view.c
@@ -0,0 +1,94 @@
+/* Copyright (C) |Meso|Star> 2016 (contact@meso-star.com)
+ *
+ * This program is free software: you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation, either version 3 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. */
+
+#include "ssf.h"
+#include "test_ssf_utils.h"
+
+#include <rsys/double3.h>
+
+int
+main(int argc, char** argv)
+{
+  const size_t NSTEPS = 10000;
+  struct ssf_bxdf* bxdf;
+  struct ssf_bxdf* bxdf2;
+  struct ssf_bsdf* bsdf;
+  struct ssf_bsdf_view* view;
+  struct mem_allocator allocator;
+  double N[3], dir[4], w[3];
+  double E, V, SE;
+  double sum, sqr_sum;
+  size_t i;
+  (void)argc, (void)argv;
+
+  mem_init_proxy_allocator(&allocator, &mem_default_allocator);
+
+  CHECK(ssf_bsdf_create(&allocator, &bsdf), RES_OK);
+  CHECK(ssf_bxdf_create(&allocator, &ssf_specular_reflection, &bxdf), RES_OK);
+  CHECK(ssf_bxdf_create(&allocator, &ssf_specular_reflection, &bxdf2), RES_OK);
+  CHECK(ssf_specular_reflection_setup(bxdf, 1.0), RES_OK);
+  CHECK(ssf_specular_reflection_setup(bxdf2, 0.6789), RES_OK);
+
+  CHECK(ssf_bsdf_view_create(NULL, NULL), RES_BAD_ARG);
+  CHECK(ssf_bsdf_view_create(bsdf, NULL), RES_BAD_ARG);
+  CHECK(ssf_bsdf_view_create(NULL, &view), RES_BAD_ARG);
+  CHECK(ssf_bsdf_view_create(bsdf, &view), RES_OK);
+
+  d3_normalize(w, d3(w, -1, -1, 0));
+  d3(N, 0.0, 1.0, 0.0);
+  CHECK(ssf_bsdf_view_sample(view, 0.5, 0.5, w, N, dir), 0.0);
+
+  CHECK(ssf_bsdf_add(bsdf, bxdf, 0.5), RES_OK);
+  CHECK(ssf_bsdf_view_sample(view, 0.5, 0.5, w, N, dir), 0.0);
+  CHECK(ssf_bsdf_view_ref_put(view), RES_OK);
+  CHECK(ssf_bsdf_view_create(bsdf, &view), RES_OK);
+  CHECK(ssf_bsdf_view_sample(view, 0.5, 0.5, w, N, dir), 1.0);
+  CHECK(ssf_specular_reflection_setup(bxdf, 0.1234), RES_OK);
+  CHECK(ssf_bsdf_view_sample(view, 0.5, 0.5, w, N, dir), 0.1234);
+
+  CHECK(ssf_bsdf_add(bsdf, bxdf2, 0.5), RES_OK);
+  FOR_EACH(i, 0, NSTEPS) {
+    CHECK(ssf_bsdf_view_sample(view, 0.5, 0.5, w, N, dir), 0.1234);
+  }
+  CHECK(ssf_bsdf_view_ref_put(view), RES_OK);
+  CHECK(ssf_bsdf_view_create(bsdf, &view), RES_OK);
+
+  sum = 0;
+  sqr_sum = 0;
+  FOR_EACH(i, 0, NSTEPS) {
+    const double weight = ssf_bsdf_view_sample(view, 0.5, 0.5, w, N, dir);
+    sum += weight;
+    sqr_sum += weight*weight;
+  }
+  E =  sum / (double)NSTEPS;
+  V = sqr_sum / (double)NSTEPS - E*E;
+  SE = sqrt(V / (double)NSTEPS);
+  CHECK(eq_eps((0.1234+0.6789) * 0.5, E, SE), 1);
+
+  CHECK(ssf_bsdf_view_ref_get(NULL), RES_BAD_ARG);
+  CHECK(ssf_bsdf_view_ref_get(view), RES_OK);
+  CHECK(ssf_bsdf_view_ref_put(NULL), RES_BAD_ARG);
+  CHECK(ssf_bsdf_view_ref_put(view), RES_OK);
+  CHECK(ssf_bsdf_view_ref_put(view), RES_OK);
+
+  CHECK(ssf_bxdf_ref_put(bxdf), RES_OK);
+  CHECK(ssf_bsdf_ref_put(bsdf), RES_OK);
+
+  check_memory_allocator(&allocator);
+  mem_shutdown_proxy_allocator(&allocator);
+  CHECK(mem_allocated_size(), 0);
+  return 0;
+}
+