myprojects-hg-reborn

Fork

#include <blitz/array.h>
using namespace blitz;
using namespace std;

// first a list of the functions using the Blitz++ library
//void rhs_eval (double x, Array<double,1> y, Array<double,1>& dydx);
void rk4_fixed (double& x, Array<double,1>& y, 
                void (*rhs_eval)(double, Array<double,1>, Array<double,1>&), 
                double h);

void rk4_adaptive (double& x, Array<double,1>& y, 
                   void (*rhs_eval)(double, Array<double,1>, Array<double,1>&),
                   double& h, double& t_err, double acc, 
                   double S, int& rept, int maxrept, 
                   double h_min, double h_max, int flag);


void rk4_fixed (double& x, Array<complex<double>,1>& y, 
                void (*rhs_eval)(double, Array<complex<double>,1>, 
                                Array<complex<double>,1>& ,double,
                            double,int ,Array<complex<double>,1>, int), 
                double h, int m,    Array<complex<double>,1> state,
                                      double U, double Jt,int periodic);

void rk4_adaptive (double& x, Array<complex<double>,1>& y, 
          void (*rhs_eval)(double, Array<complex<double>,1>, 
                                Array<complex<double>,1>& ,double,
                            double,int ,Array<complex<double>,1>, int), 
double& h, double& t_err, double acc, 
                   double S, int& rept, int maxrept, 
                   double h_min, double h_max, int flag,
                  int   m, Array<complex<double>,1> state,
     double U, double Jt, int periodic ) ;

void rhs_eval (double x, Array<complex<double>,1> y,
Array<complex<double>,1>& dydx,double U, double Jt,
int m,Array<complex<double>,1> state,int periodic);

void rhs_eval2 (double x, Array<complex<double>,1> y,
Array<complex<double>,1>& dydx,double U, double Jt,
int m,Array<complex<double>,1> state,int periodic);


void calc_lattice(double Jt, double U, Array<complex<double>,1> inistate,vector<double> & Tinteraction,vector<double> & Tkin);


void calc_cent_site(double Jt, double U, Array<complex<double>,1> inistate, vector<double> & Cinteraction,vector<double> & Ckin);

/*
void rk4_fixed2 (double& x, Array<complex<double>,1>& y, 
                void (*rhs_eval2)(double, Array<complex<double>,1>, 
                                Array<complex<double>,1>& ,double,
                            double,int ,Array<complex<double>,1>, int), 
                double h, int m,    Array<complex<double>,1> state,
                                      double U, double Jt,int periodic);

void rk4_adaptive2 (double& x, Array<complex<double>,1>& y, 
          void (*rhs_eval2)(double, Array<complex<double>,1>, 
                                Array<complex<double>,1>& ,double,
                            double,int ,Array<complex<double>,1>, int), 
double& h, double& t_err, double acc, 
                   double S, int& rept, int maxrept, 
                   double h_min, double h_max, int flag,
                  int   m, Array<complex<double>,1> state,
     double U, double Jt, int periodic );
*/



double real_integration(double f[],int NN,double length);
double complex_integration (double f[],int NN,double length);
void my_save(char fileout[80],int NN, double arr[]);
void my_save_vec(char fileout[80],int NN, vector<double> arr);
void my_density(double arr[],int NN,double dens[]);
void FFT(double data[],int NN,double V[],double delta_t,double length );
void nonlinFFT(double data[],int NN,double V[],double delta_t,double length,double G );
void imaginaryFFT(double data[],int NN,double V[],double delta_t,double length,double G );
void renormalization(double g[],int NN,double length);
void renorm_imaginaryFFT(double data[],int NN,double V[],double delta_t,double length,double G );
void imaginary_evol_FFT(double data[],int NN,double V[],double delta_t,double length,double G,double time );
void real_evol_FFT(double data[],int NN,double V[],double delta_t,double length,double G,double time );