Diffpack Documentation

---

              Search  

---

PtvGen_Type.h

Go to the documentation of this file.

#ifndef PtvDIM
#define PtvDIM 3
#endif








#if PtvDIM == 0
class Ptv(Type),1;
class Ptv(Type),2;
class Ptv(Type),3;
class Ptm(Type);
#else
class Ptv(Type),0;
#endif



class Ptv(Type)

{

#if PtvDIM == 0
  int   n;            
  Type  x[PtvDIM];    
  Type* v;            
                      
  bool  internal;     
#else
  static int n;       
  Type  v[PtvDIM]; 
#endif


#if PtvDIM == 0
  bool allocate (int n_);
  void deallocate ();
  bool pureRedim (int n_);
#else
  bool redim(int n_);
#endif

public:

#if PtvDIM == 0
  Ptv(Type) ();
  Ptv(Type)(int n_);

  Ptv(Type)(Type x_, Type y_);
  Ptv(Type)(Type x_, Type y_, Type z_);
  Ptv(Type)(Type a[], int n_, bool copy = true);
  Ptv(Type)(const Ptv(Type)& X);
  Ptv(Type)(const Ptv(Type),1& X);
  Ptv(Type)(const Ptv(Type),2& X);
  Ptv(Type)(const Ptv(Type),3& X);
  Ptv(Type)(const VecSimple(Type)& v);
 ~Ptv(Type) ();
  bool redim(int n_);
  bool redim(Type a[], int n_, bool copy = false);
#elif PtvDIM == 1
  Ptv(Type) ();
  Ptv(Type)(Type x_);
  Ptv(Type)(Type a[], int n);
  Ptv(Type)(const Ptv(Type)& X);
  Ptv(Type)(const Ptv(Type),0& X);
 ~Ptv(Type) () {}
#elif PtvDIM == 2
  Ptv(Type) ();
  Ptv(Type)(Type x_, Type y_);
  Ptv(Type)(Type a[], int n);
  Ptv(Type)(const Ptv(Type)& X);
  Ptv(Type)(const Ptv(Type),0& X);
 ~Ptv(Type) () {}
#elif PtvDIM == 3
  Ptv(Type) ();
  Ptv(Type)(Type x_, Type y_, Type z_);
  Ptv(Type)(Type a[], int n);
  Ptv(Type)(const Ptv(Type)& X);
  Ptv(Type)(const Ptv(Type),0& X);
 ~Ptv(Type) () {}
#endif

  bool ok() const;
  int  size() const;
  void fill (Type value);

#if PtvDIM == 0
  void fill (Type v1, Type v2);
  void fill (Type v1, Type v2, Type v3);
  void fill (Type a[], int n);
  void fill (const VecSimple(Type)& v_);
  void fill (const Ptv(Type)& X);
#elif PtvDIM == 2
  void fill (Type v1, Type v2);
#elif PtvDIM == 3
  void fill (Type v1, Type v2, Type v3);
#endif

  Type& operator () (int i);
  Type  operator () (int i) const;
  bool indexOk (int i) const;

  
  Ptv(Type)& operator = (const Type& value);
  Ptv(Type)& operator = (const Ptv(Type)& X);
  void                    operator = (const VecSimple(Type)& X);

  friend void setEqual (VecSimple(Type)& V,
                        const Ptv(Type)& X); 

  bool operator == (const Ptv(Type)& X) const;
  bool operator != (const Ptv(Type)& X) const;

  bool operator >  (const Ptv(Type)& y) const;
  bool operator >= (const Ptv(Type)& y) const;
  bool operator <  (const Ptv(Type)& y) const;
  bool operator <= (const Ptv(Type)& y) const;

  bool eq (const Ptv(Type)& y,
           real tolerance = comparison_tolerance) const;
  bool lt (const Ptv(Type)& y,
           real tolerance = comparison_tolerance) const;
  bool le (const Ptv(Type)& y,
           real tolerance = comparison_tolerance) const;
  bool gt (const Ptv(Type)& y,
           real tolerance = comparison_tolerance) const;
  bool ge (const Ptv(Type)& y,
           real tolerance = comparison_tolerance) const;

  Ptv(Type) operator+ (const Ptv(Type)& X) const;
  Ptv(Type) operator- (const Ptv(Type)& X) const;

  Ptv(Type) operator+ (const Type& value) const;
  Ptv(Type) operator- (const Type& value) const;
  Ptv(Type) operator* (const Type& value) const;
  Ptv(Type) operator/ (const Type& value) const;

  void add (Ptv(Type)& X,
            Ptv(Type)& Y);          
  void add (Type a, Ptv(Type)& X,
            Type b, Ptv(Type)& Y);  
  void add (const Ptv(Type)& X);    
  void add (Type a, Ptv(Type)& X);  

  void fillMin (const Ptv(Type)& X);
  void fillMax (const Ptv(Type)& X);

  Type inner (const Ptv(Type)& X) const;  

#if PtvDIM == 0 || PtvDIM == 3
  void prod (const Ptv(Type)& X,
             const Ptv(Type)& Y);   
#endif

  real distance  (const Ptv(Type)& y) const;  
  real distance2 (const Ptv(Type)& y) const;  
  real norm (Norm_type lp) const;                
  real norm () const;                            
  void mult (Type a);
  real sqr  () const;
  Type sum  () const;

  void   printAsIndex (Os os) const;
  String printAsIndex () const;
  void   print (Os os) const;
  void   scan  (Is is);

  Type* getPtr0 ();
  Type* getPtr1 ();

  friend Os& operator << (Os& os, const Ptv(Type)& x);
  friend Is& operator >> (Is& is,       Ptv(Type)& x);


#if PtvDIM == 0
  friend class Ptv(Type),1;
  friend class Ptv(Type),2;
  friend class Ptv(Type),3;
  friend class Ptm(Type);
#else
  friend class Ptv(Type),0;
#endif
};

---