---

vec.h

Go to the documentation of this file.

#ifndef _3DMANIP_H_
#define _3DMANIP_H_
#include <math.h>
#define QFLOAT float
#define XSQRT sqrtf
#define XVector Vector
#define YVector QVector
#include "xvector.h"
#undef QFLOAT
#undef XVector
#undef YVector
#undef XSQRT
#define QFLOAT double
#define XSQRT sqrt
#define XVector QVector
#define YVector Vector

#include "xvector.h"
#undef XSQRT
#undef QFLOAT
#undef XVector
#undef YVector

inline Vector QVector::operator =(const Vector &a) {
  i=a.i;
  j=a.j;
  k=a.k;
  return a;
}
inline QVector::QVector (const Vector &a) {
  i=a.i;
  j=a.j;
  k=a.k;
}
inline QVector Vector::operator = (const QVector &a) {
  i=a.i;
  j=a.j;
  k=a.k;
  return a;
}
inline Vector::Vector (const QVector &a) {
  i=a.i;
  j=a.j;
  k=a.k;
}
inline QVector Vector::Cast() const{
  return QVector (i,j,k);
}
inline Vector QVector::Cast() const{
  return Vector (i,j,k);
}

template <class XVector, class DoItType> inline void forEach3dVector (XVector first, XVector last, DoItType DoIt, XVector increment) {
    first.Set(floor(first.i/increment.i)*increment.i,
              floor(first.j/increment.j)*increment.j,
              floor(first.k/increment.k)*increment.k);
    last.Set(ceil(last.i/increment.i)*increment.i,
             ceil(last.j/increment.j)*increment.j,
             ceil(last.k/increment.k)*increment.k);
    XVector xvec (first);
    for (first.k=xvec.k;first.k<last.k;first.k+=increment.k) {
        for (first.j=xvec.j;first.j<last.j;first.j+=increment.j) {
            for (first.i=xvec.i;first.i<last.i;first.i+=increment.i) {
                DoIt(first);
            }
        }
    }
}


template <class XVector, class DoItType> inline void forEach2dVector (XVector first, XVector last, DoItType DoIt, XVector increment) {
    first.Set(floor(first.i/increment.i)*increment.i,
              floor(first.j/increment.j)*increment.j,
              first.k);
    last.Set(ceil(last.i/increment.i)*increment.i,
              ceil(last.j/increment.j)*increment.j,
              last.k);
    XVector xvec (first);
    for (first.j=xvec.j;first.j<last.j;first.j+=increment.j) {
        for (first.i=xvec.i;first.i<last.i;first.i+=increment.i) {
            DoIt(first);
        }
    }
}

#endif

---