-
Notifications
You must be signed in to change notification settings - Fork 0
/
cal_matrix.c
127 lines (114 loc) · 6.24 KB
/
cal_matrix.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
/* ************************************************************************** */
/* */
/* ::: :::::::: */
/* cal_matrix.c :+: :+: :+: */
/* +:+ +:+ +:+ */
/* By: alamy <[email protected]> +#+ +:+ +#+ */
/* +#+#+#+#+#+ +#+ */
/* Created: 2018/01/20 16:46:30 by alamy #+# #+# */
/* Updated: 2018/02/02 16:28:27 by alamy ### ########.fr */
/* */
/* ************************************************************************** */
#include "fdf.h"
t_vecteur4 ft_cal_translation(t_vecteur4 vecteur2, matrix4_t matrix_translation)
{
t_vecteur4 vec_translation;
vec_translation.x1 = (vecteur2.x1 * matrix_translation.m[0][0]) +
(vecteur2.y1 * matrix_translation.m[0][1]) + (vecteur2.z1 * matrix_translation.m[0][2])
+ (vecteur2.w1 * matrix_translation.m[0][3]);
vec_translation.y1 = (vecteur2.x1 * matrix_translation.m[1][0]) +
(vecteur2.y1 * matrix_translation.m[1][1]) + (vecteur2.z1 * matrix_translation.m[1][2])
+ (vecteur2.w1 * matrix_translation.m[1][3]);
vec_translation.z1 = (vecteur2.x1 * matrix_translation.m[2][0]) +
(vecteur2.y1 * matrix_translation.m[2][1]) + (vecteur2.z1 * matrix_translation.m[2][2])
+ (vecteur2.w1 * matrix_translation.m[2][3]);
vec_translation.w1 = (vecteur2.x1 * matrix_translation.m[3][0]) +
(vecteur2.y1 * matrix_translation.m[3][1]) + (vecteur2.z1 * matrix_translation.m[3][2])
+ (vecteur2.w1 * matrix_translation.m[3][3]);
return(vec_translation);
}
t_vecteur4 ft_cal_rotationX(t_vecteur4 vecteur, matrix4_t matrix_rotationX)
{
t_vecteur4 vec_rotationX;
vec_rotationX.x1 = (vecteur.x1 * matrix_rotationX.m[0][0]) +
(vecteur.y1 * matrix_rotationX.m[0][1]) + (vecteur.z1 * matrix_rotationX.m[0][2])
+ (vecteur.w1 * matrix_rotationX.m[0][3]);
vec_rotationX.y1 = (vecteur.x1 * matrix_rotationX.m[1][0]) +
(vecteur.y1 * matrix_rotationX.m[1][1]) + (vecteur.z1 * matrix_rotationX.m[1][2])
+ (vecteur.w1 * matrix_rotationX.m[1][3]);
vec_rotationX.z1 = (vecteur.x1 * matrix_rotationX.m[2][0]) +
(vecteur.y1 * matrix_rotationX.m[2][1]) + (vecteur.z1 * matrix_rotationX.m[2][2])
+ (vecteur.w1 * matrix_rotationX.m[2][3]);
vec_rotationX.w1 = (vecteur.x1 * matrix_rotationX.m[3][0]) +
(vecteur.y1 * matrix_rotationX.m[3][1]) + (vecteur.z1 * matrix_rotationX.m[3][2])
+ (vecteur.w1 * matrix_rotationX.m[3][3]);
return(vec_rotationX);
}
t_vecteur4 ft_cal_rotationY(t_vecteur4 vecteur, matrix4_t matrix_rotationY)
{
t_vecteur4 vec_rotationY;
vec_rotationY.x1 = (vecteur.x1 * matrix_rotationY.m[0][0]) +
(vecteur.y1 * matrix_rotationY.m[0][1]) + (vecteur.z1 * matrix_rotationY.m[0][2])
+ (vecteur.w1 * matrix_rotationY.m[0][3]);
vec_rotationY.y1 = (vecteur.x1 * matrix_rotationY.m[1][0]) +
(vecteur.y1 * matrix_rotationY.m[1][1]) + (vecteur.z1 * matrix_rotationY.m[1][2])
+ (vecteur.w1 * matrix_rotationY.m[1][3]);
vec_rotationY.z1 = (vecteur.x1 * matrix_rotationY.m[2][0]) +
(vecteur.y1 * matrix_rotationY.m[2][1]) + (vecteur.z1 * matrix_rotationY.m[2][2])
+ (vecteur.w1 * matrix_rotationY.m[2][3]);
vec_rotationY.w1 = (vecteur.x1 * matrix_rotationY.m[3][0]) +
(vecteur.y1 * matrix_rotationY.m[3][1]) + (vecteur.z1 * matrix_rotationY.m[3][2])
+ (vecteur.w1 * matrix_rotationY.m[3][3]);
return(vec_rotationY);
}
t_vecteur4 ft_cal_rotationZ(t_vecteur4 vecteur, matrix4_t matrix_rotationZ)
{
t_vecteur4 vec_rotationZ;
vec_rotationZ.x1 = (vecteur.x1 * matrix_rotationZ.m[0][0]) +
(vecteur.y1 * matrix_rotationZ.m[0][1]) + (vecteur.z1 * matrix_rotationZ.m[0][2])
+ (vecteur.w1 * matrix_rotationZ.m[0][3]);
vec_rotationZ.y1 = (vecteur.x1 * matrix_rotationZ.m[1][0]) +
(vecteur.y1 * matrix_rotationZ.m[1][1]) + (vecteur.z1 * matrix_rotationZ.m[1][2])
+ (vecteur.w1 * matrix_rotationZ.m[1][3]);
vec_rotationZ.z1 = (vecteur.x1 * matrix_rotationZ.m[2][0]) +
(vecteur.y1 * matrix_rotationZ.m[2][1]) + (vecteur.z1 * matrix_rotationZ.m[2][2])
+ (vecteur.w1 * matrix_rotationZ.m[2][3]);
vec_rotationZ.w1 = (vecteur.x1 * matrix_rotationZ.m[3][0]) +
(vecteur.y1 * matrix_rotationZ.m[3][1]) + (vecteur.z1 * matrix_rotationZ.m[3][2])
+ (vecteur.w1 * matrix_rotationZ.m[3][3]);
return(vec_rotationZ);
}
t_vecteur4 ft_cal_homothetie(t_vecteur4 vecteur, matrix4_t matrix_homothetie)
{
t_vecteur4 vec_homothetie;
vec_homothetie.x1 = (vecteur.x1 * matrix_homothetie.m[0][0]) +
(vecteur.y1 * matrix_homothetie.m[0][1]) + (vecteur.z1 * matrix_homothetie.m[0][2])
+ (vecteur.w1 * matrix_homothetie.m[0][3]);
vec_homothetie.y1 = (vecteur.x1 * matrix_homothetie.m[1][0]) +
(vecteur.y1 * matrix_homothetie.m[1][1]) + (vecteur.z1 * matrix_homothetie.m[1][2])
+ (vecteur.w1 * matrix_homothetie.m[1][3]);
vec_homothetie.z1 = (vecteur.x1 * matrix_homothetie.m[2][0]) +
(vecteur.y1 * matrix_homothetie.m[2][1]) + (vecteur.z1 * matrix_homothetie.m[2][2])
+ (vecteur.w1 * matrix_homothetie.m[2][3]);
vec_homothetie.w1 = (vecteur.x1 * matrix_homothetie.m[3][0]) +
(vecteur.y1 * matrix_homothetie.m[3][1]) + (vecteur.z1 * matrix_homothetie.m[3][2])
+ (vecteur.w1 * matrix_homothetie.m[3][3]);
return(vec_homothetie);
}
t_vecteur4 ft_cal_projection(t_vecteur4 vecteur, matrix4_t matrix_projection)
{
t_vecteur4 vec_projection;
vec_projection.x1 = (vecteur.x1 * matrix_projection.m[0][0]) +
(vecteur.y1 * matrix_projection.m[0][1]) + (vecteur.z1 * matrix_projection.m[0][2])
+ (vecteur.w1 * matrix_projection.m[0][3]);
vec_projection.y1 = (vecteur.x1 * matrix_projection.m[1][0]) +
(vecteur.y1 * matrix_projection.m[1][1]) + (vecteur.z1 * matrix_projection.m[1][2])
+ (vecteur.w1 * matrix_projection.m[1][3]);
vec_projection.z1 = (vecteur.x1 * matrix_projection.m[2][0]) +
(vecteur.y1 * matrix_projection.m[2][1]) + (vecteur.z1 * matrix_projection.m[2][2])
+ (vecteur.w1 * matrix_projection.m[2][3]);
vec_projection.w1 = (vecteur.x1 * matrix_projection.m[3][0]) +
(vecteur.y1 * matrix_projection.m[3][1]) + (vecteur.z1 * matrix_projection.m[3][2])
+ (vecteur.w1 * matrix_projection.m[3][3]);
return(vec_projection);
}