Only store u[kx,ky] with kx>=0

1 files changed, 7 insertions, 8 deletions

diff --git a/src/init.c b/src/init.c
index 4865ed9..ca693f9 100644
--- a/src/init.c
+++ b/src/init.c
@@ -26,8 +26,7 @@ int init_random (
       if (kx!=0 || ky>0){
 	x=-0.5+((double) rand())/RAND_MAX;
 	y=-0.5+((double) rand())/RAND_MAX;
-	u0[klookup(kx,ky,2*K1+1,2*K2+1)]=x+y*I;
-	u0[klookup(-kx,-ky,2*K1+1,2*K2+1)]=conj(u0[klookup(kx,ky,2*K1+1,2*K2+1)]);
+	u0[klookup_sym(kx,ky,K2)]=x+y*I;
       }
     }
   }
@@ -39,9 +38,9 @@ int init_random (
     // fix the enstrophy
     rescale=compute_enstrophy(u0, K1, K2, L);
   }
-  for(kx=-K1;kx<=K1;kx++){
+  for(kx=0;kx<=K1;kx++){
     for(ky=-K2;ky<=K2;ky++){
-      u0[klookup(kx,ky,2*K1+1,2*K2+1)]=u0[klookup(kx,ky,2*K1+1,2*K2+1)]*sqrt(init_en/rescale);
+      u0[klookup_sym(kx,ky,K2)]=u0[klookup_sym(kx,ky,K2)]*sqrt(init_en/rescale);
     }
   }
 
@@ -60,9 +59,9 @@ int init_gaussian (
   int kx,ky;
   double rescale;
 
-  for(kx=-K1;kx<=K1;kx++){
+  for(kx=0;kx<=K1;kx++){
     for(ky=-K2;ky<=K2;ky++){
-      u0[klookup(kx,ky,2*K1+1,2*K2+1)]=(kx*kx+ky*ky)*exp(-(kx*kx+ky*ky));
+      u0[klookup_sym(kx,ky,K2)]=(kx*kx+ky*ky)*exp(-(kx*kx+ky*ky));
     }
   }
 
@@ -73,9 +72,9 @@ int init_gaussian (
     // fix the enstrophy
     rescale=compute_enstrophy(u0, K1, K2, L);
   }
-  for(kx=-K1;kx<=K1;kx++){
+  for(kx=0;kx<=K1;kx++){
     for(ky=-K2;ky<=K2;ky++){
-      u0[klookup(kx,ky,2*K1+1,2*K2+1)]=u0[klookup(kx,ky,2*K1+1,2*K2+1)]*sqrt(init_en/rescale);
+      u0[klookup_sym(kx,ky,K2)]=u0[klookup_sym(kx,ky,K2)]*sqrt(init_en/rescale);
     }
   }