MPI_Type_create_subarray与MPI_File_read_all读文件乱码及段错误

• 核心逻辑混淆了MPI派生类型的两个角色
  MPI_Type_create_subarray生成的mysub类型，本来就是为「大尺寸全局数组提取子块」设计的：它的内存布局规则是，每写入连续5个整数后，自动跳过5个整数的间隔（因为全局数组每行是10个元素，只取前5个），整个类型的总内存跨度是44个int的长度（从全局数组(0,0)到(4,4)的位置）。
  但你实际分配的本地缓冲区是只有25个int的紧凑连续空间，没有预留行间隔的冗余位置，却直接把mysub同时当做了文件视图类型和内存读写类型使用——MPI会严格按照mysub的步长规则往这个小缓冲区里写数据，直接造成越界访问：越界写出来的就是随机乱码值，越界到无权限访问的内存时就触发段错误，和你观测到的输出、崩溃现象完全对应。
• 文件视图参数的匹配偏差
  MPI_File_set_view的参数逻辑是：displacement（文件偏移）、etype（文件操作的基本单位）、filetype（当前进程可访问的文件区域映射）。你虽然把mysub放到了filetype的位置，但后续读操作没有匹配对应的内存布局类型，进一步放大了读取错位的问题。

1. 拆分派生类型的用途，别拿同一个类型同时管文件布局和内存布局
   • 保留原来的mysub作为文件侧的filetype，用来标记在1010的全局文件数组中，需要提取(0,0)起点的55子块
   • 新增一个连续类型作为内存侧的布局，匹配你本地紧凑连续的缓冲区：

     MPI_Datatype memsub;
     MPI_Type_contiguous(subsize * subsize, MPI_INT, &memsub);
     MPI_Type_commit(&memsub);
     

2. 调整读操作的参数，匹配布局规则
   原来的读调用是把mysub当内存类型，步长错位，改成直接用连续缓冲区的首地址作为读写位置：

   MPI_File_read_all(fh, buf, 1, memsub, &status);
   

   这一步的逻辑是：MPI会按照filetype的定义，从文件里取出所有目标离散位置的整数，按照memsub的连续布局，依次填充到你的本地紧凑缓冲区里，不会再出现间隔跳跃的越界写问题。
3. 资源释放要对应
   用完的memsub和mysub都要通过MPI_Type_free释放，避免内存泄漏。

乱码现象的对应解释：按照错位的步长写规则，MPI前两次写的位置刚好落在你分配的25个int缓冲区范围内（第0行0-4、第1行0-4的位置），所以这部分输出是符合预期的；再往下的写位置已经超出缓冲区边界，写到的是堆内存里的随机值，和你看到的异常输出完全匹配。

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