/Users/barrand/private/dev/softinex/old/inexlib-1.2/inlib/inlib/a3d

Go to the documentation of this file.

// Copyright (C) 2010, Guy Barrand. All rights reserved.
// See the file inlib.license for terms.

#ifndef inlib_a3d
#define inlib_a3d

#include "a3"

namespace inlib {
namespace a3d {

class vec : public a3::vec<double> {
public:
  vec(): a3::vec<double>() {}
  vec(const double a_vec[3]): a3::vec<double>(a_vec) {}
  vec(const double& a0,const double& a1,const double& a2)
  : a3::vec<double>(a0,a1,a2){}
  virtual ~vec() {}
public:
  vec(const vec& a_from): a3::vec<double>(a_from){}
  vec& operator=(const vec& a_from){
    a3::vec<double>::operator=(a_from);
    return *this;
  }
public:
/*
  double length() const {
    return ::sqrt(m_vector[0] * m_vector[0] + 
                  m_vector[1] * m_vector[1] + 
                  m_vector[2] * m_vector[2]); 
  }
  double normalize() {
    double norme = length();
    if(norme==0) return 0;
    divide(norme);
    return norme;
  }
  double dot(const vec& aV) const {
    return (m_vector[0] * aV.m_vector[0] + 
            m_vector[1] * aV.m_vector[1] + 
            m_vector[2] * aV.m_vector[2]);
  }    
*/
  vec cross(const vec& aV) const {
    return vec(m_vector[1] * aV.m_vector[2] - m_vector[2] * aV.m_vector[1],
               m_vector[2] * aV.m_vector[0] - m_vector[0] * aV.m_vector[2],
               m_vector[0] * aV.m_vector[1] - m_vector[1] * aV.m_vector[0]);
  }    
};


//  1  0  0
//  0  1  0
//  0  0  1
class id : public a3::sqm<double> {
public:
  id(){
    v00(1);
    v11(1);
    v22(1);
  }
};

// Generators of rotation group :
// Rk(i,j) = epsilon(k,i,j)  i,j,k=1,2,3
//  0  0  0
//  0  0  1
//  0 -1  0
class R1 : public a3::sqm<double> {
public:
  R1(){
    v12( 1);
    v21(-1);
  }
};
//  0  0 -1
//  0  0  0
//  1  0  0
class R2 : public a3::sqm<double> {
public:
  R2(){
    v02(-1);
    v20( 1);
  }
};
//  0  1  0
// -1  0  0
//  0  0  0
class R3 : public a3::sqm<double> {
public:
  R3(){
    v01( 1);
    v10(-1);
  }
};

// Levi-Civita indicator :
class epsilon : public a3::cbm<double> {
public:
  epsilon(){
    v210(-1);
    v120( 1);
    v201( 1);
    v021(-1);
    v102(-1);
    v012( 1);
  }
};

}}

#endif