Changeset 5103

Timestamp:

Mar 8, 2015, 9:34:59 PM (9 years ago)

Author:

mmh

Message:

uq work snapshot

Location:

branches/uq

Files:

• gui/scripts/probdisteditor.tcl (modified) (4 diffs)
• gui/scripts/uq.tcl (modified) (1 diff)
• puq/puq_analyze.py (modified) (2 diffs)

branches/uq/gui/scripts/probdisteditor.tcl

@@ -326 +326 @@
         uniform {
             if {![info exists _value(min)] || ![info exists _value(max)]} {
-                puts "UNIFORM min=$min max$max"
+                # puts "UNIFORM min=$min max=$max"
                 if {$min != "" && $max != ""} {
                     set _value(min) $min
@@ -342 +342 @@
         }
         gaussian {
+            # puts "GAUSSIAN min=$min max=$max"
             if {![info exists _value(central)]} {
                 if {[info exists _value(min)] && [info exists _value(max)]} {
@@ -349 +350 @@
                 }
             }
+            # puts "central=$_value(central)"
             if {![info exists _value(stddev)]} {
                 if {[info exists _value(min)] && [info exists _value(max)]} {
@@ -355 +357 @@
                     set _value(stddev) [expr {0.3*($max - $min)}]
                 } else {
-                    set _value(stddev)  [expr _value(central) * 0.10]
-                }
-            }
-
+                    set _value(stddev)  [expr $_value(central) * 0.10]
+                }
+            }
+            if {$_value(stddev) <= 0} {
+                set _value(stddev) 1.0
+            }
+            # puts "stddev=$_value(stddev)"
             # lower bound is -3 deviations or tool min
             set trunc [expr {$_value(central) - 3*$_value(stddev)}]

branches/uq/gui/scripts/uq.tcl

@@ -53 +53 @@
 
         label $inner.type.labeltext -text "Level" -font "Arial 9"
-        Rappture::Spinint $inner.type.level -min 2 -max 20
+        Rappture::Spinint $inner.type.level -min 1 -max 20
         $inner.type.level value $_args
         bind $inner.type.level <<Value>> [itcl::code $this _adjust_level $inner]

branches/uq/puq/puq_analyze.py

@@ -1 +1 @@
 #!/usr/bin/env python
-from puq import *
+import sys
+sys.stdout = open("analyze.out", 'w')
+sys.stderr = open("analyze.err", 'w')
+
+import puq
 import csv
 import numpy as np
-import sys, os, h5py
+import os, h5py
+from puq.jpickle import unpickle
+import xml.etree.ElementTree as xml
 
 
@@ -20 +26 @@
 sw = load_from_hdf5(sys.argv[1])
 sw.analyze()
+
+h5 = h5py.File(sys.argv[1], 'r+')
+
+
+dtree = xml.parse('run_uq.xml')
+droot = dtree.getroot()
+dout = droot.find('output')
+
+uqtype = h5.attrs['UQtype']
+for v in h5[uqtype]:
+    print 'PDF', v
+    rs = unpickle(h5['/%s/%s/response' % (uqtype, v)].value)
+    pdf = rs.pdf(fit=False)
+
+    id = '%s PDF' % v
+    elem = xml.SubElement(dout, 'curve', {'id': id})
+    about = xml.SubElement(elem, 'about')
+    label = xml.SubElement(about, 'label')
+    label.text = id
+    desc = xml.SubElement(about, 'description')
+    desc.text = "PDF for %s" % id
+    yaxis = xml.SubElement(elem, 'yaxis')
+    label = xml.SubElement(yaxis, 'label')
+    label.text = 'Probability'
+
+    component = xml.SubElement(elem, 'component')
+    xy = xml.SubElement(component, 'xy')
+    pts = "%s 0\n" % pdf.x[0]
+    for x, y in zip(pdf.x, pdf.y):
+        pts += "%s %s\n" % (x, y)
+    pts += "%s 0\n" % pdf.x[-1]
+    xy.text = pts
+
+with open('run_uq.xml', 'w') as f:
+    f.write("<?xml version=\"1.0\"?>\n")
+    dtree.write(f)
+
+
+h5.close()