star-sf

commit f0830178e6791f9e21d5c4833ecd410839a1d668
parent 3e2fe000de4c03b4a7c2be7eaa0ada2270ddb689
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri,  5 Feb 2021 13:19:39 +0100

Begin to test the RDGFA builtin phase function

Diffstat:
M
cmake/CMakeLists.txt
 | 
1
+
A
src/test_ssf_phase_rdgfa.c
 | 
70
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

2 files changed, 71 insertions(+), 0 deletions(-)
diff --git a/cmake/CMakeLists.txt b/cmake/CMakeLists.txt
@@ -120,6 +120,7 @@ if(NOT NO_TEST)
   new_test(test_ssf_phase)
   new_test(test_ssf_phase_hg)
   new_test(test_ssf_phase_rayleigh)
+  new_test(test_ssf_phase_rdgfa)
   new_test(test_ssf_specular_dielectric_dielectric_reflection)
   new_test(test_ssf_specular_reflection)
   new_test(test_ssf_thin_specular_dielectric)
diff --git a/src/test_ssf_phase_rdgfa.c b/src/test_ssf_phase_rdgfa.c
@@ -0,0 +1,70 @@
+/* Copyright (C) 2016-2018, 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/>. */
+
+#include "ssf.h"
+#include "test_ssf_utils.h"
+
+int
+main(int argc, char** argv)
+{
+  struct mem_allocator allocator;
+  struct ssf_phase_rdgfa_setup_args args = SSF_PHASE_RDGFA_SETUP_ARGS_DEFAULT;
+  struct ssf_phase_rdgfa_desc desc = SSF_PHASE_RDGFA_DESC_NULL;
+  struct ssf_phase* phase;
+  struct ssf_phase* dummy;
+  struct ssp_rng* rng;
+
+  /* Reference values computed from a separated code */
+  const double ref_theta_limit = 0.38380861601426791;
+  const double ref_norm_factor = 1.0412632578227441;
+
+  (void)argc, (void)argv;
+
+  mem_init_proxy_allocator(&allocator, &mem_default_allocator);
+  CHK(ssp_rng_create(&allocator, &ssp_rng_mt19937_64, &rng) == RES_OK);
+
+  CHK(ssf_phase_create(&allocator, &ssf_phase_rdgfa, &phase) == RES_OK);
+  CHK(ssf_phase_create(&allocator, &phase_dummy, &dummy) == RES_OK);
+
+  args.wavelength = 633;
+  args.fractal_dimension = 1.75;
+  args.gyration_radius =  385.17247758644652;
+  CHK(ssf_phase_rdgfa_setup(NULL, &args) == RES_BAD_ARG);
+  CHK(ssf_phase_rdgfa_setup(phase, NULL) == RES_BAD_ARG);
+  CHK(ssf_phase_rdgfa_setup(dummy, &args) == RES_BAD_ARG);
+  CHK(ssf_phase_rdgfa_setup(phase, &args) == RES_OK);
+
+  CHK(ssf_phase_rdgfa_get_desc(NULL, &desc) == RES_BAD_ARG);
+  CHK(ssf_phase_rdgfa_get_desc(phase, NULL) == RES_BAD_ARG);
+  CHK(ssf_phase_rdgfa_get_desc(dummy, &desc) == RES_BAD_ARG);
+  CHK(ssf_phase_rdgfa_get_desc(phase, &desc) == RES_OK);
+
+  CHK(desc.wavelength == args.wavelength);
+  CHK(desc.fractal_dimension == args.fractal_dimension);
+  CHK(desc.gyration_radius == args.gyration_radius);
+  CHK(desc.nintervals = args.nintervals[0] + args.nintervals[1]);
+
+  CHK(eq_eps(desc.theta_limit, ref_theta_limit, ref_theta_limit*1.e-6));
+  CHK(eq_eps(desc.normalization_factor, ref_norm_factor, ref_norm_factor*1.e-6));
+
+  CHK(ssf_phase_ref_put(phase) == RES_OK);
+  CHK(ssf_phase_ref_put(dummy) == RES_OK);
+  CHK(ssp_rng_ref_put(rng) == RES_OK);
+
+  check_memory_allocator(&allocator);
+  mem_shutdown_proxy_allocator(&allocator);
+  CHK(mem_allocated_size() == 0);
+  return 0;
+}