ValueError: фигуры (6,108) и (36,) не выровнены: 108 (тусклый 1)!= 36 (тусклый 0)

Question

Я вижу некоторые ошибки измерения от метода линеаризации. В частности,

File "/usr/local/lib/python2.7/site-packages/openmdao/core/system.py", line 726, in _apply_linear_jac
    arg_vec[param] += J.T.dot(dresids._flat(unknown)).reshape(shape)
ValueError: shapes (6,108) and (36,) not aligned: 108 (dim 1) != 36 (dim 0)

У меня есть несколько 2d массивов, которые определяют часть якобиана для одного выходного вектора со всеми входными векторами. Я не смог точно определить точную причину ошибки, так как есть несколько случаев, которые соответствуют заданным формам. Мне интересно, если проблема может быть решена путем выравнивания массива самостоятельно.

Любые другие мысли о том, как устранить эту ошибку, будут оценены. Эти градиенты отлично работали в до-1.0.

Полное сообщение об ошибке ниже:

Traceback (most recent call last):
  File "/Users/user/directory/exampleOptimization.py", line 82, in <module>
    prob.run()
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/problem.py", line 789, in run
    self.driver.run(self)
  File "/usr/local/lib/python2.7/site-packages/openmdao/drivers/scipy_optimizer.py", line 191, in run
    options=self.opt_settings)
  File "/usr/local/lib/python2.7/site-packages/scipy/optimize/_minimize.py", line 452, in minimize
    constraints, callback=callback, **options)
  File "/usr/local/lib/python2.7/site-packages/scipy/optimize/slsqp.py", line 383, in _minimize_slsqp
    g = append(fprime(x),0.0)
  File "/usr/local/lib/python2.7/site-packages/scipy/optimize/optimize.py", line 285, in function_wrapper
    return function(*(wrapper_args + args))
  File "/usr/local/lib/python2.7/site-packages/openmdao/drivers/scipy_optimizer.py", line 306, in _gradfunc
    return_format='array')
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/driver.py", line 666, in calc_gradient
    sparsity=sparsity)
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/problem.py", line 899, in calc_gradient
    sparsity=sparsity)
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/problem.py", line 1261, in _calc_gradient_ln_solver
    dx_mat = root.ln_solver.solve(rhs, root, mode)
  File "/usr/local/lib/python2.7/site-packages/openmdao/solvers/scipy_gmres.py", line 104, in solve
    callback=self.monitor)
  File "<string>", line 2, in gmres
  File "/usr/local/lib/python2.7/site-packages/scipy/sparse/linalg/isolve/iterative.py", line 85, in non_reentrant
    return func(*a, **kw)
  File "/usr/local/lib/python2.7/site-packages/scipy/sparse/linalg/isolve/iterative.py", line 478, in gmres
    work[slice2] += sclr1*matvec(work[slice1])
  File "/usr/local/lib/python2.7/site-packages/scipy/sparse/linalg/interface.py", line 220, in matvec
    y = self._matvec(x)
  File "/usr/local/lib/python2.7/site-packages/scipy/sparse/linalg/interface.py", line 460, in _matvec
    return self.__matvec_impl(x)
  File "/usr/local/lib/python2.7/site-packages/openmdao/solvers/scipy_gmres.py", line 162, in mult
    system._sys_apply_linear(mode, ls_inputs=self.system._ls_inputs, vois=(voi,))
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/group.py", line 680, in _sys_apply_linear
    gs_outputs=gs_outputs)
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/group.py", line 680, in _sys_apply_linear
    gs_outputs=gs_outputs)
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/system.py", line 619, in _sys_apply_linear
    self.apply_linear(self.params, self.unknowns, dparams, dunknowns, dresids, mode)
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/component.py", line 512, in apply_linear
    mode)
  File "/usr/local/lib/python2.7/site-packages/openmdao/core/system.py", line 726, in _apply_linear_jac
    arg_vec[param] += J.T.dot(dresids._flat(unknown)).reshape(shape)
ValueError: shapes (6,108) and (36,) not aligned: 108 (dim 1) != 36 (dim 0)

Answer 1

Просто добавили исправление в dev, которое улучшает сообщения об ошибках для якобианов неправильного размера. Теперь они будут выглядеть так:

"In component 'comp', the derivative of 'y2' wrt 'x1' should have shape '(3, 2)' but has shape '(3, 3)' instead."

Если вы выберете последнюю версию и попробуете запустить ее снова, я думаю, что новое сообщение поможет точно определить проблему.