AMROC Main     Blockstructured Adaptive Mesh Refinement in object-oriented C++


Main Page   Class Hierarchy   Compound List   File List  

Vector17.h

Go to the documentation of this file.
00001 #ifndef AMROC_VECTOR17_H
00002 #define AMROC_VECTOR17_H
00003 
00011 #include "CompConf.h"
00012 #include "generic.h"
00013 #include <iostream.h>
00014 #include <assert.h>
00015  
00016 #ifndef VectorName
00017 #define Vector(dim)      name2(Vector,dim)
00018 #define VectorName
00019 #endif
00020 
00028 template <class DataType>
00029 class Vector(17) {
00030   typedef Vector(17)<DataType> TVector;
00031 
00032 public:
00033   typedef DataType InternalDataType;
00034 
00035   Vector(17)() {}
00036   Vector(17)(const TVector& x) { 
00037     memcpy((void *) data_, (void *) x.data_, sizeof(TVector));
00038   }
00039   Vector(17)(const DataType& v) { 
00040     data_[0] = v; 
00041     data_[1] = v; 
00042     data_[2] = v; 
00043     data_[3] = v; 
00044     data_[4] = v; 
00045     data_[5] = v; 
00046     data_[6] = v; 
00047     data_[7] = v; 
00048     data_[8] = v; 
00049     data_[9] = v; 
00050     data_[10] = v; 
00051     data_[11] = v; 
00052     data_[12] = v; 
00053     data_[13] = v; 
00054     data_[14] = v; 
00055     data_[15] = v; 
00056     data_[16] = v; 
00057   }
00058 
00059   ~Vector(17)() {}
00060   
00061   DataType * data() { return data_; }
00062   const DataType * data() const { return data_; }
00063   int length() const { return 17; }
00064   static int Length() { return 17; }
00065 
00066   inline DataType operator()(int i) const {
00067 #ifdef DEBUG_PRINT
00068     assert((i >= 0) && (i < 17)); 
00069 #endif
00070     return data_[i];
00071   }
00072   inline DataType& operator()(int i) { 
00073 #ifdef DEBUG_PRINT
00074     assert((i >= 0) && (i < 17)); 
00075 #endif
00076     return data_[i];
00077   }
00078 
00079   inline DataType operator[](int i) const {
00080 #ifdef DEBUG_PRINT
00081     assert((i >= 0) && (i < 17)); 
00082 #endif
00083     return data_[i];
00084   }
00085   inline DataType& operator[](int i) { 
00086 #ifdef DEBUG_PRINT
00087     assert((i >= 0) && (i < 17)); 
00088 #endif
00089     return data_[i];
00090   }
00091 
00092   inline TVector& operator=(const TVector& x) { 
00093     memcpy((void *) data_, (void *) x.data_, sizeof(TVector));
00094     return *this; 
00095   }
00096   inline TVector& operator=(const DataType& v) { 
00097     data_[0] = v; 
00098     data_[1] = v; 
00099     data_[2] = v; 
00100     data_[3] = v; 
00101     data_[4] = v; 
00102     data_[5] = v; 
00103     data_[6] = v; 
00104     data_[7] = v; 
00105     data_[8] = v; 
00106     data_[9] = v; 
00107     data_[10] = v; 
00108     data_[11] = v; 
00109     data_[12] = v; 
00110     data_[13] = v; 
00111     data_[14] = v; 
00112     data_[15] = v; 
00113     data_[16] = v; 
00114     return *this; 
00115   }
00116 
00117 #define Vector_Vector_Operator(ope,op)                          \
00118    inline TVector& operator ope (const TVector &x)              \
00119      {                                                          \
00120       data_[0] ope x.data_[0];                                  \
00121       data_[1] ope x.data_[1];                                  \
00122       data_[2] ope x.data_[2];                                  \
00123       data_[3] ope x.data_[3];                                  \
00124       data_[4] ope x.data_[4];                                  \
00125       data_[5] ope x.data_[5];                                  \
00126       data_[6] ope x.data_[6];                                  \
00127       data_[7] ope x.data_[7];                                  \
00128       data_[8] ope x.data_[8];                                  \
00129       data_[9] ope x.data_[9];                                  \
00130       data_[10] ope x.data_[10];                                \
00131       data_[11] ope x.data_[11];                                \
00132       data_[12] ope x.data_[12];                                \
00133       data_[13] ope x.data_[13];                                \
00134       data_[14] ope x.data_[14];                                \
00135       data_[15] ope x.data_[15];                                \
00136       data_[16] ope x.data_[16];                                \
00137       return(*this);                                            \
00138      }                                                          \
00139    inline TVector operator op (const TVector &x) const          \
00140      {                                                          \
00141       TVector ret(*this);                                       \
00142       ret.data_[0] ope x.data_[0];                              \
00143       ret.data_[1] ope x.data_[1];                              \
00144       ret.data_[2] ope x.data_[2];                              \
00145       ret.data_[3] ope x.data_[3];                              \
00146       ret.data_[4] ope x.data_[4];                              \
00147       ret.data_[5] ope x.data_[5];                              \
00148       ret.data_[6] ope x.data_[6];                              \
00149       ret.data_[7] ope x.data_[7];                              \
00150       ret.data_[8] ope x.data_[8];                              \
00151       ret.data_[9] ope x.data_[9];                              \
00152       ret.data_[10] ope x.data_[10];                            \
00153       ret.data_[11] ope x.data_[11];                            \
00154       ret.data_[12] ope x.data_[12];                            \
00155       ret.data_[13] ope x.data_[13];                            \
00156       ret.data_[14] ope x.data_[14];                            \
00157       ret.data_[15] ope x.data_[15];                            \
00158       ret.data_[16] ope x.data_[16];                            \
00159       return(ret);                                              \
00160      }
00161    Vector_Vector_Operator(+=,+)
00162    Vector_Vector_Operator(-=,-)
00163    Vector_Vector_Operator(*=,*)
00164    Vector_Vector_Operator(/=,/)
00165 #undef Vector_Vector_Operator
00166 
00167 #define Vector_Scalar_Operator(ope,op)                          \
00168    inline TVector& operator ope (const DataType &v)             \
00169      {                                                          \
00170       data_[0] ope v;                                           \
00171       data_[1] ope v;                                           \
00172       data_[2] ope v;                                           \
00173       data_[3] ope v;                                           \
00174       data_[4] ope v;                                           \
00175       data_[5] ope v;                                           \
00176       data_[6] ope v;                                           \
00177       data_[7] ope v;                                           \
00178       data_[8] ope v;                                           \
00179       data_[9] ope v;                                           \
00180       data_[10] ope v;                                          \
00181       data_[11] ope v;                                          \
00182       data_[12] ope v;                                          \
00183       data_[13] ope v;                                          \
00184       data_[14] ope v;                                          \
00185       data_[15] ope v;                                          \
00186       data_[16] ope v;                                          \
00187       return(*this);                                            \
00188      }                                                          \
00189    inline TVector operator op (const DataType &v) const         \
00190      {                                                          \
00191       TVector ret(*this);                                       \
00192       ret.data_[0] ope v;                                       \
00193       ret.data_[1] ope v;                                       \
00194       ret.data_[2] ope v;                                       \
00195       ret.data_[3] ope v;                                       \
00196       ret.data_[4] ope v;                                       \
00197       ret.data_[5] ope v;                                       \
00198       ret.data_[6] ope v;                                       \
00199       ret.data_[7] ope v;                                       \
00200       ret.data_[8] ope v;                                       \
00201       ret.data_[9] ope v;                                       \
00202       ret.data_[10] ope v;                                      \
00203       ret.data_[11] ope v;                                      \
00204       ret.data_[12] ope v;                                      \
00205       ret.data_[13] ope v;                                      \
00206       ret.data_[14] ope v;                                      \
00207       ret.data_[15] ope v;                                      \
00208       ret.data_[16] ope v;                                      \
00209      return(ret);                                               \
00210      }
00211    Vector_Scalar_Operator(+=,+)
00212    Vector_Scalar_Operator(-=,-)
00213    Vector_Scalar_Operator(*=,*)
00214    Vector_Scalar_Operator(/=,/)
00215 #undef Vector_Scalar_Operator
00216 
00217    /* Define the negation operator */
00218    inline TVector operator - () const
00219      { 
00220        TVector ret;
00221        ret.data_[0] = -data_[0];
00222        ret.data_[1] = -data_[1];
00223        ret.data_[2] = -data_[2];
00224        ret.data_[3] = -data_[3];
00225        ret.data_[4] = -data_[4];
00226        ret.data_[5] = -data_[5];
00227        ret.data_[6] = -data_[6];
00228        ret.data_[7] = -data_[7];
00229        ret.data_[8] = -data_[8];
00230        ret.data_[9] = -data_[9];
00231        ret.data_[10] = -data_[10];
00232        ret.data_[11] = -data_[11];
00233        ret.data_[12] = -data_[12];
00234        ret.data_[13] = -data_[13];
00235        ret.data_[14] = -data_[14];
00236        ret.data_[15] = -data_[15];
00237        ret.data_[16] = -data_[16];
00238        return(ret);
00239      }
00240 
00241 #define Vector_Vector_RelOperator(op)                           \
00242    inline int operator op (const TVector &x) const              \
00243      { return(                                                  \
00244                (data_[0] op x.data_[0])                         \
00245             && (data_[1] op x.data_[1])                         \
00246             && (data_[2] op x.data_[2])                         \
00247             && (data_[3] op x.data_[3])                         \
00248             && (data_[4] op x.data_[4])                         \
00249             && (data_[5] op x.data_[5])                         \
00250             && (data_[6] op x.data_[6])                         \
00251             && (data_[7] op x.data_[7])                         \
00252             && (data_[8] op x.data_[8])                         \
00253             && (data_[9] op x.data_[9])                         \
00254             && (data_[10] op x.data_[10])                       \
00255             && (data_[11] op x.data_[11])                       \
00256             && (data_[12] op x.data_[12])                       \
00257             && (data_[13] op x.data_[13])                       \
00258             && (data_[14] op x.data_[14])                       \
00259             && (data_[15] op x.data_[15])                       \
00260             && (data_[16] op x.data_[16])                       \
00261              ) ; }                     
00262 
00263    Vector_Vector_RelOperator(==)
00264    Vector_Vector_RelOperator(!=)
00265 #undef Vector_Vector_RelOperator
00266 
00267 #define Vector_Vector_RelOperator(op)                           \
00268    inline int operator op (const TVector &x) const              \
00269      { return(                                                  \
00270                (data_[0] op x.data_[0])                         \
00271             || (data_[1] op x.data_[1])                         \
00272             || (data_[2] op x.data_[2])                         \
00273             || (data_[3] op x.data_[3])                         \
00274             || (data_[4] op x.data_[4])                         \
00275             || (data_[5] op x.data_[5])                         \
00276             || (data_[6] op x.data_[6])                         \
00277             || (data_[7] op x.data_[7])                         \
00278             || (data_[8] op x.data_[8])                         \
00279             || (data_[9] op x.data_[9])                         \
00280             || (data_[10] op x.data_[10])                       \
00281             || (data_[11] op x.data_[11])                       \
00282             || (data_[12] op x.data_[12])                       \
00283             || (data_[13] op x.data_[13])                       \
00284             || (data_[14] op x.data_[14])                       \
00285             || (data_[15] op x.data_[15])                       \
00286             || (data_[16] op x.data_[16])                       \
00287              ) ; }                     
00288 
00289    Vector_Vector_RelOperator(>)
00290    Vector_Vector_RelOperator(<)
00291    Vector_Vector_RelOperator(>=)
00292    Vector_Vector_RelOperator(<=)
00293 #undef Vector_Vector_RelOperator
00294 
00295 #define Vector_Scalar_RelOperator(op)                           \
00296    inline int operator op (const DataType &v) const             \
00297      { return(                                                  \
00298                (data_[0] op v)                                  \
00299             && (data_[1] op v)                                  \
00300             && (data_[2] op v)                                  \
00301             && (data_[3] op v)                                  \
00302             && (data_[4] op v)                                  \
00303             && (data_[5] op v)                                  \
00304             && (data_[6] op v)                                  \
00305             && (data_[7] op v)                                  \
00306             && (data_[8] op v)                                  \
00307             && (data_[9] op v)                                  \
00308             && (data_[10] op v)                                 \
00309             && (data_[11] op v)                                 \
00310             && (data_[12] op v)                                 \
00311             && (data_[13] op v)                                 \
00312             && (data_[14] op v)                                 \
00313             && (data_[15] op v)                                 \
00314             && (data_[16] op v)                                 \
00315              ); }                     
00316 
00317    Vector_Scalar_RelOperator(==)
00318    Vector_Scalar_RelOperator(!=)
00319 #undef Vector_Scalar_RelOperator
00320 
00321 #define Vector_Scalar_RelOperator(op)                           \
00322    inline int operator op (const DataType &v) const             \
00323      { return(                                                  \
00324                (data_[0] op v)                                  \
00325             || (data_[1] op v)                                  \
00326             || (data_[2] op v)                                  \
00327             || (data_[3] op v)                                  \
00328             || (data_[4] op v)                                  \
00329             || (data_[5] op v)                                  \
00330             || (data_[6] op v)                                  \
00331             || (data_[7] op v)                                  \
00332             || (data_[8] op v)                                  \
00333             || (data_[9] op v)                                  \
00334             || (data_[10] op v)                                 \
00335             || (data_[11] op v)                                 \
00336             || (data_[12] op v)                                 \
00337             || (data_[13] op v)                                 \
00338             || (data_[14] op v)                                 \
00339             || (data_[15] op v)                                 \
00340             || (data_[16] op v)                                 \
00341               ); }                     
00342 
00343    Vector_Scalar_RelOperator(>)
00344    Vector_Scalar_RelOperator(<)
00345    Vector_Scalar_RelOperator(>=)
00346    Vector_Scalar_RelOperator(<=)
00347 #undef Vector_Scalar_RelOperator
00348 
00349    /* Define elementwise minimum and maximum functions for Vectors */
00350    inline void min(const TVector &x) { 
00351      data_[0] = (data_[0] < x.data_[0]) ? data_[0] : x.data_[0]; 
00352      data_[1] = (data_[1] < x.data_[1]) ? data_[1] : x.data_[1]; 
00353      data_[2] = (data_[2] < x.data_[2]) ? data_[2] : x.data_[2]; 
00354      data_[3] = (data_[3] < x.data_[3]) ? data_[3] : x.data_[3]; 
00355      data_[4] = (data_[4] < x.data_[4]) ? data_[4] : x.data_[4]; 
00356      data_[5] = (data_[5] < x.data_[5]) ? data_[5] : x.data_[5]; 
00357      data_[6] = (data_[6] < x.data_[6]) ? data_[6] : x.data_[6]; 
00358      data_[7] = (data_[7] < x.data_[7]) ? data_[7] : x.data_[7]; 
00359      data_[8] = (data_[8] < x.data_[8]) ? data_[8] : x.data_[8]; 
00360      data_[9] = (data_[9] < x.data_[9]) ? data_[9] : x.data_[9]; 
00361      data_[10] = (data_[10] < x.data_[10]) ? data_[10] : x.data_[10]; 
00362      data_[11] = (data_[11] < x.data_[11]) ? data_[11] : x.data_[11]; 
00363      data_[12] = (data_[12] < x.data_[12]) ? data_[12] : x.data_[12]; 
00364      data_[13] = (data_[13] < x.data_[13]) ? data_[13] : x.data_[13]; 
00365      data_[14] = (data_[14] < x.data_[14]) ? data_[14] : x.data_[14]; 
00366      data_[15] = (data_[15] < x.data_[15]) ? data_[15] : x.data_[15]; 
00367      data_[16] = (data_[16] < x.data_[16]) ? data_[16] : x.data_[16]; 
00368    }
00369    inline void max(const TVector &x) {
00370      data_[0] = (data_[0] > x.data_[0]) ? data_[0] : x.data_[0]; 
00371      data_[1] = (data_[1] > x.data_[1]) ? data_[1] : x.data_[1]; 
00372      data_[2] = (data_[2] > x.data_[2]) ? data_[2] : x.data_[2]; 
00373      data_[3] = (data_[3] > x.data_[3]) ? data_[3] : x.data_[3]; 
00374      data_[4] = (data_[4] > x.data_[4]) ? data_[4] : x.data_[4]; 
00375      data_[5] = (data_[5] > x.data_[5]) ? data_[5] : x.data_[5]; 
00376      data_[6] = (data_[6] > x.data_[6]) ? data_[6] : x.data_[6]; 
00377      data_[7] = (data_[7] > x.data_[7]) ? data_[7] : x.data_[7]; 
00378      data_[8] = (data_[8] > x.data_[8]) ? data_[8] : x.data_[8]; 
00379      data_[9] = (data_[9] > x.data_[9]) ? data_[9] : x.data_[9]; 
00380      data_[10] = (data_[10] > x.data_[10]) ? data_[10] : x.data_[10]; 
00381      data_[11] = (data_[11] > x.data_[11]) ? data_[11] : x.data_[11]; 
00382      data_[12] = (data_[12] > x.data_[12]) ? data_[12] : x.data_[12]; 
00383      data_[13] = (data_[13] > x.data_[13]) ? data_[13] : x.data_[13]; 
00384      data_[14] = (data_[14] > x.data_[14]) ? data_[14] : x.data_[14]; 
00385      data_[15] = (data_[15] > x.data_[15]) ? data_[15] : x.data_[15]; 
00386      data_[16] = (data_[16] > x.data_[16]) ? data_[16] : x.data_[16]; 
00387    }
00388 
00389    /* Define minimum and maximum of components */
00390    inline DataType mincomp() const { 
00391      DataType min;
00392      min = data_[0];
00393      if (data_[1] < min) min = data_[1]; 
00394      if (data_[2] < min) min = data_[2]; 
00395      if (data_[3] < min) min = data_[3]; 
00396      if (data_[4] < min) min = data_[4]; 
00397      if (data_[5] < min) min = data_[5]; 
00398      if (data_[6] < min) min = data_[6]; 
00399      if (data_[7] < min) min = data_[7]; 
00400      if (data_[8] < min) min = data_[8]; 
00401      if (data_[9] < min) min = data_[9]; 
00402      if (data_[10] < min) min = data_[10]; 
00403      if (data_[11] < min) min = data_[11]; 
00404      if (data_[12] < min) min = data_[12]; 
00405      if (data_[13] < min) min = data_[13]; 
00406      if (data_[14] < min) min = data_[14]; 
00407      if (data_[15] < min) min = data_[15]; 
00408      if (data_[16] < min) min = data_[16]; 
00409      return min;
00410    }
00411    inline DataType maxcomp() const { 
00412      DataType max;
00413      max = data_[0];
00414      if (data_[1] > max) max = data_[1]; 
00415      if (data_[2] > max) max = data_[2]; 
00416      if (data_[3] > max) max = data_[3]; 
00417      if (data_[4] > max) max = data_[4]; 
00418      if (data_[5] > max) max = data_[5]; 
00419      if (data_[6] > max) max = data_[6]; 
00420      if (data_[7] > max) max = data_[7]; 
00421      if (data_[8] > max) max = data_[8]; 
00422      if (data_[9] > max) max = data_[9]; 
00423      if (data_[10] > max) max = data_[10]; 
00424      if (data_[11] > max) max = data_[11]; 
00425      if (data_[12] > max) max = data_[12]; 
00426      if (data_[13] > max) max = data_[13]; 
00427      if (data_[14] > max) max = data_[14]; 
00428      if (data_[15] > max) max = data_[15]; 
00429      if (data_[16] > max) max = data_[16]; 
00430      return max;
00431    }
00432 
00433    inline double abs() const {
00434      DataType s = data_[0]*data_[0];
00435      s += data_[1]*data_[1]; 
00436      s += data_[2]*data_[2]; 
00437      s += data_[3]*data_[3]; 
00438      s += data_[4]*data_[4]; 
00439      s += data_[5]*data_[5]; 
00440      s += data_[6]*data_[6]; 
00441      s += data_[7]*data_[7]; 
00442      s += data_[8]*data_[8]; 
00443      s += data_[9]*data_[9];
00444      s += data_[10]*data_[10];
00445      s += data_[11]*data_[11];
00446      s += data_[12]*data_[12];
00447      s += data_[13]*data_[13];
00448      s += data_[14]*data_[14];
00449      s += data_[15]*data_[15];
00450      s += data_[16]*data_[16];
00451      return (sqrt(double(s)));
00452    }
00453 
00454    inline TVector getmin(const TVector &x) const
00455      { TVector ret(*this); ret.min(x); return(ret); }
00456    inline TVector getmax(const TVector &x) const
00457      { TVector ret(*this); ret.max(x); return(ret); }
00458 
00459 public:
00460    DataType data_[17];
00461 };
00462 
00463 
00464 template <class DataType>
00465 inline ostream&  operator << (ostream& os, const Vector(17)<DataType> &v) {
00466   os << "(" 
00467      << v.data_[0] 
00468      << "," << v.data_[1] 
00469      << "," << v.data_[2] 
00470      << "," << v.data_[3] 
00471      << "," << v.data_[4] 
00472      << "," << v.data_[5] 
00473      << "," << v.data_[6] 
00474      << "," << v.data_[7] 
00475      << "," << v.data_[8] 
00476      << "," << v.data_[9] 
00477      << "," << v.data_[10] 
00478      << "," << v.data_[11] 
00479      << "," << v.data_[12] 
00480      << "," << v.data_[13] 
00481      << "," << v.data_[14] 
00482      << "," << v.data_[15] 
00483      << "," << v.data_[16] 
00484      << ")";
00485   return os;
00486 }
00487 
00488 template <class DataType>
00489 inline Vector(17)<DataType> Min(const Vector(17)<DataType> &x, 
00490                                    const Vector(17)<DataType> &y) 
00491   { Vector(17)<DataType> ret(x); ret.min(y); return(ret); }
00492 
00493 template <class DataType>
00494 inline Vector(17)<DataType> Max(const Vector(17)<DataType> &x, 
00495                                    const Vector(17)<DataType> &y) 
00496   { Vector(17)<DataType> ret(x); ret.max(y); return(ret); }
00497 
00498 template <class DataType>
00499 inline double abs(const Vector(17)<DataType> &x) 
00500   { return (x.abs()); }
00501 
00502 template <class DataType>
00503 inline DataType mincomp(const Vector(17)<DataType> &x) 
00504   { return(x.mincomp()); }
00505 
00506 template <class DataType>
00507 inline DataType maxcomp(const Vector(17)<DataType> &x) 
00508   { return(x.maxcomp()); }
00509 
00510 
00511 #define Global_Vector_Scalar_Operator(op)                                  \
00512 template <class DataType>                                                  \
00513 inline Vector(17)<DataType> operator op (const DataType &v,                 \
00514    const Vector(17)<DataType> &x)                                           \
00515    { return(x op v); }
00516 
00517 Global_Vector_Scalar_Operator(+)
00518 Global_Vector_Scalar_Operator(-)
00519 Global_Vector_Scalar_Operator(*)
00520 Global_Vector_Scalar_Operator(/)
00521 #undef Global_Vector_Scalar_Operator
00522 
00523 
00524 #define Global_Vector_Function(fct)                                       \
00525 template <class DataType>                                                 \
00526 inline Vector(17)<DataType> fct(const Vector(17)<DataType> &x)              \
00527    {                                                                      \
00528       Vector(17)<DataType> ret;                                            \
00529       ret[0] = fct(x[0]);                                                 \
00530       ret[1] = fct(x[1]);                                                 \
00531       ret[2] = fct(x[2]);                                                 \
00532       ret[3] = fct(x[3]);                                                 \
00533       ret[4] = fct(x[4]);                                                 \
00534       ret[5] = fct(x[5]);                                                 \
00535       ret[6] = fct(x[6]);                                                 \
00536       ret[7] = fct(x[7]);                                                 \
00537       ret[8] = fct(x[8]);                                                 \
00538       ret[9] = fct(x[9]);                                                 \
00539       ret[10] = fct(x[10]);                                               \
00540       ret[11] = fct(x[11]);                                               \
00541       ret[12] = fct(x[12]);                                               \
00542       ret[13] = fct(x[13]);                                               \
00543       ret[14] = fct(x[14]);                                               \
00544       ret[15] = fct(x[15]);                                               \
00545       ret[16] = fct(x[16]);                                               \
00546       return(ret);                                                        \
00547      }
00548 
00549 Global_Vector_Function(fabs)
00550 Global_Vector_Function(sqrt)
00551 #undef Global_Vector_Function
00552 
00553 
00554 #define Vector_Vector_Functions(NameTo,NameFrom,ope)                      \
00555   template <class DataType, class VectorType>                             \
00556   inline void NameTo (Vector(17)<DataType> &x, const VectorType &y)        \
00557     {                                                                     \
00558       x(0) ope y(0);                                                      \
00559       x(1) ope y(1);                                                      \
00560       x(2) ope y(2);                                                      \
00561       x(3) ope y(3);                                                      \
00562       x(4) ope y(4);                                                      \
00563       x(5) ope y(5);                                                      \
00564       x(6) ope y(6);                                                      \
00565       x(7) ope y(7);                                                      \
00566       x(8) ope y(8);                                                      \
00567       x(9) ope y(9);                                                      \
00568       x(10) ope y(10);                                                    \
00569       x(11) ope y(11);                                                    \
00570       x(12) ope y(12);                                                    \
00571       x(13) ope y(13);                                                    \
00572       x(14) ope y(14);                                                    \
00573       x(15) ope y(15);                                                    \
00574       x(16) ope y(16);                                                    \
00575     }                                                                     \
00576   template <class DataType, class VectorType>                             \
00577   inline void NameFrom (VectorType &x, const Vector(17)<DataType> &y)      \
00578     {                                                                     \
00579       x(0) ope y(0);                                                      \
00580       x(1) ope y(1);                                                      \
00581       x(2) ope y(2);                                                      \
00582       x(3) ope y(3);                                                      \
00583       x(4) ope y(4);                                                      \
00584       x(5) ope y(5);                                                      \
00585       x(6) ope y(6);                                                      \
00586       x(7) ope y(7);                                                      \
00587       x(8) ope y(8);                                                      \
00588       x(9) ope y(9);                                                      \
00589       x(10) ope y(10);                                                    \
00590       x(11) ope y(11);                                                    \
00591       x(12) ope y(12);                                                    \
00592       x(13) ope y(13);                                                    \
00593       x(14) ope y(14);                                                    \
00594       x(15) ope y(15);                                                    \
00595       x(16) ope y(16);                                                    \
00596     }
00597 
00598    Vector_Vector_Functions(equals_to,equals_from,=)
00599    Vector_Vector_Functions(plus_to,plus_from,+=)
00600    Vector_Vector_Functions(minus_to,minus_from,-=)
00601    Vector_Vector_Functions(multiply_to,multiply_from,*=)
00602    Vector_Vector_Functions(divide_to,divide_from,/=)
00603 #undef Vector_Vector_Functions
00604 
00605 
00606 #endif


Quickstart     Users Guide     Programmers Reference     Installation      Examples     Download



AMROC Main      Home      Contact
last update: 06/01/04