for j, const in enumerate(constraints):
dg_j = ["ssa.dG[%d][%d]*%s" % (j, i, x)
for i, x in enumerate(delta_z)]
constraint_parts = ["ssa.G[%d]" % j]
constraint_parts.extend(dg_j)
lin_constraint = "%s < 0" % "+".join(constraint_parts)
sysopt.add_constraint(lin_constraint)
#self.parent.driver.workflow.add("mda")
self.parent.driver.workflow.add(sa_s)
if global_dvs:
self.parent.driver.workflow.add("ssa")
self.parent.driver.workflow.add(bbopts)
if global_dvs:
self.parent.driver.workflow.add("sysopt")
if __name__ == "__main__":
from openmdao.lib.optproblems.api import SellarProblem
#from openmdao.main.api import ArchitectureAssembly
sp = SellarProblem()
sp.architecture = BLISS()
sp.run()
for k, v in sp.check_solution().iteritems():
print " ", k, ": ", v
if __name__ == "__main__": # pragma: no cover
import time
from openmdao.main.api import set_as_top
prob = SellarProblem()
prob.architecture = MDF()
prob.configure()
prob.driver.recorders = [DBCaseRecorder()]
tt = time.time()
prob.run()
error = prob.check_solution()
print "\nUsing MDF Architecture"
print "Minimum found at (%f, %f, %f)" % (prob.dis1.z1, prob.dis1.z2,
prob.dis1.x1)
print "Minimum differs from expected by (%f, %f, %f)" % (
error["z1"], error["z2"], error['dis1.x1'])
print "Couping vars: %f, %f" % (prob.dis1.y1, prob.dis2.y2)
print "Minimum objective: ", prob.solution['obj1']
print "Elapsed time: ", time.time() - tt, "seconds"
print "\n"
prob = SellarProblem()
prob.architecture = BLISS()
prob.configure()
for j, const in enumerate(constraints):
dg_j = [
"ssa.dG[%d][%d]*%s" % (j, i, x) for i, x in enumerate(delta_z)
]
constraint_parts = ["ssa.G[%d]" % j]
constraint_parts.extend(dg_j)
lin_constraint = "%s < 0" % "+".join(constraint_parts)
sysopt.add_constraint(lin_constraint)
#self.parent.driver.workflow.add("mda")
self.parent.driver.workflow.add(sa_s)
if global_dvs:
self.parent.driver.workflow.add("ssa")
self.parent.driver.workflow.add(bbopts)
if global_dvs:
self.parent.driver.workflow.add("sysopt")
if __name__ == "__main__":
from openmdao.lib.optproblems.api import SellarProblem
#from openmdao.main.api import ArchitectureAssembly
sp = SellarProblem()
sp.architecture = BLISS()
sp.run()
for k, v in sp.check_solution().iteritems():
print " ", k, ": ", v
from openmdao.lib.optproblems.api import SellarProblem
if __name__ == "__main__": # pragma: no cover
import time
prob = SellarProblem()
prob.architecture = MDF()
prob.configure()
prob.driver.recorders = [DBCaseRecorder()]
tt = time.time()
prob.run()
error = prob.check_solution()
print "\nUsing MDF Architecture"
print "Minimum found at (%f, %f, %f)" % (prob.dis1.z1,
prob.dis1.z2,
prob.dis1.x1)
print "Minimum differs from expected by (%f, %f, %f)" % (error["z1"],
error["z2"],
error['dis1.x1'])
print "Couping vars: %f, %f" % (prob.dis1.y1, prob.dis2.y2)
print "Minimum objective: ", prob.get_objectives()['obj1'].evaluate()
print "Elapsed time: ", time.time()-tt, "seconds"
print "\n"
prob = SellarProblem()
