Implemented starting configuration.

[phys/u1casc.git] / u1casc-ordinary / sim.hpp

diff --git a/u1casc-ordinary/sim.hpp b/u1casc-ordinary/sim.hpp
index f88e650d6c2bc407b5aa44565a431040ed5e7591..986ae371fe10f54b8a041848e8e401e4b735745e 100644 (file)
--- a/u1casc-ordinary/sim.hpp
+++ b/u1casc-ordinary/sim.hpp
@@ -18,10 +18,11 @@ public:
   neigh *nb;
   complex<double> *U, *phi;
   double kappa[2], lambda[2], beta;
   neigh *nb;
   complex<double> *U, *phi;
   double kappa[2], lambda[2], beta;

-
+  void _newParas();
+  

 private:
 private:

+  void _resetConfig();

   void _makeSweep();
   void _makeSweep();

-  void _newParas();

   gsl_rng* rangsl;
   double rhoPhi(const int& iphi, const int& x0, const complex<double>& candPhi);
   double rhoU(const int& x0, const int& nu0, const complex<double>& candU);
   gsl_rng* rangsl;
   double rhoPhi(const int& iphi, const int& x0, const complex<double>& candPhi);
   double rhoU(const int& x0, const int& nu0, const complex<double>& candU);


@@ -60,7 +61,9 @@ void sim::_newParas() {
   lambda[0] = (*O815->paraQ)["lambdaone"];
   lambda[1] = (*O815->paraQ)["lambdatwo"];
   beta = (*O815->paraQ)["beta"];
   lambda[0] = (*O815->paraQ)["lambdaone"];
   lambda[1] = (*O815->paraQ)["lambdatwo"];
   beta = (*O815->paraQ)["beta"];

+}

 
 

+void sim::_resetConfig() {

   for(int ix=0; ix<lsize4; ix++) {
     for(int i=0; i<2; i++) phi[ i*lsize4 + ix ] = 0;
     for(int nu=0; nu<4; nu++) U[ ix*4 + nu ] = 1;
   for(int ix=0; ix<lsize4; ix++) {
     for(int i=0; i<2; i++) phi[ i*lsize4 + ix ] = 0;
     for(int nu=0; nu<4; nu++) U[ ix*4 + nu ] = 1;