【刘庆源码共享】稀疏线性系统求解算法MGMRES(m) 之 MGMRES类定义（C++）

最新推荐文章于 2023-03-19 18:53:27 发布

lq_0402

最新推荐文章于 2023-03-19 18:53:27 发布

阅读量1k

点赞数

文章标签：算法 c++ matrix delete null c

本文链接：https://blog.csdn.net/lq_0402/article/details/6007777

版权

本文档展示了C++中MGMRES类的定义，用于求解稀疏线性系统。类包括构造函数、析构函数、设置矩阵值、执行求解算法和验证结果的方法。通过读取文件设置系数矩阵和右端向量，并使用迭代法计算解，监控误差以判断收敛性。

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

/*
* Copyright (c) 2009 湖南师范大学数计院一心飞翔项目组
* All Right Reserved
*
* 文件名：mgmres.cpp
* 摘要：MGMRES类的函数定义
*
* 作者：刘庆
* 完成日期：2009年4月28日
*
*/
#include "stdafx.h"

#include <math.h>
#include <iostream>
#include <time.h>
#include "mgmres.h"
#include "EvalFile.h"

using namespace std;

MGMRES::MGMRES()
{
max_outer_itr = 1;
resultOk = 0;
}
MGMRES::~MGMRES()
{
if(resultOk!=NULL)
delete[] resultOk;
resultOk = NULL;
}

/* 替算法所涉及的矩阵赋初始值，所有数据在fileName的文件中 */
void MGMRES::SetMatrixsValue(const char* fileName)
{
EvalFile file = EvalFile(fileName);
file.SetWithSub(A); /* 替系数矩阵赋值 */
file.SetWithoutSub(b); /* 替右端向量赋值 n*1 */

if(A.GetTotal_col()<A.GetTotal_ln())
{
  A.SetTotal_col(A.GetTotal_ln());
  b.PosSetValue(A.GetTotal_col()-1, 0, 0);
}
if(A.GetTotal_col()>A.GetTotal_ln())
{
  A.PosSetValue(A.GetTotal_col()-1,0,0);
  b.PosSetValue(A.GetTotal_col()-1,0,0);
}
x = Matrix::GetE(A.GetTotal_col());   /* 将x的各个值初始化为0 */

// A.DisplayAsMatrix("A", 1);
}

/* A*B、A+B运算中A都会发生变化，B不会 */
Matrix& MGMRES::MGMRESCode()
{
Double mu=0, residue=0, h=0;
Double *c = new Double[A.GetTotal_ln()+1], *s =