Question

是否可以计算具有复数双精度的双打uBLAS向量的元素积？以下代码无法编译，因为它找不到重载的运算符*。我希望它能够起作用，因为将复数双倍乘以一个很好的定义。

#include <complex>

#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/io.hpp>

int main(int argc, char **argv)
{
    using namespace boost::numeric::ublas;

    vector<double> v(3);
    for (unsigned i = 0; i < v.size(); ++i)
    {
        v (i) = i;
    }

    vector<std::complex<double> > w = v * std::complex<double>(3.0, -1.0);

    return 0;
}

使用GCC 4.6和Boost 1.55.0进行编译会产生以下结果：

error: no match for ‘operator*’ (operand types are ‘boost::numeric::ublas::vector<double>’ and ‘std::complex<double>’)

Answer 1

通过查看vector_expression.hpp中的重载*运算符：

// (t * v) [i] = t * v [i]
template<class T1, class E2>
BOOST_UBLAS_INLINE
typename enable_if< is_convertible<T1, typename E2::value_type >,    
typename vector_binary_scalar1_traits<const T1, E2, scalar_multiplies<T1, typename E2::value_type> >::result_type
>::type
operator * (const T1 &e1,
            const vector_expression<E2> &e2) {
    typedef typename vector_binary_scalar1_traits<const T1, E2, scalar_multiplies<T1, typename E2::value_type> >::expression_type expression_type;
    return expression_type (e1, e2 ());
}

看起来问题是＆＃34; is_convertible＆＃34;的输出，它没有双向，所以在这种情况下没有从std :: complex到double的转换，它没有＆＃39;工作。我添加了一个新的定义，只交换了这个模板中参数的顺序，它可以工作......

抱歉英语不好

如何将双倍的Boost uBLAS向量与复数双因子相乘？

1 个答案: