Changeset 172 for branches/UQ/TEST.py

Timestamp:

08/07/09 11:28:29 (7 years ago)

Author:

mmckerns

Message:

bugfix in termination condition choice for UQ & parameter sensitivity
added explicit parameter sensitivity case (TEST2b.py)

File:

• branches/UQ/TEST.py (modified) (6 diffs)

branches/UQ/TEST.py

@@ -9 +9 @@
 #XXX: <mpi config goes here>
 
-npop = 100
-maxiter = 2 #XXX: maxiter too small?
-maxfun = 10 #XXX: maxfun too small?
-convergence_tol = 0.05
-crossover = 0.5
-percent_change = 1.0 #XXX: is 0.9 better?
-
+npop = 20
+maxiter = 500
+maxfun = 1e+6
+convergence_tol = 1e-4
+crossover = 0.9
+percent_change = 0.9
 
 #######################################################################
@@ -28 +27 @@
   def model(t):
     """a simple time-dependent model function
-  f(t) = x1*x2 + x3 + 2.0  [yes, it's actually time-independent...]
+  f(t) = x1*x2 + x3 + 2.0
 
   Input:
@@ -84 +83 @@
 def dakota(cost,lb,ub):
   from mystic.differential_evolution import DifferentialEvolutionSolver2
-  from mystic.termination import VTR
+  from mystic.termination import CandidateRelativeTolerance as CRT
   from mystic.strategy import Best1Exp
   from mystic import getch, random_seed, VerboseSow, Sow
@@ -90 +89 @@
   random_seed(123)
 
-  stepmon = VerboseSow()
+ #stepmon = VerboseSow(100)
+  stepmon = Sow()
   evalmon = Sow()
 
@@ -100 +100 @@
   solver.SetEvaluationLimits(maxiter,maxfun)
 
-  solver.Solve(cost,termination=VTR(convergence_tol),strategy=Best1Exp, \
+  tol = convergence_tol
+  solver.Solve(cost,termination=CRT(tol,tol),strategy=Best1Exp, \
                CrossProbability=crossover,ScalingFactor=percent_change, \
                StepMonitor=stepmon, EvaluationMonitor=evalmon)
 
+  print solver.Solution()
   diameter = -solver.bestEnergy / scale
   func_evals = len(evalmon.y)
@@ -152 +154 @@
   print " model: f(t) = x1*x2 + x3 + 2.0"
   print " parameters: ['x1', 'x2', 'x3']"
+  print " lower bounds: %s" % lower_bounds
+  print " upper bounds: %s" % upper_bounds
   print " ..."