C语言中如何实现double类型指针数组的拷贝？

现有代码问题排查 #

strtod赋值逻辑问题 #

• 你所担心的「将strtod输出作为地址赋值」的问题不存在：strtod 本身返回double类型的数值，赋值给positions_particles_source[i][j]（一个double类型的变量）的语法逻辑完全合法，不存在将数值当作地址使用的问题。
• 该代码的实际错误：内层j循环始终读取同一个字符串initial_structure_file_entries[i]，导致同一行的所有粒子坐标值完全相同。需要每次调用strtod后推进字符串指针到下一个数值的起始位置，才能正确读取所有粒子的坐标。

拷贝失败原因 #

• 你当前代码中逐行memcpy的写法是深拷贝的正确实现：每个维度单独为目标数组分配内存，再拷贝对应行的所有double数值，拷贝完成后两个数组的指针指向独立内存，修改其中一个不会影响另一个。
• 如果你之前尝试直接对数组整体赋值，或者直接memcpy(positions_particles_destination, positions_particles_source, 3*sizeof(double*))，这两种属于浅拷贝，仅复制了3个指针的地址，两个数组的指针会指向同一块内存，自然无法得到独立的拷贝结果。
• 额外提示：C语言中calloc返回的void*可以隐式转换为任意指针类型，不需要强制类型转换，多余的强制转换反而可能隐藏潜在的编译错误。

更优的矩阵存储与操作方案 #

你当前使用的「指针数组+逐行分配内存」是C语言动态矩阵的常见实现方式，但对于固定列数为3的粒子坐标场景，更推荐使用连续内存分配方案，优势如下：

• 内存连续，CPU缓存命中率更高，访问速度远快于非连续的逐行分配内存
• 整个矩阵的拷贝、释放都只需要一次操作，代码更简洁，不易出现内存泄漏

连续内存方案的声明方式如下，访问语法仍兼容[粒子索引][坐标维度]的形式：

// 声明维度为 number_particles_total × 3 的连续内存矩阵
double (*positions_particles_source)[3] = calloc(number_particles_total, sizeof(*positions_particles_source));

修正后的代码参考 #

原有指针数组写法（修正读取逻辑） #

double *positions_particles_source[3];
double *positions_particles_destination[3];
char *initial_structure_file_entries[3]; // 已预先存储三行坐标文本

for (int i = 0; i < 3; i++) {
    positions_particles_source[i] = calloc(number_particles_total, sizeof(double));
    char *cur_ptr = initial_structure_file_entries[i];
    char *end_ptr;
    for (int j = 0; j < number_particles_total; j++) {
        positions_particles_source[i][j] = strtod(cur_ptr, &end_ptr);
        cur_ptr = end_ptr; // 推进指针到下一个数值起始位置
    }
    // 深拷贝到目标数组
    positions_particles_destination[i] = calloc(number_particles_total, sizeof(double));
    memcpy(positions_particles_destination[i], positions_particles_source[i], number_particles_total * sizeof(double));
}

连续内存矩阵写法 #

// 分配连续内存的矩阵，索引形式为 [粒子序号][x/y/z维度]
double (*positions_source)[3] = calloc(number_particles_total, sizeof(*positions_source));
double (*positions_dest)[3] = calloc(number_particles_total, sizeof(*positions_dest));

// 读取坐标数据
for (int dim = 0; dim < 3; dim++) {
    char *cur_ptr = initial_structure_file_entries[dim];
    char *end_ptr;
    for (int p_idx = 0; p_idx < number_particles_total; p_idx++) {
        positions_source[p_idx][dim] = strtod(cur_ptr, &end_ptr);
        cur_ptr = end_ptr;
    }
}

// 一次操作完成整个矩阵深拷贝
memcpy(positions_dest, positions_source, number_particles_total * sizeof(*positions_source));

// 使用完释放内存，仅需一次free操作
free(positions_source);
free(positions_dest);

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