Changeset 373 for branches/alta/mystic-0.2a1/test_constraints_cvxopt_qp_allsolvers.py

Timestamp:

08/04/10 16:07:50 (6 years ago)

Author:

altafang

Message:

Fixing minor things, adjusting termination conditions on tests

File:

• branches/alta/mystic-0.2a1/test_constraints_cvxopt_qp_allsolvers.py (modified) (4 diffs)

branches/alta/mystic-0.2a1/test_constraints_cvxopt_qp_allsolvers.py

@@ -232 +232 @@
 # Additional parameters
 ndim = 2        # Number of variables in the x vector.
-x0 = [0.5, 0.5]   # Infeasible initial point is ok!
-#x0 = [0.8, 0.3]
-#x0 = [1., 1.]
+#x0 = [0.5, 0.5]   # Feasible
+#x0 = [0.8, 0.3] # Infeasible
+x0 = [1., 1.]  # Infeasible
 npop = 25       # Required for differential evolution solver
 
@@ -589 +589 @@
     from scipy_bfgs import BFGSSolver
     from mystic.termination import VTR
+    from mystic.termination import GradientTermination
+    from mystic.termination import ChangeOverGeneration as COG
     solver = BFGSSolver(ndim)
       
     solver.SetInitialPoints(x0)
     solver.enable_signal_handler()
-    term = VTR()
+    term = VTR() # works well
+    #term = GradientTermination() # works well
+    #term= COG() # works well
     solver.Solve(costfunc, term, constraints=constraints_string, \
                 constraints_method='auglag')
@@ -619 +623 @@
     from scipy_ncg import NCGSolver
     from mystic.termination import VTR
+    from mystic.termination import SolutionImprovement
+    from mystic.termination import ChangeOverGeneration as COG
     solver = NCGSolver(ndim)
       
     solver.SetInitialPoints(x0)
     solver.enable_signal_handler()
-    term = VTR()
-    solver.Solve(costfunc, term, constraints=constraints_string, \
-                constraints_method='auglag')
+    term = VTR() # very slow with [0.8, 0.3]
+    #term = SolutionImprovement() # very slow with [0.8, 0.3]
+    #term = COG() # very slow with [0.8, 0.3]
+    solver.Solve(costfunc, term, constraints=constraints_string, \
+                constraints_method='auglag', disp=True)
     soln = solver.Solution()
 
@@ -1135 +1143 @@
     test_powelldirectional_auglag()
     test_bfgs_auglag() 
-    test_ncg_auglag() # little off. gets stuck with x0=[0.8, 0.3]
+    test_ncg_auglag() # very slow, esp. in the first SUMT iteration. probably
+                       # relies on the hack to not get completely stuck.
     test_cg_auglag()