Changeset 24313 for issm/trunk/src/m/qmu/setupdesign/QmuSetupVariables.py

issm/trunk

Property svn:mergeinfo changed

merged: 23396,23399-23402,23406,23408-23435,23437-23446,23448,23451,23453,23455,23457,23459-23563,23565-23566,23568-23914,23919-23926,23935-24306,24309-24310

issm/trunk/src

Property svn:mergeinfo changed

/issm/trunk-jpl/src

merged: 23399-23402,23422,23430-23434,23437,23439-23440,23448,23459-23463,23465-23469,23471-23475,23477,23481-23490,23493-23495,23499-23508,23513-23516,23518-23524,23526,23528-23529,23531-23534,23536-23538,23540,23542,23544,23546-23548,23551,23556-23563,23568-23576,23578-23593,23597-23607,23609-23610,23612-23615,23617-23618,23620-23621,23624-23625,23627-23655,23657-23658,23661,23664-23670,23672-23674,23677-23678,23689-23691,23695-23697,23699,23701-23705,23708-23709,23711,23713-23723,23729-23730,23735,23737-23740,23748,23750-23751,23754-23756,23758-23761,23764-23777,23779-23784,23787-23788,23795,23797,23799-23815,23817,23819,23822,23824-23828,23831,23833-23834,23837-23847,23849,23853-23857,23859,23861-23864,23866-23868,23870-23876,23880,23882-23883,23885,23887-23889,23891,23893-23897,23900-23907,23909-23914,23920-23926,23935-23937,23941-23948,23950-23955,23957-23975,23977,23980-23985,23987,23990-24002,24010-24011,24013-24014,24016-24024,24027-24032,24034-24042,24044-24045,24047-24049,24053-24064,24066-24070,24072-24073,24077-24080,24082-24083,24086,24088-24101,24103-24104,24106-24120,24126,24130-24137,24139-24141,24145-24146,24151-24158,24160-24163,24165-24173,24175-24179,24182,24184-24186,24188-24189,24191-24192,24194-24195,24197,24199,24202-24205,24207-24209,24211,24213,24216,24220-24222,24227-24228,24240-24241,24250-24256,24258-24262,24269-24270,24275,24279-24280,24282,24285-24295,24301-24302,24304-24305,24309-24310

issm/trunk/src/m/qmu/setupdesign/QmuSetupVariables.py

-              r23189
+              r24313
 from MatlabFuncs import *
 from uniform_uncertain import*
+from uniform_uncertain import *
 from normal_uncertain import *
 from copy import deepcopy
