70 lines
2.1 KiB
C
70 lines
2.1 KiB
C
/* ************************************************************************** */
|
|
/* */
|
|
/* ::: :::::::: */
|
|
/* maths_matrix.c :+: :+: :+: */
|
|
/* +:+ +:+ +:+ */
|
|
/* By: ntrahy <marvin@42.fr> +#+ +:+ +#+ */
|
|
/* +#+#+#+#+#+ +#+ */
|
|
/* Created: 2017/01/03 00:43:22 by ntrahy #+# #+# */
|
|
/* Updated: 2017/01/03 00:43:24 by ntrahy ### ########.fr */
|
|
/* */
|
|
/* ************************************************************************** */
|
|
|
|
#include "rt.h"
|
|
#include <math.h>
|
|
|
|
static t_matrix transposed_matrix(t_matrix m)
|
|
{
|
|
t_matrix matrix;
|
|
|
|
matrix.type = m.type;
|
|
matrix.m11 = m.m11;
|
|
matrix.m12 = m.m21;
|
|
matrix.m13 = m.m31;
|
|
matrix.m14 = m.m41;
|
|
matrix.m21 = m.m12;
|
|
matrix.m22 = m.m22;
|
|
matrix.m23 = m.m32;
|
|
matrix.m24 = m.m42;
|
|
matrix.m31 = m.m13;
|
|
matrix.m32 = m.m23;
|
|
matrix.m33 = m.m33;
|
|
matrix.m34 = m.m43;
|
|
matrix.m41 = m.m14;
|
|
matrix.m42 = m.m24;
|
|
matrix.m43 = m.m34;
|
|
matrix.m44 = m.m44;
|
|
return (matrix);
|
|
}
|
|
|
|
void apply_matrix(t_matrix matrix, t_vect *v)
|
|
{
|
|
double x;
|
|
double y;
|
|
double z;
|
|
double w;
|
|
|
|
x = v->x;
|
|
y = v->y;
|
|
z = v->z;
|
|
w = v->w;
|
|
v->x = matrix.m11 * x + matrix.m12 * y + matrix.m13 * z + matrix.m14 * w;
|
|
v->y = matrix.m21 * x + matrix.m22 * y + matrix.m23 * z + matrix.m24 * w;
|
|
v->z = matrix.m31 * x + matrix.m32 * y + matrix.m33 * z + matrix.m34 * w;
|
|
v->w = matrix.m41 * x + matrix.m42 * y + matrix.m43 * z + matrix.m44 * w;
|
|
}
|
|
|
|
t_matrix get_inverse_matrix(t_matrix m)
|
|
{
|
|
t_matrix inverse;
|
|
|
|
if (m.type == ROTATION_Y || m.type == ROTATION_X
|
|
|| m.type == ROTATION_Z || m.type == ROTATION_XYZ)
|
|
inverse = transposed_matrix(m);
|
|
else if (m.type == SCALE)
|
|
inverse = init_scale_matrix(1 / m.m11, 1 / m.m22, 1 / m.m33);
|
|
else
|
|
inverse = init_translation_matrix(-m.m14, -m.m24, -m.m34);
|
|
return (inverse);
|
|
}
|