test_ssf_fresnel_dielectric_dielectric.c (2137B)
1 /* Copyright (C) 2016-2018, 2021-2025 |Méso|Star> (contact@meso-star.com) 2 * 3 * This program is free software: you can redistribute it and/or modify 4 * it under the terms of the GNU General Public License as published by 5 * the Free Software Foundation, either version 3 of the License, or 6 * (at your option) any later version. 7 * 8 * This program is distributed in the hope that it will be useful, 9 * but WITHOUT ANY WARRANTY; without even the implied warranty of 10 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 11 * GNU General Public License for more details. 12 * 13 * You should have received a copy of the GNU General Public License 14 * along with this program. If not, see <http://www.gnu.org/licenses/>. */ 15 16 #include "ssf.h" 17 #include "test_ssf_utils.h" 18 19 #include <rsys/math.h> 20 21 int 22 main(int argc, char** argv) 23 { 24 const size_t nsteps = 100; 25 const double step = (PI/2.0) / (double)nsteps; 26 size_t i; 27 struct mem_allocator allocator; 28 struct ssf_fresnel* fresnel; 29 struct ssf_fresnel* dummy; 30 double val; 31 (void)argc, (void)argv; 32 33 mem_init_proxy_allocator(&allocator, &mem_default_allocator); 34 CHK(ssf_fresnel_create(&allocator, &fresnel_dummy, &dummy) == RES_OK); 35 CHK(ssf_fresnel_create 36 (&allocator, &ssf_fresnel_dielectric_dielectric, &fresnel) == RES_OK); 37 38 CHK(ssf_fresnel_dielectric_dielectric_setup(NULL, -1.0, -1.0) == RES_BAD_ARG); 39 CHK(ssf_fresnel_dielectric_dielectric_setup(fresnel, -1.0, -1.0) == RES_OK); 40 CHK(ssf_fresnel_dielectric_dielectric_setup(dummy,-1.0, -1.0) == RES_BAD_ARG); 41 CHK(ssf_fresnel_dielectric_dielectric_setup(fresnel,1.000277, 1.5) == RES_OK); 42 43 val = ssf_fresnel_eval(fresnel, 1); 44 CHK(eq_eps(val, 0.04, 1.e-4) == 1); 45 46 FOR_EACH(i, 0, nsteps+1) { 47 const double theta_i = MMIN((double)i*step, PI/2); 48 const double cos_theta_i = cos(theta_i); 49 val = ssf_fresnel_eval(fresnel, cos_theta_i); 50 printf("%g %g\n", theta_i, val); 51 } 52 53 CHK(ssf_fresnel_ref_put(fresnel) == RES_OK); 54 CHK(ssf_fresnel_ref_put(dummy) == RES_OK); 55 56 check_memory_allocator(&allocator); 57 mem_shutdown_proxy_allocator(&allocator); 58 CHK(mem_allocated_size() == 0); 59 return 0; 60 }