/* Mathematiques de l'informatique graphique    */
/* Direction en 3D                              */
/*                                              */
/* Mars 2021                                    */

#include <math.h>

#include "CH3D.h"
#include "Dir3D.h"
#include "Pos3D.h"

/* Constructeurs                                */

/* Construit le vecteur z                       */

Dir3D::Dir3D(void):CH3D(0.0F,0.0F,1.0F,0.0F) {
}

/* Construit la direction (xp,yp,zp)               */

Dir3D::Dir3D(float xp,float yp,float zp):CH3D(xp,yp,zp,0.0F) {
}

/* Construit un clone de la Dir3D d             */

Dir3D::Dir3D(Dir3D *d):CH3D(d) {
}

/* Construit la direction p1 p2                 */

Dir3D::Dir3D(const Pos3D &p1,const Pos3D &p2) : CH3D(p2.x-p1.x,p2.y-p1.y,p2.z-p1.z,0.0f) {
}

/* Destructeur                                  */

Dir3D::~Dir3D(void) {
}

/* Methode de calcul de la norme de this        */

float Dir3D::norme(void) {
  return((float) sqrt(x*x+y*y+z*z));
}

/* Calcul de la direction normalisée de this    */

Dir3D Dir3D::normalize(void) {
  float norm = norme();
  Dir3D res;
  if (norm != 0.0F) {
    res.x = x / norm;
    res.y = y / norm;
    res.z = z / norm; }
  return res;
}

/* Calcul du produit scalaire                   */
/* de this et de la direction dir               */

float Dir3D::operator*(const Dir3D& dir) {
  return this->x*dir.x + this->y*dir.y + this->z * dir.z;
}

/* Calcul du produit vectoriel                  */
/* de this et de la direction dir               */

Dir3D Dir3D::operator^(const Dir3D& dir) {
  Dir3D res;
  res.x = this->y * dir.z - dir.y * this->z;
  res.y = this->z * dir.x - dir.z * this->x;
  res.z = this->x * dir.y - dir.x * this->y;
  return res;
}