0

これは、値の繰り返しとopenMDAOの問題に関する以前の質問への追加です。top.run() で値を計算できるようにしたい。次に、変更を調べます。ただし、以前のアドバイスは機能していないようです。

問題が変更または実行されても、prob.root.unknowns は実際には変更されないことがわかりました。

class TestObj(Component):

    def __init__(self):

        Component.__init__(self)

        self.add_param('x',5.0)
        self.add_output('y',5.0)

    def solve_nonlinear(self,params,unknowns,resids):

        unknowns['y']=params['x']

top = Problem()
root = Group()
top.root=root

root.add('split',IndepVarComp('x',2.0),promotes=['*'])
root.add('test',TestObj(),promotes=['*'])

top.setup()
top.root.list_connections()

print 'starting'
print top['x']
print top['y']
for k,v in top.root.unknowns.iteritems():
    print k
    print v['val']
top.run()

print 'after first solve'
print top['x']
print top['y']
for k,v in top.root.unknowns.iteritems():
    print k
    print v['val']

top['x']=33.0

print 'after second solve'
print top['x']
print top['y']
for k,v in top.root.unknowns.iteritems():
    print k
    print v['val']

top.run()

print 'after third solve'
print top['x']
print top['y']
for k,v in top.root.unknowns.iteritems():
    print k
    print v['val']
4

1 に答える 1

1

vec_wrapperのイテレータは、キーとメタデータ ディクショナリのイテレータを返します。メタデータでアクセスしている「val」属性は、実際には初期値です。

代わりにこれを試してみてください。変更されていることがわかります。

from openmdao.api import Component, Problem, Group, IndepVarComp

class TestObj(Component):

    def __init__(self):

        Component.__init__(self)

        self.add_param('x',5.0)
        self.add_output('y',5.0)

    def solve_nonlinear(self,params,unknowns,resids):

        unknowns['y']=params['x']

top = Problem()
root = Group()
top.root=root

root.add('split',IndepVarComp('x',2.0),promotes=['*'])
root.add('test',TestObj(),promotes=['*'])

top.setup()
top.root.list_connections()

print 'starting'
print top['x']
print top['y']
for k,v in root.unknowns.iteritems():
    print k
    print top.root.unknowns[k]
top.run()

print 'after first solve'
print top['x']
print top['y']
for k,v in root.unknowns.iteritems():
    print k
    print top.root.unknowns[k]

top['x']=33.0

print 'after second solve'
print top['x']
print top['y']
for k,v in top.root.unknowns.iteritems():
    print k
    print top.root.unknowns[k]

top.run()

print 'after third solve'
print top['x']
print top['y']
for k,v in top.root.unknowns.iteritems():
    print k
    print top.root.unknowns[k]
于 2016-05-22T17:57:49.167 に答える