【无标题】

已于 2023-06-21 20:01:56 修改
阅读量49 收藏
点赞数
文章标签： c++ 算法 开发语言
于 2023-06-21 19:31:24 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_51469755/article/details/131333362
版权
这段代码实现了一个使用OpenMP进行并行化的Gauss-Seidel迭代法，用于求解稀疏矩阵的线性方程组。程序首先检查B矩阵是否有全零列，然后进行迭代计算，每一步都考虑了矩阵元素的稀疏性以优化性能。在每次迭代中，它会计算误差范数并检查是否满足预设的精度要求。如果解超出指定范围或超过最大迭代次数，函数将返回相应错误码。
摘要由CSDN通过智能技术生成 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <vector>
#include <iostream>
#include <cstring>
#include <fstream>
#include <sstream>
#include <time.h>
#include "omp.h"
using namespace std;
#define MAX_SIZE 2500 // 矩阵最大维数
#define bound pow(2,127) // 判断迭代发散边界
#define ZERO 0.000001 //当一个数字小于0.000001时认为其为0

// solve the problem with OpenMP
int Gauss_Seidel ( int n, float a[][MAX_SIZE], float b[], float x[], float TOL, int MAXN )
{
    //Turning "ax=b" into "x=Bx+f"
    float B[n][n]; //B矩阵为迭代矩阵
    int i,j;
    // 减少稀疏矩阵的计算开销
    int count[n]; //计算B矩阵的某一指定列中元素为0的个数
    for(i=0;i<n;i++)
    {
        B[i][i] = 0;
        for(j=0;j<n;j++)
        {
            if(j!=i)
                B[i][j] = -a[i][j]/a[i][i];
        }
    }
    //若B矩阵中有一列的元素值全为0,则Gauss_Seidel函数返回-1
    for(i=0;i<n;i++)
    {
        count[i] = 0;
        for(j=0;j<n;j++)
        {
            if(B[i][j] == 0)
                count[i]++;
        }
    }
    if(count[i]==n)
        return -1;
    //Iteration
    int k = 0; //定义并初始化迭代次数k=0
    float x_kp1[n]; //创建新的数组用于存储第k+1次后x的迭代值
    float sum;
    float max = 0; //用于存储最大范数
#pragma omp parallel for num_threads(thread_count)
    for(int k = 0; k < MAXN; ++k)
    {
        for(i=0;i<n;i++)
        {
            sum = 0;
            for(j=0;j<n;j++)
            {
                if(j<i)
                    sum += a[i][j]*x_kp1[j];
                if(j>i)
                    sum += a[i][j]*x[j];
            }
            x_kp1[i] = (b[i]-sum)/a[i][i];
        }
        //求第k+1次迭代所得的x与第k次迭代所得的x之间差的无穷范数
        max = fabs(x_kp1[0]-x[0]);
        for(i=1;i<n;i++)
        {
            if(max<fabs(x_kp1[i]-x[i]))
                max = fabs(x_kp1[i]-x[i]);
        }
        //若第k次迭代中x的任意某个分量值超出了最大区间[-bound,bound]限制
        //则Gauss_Seidel函数返回-2
        for(i=0;i<n;i++)
        {
            if(fabs(x_kp1[i])>bound)
               return -2;
            x[i] = x_kp1[i]; //重复利用x向量，用于存储第k次迭代的x值
        }

        for(i=0;i<n;i++)
            x[i] = x_kp1[i]; //重复利用x向量，用于存储第k次迭代的x值
        //若经过k次叠代达到了精度TOL的要求
        //则Gauss_Seidel函数返回k值（求得的解存储在x中）
        if(max<TOL)
            return k;
    }
    //若算法经过MAXN次叠代后还没有达到精度TOL要求
    //则Gauss_Seidel函数返回0
    return 0;
}

// the size of the matrix?
// 读入稀疏矩阵
int main(int argc, char ** argv){ // row_size = 2500, col_size = 2500
    int thread_count = strtol(argv[1],NULL,10);
    int i, j, k, MAXN,n;// MAXN 最大迭代次数，n 计算稀疏矩阵特性的参数
    double TOL = ZERO; // TOL 精度要求
    float A[MAX_SIZE][MAX_SIZE], b[MAX_SIZE],x[MAX_SIZE];
    ifstream inFileA("./A5.txt");
    ifstream inFileb("./b.txt");
    float temp1;
    for(i = 0;i < MAX_SIZE;i++){
        for(j = 0;j < MAX_SIZE;j++){
            inFileA >> temp1;
            A[i][j] = temp1;
        }
    }
    int temp2;
    for(j = 0;j < MAX_SIZE;j++){
        inFileb >> temp2;
        b[j] = temp2;
        x[j] = 1;
    }
    n = MAX_SIZE;
    MAXN = 10; //最多迭代MAXN次，如果达到精度要求，则返回x,如果没有达到精度要求，则返回0
    Gauss_Seidel(n,A,b,x,TOL,MAXN);
}
/*
    vector<float> A;
    vector<int> b;
    // int row_size = 2500;
    // int col_size = 2500;
    ifstream inFile_A("./A5.txt");
    ifstream inFile_b("./b.txt");
    // 默认为char类型
    float temp1;
    while(!inFile_A.eof()){
        inFile_A >> temp1;
        A.push_back(temp1);
    }
    int temp2;
    while(!inFile_b.eof()){
        inFile_b >> temp2;
        b.push_back(temp2);
    }
*/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <vector>
#include <iostream>
#include <cstring>
#include <fstream>
#include <sstream>
#include <time.h>
#include "omp.h"
using namespace std;
#define MAX_SIZE 2500 // 矩阵最大维数
#define bound 1024 // 判断迭代发散边界
#define ZERO 0.000001 //当一个数字小于0.000001时认为其为0

// solve the problem with OpenMP
int Gauss_Seidel ( int n, float a[][MAX_SIZE], float b[], float x[], float TOL, int MAXN,int thread_count)
{
    //Turning "ax=b" into "x=Bx+f"
    float B[n][n]; //B矩阵为迭代矩阵
    int i,j;
    // 减少稀疏矩阵的计算开销
    int count[n]; //计算B矩阵的某一指定列中元素为0的个数
    for(i=0;i<n;i++)
    {
        B[i][i] = 0;
        for(j=0;j<n;j++)
        {
            if(j!=i)
                B[i][j] = -a[i][j]/a[i][i];
        }
    }
    //若B矩阵中有一列的元素值全为0,则Gauss_Seidel函数返回-1
    for(i=0;i<n;i++)
    {
        count[i] = 0;
        for(j=0;j<n;j++)
        {
            if(B[i][j] == 0)
                count[i]++;
        }
    }
    if(count[i]==n)
        return -1;
    //Iteration
    int k = 0; //定义并初始化迭代次数k=0
    float x_kp1[n]; //创建新的数组用于存储第k+1次后x的迭代值
    float sum;
    float max = 0; //用于存储最大范数

    for(int k = 0; k < MAXN; ++k)
    {
        #pragma omp parallel for num_threads(thread_count)
        for(i=0;i<n;i++)
        {
            sum = 0;
            for(j=0;j<n;j++)
            {
                if(j<i)
                    sum += a[i][j]*x_kp1[j];
                if(j>i)
                    sum += a[i][j]*x[j];
            }
            x_kp1[i] = (b[i]-sum)/a[i][i];
        }
        //求第k+1次迭代所得的x与第k次迭代所得的x之间差的无穷范数
        max = fabs(x_kp1[0]-x[0]);
        for(i=1;i<n;i++)
        {
            if(max<fabs(x_kp1[i]-x[i]))
                max = fabs(x_kp1[i]-x[i]);
        }
        //若第k次迭代中x的任意某个分量值超出了最大区间[-bound,bound]限制
        //则Gauss_Seidel函数返回-2
        for(i=0;i<n;i++)
        {
            if(fabs(x_kp1[i])>bound)
               return -2;
            x[i] = x_kp1[i]; //重复利用x向量，用于存储第k次迭代的x值
        }

        for(i=0;i<n;i++)
            x[i] = x_kp1[i]; //重复利用x向量，用于存储第k次迭代的x值
        //若经过k次叠代达到了精度TOL的要求
        //则Gauss_Seidel函数返回k值（求得的解存储在x中）
        if(max<TOL)
            return k;
    }
    //若算法经过MAXN次叠代后还没有达到精度TOL要求
    //则Gauss_Seidel函数返回0
    return 0;
}

