使用Eigen库计算向量lpNorm时循环调用报错的解决方法

尝试用Eigen库计算向量的lpNorm，显式指定模板参数（比如v.lpNorm<1>()）时能正常运行，但在for循环里用变量p作为模板参数调用v.lpNorm<p>()就报错，求修复方法。

代码示例：

#include <iostream>
#include <Eigen/Dense>

using Eigen::VectorXd;

int main()
{

    int sizev = 3;
    int p;

    Eigen::VectorXd v(sizev);

    v(0) = 3.;
    v(1) = 2.;
    v(2) = 1.;

    // test 1, passes

    std::cout << "||v||_" << 1 << " = " << v.lpNorm<1>() << std::endl;
    std::cout << "||v||_" << 2 << " = " << v.lpNorm<2>() << std::endl;
    std::cout << "||v||_" << 3 << " = " << v.lpNorm<3>() << std::endl;
    std::cout << "||v||_" << 4 << " = " << v.lpNorm<4>() << std::endl;
    std::cout << "||v||_inf = "         << v.lpNorm<Eigen::Infinity>() << std::endl;

    // test 2, fails

    for (int p=1; p<5; p++)
    {
        std::cout << "||v||_" << p << " = " << v.lpNorm<p>() << std::endl;
    }
    
    return 0;

}

编译报错信息：

error: no matching member function for call to 'lpNorm'
        std::cout << "||v||_" << p << " = " << v.lpNorm<p>() << std::endl;
                                               ~~^~~~~~~~~

note: candidate template ignored: invalid explicitly-specified argument for template parameter 'p'
    template<int p> EIGEN_DEVICE_FUNC RealScalar lpNorm() const;
                                                 ^
1 error generated.

问题原因 #

C++模板参数必须是编译期常量，而for循环里的变量p是运行时才确定值的变量，无法作为模板参数传递给lpNorm。

修复方法 #

方法1：使用运行时参数版本的lpNorm函数 #

Eigen 3.3及以上版本提供了接受运行时整数参数的lpNorm重载，直接将p作为函数参数传入即可：

for (int p=1; p<5; p++)
{
    std::cout << "||v||_" << p << " = " << v.lpNorm(p) << std::endl;
}

这个版本会在运行时处理不同的p值，适配循环场景。

方法2：编译期枚举固定范围的p值（追求性能时用） #

如果p的取值范围很小且固定，可通过编译期分支枚举所有可能的p，保留模板版本的编译期优化优势（需C++17及以上）：

template<int P>
void printNorm(const Eigen::VectorXd& v) {
    std::cout << "||v||_" << P << " = " << v.lpNorm<P>() << std::endl;
}

int main() {
    // 向量初始化代码...

    for (int p=1; p<5; p++) {
        switch(p) {
            case 1: printNorm<1>(v); break;
            case 2: printNorm<2>(v); break;
            case 3: printNorm<3>(v); break;
            case 4: printNorm<4>(v); break;
        }
    }
}

内容的提问来源于stack exchange，提问作者gypsophila