src/HEXABLOCK/HexVector.hxx

   1
   2 // class : Les vecteurs
   3
   4 // Copyright (C) 2009-2013  CEA/DEN, EDF R&D
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 #ifndef __VECTOR_H
  23 #define __VECTOR_H
  24
  25 #include "HexEltBase.hxx"
  26
  27 #include <cmath>
  28
  29 BEGIN_NAMESPACE_HEXA
  30
  31 class Vector : public EltBase
  32 {
  33 public:
  34     double getDx ()    { return v_dx ; }
  35     double getDy ()    { return v_dy ; }
  36     double getDz ()    { return v_dz ; }
  37     double getNorm   ()    { return sqrt (v_dx*v_dx + v_dy*v_dy + v_dz*v_dz); }
  38     double getAngleX ();
  39
  40
  41 public:
  42     Vector (Document* doc, double dx=0, double dy=0, double dz=0);
  43     Vector (Vector* lautre);
  44     virtual ~Vector () {}
  45     virtual void dump ();
  46     double  getCoord (int dir);
  47     double* getCoord (double coord[]);
  48
  49     virtual void saveXml (XmlWriter* xml);
  50     int    renormer  ();
  51     void   vectoriel  (Vector* a, Vector*b);
  52     void   multiplier (double scal);
  53     int    calculNormale (double& nx, double& ny, double& nz);
  54
  55 private:
  56     double v_dx;
  57     double v_dy;
  58     double v_dz;
  59 };
  60 // ------------------------------------------- Inlining
  61 // ===================================================== getCoord
  62 inline double Vector::getCoord (int dir)
  63 {
  64    double val = 0;
  65    switch (dir)
  66           {
  67           case dir_x : val = v_dx;
  68                break;
  69           case dir_y : val = v_dy;
  70                break;
  71           case dir_z : val = v_dz;
  72                break;
  73           }
  74    return val;
  75 }
  76 // ===================================================== renormer
  77 inline int Vector::renormer()
  78 {
  79    double dn = getNorm ();
  80    if (dn < 1e-30)
  81       return HERR;
  82
  83    v_dx /= dn;
  84    v_dy /= dn;
  85    v_dz /= dn;
  86
  87    return HOK;
  88 }
  89 // ===================================================== vectoriel
  90 inline void Vector::vectoriel (Vector* a, Vector*b)
  91 {
  92    v_dx =   a->v_dy * b->v_dz - b->v_dy * a->v_dz;
  93    v_dy =   a->v_dz * b->v_dx - b->v_dz * a->v_dx;
  94    v_dz =   a->v_dx * b->v_dy - b->v_dx * a->v_dy;
  95 }
  96
  97 // ===================================================== getCoord
  98 inline double* Vector::getCoord (double coord[])
  99 {
 100    coord [dir_x] = v_dx;
 101    coord [dir_y] = v_dy;
 102    coord [dir_z] = v_dz;
 103    return coord;
 104 }
 105 // ===================================================== multiplier
 106 inline void Vector::multiplier (double scalaire)
 107 {
 108    v_dx *= scalaire;
 109    v_dy *= scalaire;
 110    v_dz *= scalaire;
 111 }
 112 // ===================================================== dump
 113 inline void Vector::dump ()
 114 {
 115    printf ("Vector ");
 116    printName (" = ");
 117    if (NOT isHere ())
 118       {
 119       printf ("(*** deleted ***)\n");
 120       return;
 121       }
 122
 123    printf ("(%g, %g, %g)\n", v_dx,v_dy,v_dz);
 124 }
 125
 126 END_NAMESPACE_HEXA
 127 #endif