commit 30ba1221268c7ecf0d4d7d92909dd1fbe2c59c44
parent d3e231310f17a79e3c83b9c2b992d97b2f2c683b
Author: Vincent Forest <vincent.forest@meso-star.com>
Date: Tue, 27 Jul 2021 08:44:31 +0200
Begin the tests of the this lens model
Diffstat:
2 files changed, 137 insertions(+), 0 deletions(-)
diff --git a/cmake/CMakeLists.txt b/cmake/CMakeLists.txt
@@ -82,6 +82,7 @@ if(NOT NO_TEST)
endfunction()
new_test(test_scam_pinhole)
+ new_test(test_scam_thin_lens)
if(Star3D_FOUND)
new_test(test_scam_cbox)
diff --git a/src/test_scam_thin_lens.c b/src/test_scam_thin_lens.c
@@ -0,0 +1,136 @@
+/* Copyright (C) 2021 |Meso|Star> (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/>. */
+
+#define _POSIX_C_SOURCE 200112L /* nextafterf */
+
+#include "scam.h"
+#include "test_scam_utils.h"
+
+#include <rsys/double3.h>
+#include <rsys/math.h>
+#include <rsys/mem_allocator.h>
+
+int
+main(int argc, char** argv)
+{
+ struct scam* cam = NULL;
+ struct scam_perspective_args args = SCAM_PERSPECTIVE_ARGS_DEFAULT;
+ struct scam_sample sample = SCAM_SAMPLE_NULL;
+ struct scam_ray ray = SCAM_RAY_NULL;
+ const size_t nsamps = 1000;
+ size_t i;
+ double axis_x[3];
+ double axis_y[3];
+ double axis_z[3];
+ double hori_hfov; /* horizontal half field of view */
+ double sin_hori_hfov; /* Sinus of the horizontal half field of view */
+ double sin_vert_hfov; /* Sinos of the vertical half field of view */
+ double depth;
+ double fov;
+ double length;
+ double radius;
+ double ref;
+ (void)argc, (void)argv;
+
+ args.lens_radius = -1;
+ CHK(scam_create_perspective(NULL, NULL, 1, &args, &cam) == RES_BAD_ARG);
+ args.lens_radius = 1;
+ args.focal_distance = -1;
+ CHK(scam_create_perspective(NULL, NULL, 1, &args, &cam) == RES_BAD_ARG);
+ args.focal_distance = 1;
+ CHK(scam_create_perspective(NULL, NULL, 1, &args, &cam) == RES_OK);
+ CHK(scam_ref_put(cam) == RES_OK);
+
+ length = 1.2;
+ radius = 0.01;
+ CHK(scam_focal_length_to_field_of_view(-1, length, &fov) == RES_BAD_ARG);
+ CHK(scam_focal_length_to_field_of_view(radius, 0, &fov) == RES_BAD_ARG);
+ CHK(scam_focal_length_to_field_of_view(radius, length, &fov) == RES_OK);
+
+ ref = 2*atan(radius/length);
+ CHK(eq_eps(fov, ref, ref*1.e-6));
+
+ ref = length;
+ CHK(scam_field_of_view_to_focal_length(-1, fov, &length) == RES_BAD_ARG);
+ CHK(scam_field_of_view_to_focal_length(radius, 0, &length) == RES_BAD_ARG);
+ CHK(scam_field_of_view_to_focal_length(radius, PI, &length) == RES_BAD_ARG);
+ CHK(scam_field_of_view_to_focal_length(radius, fov, &length) == RES_OK);
+
+ CHK(eq_eps(length, ref, ref*1.e-6));
+
+ args.position[0] = rand_canonical();
+ args.position[1] = rand_canonical();
+ args.position[2] = rand_canonical();
+ args.target[0] = rand_canonical();
+ args.target[1] = rand_canonical();
+ args.target[2] = rand_canonical();
+ args.up[0] = rand_canonical();
+ args.up[1] = rand_canonical();
+ args.up[2] = rand_canonical();
+ args.field_of_view = PI/2.0;
+ args.aspect_ratio = 4.0/3.0;
+ args.lens_radius = 0.01;
+ args.focal_distance = 12.34;
+
+ CHK(scam_create_perspective(NULL, NULL, 1, &args, &cam) == RES_OK);
+
+ /* Precompute some view frustum constants */
+ d3_normalize(axis_z, d3_sub(axis_z, args.target, args.position));
+ d3_normalize(axis_x, d3_cross(axis_x, axis_z, args.up));
+ d3_normalize(axis_y, d3_cross(axis_y, axis_z, axis_x));
+ depth = 1.0/tan(args.field_of_view*0.5);
+ hori_hfov = atan(args.aspect_ratio / depth);
+ sin_hori_hfov = sin(hori_hfov);
+ sin_vert_hfov = sin(args.field_of_view*0.5);
+
+ sample.lens[0] = nextafterf(0, -1);
+ CHK(scam_generate_ray(cam, &sample, &ray) == RES_BAD_ARG);
+ sample.lens[0] = 1;
+ CHK(scam_generate_ray(cam, &sample, &ray) == RES_BAD_ARG);
+ sample.lens[0] = 0;
+ sample.lens[1] = nextafterf(0, -1);
+ CHK(scam_generate_ray(cam, &sample, &ray) == RES_BAD_ARG);
+ sample.lens[1] = 1;
+ CHK(scam_generate_ray(cam, &sample, &ray) == RES_BAD_ARG);
+ sample.film[0] = sample.lens[0] = 0;
+ sample.film[1] = sample.lens[1] = nextafterf(1, 0);
+ CHK(scam_generate_ray(cam, &sample, &ray) == RES_OK);
+
+ FOR_EACH(i, 0, nsamps) {
+ double cos_ray_axis_x;
+ double cos_ray_axis_y;
+ double cos_ray_axis_z;
+ sample.film[0] = rand_canonical();
+ sample.film[1] = rand_canonical();
+ sample.lens[0] = rand_canonical();
+ sample.lens[1] = rand_canonical();
+ CHK(scam_generate_ray(cam, &sample, &ray) == RES_OK);
+
+ /* Check that the ray origin lies on the lens */
+ /*CHK(d3_eq(ray.org, args.position));*/
+
+ /* Check that the generated ray is in the view frustum */
+ /*cos_ray_axis_x = d3_dot(ray.dir, axis_x);
+ cos_ray_axis_y = d3_dot(ray.dir, axis_y);
+ cos_ray_axis_z = d3_dot(ray.dir, axis_z);
+ CHK(cos_ray_axis_z >= 0);
+ CHK(cos_ray_axis_y <= sin_vert_hfov);
+ CHK(cos_ray_axis_x <= sin_hori_hfov);*/
+ }
+
+ CHK(scam_ref_put(cam) == RES_OK);
+ CHK(mem_allocated_size() == 0);
+ return 0;
+}
