# HG changeset patch
# User nikcleju
# Date 1320797607 0
# Node ID c02eb33d2c54db014b2e963f8e7725ef73d6f1b4
# Parent  2dd78e37b23a10565c6e574544f98865fafbf34b
Prepare to run whole script

diff -r 2dd78e37b23a -r c02eb33d2c54 scripts/ABSapprox.py
--- a/scripts/ABSapprox.py	Wed Nov 09 00:11:14 2011 +0000
+++ b/scripts/ABSapprox.py	Wed Nov 09 00:13:27 2011 +0000
@@ -57,15 +57,15 @@
   #Set up experiment parameters
   d = 50;
   sigma = 2.0
-  #deltas = np.arange(0.05,0.95,0.05)
-  #rhos = np.arange(0.05,0.95,0.05)
-  deltas = np.array([0.05,0.95])
-  rhos = np.array([0.05,0.95])
+  deltas = np.arange(0.05,0.95,0.05)
+  rhos = np.arange(0.05,0.95,0.05)
+  #deltas = np.array([0.05,0.95])
+  #rhos = np.array([0.05,0.95])
   #deltas = np.array([0.05])
   #rhos = np.array([0.05])
   #delta = 0.8;
   #rho   = 0.15;
-  numvects = 10; # Number of vectors to generate
+  numvects = 100; # Number of vectors to generate
   SNRdb = 20.;    # This is norm(signal)/norm(noise), so power, not energy
   # Values for lambda
   #lambdas = [0 10.^linspace(-5, 4, 10)];
@@ -119,14 +119,14 @@
       print "Oops, Type Error"
       raise    
   # Show
-  for algotuple in algosN:
-    plt.figure()
-    plt.imshow(meanmatrix[algotuple[1]], cmap=cm.gray, interpolation='nearest')
-  for algotuple in algosL:
-    for ilbd in np.arange(lambdas.size):
-      plt.figure()
-      plt.imshow(meanmatrix[algotuple[1]][ilbd], cmap=cm.gray, interpolation='nearest')
-  plt.show()
+  #  for algotuple in algosN:
+  #    plt.figure()
+  #    plt.imshow(meanmatrix[algotuple[1]], cmap=cm.gray, interpolation='nearest')
+  #  for algotuple in algosL:
+  #    for ilbd in np.arange(lambdas.size):
+  #      plt.figure()
+  #      plt.imshow(meanmatrix[algotuple[1]][ilbd], cmap=cm.gray, interpolation='nearest')
+  #  plt.show()
   print "Finished."
   
 def genData(d,sigma,delta,rho,numvects,SNRdb):