QMCHarmonicOscillator.h

Go to the documentation of this file.

//            QMcBeaver
//
//         Constructed by 
//
//     Michael Todd Feldmann 
//              and 
//   David Randall "Chip" Kent IV
//
// Copyright 2000.  All rights reserved.
//
// drkent@users.sourceforge.net mtfeldmann@users.sourceforge.net

#ifndef QMCHarmonicOscillator_H
#define QMCHarmonicOscillator_H

#include <iostream>
#include <sstream>

#include "QMCInput.h"
#include "QMCDouble.h"
#include "QMCConfigIO.h"
#include "IeeeMath.h"
#include "QMCWalkerData.h"
#include "QMCFunctions.h"

using namespace std;

class QMCHarmonicOscillator : public QMCFunctions
{
public:
  QMCHarmonicOscillator() {};

  QMCHarmonicOscillator(QMCInput *input);

  QMCHarmonicOscillator(const QMCHarmonicOscillator & rhs )
    {
      *this = rhs;
    }

  ~QMCHarmonicOscillator() {};

  void evaluate(Array1D<QMCWalkerData *> &walkerData,
                Array1D<Array2D<double> * > &xData,
                int num);

  void evaluate(Array2D<double> &X, QMCWalkerData & data)
    {
      Array1D< Array2D<double>* > tempX;
      tempX.allocate(1);
      tempX(0) = &X;
      Array1D< QMCWalkerData * > tempWD;
      tempWD.allocate(1);
      tempWD(0) = &data;
      evaluate(tempWD,tempX,1);
    }

  void operator=(const QMCHarmonicOscillator & rhs )
    {
      Input = rhs.Input;
    }

 private:
  double a, w, w2div2, a2;

};

#endif

---