Changeset 561

Timestamp:

09/17/12 10:29:41 (4 years ago)

Author:

mmckerns

Message:

added option to scale color intensity for weights in hypercube plot

File:

• mystic/scripts/support_hypercube_measures.py (modified) (4 diffs)

mystic/scripts/support_hypercube_measures.py

@@ -59 +59 @@
                     metavar="INT",default=None,
                     help="id # of the nth simultaneous points to plot")
+  parser.add_option("-s","--scale",action="store",dest="scale",\
+                    metavar="FLOAT",default=1.0,
+                    help="multiplier for scaling color of weights in plot")
   parser.add_option("-f","--flat",action="store_true",dest="flatten",\
                     default=False,help="show selected iterations in a single plot")
@@ -114 +117 @@
   except:
     id = None # i.e. 'all' **or** use id=0, which should be 'best' energy ?
+
+  try: # scale the color in plotting the weights
+    scale = float(parsed_opts.scale)
+  except:
+    scale = 1.0 # color = color**scale
 
   # ensure all terms of bounds and xyz are tuples
@@ -219 +227 @@
     steps = [list(flatten(steps))]
 
+  # adjust for logarithmic scaling of intensity
+  from numpy import e
+  scale = e**(scale - 1.0)
+
   # dot color is based on a product of weights
   t = []
@@ -227 +239 @@
         for j in eval("[params[q][%s] for q in wxyz[1]]" % s):
           for k in eval("[params[q][%s] for q in wxyz[2]]" % s):
-            t[v].append([str(1.0 - i[q]*j[q]*k[q]) for q in range(len(i))])
+            t[v].append([str((1.0 - i[q]*j[q]*k[q])**scale) for q in range(len(i))])
             if float(t[v][-1][-1]) > 1.0 or float(t[v][-1][-1]) < 0.0:
               raise ValueError, "Weights must be in range 0-1. Check normalization and/or assignment."