-def QmuSetupVariables(md,dvar,variables):
+        #get descriptor
+        descriptor=variables.descriptor
+def QmuSetupVariables(md, dvar, variables):
+        #decide whether this is a distributed variable, which will drive whether we expand it into npart values,
+        #or if we just carry it forward as is.
+    #get descriptor
+    descriptor = variables.descriptor
+        #ok, key off according to type of descriptor:
+        if strncmp(descriptor,'scaled_',7):
+                #we have a scaled variable, expand it over the partition.
+    #decide whether this is a distributed variable, which will drive whether we expand it into npart values,
+    #or if we just carry it forward as is.
+                if isinstance(variables,uniform_uncertain):
+                        if ((type(variables.lower) in [list,np.ndarray] and len(variables.lower) > md.qmu.numberofpartitions) or (type(variables.upper) in [list,np.ndarray] and len(variables.upper) > md.qmu.numberofpartitions)):
+                                raise RuntimeError('QmuSetupDesign error message: upper and lower should be either a scalar or a "npart" length vector')
+                elif isinstance(variables,normal_uncertain):
+                        if type(variables.stddev) in [list,np.ndarray] and len(variables.stddev) > md.qmu.numberofpartitions:
+                                raise RuntimeError('QmuSetupDesign error message: stddev should be either a scalar or a "npart" length vector')
+    #ok, key off according to type of descriptor:
+    if strncmp(descriptor, 'scaled_', 7):
+        #we have a scaled variable, expand it over the partition.
+        if isinstance(variables, uniform_uncertain):
+            if ((type(variables.lower) in [list, np.ndarray] and len(variables.lower) > md.qmu.numberofpartitions) or (type(variables.upper) in [list, np.ndarray] and len(variables.upper) > md.qmu.numberofpartitions)):
+                raise RuntimeError('QmuSetupDesign error message: upper and lower should be either a scalar or a "npart" length vector')
+                #ok, dealing with semi-discrete distributed variable. Distribute according to how many
+                #partitions we want
+        elif isinstance(variables, normal_uncertain):
+            if type(variables.stddev) in [list, np.ndarray] and len(variables.stddev) > md.qmu.numberofpartitions:
+                raise RuntimeError('QmuSetupDesign error message: stddev should be either a scalar or a "npart" length vector')
+                for j in range(md.qmu.numberofpartitions):
+                        dvar.append(deepcopy(variables))
+                        # "'" is because qmu.in files need for strings to be in actual ''
+                        # must also account for whether we are given 1 instance or an array of instances
+        #ok, dealing with semi - discrete distributed variable. Distribute according to how many
+        #partitions we want
+        for j in range(md.qmu.numberofpartitions):
+            dvar.append(deepcopy(variables))
+            # "'" is because qmu.in files need for strings to be in actual ''
+            # must also account for whether we are given 1 instance or an array of instances
+            # handle descriptors for everything
+            if type(dvar[-1].descriptor) in [list, np.ndarray] and len(variables.descriptor) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                if type(variables.descriptor) == np.ndarray:
+                    dvar[-1].descriptor = np.append(dvar[-1].descriptor, "'" + str(variables.descriptor) + '_' + str(j + 1) + "'")
+                else:
+                    dvar[-1].descriptor.append("'" + str(variables.descriptor) + '_' + str(j + 1) + "'")
+            else:
+                dvar[-1].descriptor = "'" + str(variables.descriptor) + '_' + str(j + 1) + "'"
+                        # handle descriptors for everything
+                        if type(dvar[-1].descriptor) in [list,np.ndarray] and len(variables.descriptor) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                                if type(variables.descriptor) == np.ndarray:
+                                        dvar[-1].descriptor = np.append(dvar[-1].descriptor,"'"+str(variables.descriptor)+'_'+str(j+1)+"'")
+                                else:
+                                        dvar[-1].descriptor.append("'"+str(variables.descriptor)+'_'+str(j+1)+"'")
+                        else:
+                                dvar[-1].descriptor = "'"+str(variables.descriptor)+'_'+str(j+1)+"'"
+    # handle uniform_uncertain
+            if isinstance(variables, uniform_uncertain):
+                if type(variables.lower) in [list, np.ndarray] and len(variables.lower) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                    if type(variables.lower) == np.ndarray:
+                        dvar[-1].lower = np.append(dvar[-1].lower, variables.lower[j])
+                    else:
+                        dvar[-1].lower.append(variables.lower[j])
+                else:
+                    dvar[-1].lower = variables.lower
+                        # handle uniform_uncertain
+                        if isinstance(variables,uniform_uncertain):
+                                if type(variables.lower) in [list,np.ndarray] and len(variables.lower) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                                        if type(variables.lower) == np.ndarray:
+                                                dvar[-1].lower = np.append(dvar[-1].lower, variables.lower[j])
+                                        else:
+                                                dvar[-1].lower.append(variables.lower[j])
+                                else:
+                                        dvar[-1].lower = variables.lower
+                if type(variables.upper) in [list, np.ndarray] and len(variables.upper) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                    if type(variables.upper) == np.ndarray:
+                        dvar[-1].upper = np.append(dvar[-1].upper, variables.upper[j])
+                    else:
+                        dvar[-1].upper.append(variables.upper[j])
+                else:
+                    dvar[-1].upper = variables.upper
+                                if type(variables.upper) in [list,np.ndarray] and len(variables.upper) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                                        if type(variables.upper) == np.ndarray:
+                                                dvar[-1].upper = np.append(dvar[-1].upper, variables.upper[j])
+                                        else:
+                                                dvar[-1].upper.append(variables.upper[j])
+                                else:
+                                        dvar[-1].upper = variables.upper
+    # handle normal_uncertain
+            elif isinstance(variables, normal_uncertain):
+                if type(variables.stddev) in [list, np.ndarray] and len(variables.stddev) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                    if type(variables.stddev) == np.ndarray:
+                        dvar[-1].stddev = np.append(dvar[-1].stddev, variables.stddev[j])
+                    else:
+                        dvar[-1].stddev.append(variables.stddev[j])
+                else:
+                    dvar[-1].stddev = variables.stddev
+                        # handle normal_uncertain
+                        elif isinstance(variables,normal_uncertain):
+                                if type(variables.stddev) in [list,np.ndarray] and len(variables.stddev) > 1 and len(variables.upper) != md.qmu.numberofpartitions:
+                                        if type(variables.stddev) == np.ndarray:
+                                                dvar[-1].stddev = np.append(dvar[-1].stddev, variables.stddev[j])
+                                        else:
+                                                dvar[-1].stddev.append(variables.stddev[j])
+                                else:
+                                        dvar[-1].stddev = variables.stddev
+    # running with a single instance, and therefore length 1 arrays of qmu classes
+    else:
+        dvar.append(variables)
+    # text parsing in dakota requires literal "'identifier'" not just "identifier"
+        for v in dvar:
+            v.descriptor = "'" + str(v.descriptor) + "'"
+        # running with a single instance, and therefore length 1 arrays of qmu classes
+        else:
+                dvar.append(variables)
+                # text parsing in dakota requires literal "'identifier'" not just "identifier"
+                for v in dvar:
+                        v.descriptor = "'"+str(v.descriptor)+"'"
+        return dvar
+    return dvar

Context Navigation

Legend:

issm/trunk

issm/trunk/src

issm/trunk/src/m/qmu/setupdesign/QmuSetupVariables.py

Download in other formats: