Changeset 5931
- Timestamp:
- Nov 3, 2015, 9:00:05 AM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/1.4/puq/puq_analyze.py
r5829 r5931 52 52 yp = np.empty(len(xvals[vname])) 53 53 ym = np.empty(len(xvals[vname])) 54 54 55 for vindex in sorted(xvals[vname].keys()): 55 56 if label is None: … … 130 131 131 132 def plot_pdf(io, v, pdf, desc): 132 p = io['output.curve(%s)' % desc] 133 print("plot_pdf %s desc=%s" % (v, desc)) 134 p = io['output.curve(pdf-%s)' % v] 133 135 p['about.label'] = desc 134 136 p['about.uqtype'] = "PDF" … … 145 147 uqtype = h5.attrs['UQtype'] 146 148 for v in h5[uqtype]: 149 print("write_responses", v) 147 150 if '[' in v: 148 151 # It is a curve. Ignore. … … 158 161 159 162 rsp = h5['/%s/%s/response' % (uqtype, v)].value 160 161 rout = io['output.response(%s)' % label] 163 rout = io['output.response(%s)' % v] 162 164 rout['value'] = rsp 163 165 rout['about.description'] = desc … … 272 274 uqtype = h5.attrs['UQtype'] 273 275 for v in h5[uqtype]: 276 print('v=', v) 274 277 rsp = h5['/%s/%s/response' % (uqtype, v)].value 275 278 rs = unpickle(rsp) … … 281 284 matches = reg1.findall(v) 282 285 vname, vindex = matches[0] 286 print('CURVE: vname=%s vindex=%s' % (vname, vindex)) 283 287 vindex = int(vindex) 284 288 if vname not in pcurves: … … 316 320 write_responses(io, h5) 317 321 write_summary(io, h5) 322 318 323 io.close() 319 320 324 h5.close()
Note: See TracChangeset
for help on using the changeset viewer.