#include <string.h>
+#include <math.h>
+
#include "math.h"
int v3_add(t_3dvec *sum,t_3dvec *a,t_3dvec *b) {
return(memcmp(a,b,sizeof(t_3dvec)));
}
+double v3_absolute_square(t_3dvec *a) {
+
+ return(a->x*a->x+a->y*a->y+a->z*a->z);
+}
+
+double v3_norm(t_3dvec *a) {
+
+ return(sqrt(v3_absolute_square(a)));
+}
+
+int v3_per_bound(t_3dvec *a,t_3dvec *dim) {
+
+ double x,y,z;
+
+ x=0.5*dim->x;
+ y=0.5*dim->y;
+ z=0.5*dim->z;
+
+ if(a->x>=x) a->x-=dim->x;
+ else if(-a->x>x) a->x+=dim->x;
+ if(a->y>=y) a->y-=dim->y;
+ else if(-a->y>y) a->y+=dim->y;
+ if(a->z>=z) a->z-=dim->z;
+ else if(-a->z>z) a->z+=dim->z;
+
+ return 0;
+}