C ++代码中带有MPI的每个进程具有相同的随机数

如何解决C ++代码中带有MPI的每个进程具有相同的随机数

|| 我有可以使用的C ++ MPI代码,因为它可以编译并确实在指定数量的处理器(n)上启动。问题在于它只是简单地执行了n次相同的计算,而不是进行了n次较快的计算。 我已经破解了在许多站点上找到的许多示例,并且似乎缺少对MPI_Send和MPI_Receive的正确使用,但是我找不到这些将功能作为输入的命令的实例(这些MPI命令为何对功能以外的其他功能有用)。 我的代码如下。本质上,它调用了我编写的C ++函数以获取Fisher的Exact Test p值。我只是为了测试速度而使用了随机数位。 该程序要做的是将Fisher.TwoTailed与每个随机变量集(即A,B,C和D)拖到不同的处理器上,而不是在多个处理器上进行完全相同的计算。预先感谢您的任何见解-欢呼! 这是代码:
int 
main (int argc,char* argv[]) 
{ 

  int id;
  int p;
//
//  Initialize MPI.
//
  MPI::Init ( argc,argv );
//
//  Get the number of processors.
//
  p = MPI::COMM_WORLD.Get_size ( );
//
//  Get the rank of this processor.
//
 id = MPI::COMM_WORLD.Get_rank ( );

  FishersExactTest Fishers;

  int i = 0;
  while (i < 10) {
    int A = 0 + rand() % (100 - 0);
    int B = 0 + rand() % (100 - 0);
    int C = 0 + rand() % (100 - 0);
    int D = 0 + rand() % (100 - 0); 
    cout << Fishers.TwoTailed(A,B,C,D) << endl;
    i += 1;
  }

  MPI::Finalize ( );
  return 0;
}
    

解决方法

        您应该研究一些有关并行计算和MPI的基础培训。国家基础超级计算应用中心(NCSA)提供了一套免费的在线课程,这些资源可以教我一些基础知识。     ,        您必须告诉MPI如何并行处理代码-它不会自动执行。 换句话说,您不能在所有系统上初始化MPI,然后将它们传递给同一循环。您想使用每个处理器的ѭ1来确定它将在循环的哪一部分工作。然后,您需要他们将其结果全部传递回ID 0。     ,        以上所有答案都是完全正确的。让我补充一点: 在这里,由于您似乎只是在进行随机抽样,因此您需要做的就是让不同的处理器生成不同的随机数以提供给Fishers。TwoTailed是确保它们对PRNG都有不同的种子:
int 
main (int argc,char* argv[]) 
{ 

  int id;
  int p;
//
//  Initialize MPI.
//
  MPI::Init ( argc,argv );
//
//  Get the number of processors.
//
  p = MPI::COMM_WORLD.Get_size ( );
//
//  Get the rank of this processor.
//
 id = MPI::COMM_WORLD.Get_rank ( );

  FishersExactTest Fishers;

  srand(id);  // <--- each rank gets a different seed

  int i = 0;
  while (i < 10) {
    int A = 0 + rand() % (100 - 0);
    int B = 0 + rand() % (100 - 0);
    int C = 0 + rand() % (100 - 0);
    int D = 0 + rand() % (100 - 0); 
    cout << Fishers.TwoTailed(A,B,C,D) << endl;
    i += 1;
  }

  MPI::Finalize ( );
  return 0;
}
由于循环是从1..10开始的,因此您仍将使每个进程进行10个采样。如果希望他们总共做10次,可以将10除以p并做一些分配余数的操作:例如
int niters = (10+id)/p;
int i=0;
while (i < niters) {
  ...
}
    ,        那么,当您运行MPI作业时会收到什么消息?只是要重申其他人说的话,您将必须明确定义每个处理器的工作...例如。.如果您的等级为0(默认),则执行此操作...如果您的等级为1,依此类推。 ..(或其他语法)定义每个等级的角色。然后,根据您的代码结构,可以使节点具有Send / Recv,Gather,Scatter等。     

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 &lt;select id=&quot;xxx&quot;&gt; SELECT di.id, di.name, di.work_type, di.updated... &lt;where&gt; &lt;if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 &lt;property name=&quot;dynamic.classpath&quot; value=&quot;tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams[&#39;font.sans-serif&#39;] = [&#39;SimHei&#39;] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -&gt; systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping(&quot;/hires&quot;) public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate&lt;String
使用vite构建项目报错 C:\Users\ychen\work&gt;npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-