使用多维数组作为输入时，OpenMDAO ListGenerator DOE 失败答案

【问题标题】：OpenMDAO ListGenerator DOEs fail when using multidimensional arrays as inputs使用多维数组作为输入时，OpenMDAO ListGenerator DOE 失败
【发布时间】：2021-03-25 17:51:53
【问题描述】：

我在使用 OpenMDAO 时遇到问题，当输入是多维数组时，使用 ListGenerator DOE 驱动程序创建的 DOE 将失败。看起来数组在某个时候被展平，然后在分配输入数组时导致尺寸不匹配错误。我之前成功地使用了带有标量输入的 ListGenerator。

这里有一小段代码应该可以演示这个问题：

import openmdao.api as om
import numpy as np

# Create a subsystem with multidimensional array inputs
test_system = om.ExecComp('z = matmul(x,y)',
                          x=np.ones((3,3)),y=np.ones((3,3)),z=np.ones((3,3)))

# Single execution test
prob = om.Problem()
model = prob.model
model.add_subsystem('testSystem',test_system,promotes=['*'])
prob.setup()
prob['x'] = np.array([[1,2,3],
                      [4,5,6],
                      [7,8,9]])
prob['y'] = np.array([[9,8,7],
                      [6,5,4],
                      [3,2,1]])
prob.run_model()
print(prob['z'])

# DOE test
prob2 = om.Problem()
model2 = prob2.model
model2.add_subsystem('testSystem',test_system,promotes=['*'])
model2.add_design_var('x')
model2.add_design_var('y')
model2.add_objective('z')

prob2.setup()

case_list = [
    [['x',np.array([[1,2,3],[4,5,6],[7,8,9]])],['y',np.array([[9,8,7],[6,5,4],[3,2,1]])]]
    ]

prob2.driver = om.DOEDriver(case_list)
prob2.driver.add_recorder(om.SqliteRecorder("cases.sql"))

prob2.run_driver()
prob2.cleanup()

第一个问题（“prob”）按预期工作。第二个问题（“prob2”）应该通过使用单一案例 DOE 来完成同样的事情，但失败并出现以下错误：

Traceback (most recent call last):

  File "/Users/josephrobinson/Documents/GitHub/thesis/work/python/testing/doe_array_test.py", line 46, in <module>
    prob2.run_driver()

  File "/Users/josephrobinson/opt/anaconda3/envs/thesis/lib/python3.8/site-packages/openmdao/core/problem.py", line 665, in run_driver
    return self.driver.run()

  File "/Users/josephrobinson/opt/anaconda3/envs/thesis/lib/python3.8/site-packages/openmdao/drivers/doe_driver.py", line 171, in run
    self._run_case(case)

  File "/Users/josephrobinson/opt/anaconda3/envs/thesis/lib/python3.8/site-packages/openmdao/drivers/doe_driver.py", line 195, in _run_case
    raise(ValueError(msg))

ValueError: Error assigning x = [[1 2 3]
 [4 5 6]
 [7 8 9]]: could not broadcast input array from shape (3,3) into shape (9)

我将 OpenMDAO 3.8.0 与 Python 3.8.5 一起使用，尽管在更新之前我在 OpenMDAO 3.4.0 上看到了同样的错误。对于更大的三维数组，我也看到了同样的问题。

有人知道我该如何解决这个问题吗？

提前致谢。

【问题讨论】：

标签： python arrays numpy openmdao

【解决方案1】：

这似乎是 OpenMDAO 3.8.1 及更早版本中的一个错误。看起来我们需要先将提供给 ListGenerator 的设计变量值展平，然后再将它们设置到模型中。我已经在github上提交了一个问题：

https://github.com/OpenMDAO/OpenMDAO/issues/1983

您可以通过该链接监控进度。

作为目前的解决方法，如果您展平列表中的设计变量，您提供的示例将毫无例外地运行：

  case_list = [
      [['x',np.array([[1,2,3],[4,5,6],[7,8,9]]).flatten()],
       ['y',np.array([[9,8,7],[6,5,4],[3,2,1]]).flatten()]]
      ]

【讨论】：