intercept.h

/*	intercept.h
 *
 * Copyright (C) 2014 - 2018 Ivo Alxneit, Paul Scherrer Institute
 *
 * This file is part of rt
 *
 * rt 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.
 *
 * rt 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 rt. If not, see <http://www.gnu.org/licenses/>.
 *
 */

#ifndef __INTERCEPT_H__
#define __INTERCEPT_H__

#include <pthread.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <unistd.h>

#include <gsl/gsl_math.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_spline.h>

#include <libconfig.h>

#include "ray.h"
#include "math_utils.h"
#include "reflect.h"
#include "targets.h"


/*
 * surface normals of non-planar targets (pointing away from convex side
 */
extern void cone_surf_normal(double *const intercept, const double tan2_a,
			     const double H, double *const normal);
extern void cyl_surf_normal(double *const icpt, const double *C,
			    const double *a, const double r,
			    double *const normal);
extern void ell_surf_normal(const double *point, const double *axes,
			    double *const normal);
extern void par_surf_normal(const double *point, const double foc2,
			    double *const normal);
extern void sph_surf_normal(const double *point, double *normal);



/*
 * functions to return intercepts of ray with target
 */
double *intercept_cone(const ray_t * ray, const double *M,
		       const double origin[3], const double tan2_a,
		       const double H, const double z_max,
		       int *hits_outside);
extern double *intercept_cylinder(const ray_t * ray, const double *c,
				  const double *a, const double r,
				  const double l, int *hits_outside);
extern double *intercept_disk(const ray_t * ray, const double *origin,
			      const double *M, const double R2,
			      int *hits_front);
extern double *intercept_ellipsoid(const ray_t * ray, const double *M,
				   const double center[3],
				   const double axes[3],
				   const double z_min, const double z_max,
				   int *hits_outside);
extern double *intercept_paraboloid(const ray_t * ray, const double *M,
				    const double vertex[3],
				    const double foc2, const double foc4,
				    const double z_min, const double z_max,
				    int *hits_outside);
extern double *intercept_plane(const ray_t * ray,
			       const double *plane_normal,
			       const double *plane_point, int *hits_front);
extern double *intercept_sphere(const ray_t * ray, const double *M,
				const double *center, const double R2,
				const double z_min, const double z_max,
				int *hits_outside);

#endif				/* __INTERCEPT_H__ */