C++多线程实现矩阵对象点积的代码技术咨询

嘿，我看到你正在尝试用C++多线程实现矩阵点积，先帮你梳理下当前代码里的问题，再给出完整的可运行实现方案：

首先修正核心计算逻辑的错误 #

你的dot_t函数里把点积的元素相乘写成了元素相加，这是点积的核心错误，先修正这个：

double mat_cont::dot_t(mat_cont& other, const int offset, const int upper_lim) {
    double local_sum = 0;
    // 点积是对应元素相乘后累加，修正运算符
    for (int i = offset; i < offset + upper_lim; ++i) {
        local_sum += this->mat[i] * other.mat[i];
    }
    return local_sum;
}

完善多线程调度的dot函数 #

你之前的dot函数只声明了future数组但没启动线程，这里补充完整的多线程逻辑，同时加入维度校验、任务分片和剩余元素处理：

#include <future>
#include <stdexcept>

double mat_cont::dot(mat_cont& other) {
    // 先校验两个矩阵维度是否一致，避免越界访问
    if (this->dim != other.dim) {
        throw std::invalid_argument("Matrices must have identical dimensions for dot product calculation");
    }

    const int total_elements = dim * dim;
    const int num_threads = 4; // 对应你之前的4分片逻辑
    std::future<double> threads[num_threads - 1]; // 3个线程处理前3份
    int part_size = total_elements / num_threads;
    int remaining_elements = total_elements % num_threads;

    // 启动前3个工作线程，分配各自的元素区间
    for (int i = 0; i < num_threads - 1; ++i) {
        int start_offset = i * part_size;
        // 用std::launch::async强制创建新线程，传递this和引用参数
        threads[i] = std::async(std::launch::async, &mat_cont::dot_t, this, std::ref(other), start_offset, part_size);
    }

    // 主线程处理最后一份，包含剩余的元素（如果有的话）
    double total_sum = dot_t(other, (num_threads - 1) * part_size, part_size + remaining_elements);

    // 等待所有线程完成，累加各自的计算结果
    for (auto& thread_future : threads) {
        total_sum += thread_future.get();
    }

    return total_sum;
}

关键细节说明 #

• 维度校验：必须确保两个矩阵维度相同，否则点积无意义，还会导致数组越界，这里用异常抛出明确的错误信息。
• 任务分片：将总元素数分成4份，前3份交给独立线程处理，最后一份由主线程处理，同时处理不能整除时的剩余元素。
• std::async的使用：std::launch::async参数确保系统创建新线程执行任务（默认策略可能会用同步执行），传递other时要用std::ref，否则会传递副本，浪费内存且不符合逻辑。
• 结果收集：通过future.get()等待线程完成并获取局部和，这里如果线程抛出异常，get()会重新抛出，你可以根据需要添加异常处理逻辑。

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