Changeset 5931 for branches/1.4

Timestamp:

Nov 3, 2015 9:00:05 AM (8 years ago)

Author:

mmh

Message:

use varnames as node names instead of labels for pdfs

File:

• branches/1.4/puq/puq_analyze.py (modified) (7 diffs)

branches/1.4/puq/puq_analyze.py

@@ -52 +52 @@
     yp = np.empty(len(xvals[vname]))
     ym = np.empty(len(xvals[vname]))
+
     for vindex in sorted(xvals[vname].keys()):
         if label is None:
@@ -130 +131 @@
 
 def plot_pdf(io, v, pdf, desc):
-    p = io['output.curve(%s)' % desc]
+    print("plot_pdf %s desc=%s" % (v, desc))
+    p = io['output.curve(pdf-%s)' % v]
     p['about.label'] = desc
     p['about.uqtype'] = "PDF"
@@ -145 +147 @@
     uqtype = h5.attrs['UQtype']
     for v in h5[uqtype]:
+        print("write_responses", v)
         if '[' in v:
             # It is a curve. Ignore.
@@ -158 +161 @@
 
         rsp = h5['/%s/%s/response' % (uqtype, v)].value
-
-        rout = io['output.response(%s)' % label]
+        rout = io['output.response(%s)' % v]
         rout['value'] = rsp
         rout['about.description'] = desc
@@ -272 +274 @@
 uqtype = h5.attrs['UQtype']
 for v in h5[uqtype]:
+    print('v=', v)
     rsp = h5['/%s/%s/response' % (uqtype, v)].value
     rs = unpickle(rsp)
@@ -281 +284 @@
         matches = reg1.findall(v)
         vname, vindex = matches[0]
+        print('CURVE: vname=%s   vindex=%s' % (vname, vindex))
         vindex = int(vindex)
         if vname not in pcurves:
@@ -316 +320 @@
 write_responses(io, h5)
 write_summary(io, h5)
+
 io.close()
-
 h5.close()