LU 分解，采用行连续划分方式下的 MPI 实现

koastal

于 2016-01-20 15:30:57 发布

阅读量2.5k

点赞数

分类专栏：并行程序设计

本文链接：https://blog.youkuaiyun.com/koastal/article/details/50548800

版权

这篇博客详细介绍了如何使用 MPI 并行计算实现 LU 分解，通过行连续划分的方式，将一个 MxM 矩阵分解为 L 和 U 两个矩阵。程序首先读取矩阵数据，然后通过 MPI 进程间通信进行矩阵的分布式存储和运算，最后输出分解后的 L 和 U 矩阵。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include "stdio.h"
#include "stdlib.h"
#include "mpi.h"
#define a(x,y) a[x*M+y]
/*A为M*M矩阵*/
#define A(x,y) A[x*M+y]
#define l(x,y) l[x*M+y]
#define u(x,y) u[x*M+y]
#define floatsize sizeof(float)
#define intsize sizeof(int)

int M,N;
int m;
float *A;
int my_rank;
int p;
MPI_Status status;

void fatal(char *message)
{
    printf("%s\n",message);
    exit(1);
}

void Environment_Finalize(float *a,float *f)
{
    free(a);
    free(f);
}

int main(int argc, char **argv)
{
    int i,j,k,my_rank,group_size;
    int i1,i2;
    int v,w;
    float *a,*f,*l,*u;
    FILE *fdA;

    MPI_Init(&argc,&argv);
    MPI_Comm_size(MPI_COMM_WORLD,&group_size);
    MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);

    p=group_size;

    if (my_rank==0)
    {
        fdA=fopen("dataIn.txt","r");
        fscanf(fdA,"%d %d", &M, &N);
        if(M != N)
        {