Vector类hadamard_product函数触发segmentation fault的原因及修复

错误根源 #

触发段错误的直接原因是未初始化的野指针：

• 在main函数中，你声明了Vector ** returns;但没有给它赋值，这个指针指向随机的内存地址（野指针）。
• 当hadamard_product函数执行*returns = new Vector[this->length];时，尝试对野指针解引用，往非法内存地址写入数据，直接触发segmentation fault。

修复方案 #

推荐两种修复方式，优先选择更安全的C++现代风格实现：

方案1：改用std::vector返回结果（最安全） #

完全抛弃手动内存管理，让标准容器自动处理内存分配与释放，从根源上避免指针错误：

1. 先添加必要的头文件，再修改hadamard_product函数定义：

#include <vector>  // 添加这个头文件到代码顶部

// ... 其他类成员代码保持不变 ...

std::vector<Vector> hadamard_product(const Vector& rhs) const
{
    if (!(this->is_column && !rhs.is_column_vector()))
    {
        throw std::runtime_error("Error: rhs must be a row-vector");
    }

    if (this->length != rhs.length)
    {
        throw std::runtime_error("Error: Vectors must have the same length to perform Hadamard product");
    }

    std::vector<Vector> result;
    result.reserve(this->length);  // 预分配内存提升性能
    for (int i = 0; i < this->length; i++)
    {
        Vector temp = rhs.Clone();
        float val = this->elements[i];
        temp.scalar_product(val);
        result.push_back(temp);
    }
    return result;
}

2. 修改main中的调用逻辑：

try 
{
    std::vector<Vector> returns = u.hadamard_product(v);
    // 遍历输出结果
    for (auto& vec : returns)
    {
        vec.display();
    }
    // 无需手动释放内存，std::vector会自动处理
} 
catch (std::exception& e) 
{
    std::cerr << "Error: " << e.what() << std::endl;
}

方案2：修正指针传递逻辑（兼容原有设计） #

如果一定要保留指针输出的方式，需要确保传入的指针是有效的，或者改用指针的引用传递：

1. 修改hadamard_product的函数签名：

void hadamard_product(const Vector& rhs, Vector**& returns)
{
    // 函数内部代码保持不变
}

2. 在main中初始化指针后再传入：

try 
{
    Vector** returns = nullptr;
    u.hadamard_product(v, returns);
    // 遍历输出结果
    for (int i=0 ; i<u.get_length() ; i++)
    {
        returns[i]->display();
    }
    delete[] returns; // 正确释放整个Vector数组，不要单独delete每个元素
} 
catch (std::exception& e) 
{
    std::cerr << "Error: " << e.what() << std::endl;
}

注意：原注释中的delete returns[i];是错误的，因为new Vector[this->length]分配的是连续数组，必须用delete[] returns一次性释放，单独delete会破坏内存布局导致未定义行为。

额外优化建议 #

• 移除dot_product函数中throw语句后的冗余return 0.0;，因为throw会直接终止函数执行，永远不会走到return。
• 给Vector添加带长度参数的构造函数，简化对象初始化：

Vector(int length, bool is_column = true) : length(length), is_column(is_column)
{
    elements = new float[length];
}

• 考虑用std::vector<float>代替float* elements，彻底避免手动内存管理的风险。

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