Векторные операции над собственным массивом векторов ⇐ C++

[Anonymous]

У меня есть двумерный собственный массив, где каждый элемент массива представляет собой трехэлементный собственный вектор (например, поле скорости над поверхностью).

Я хочу умножить каждый элемент двумерного массива на трехэлементный собственный вектор строк, эффективно взяв скалярное произведение.

Код: Выделить всё

    Eigen::Array velField(5, 5);
Eigen::Vector3d n;
// ... initialisation of n and velField not shown

Eigen::Array result(5, 5);
result = n.transpose() * velField;

Это выдает ошибку компиляции. ВЫ СМЕШАЛИ РАЗНЫЕ ЧИСЛОВЫЕ ТИПЫ. Однако если я не присвою результат, а просто вычислю его:

Код: Выделить всё

    n.transpose() * velField;

он ​​компилируется. Каков правильный тип возвращаемого значения для «результата» или как еще это исправить?

РЕДАКТИРОВАТЬ

Ту же проблему можно наблюдать и в этом более простом случае (умножение на двойное):

Код: Выделить всё

Eigen::Array velField(5, 5);
// ... initialisation of velField not shown
Eigen::Array result(5, 5);
result = 3.0 * velField;

Однако тот же код работает при использовании std::complex вместо Eigen::Vector3d:

Код: Выделить всё

Eigen::Array velField(5, 5);
// ... initialisation of velField not shown
Eigen::Array result(5, 5);
result = 3.0 * velField;

Интересно, почему это работает для std::complex, но не с Eigen::Vector3d. Оба типа определяют оператор * с двойным .

Подробнее здесь: https://stackoverflow.com/questions/225 ... of-vectors