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()