TDSE.h

#ifndef qprop_h
#define qprop_h qprop_h
#include<iostream>
#include<complex>
#include<cmath>
#include<fstream>
#include<string>


using namespace std;

class wavefunction;
class fluid;
class grid;
class hamop;
double alpha_y(double time, int me);
double alpha_x(double time, int me);
double vecpot_x(double time, int me);
double vecpot_y(double time, int me);
double vecpot_z(double time, int me);
double scalarpotx(double x, double y, double z, double time, int me);
double scalarpoty(double x, double y, double z, double time, int me);
double scalarpotz(double x, double y, double z, double time, int me);
double interactionpotxy(double x, double y, double z, double time, int me);
double imagpotx(long xindex, long yindex, long zindex, double time, grid g);
double imagpoty(long xindex, long yindex, long zindex, double time, grid g);
double field(double time, int me);
double dftpot(grid g, double x, double y, double z, double time, int me, 
	      const fluid &wf_one, const wavefunction &wf_two);
double scalarpotxtwo(double x, double y, double z, double time, int me);
double scalarpotytwo(double x, double y, double z, double time, int me);
double scalarpotztwo(double x, double y, double z, double time, int me);
double interactionpotxytwo(double x, double y, double z, double time, int me);
double imagpotxtwo(long xindex, long yindex, long zindex, double time, grid g);
double imagpotytwo(long xindex, long yindex, long zindex, double time, grid g);

double imagpotxbig(long xindex, long yindex, long zindex, double time, grid gbig);
double imagpotybig(long xindex, long yindex, long zindex, double time, grid gbig);




double dfthartree(grid g, double x, double y, double z, double time, int me,const wavefunction & wf, double eps, long box);

#endif // qprop_h