QMCJastrowElectronElectron.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 QMCJastrowElectronElectron_H
#define QMCJastrowElectronElectron_H

#include "Array1D.h"
#include "Array2D.h"
#include "QMCInput.h"
#include "QMCJastrowParameters.h"
#include "QMCWalkerData.h"

class QMCJastrowElectronElectron
{
public:
  QMCJastrowElectronElectron();

  ~QMCJastrowElectronElectron();
  
  void initialize(QMCInput * input);

  void evaluate( QMCJastrowParameters & JP,
                 QMCWalkerData * wData,
                 Array2D<double> & X);

  double get_p_a_ln(int ai);

  Array2D<double> * get_p2_xa_ln(int ai);

  double get_p3_xxa_ln(int ai);

  double jastrowOnGrid(QMCJastrowParameters & JP,
                       int E,
                       Array2D<double> & R,
                       Array2D<double> & grid,
                       Array1D<double> & integrand);

protected:
  Array1D<double>  p_a;
  Array1D< Array2D<double> > p2_xa;
  Array1D<double> p3_xxa;

private:  
  QMCWalkerData * wd;

  void updateOne( QMCJastrowParameters & JP,
                  Array2D<double> & X);

  void updateAll( QMCJastrowParameters & JP,
                  Array2D<double> & X);
  
  void collectForPair(int el1, int el2,
                      QMCCorrelationFunction *U_Function,
                      Array2D<double> & X,
                      int index, int numP);
};

#endif

---