MPI 梯形积分法（Fortran）

最新推荐文章于 2024-08-07 14:39:08 发布

「已注销」

最新推荐文章于 2024-08-07 14:39:08 发布

阅读量4.1k

点赞数 1

分类专栏： Fortran MPI 文章标签： fortran mpi 并行计算

本文链接：https://blog.csdn.net/vaf714/article/details/88414913

版权

本文介绍了如何使用 MPI 并行编程接口在 Fortran 中实现梯形积分法，通过并行计算提升数值积分的效率。内容涵盖了 MPI 的基本概念，以及将并行思想应用到梯形法则的具体步骤，帮助读者理解并行计算在科学计算中的应用。

摘要由CSDN通过智能技术生成

MPI 梯形积分法（Fortran）

使用MPI分布式接口，Fortran语言实现，梯形积分法实现计算函数积分，下面的代码中示例函数是 y=x^2，可以在 f(x)中修改成其他函数。

program main
    use mpi
    implicit none

    ! MPI 相关变量
    integer :: rank, size, ierr, stat(MPI_STATUS_SIZE), i
    ! 运算相关变量
    real :: a, b, h, local_a, local_b, local_n, local_total, total
    integer :: n
    ! 计算运算时间相关变量
    real(8) :: start_time, finish_time
    real :: local_elapsed, elapsed
    ! 函数声明
    real, external :: trap
    
    ! 初始化
    call MPI_INIT(ierr)
    call MPI_COMM_RANK(MPI_COMM_WORLD, rank, ierr)
    call MPI_COMM_SIZE(MPI_COMM_WORLD, size, ierr)

    ! 获取输入
    call get_input(rank, size, stat, ierr, a, b, n)
    
    ! 同步各进程，各进程开始计时
    call MPI_BARRIER(MPI_COMM_WORLD, ierr)
    start_time = MPI_WTIME()

    ! 各进程计算
    h = (b - a) / n
    local_n = n / size
    local_a = a + rank*local_n*h
    local_b = local_a + local_n*h
    local_total = trap(local_a, local_b, local_n, h)
    ! 求和并发送到0号进程
    call MPI_REDUCE(local_total, total, 1, MPI_REAL, MPI_SUM,