Changeset 25031
- Timestamp:
- 06/12/20 09:53:07 (5 years ago)
- Location:
- issm/trunk-jpl/src/m/qmu
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/qmu/expandresponses.py
r24873 r25031 1 from helpers import * 1 2 from QmuSetupResponses import * 2 from helpers import *3 3 4 4 5 5 def expandresponses(md, responses): 6 #EXPANDRESPONSES - expand responses 7 6 8 7 fnames = fieldnames(responses) 9 8 … … 11 10 dresp = OrderedStruct() 12 11 13 for k in fnames:14 v = getattr(responses, k)15 exec('dresp.{} = type(v )()'.format(k))16 for j in range(len(v )):12 for key in fnames: 13 value = getattr(responses, key) 14 exec('dresp.{} = type(value)()'.format(key)) 15 for j in range(len(value)): 17 16 #call setupdesign 18 exec('dresp.{}=QmuSetupResponses(md, dresp.{}, v [j])'.format(k, k))17 exec('dresp.{}=QmuSetupResponses(md, dresp.{}, value[j])'.format(key, key)) 19 18 20 19 return dresp -
issm/trunk-jpl/src/m/qmu/expandvariables.py
r24870 r25031 1 from QmuSetupVariables import *2 1 from helpers import * 3 2 from qmu_classes import * 3 from QmuSetupVariables import * 4 4 5 5 … … 11 11 dvar = OrderedStruct() 12 12 13 for k in fnames:14 v = getattr(variables, k)13 for key in fnames: 14 value = getattr(variables, key) 15 15 16 16 # for linear constraints, just copy 17 if isinstance(v , linear_inequality_constraint) or isinstance(v, linear_equality_constraint):18 exec('dvar.{} = v '.format(k))17 if isinstance(value, linear_inequality_constraint) or isinstance(value, linear_equality_constraint): 18 exec('dvar.{} = value'.format(key)) 19 19 20 20 # for variables, call the setup function 21 21 else: 22 exec('dvar.{} = type(v )()'.format(k))23 for j in range(len(v )):22 exec('dvar.{} = type(value)()'.format(key)) 23 for j in range(len(value)): 24 24 #call setupdesign 25 exec('dvar.{}=QmuSetupVariables(md, v [j])'.format(k, k))25 exec('dvar.{}=QmuSetupVariables(md, value[j])'.format(key, key)) 26 26 return dvar -
issm/trunk-jpl/src/m/qmu/setupdesign/QmuSetupVariables.py
r25028 r25031 32 32 nlower = variables.lower.shape[1] 33 33 nupper = variables.upper.shape[1] 34 if nlower != n part or nupper != npart:34 if nlower != nt or nupper != nt: 35 35 raise RuntimeError('QmuSetupVariables error message: upper and lower fields should have the same number of cols as the number of time steps') 36 36 elif isinstance(variables, normal_uncertain):
Note:
See TracChangeset
for help on using the changeset viewer.