Still a problem with initializing

XanaduAI · Jan 14, 2020 · c66f3c1 · c66f3c1
1 parent 9f0fb80
commit c66f3c1
Show file tree

Hide file tree

Showing 2 changed files with 26 additions and 12 deletions.
diff --git a/examples/Makefile b/examples/Makefile
@@ -25,5 +25,11 @@ example.o: example.cpp
 example-cpp: example.o
 	$(CC) $^ $(LFLAGS) -o $@
 
+example_hermite.o: example_hermite.cpp
+	$(CC) $^ $(CFLAGS) -c
+
+example_hermite.out: example_hermite.o
+	$(CC) $^ $(LFLAGS) -o $@
+
 clean:
 	rm -rf *~ *.out *.o *.so *.pyc *.mod
diff --git a/include/hermite_multidimensional.hpp b/include/hermite_multidimensional.hpp
@@ -105,13 +105,12 @@ inline T* hermite_multidimensional_cpp(const std::vector<T> &R, const std::vecto
     ullint Hdim = pow(resolution, dim);
     T *H;
     H = (T*) malloc(sizeof(T)*Hdim);
-
-
+    T val = 0;
+    memset(H,sizeof(val),sizeof(H));
     H[0] = 1;
 
     std::vector<int> nextPos(dim, 1);
     std::vector<int> jumpFrom(dim, 1);
-    std::vector<double> factors(resolution+1, 0);
     int jump = 0;
     int k;
     ullint nextCoordinate, fromCoordinate;
@@ -164,15 +163,18 @@ inline T*  renorm_hermite_multidimensional_cpp(const std::vector<T> &R, const st
     ullint Hdim = pow(resolution, dim);
     T *H;
     H = (T*) malloc(sizeof(T)*Hdim);
-
+    T val = 0;
+    memset(H,sizeof(val),sizeof(H));	//std::cout<<H[0];
+    //std::cout<<H[1];
+    //std::cout<<H[2];
+
     H[0] = 1;
     std::vector<double> intsqrt(resolution+1, 0);
     for (int ii = 0; ii<=resolution; ii++) {
         intsqrt[ii] = std::sqrt((static_cast<double>(ii)));
     }
     std::vector<int> nextPos(dim, 1);
     std::vector<int> jumpFrom(dim, 1);
-    std::vector<double> factors(resolution+1, 0);
     int jump = 0;
     int k;
     ullint nextCoordinate, fromCoordinate;
@@ -216,15 +218,15 @@ inline T* interferometer_cpp(const std::vector<T> &R, const int &resolution) {
     ullint Hdim = pow(resolution, dim);
     T *H;
     H = (T*) malloc(sizeof(T)*Hdim);
-
+    T val = 0;
+    memset(H,sizeof(val),sizeof(H));
     H[0] = 1;
     std::vector<double> intsqrt(resolution+1, 0);
     for (int ii = 0; ii<=resolution; ii++) {
         intsqrt[ii] = std::sqrt((static_cast<double>(ii)));
     }
     std::vector<int> nextPos(dim, 1);
     std::vector<int> jumpFrom(dim, 1);
-    std::vector<double> factors(resolution+1, 0);
     int jump = 0;
     int k;
     ullint nextCoordinate, fromCoordinate;
@@ -285,7 +287,8 @@ inline T* squeezing_cpp(const std::vector<T> &R, const int &resolution) {
     ullint Hdim = pow(resolution, dim);
     T *H;
     H = (T*) malloc(sizeof(T)*Hdim);
-
+    T val = 0;
+    memset(H,sizeof(val),sizeof(H));
     H[0] = std::sqrt(-R[1]);
     std::vector<double> intsqrt(resolution+1, 0);
     for (int ii = 0; ii<=resolution; ii++) {
@@ -340,15 +343,15 @@ inline T* displacement_cpp(const std::vector<T> &y, const int &resolution) {
     ullint Hdim = pow(resolution, dim);
     T *H;
     H = (T*) malloc(sizeof(T)*Hdim);
-
+    T val = 0;
+    memset(H,sizeof(val),sizeof(H));
     H[0] = std::exp(0.5*y[0]*y[1]);
     std::vector<double> intsqrt(resolution+1, 0);
     for (int ii = 0; ii<=resolution; ii++) {
         intsqrt[ii] = std::sqrt((static_cast<double>(ii)));
     }
     std::vector<int> nextPos(dim, 1);
     std::vector<int> jumpFrom(dim, 1);
-    std::vector<double> factors(resolution+1, 0);
     int jump = 0;
     int k;
     ullint nextCoordinate, fromCoordinate;
@@ -399,6 +402,13 @@ inline T* two_mode_squeezing_cpp(const std::vector<T> &R, const int &resolution)
     ullint Hdim = pow(resolution, dim);
     T *H;
     H = (T*) malloc(sizeof(T)*Hdim);
+    T val = 0;
+    memset(H,sizeof(val),sizeof(H));
+    //memset(H, 0x00, sizeof(H));
+    for (ullint jj = 0; jj < Hdim; jj++){
+    	H[jj] = 0.0;
+    	std::cout << H[jj];
+    }
 
     H[0] = -R[2];
     std::vector<double> intsqrt(resolution+1, 0);
@@ -407,8 +417,6 @@ inline T* two_mode_squeezing_cpp(const std::vector<T> &R, const int &resolution)
     }
     std::vector<int> nextPos(dim, 1);
     std::vector<int> jumpFrom(dim, 1);
-    std::vector<int> ek(dim, 0);
-    std::vector<double> factors(resolution+1, 0);
     int jump = 0;
     int k;
     ullint nextCoordinate, fromCoordinate;