// the size of the matrix?
// 读入稀疏矩阵
int main(int argc, char ** argv){ // row_size = 2500, col_size = 2500
    int thread_count = strtol(argv[1],NULL,10);
    int i, j, k, MAXN,n;// MAXN 最大迭代次数，n 计算稀疏矩阵特性的参数
    double TOL = ZERO; // TOL 精度要求
    float A[MAX_SIZE][MAX_SIZE], b[MAX_SIZE],x[MAX_SIZE];
    ifstream inFileA("./A5.txt");
    ifstream inFileb("./b.txt");
    float temp1;
    for(i = 0;i < MAX_SIZE;i++){
        for(j = 0;j < MAX_SIZE;j++){
            inFileA >> temp1;
            A[i][j] = temp1;
        }
    }
    float temp2;
    for(j = 0;j < MAX_SIZE;j++){
        inFileb >> temp2;
        b[j] = temp2;
        x[j] = 1;
    }
    n = MAX_SIZE;
    MAXN = 10; //最多迭代MAXN次，如果达到精度要求，则返回x,如果没有达到精度要求，则返回0
    Gauss_Seidel(n,A,b,x,TOL,MAXN,thread_count);
}
/*
    vector<float> A;
    vector<int> b;
    // int row_size = 2500;
    // int col_size = 2500;
    ifstream inFile_A("./A5.txt");
    ifstream inFile_b("./b.txt");
    // 默认为char类型
    float temp1;
    while(!inFile_A.eof()){
        inFile_A >> temp1;
        A.push_back(temp1);
    }
    int temp2;
    while(!inFile_b.eof()){
        inFile_b >> temp2;
        b.push_back(temp2);
    }
*/
沈离尘
关注
C# winform中无标题栏窗口如何实现鼠标拖动？
汽车电子软件领域知识分享
08-23 1347
在C# winform中，可以通过重写窗体的鼠标事件来实现无标题栏窗体的拖动。
Qt-无标题窗口
学而思
05-21 7665
        今天学习了下无标题窗口，于是来记录下我的学习过程首先建立一个Qt Widgets Application继承QWidget类创建完后直接运行就得到了如下效果那么，既然是无标题窗口就要去掉窗口栏啊，于是加上setWindowFlags(Qt::FramelessWindowHint | Qt::WindowMinimizeButtonHint);去掉标题栏，但是这样整个窗口就不能移动...
