kvol.h

Go to the documentation of this file.

// This file os part of FVM
// Copyright (c) 2012 FVM Authors
// See LICENSE file for terms.

#ifndef _KVOL_H_
#define _KVOL_H_

#include <vector>
#include "Vector.h"
#include "Field.h"
#include "pmode.h"

template<class T>
class kvol
{
  
 public:

  typedef kvol<T> Tkvol;
  typedef typename NumTypeTraits<T>::T_Scalar T_Scalar;
  typedef Vector<T_Scalar,3> Tvec;
  typedef pmode<T> Tmode;
  typedef shared_ptr<Tmode> Tmodeptr;
  typedef vector<Tmodeptr> Modes;

 kvol(Tmodeptr mode,T dk3):     //used in gray approx.
  _dk3(dk3),
    _modenum(1),
    _modes(1,mode)
      {}
  
  kvol()
    {}
  
 kvol(const int modes):
  _modenum(modes)
    {}
  
  Tvec getkvec() {return _Kvector;}
  void setkvec(Tvec K) {_Kvector=K;}
  void setdk3(T dk3) {_dk3=dk3;}
  T getdk3() {return _dk3;}
  int getmodenum() {return _modenum;}
  Tmode& getmode(int n) const {return *_modes[n];}
  Modes& getModes() {return _modes;}
  Tkvol& operator=(Tkvol& o)
    {
      const int m=o.getmodenum();
      _modes.clear();
      _modes.resize(m);
      for(int i=0;i<m;i++)
        {
          Tmodeptr newPmode=Tmodeptr(new Tmode());
          _modes[i]=newPmode;
          (*(_modes[i]))=o.getmode(i);
        }
      return *this;
    }
  void copyKvol(Tkvol& inKvol)
  {
    const int m=inKvol.getmodenum();
    _modenum=m;
    _dk3=inKvol.getdk3();
    _modes.clear();
    for(int i=0;i<m;i++)
      {
        Tmodeptr newPmode=Tmodeptr(new Tmode());
        newPmode->copyPmode(inKvol.getmode(i));
        _modes.push_back(newPmode);
      }
  }
  
 private:

  kvol(const kvol&);

  //weight factor (1/m3)
  T _dk3;

  //K vector
  Tvec _Kvector;

  //number of modes
  int _modenum;

  //phonon modes
  Modes _modes;
  
};

#endif