53 lines
1.7 KiB
C
53 lines
1.7 KiB
C
/* ************************************************************************** */
|
|
/* */
|
|
/* ::: :::::::: */
|
|
/* maths_obj_sphere.c :+: :+: :+: */
|
|
/* +:+ +:+ +:+ */
|
|
/* By: ntrahy <marvin@42.fr> +#+ +:+ +#+ */
|
|
/* +#+#+#+#+#+ +#+ */
|
|
/* Created: 2017/01/03 00:45:18 by ntrahy #+# #+# */
|
|
/* Updated: 2017/01/03 13:54:11 by scebula ### ########.fr */
|
|
/* */
|
|
/* ************************************************************************** */
|
|
|
|
#include "rt.h"
|
|
#include <math.h>
|
|
|
|
static void check_for_solutions_sphere(double *dist)
|
|
{
|
|
if (dist[0] > dist[1])
|
|
swap_double(&dist[0], &dist[1]);
|
|
if (dist[0] < 0 && dist[1] > 0)
|
|
dist[2] = dist[1];
|
|
else if (dist[0] > 0 && dist[1] > 0)
|
|
dist[2] = dist[0];
|
|
else
|
|
{
|
|
dist[0] = 0;
|
|
dist[1] = 0;
|
|
}
|
|
}
|
|
|
|
void get_sphere_n(t_ray ray, t_obj *obj, t_thread *thread,
|
|
double *dist)
|
|
{
|
|
double a;
|
|
double b;
|
|
double c;
|
|
double delta;
|
|
|
|
(void)thread;
|
|
(void)obj;
|
|
a = pow(ray.dir.x, 2) + pow(ray.dir.y, 2) + pow(ray.dir.z, 2);
|
|
b = (2 * ray.ori.x * ray.dir.x) + (2 * ray.ori.y * ray.dir.y) +
|
|
(2 * ray.ori.z * ray.dir.z);
|
|
c = pow(ray.ori.x, 2) + pow(ray.ori.y, 2) + pow(ray.ori.z, 2) - 1;
|
|
delta = b * b - 4 * a * c;
|
|
if (delta > 0)
|
|
{
|
|
dist[0] = (((-1 * b - sqrt(delta))) / (a + a)) - ACC;
|
|
dist[1] = (((-1 * b + sqrt(delta))) / (a + a)) - ACC;
|
|
check_for_solutions_sphere(dist);
|
|
}
|
|
}
|