谷歌浏览器打开无标题
mp2020xuexi的博客
01-16 1017
谷歌浏览器
C#中实现无标题栏窗体拖动的代码
热门推荐
百锦再的博客
01-18 1万+
事件来实现拖动功能。当鼠标按下时，我们标记鼠标正在拖动，并保存当前坐标。在鼠标移动时，如果鼠标正在拖动，我们更新窗体的位置，使其跟随鼠标移动。当鼠标释放时，我们取消拖动标记。请注意，在窗体上添加了对应事件处理程序的窗体控件属性应该设置为对应的事件处理程序方法。在窗体设计视图，选择属性窗格，在事件选项卡中为。在C#中实现无标题栏窗体拖动的代码可以通过处理鼠标事件来实现。这样，当你在运行窗体时，你就可以通过鼠标拖动窗体来改变其位置。事件分别选择对应的处理程序方法。在这个示例代码中，我们通过处理。
C# 移动无标题栏无边框窗体的3种方法
SabreWulf
03-13 2282
目录 第一种：手工移动 第二种：调用系统API 1.引入下面代码 前提需要引入命名空间using System.Runtime.InteropServices 2.增加鼠标按下事件发送消息，让系统误以为按下是标题栏 第三种，重写 WndProc 第一种：手工移动 直接通过修改窗体位置从而达到移动窗体的效果 //定义一个位置信息Point用于存储鼠标位置 private Point mPoint; /// <summary>
Qt实现改变无标题栏窗体大小及移动窗体
danshiming的博客
11-10 3601
Qt实现鼠标左键拖动无标题栏窗体右下角改变窗体大小及移动窗体
Chrome显示无标题错误代码STATUS_INVALID_IMAGE_HASH
precy_wm的博客
01-03 1763
Chrome显示无标题错误代码STATUS_INVALID_IMAGE_HASH
WIN32无标题栏窗口移动方法种种
ZDT_zdh的博客
11-02 2220
WIN32无标题栏窗口移动方法种种 转载来源：http://blog.csdn.net/fuyun_cloud/article/details/8008197     首先，看看在正常情况下系统是怎样来移动程序窗口的。当使用者在程序窗口标题栏区域（非工作区）内，按下鼠标左键时将会发生下列事情：  ◆ 系统向该窗口过程函数发送WM_NCLBUTTONDOWN消息。  ◆ WM_NC...
Chrome 无标题 哭脸 已崩溃
wkj52719的博客
05-24 4765
谷歌浏览器chrome突然打开后任何页面均显示为无标题，哭脸。鼠标选中sheet页面时，提示已崩溃。 无法打开设置、开发工具箱 经过检查，是chrome的Userdata数据出错，导致页面加载无法正常展示页面。 修复办法： 打开文件夹C:\Users\{privateUserId}\AppData\Local\Google\Chrome，删除User data文件夹，如果需要备份UserData,请备份User Data到其他地方。 ...
c#实现无标题栏窗口的拖动
09-05
本篇文章将详细讲解如何在C#中实现这种无标题栏窗口的拖动以及其他基本操作。 首先，为了实现无标题栏窗口的拖动，我们需要处理窗体的鼠标事件。窗体的`FormBorderStyle`属性应设置为`None`，这会去除窗体的边框和...
Qt无标题窗口
05-21
在本文中，我们将深入探讨如何使用Qt框架创建一个无标题窗口，这在界面设计和自定义标题栏美化中非常有用。我们将关注`framelesswidget.cpp`、`main.cpp`、`framelesswidget.h`等核心文件，理解它们在实现无标题窗口...
C#创建无标题栏窗体
03-16
摘要:C#源码,菜单窗体,无标题栏窗体　C#创建无标题栏窗体源码，实际上是动态显示或隐藏窗体的标题栏，当隐藏的时候就类似窗体没有了标题栏，当显示标题栏的时候，鼠标按住标题栏即可拖动窗体，以前记得网友需要此...
Delphi实现拖动按钮移动无标题栏窗口
05-11
摘要:Delphi源码,界面编程,窗体拖动,无标题栏 无标题栏的窗体的拖动功能实现，Delphi添加一个可拖动窗体的按钮，通过对此按钮的控制可移动窗体，实现按住标题栏移动窗口的功能，无标题栏也就不能显示最大化、最小化...
VC去掉'无标题'字样的方法.zip
11-21
在使用Visual C++进行编程时，我们经常会在创建新的窗口或者对话框时默认看到标题栏显示为“无标题”。这可能会对用户界面的美观性和专业性造成一定的影响。本资料包提供了一种方法来帮助开发者去掉这个“无标题”...
C语言 | 第十章 | 函数 作用域
ZJC744575的博客
10-05 1212
为完成某一功能的程序指令(语句)的集合,称为函数。在C语言中,函数分为:自定义函数、系统函数(查看C语言函数手册)函数还有其它叫法，比如方法等，在本视频课程中，我们统一称为 函数。返回类型 函数名（形参列表）{执行语句...;// 函数体return 返回值;// 可选形参列表：表示函数的输入函数中的语句：表示为了实现某一功能代码块函数可以有返回值,也可以没有, 如果没有返回值，返回类型 声明为 void。
使用socket编程来实现一个简单的C/S模型（TCP协议）
caiji0169的博客
10-02 1723
下图是基于TCP协议的客户端/服务器程序的一般流程：服务器调用socket()、bind()、listen()完成初始化后，调用accept()阻塞等待，处于监听端口的状态，客户端调用socket()初始化后，调用connect()发出SYN段并阻塞等待服务器应答，服务器应答一个SYN-ACK段，客户端收到后从connect()返回，同时应答一个ACK段，服务器收到后从accept()返回。数据传输的过程：建立连接后，TCP协议提供全双工的通信服务，但是一般的客户端/服务器程序的流程是由客户端主动发起请求，
大话C++：第15篇 友元
最新发布
whccf的博客
10-09 756
友元（Friend）是C++中的一个特殊概念，它允许一个类或函数访问另一个类的私有和保护成员。友元可以是一个函数，也可以是一个类。
windows C++-创建数据流代理(三)
苏洛的博客
10-08 289
控制流指的是程序中操作的执行顺序。 控制流是通过使用条件语句、循环等控制结构调节的。 此外，数据流指的是一种编程模型，在这种模型中，只有当所有需要的数据都可用时，才会进行计算。 数据流编程模型与消息传递这一概念相关，其中程序的独立组件通过发送消息相互通信。 异步代理同时支持控制流和数据流编程模型。 尽管控制流模型适用于许多情况，但数据流模型也适用于其他一些情况，例如，当代理接收数据并执行基于该数据的有效负载的操作时。
设置无标题 android
05-02
在Android应用中，设置无标题是一种简洁和美观的设计方式，可以为用户提供更流畅的使用体验。以下是一些实现设置无标题的方法： 1. 在Android的Manifest文件中，找到对应Activity的标签，并添加”android:theme="@...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值