чмы Холесский

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
using System;
using System.Collections.Generic;
using System.IO;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
namespace Holessky
{
class Program
{
public static double[][] CholeskyDecomposition(double[][] A)
{
double[][] L = new double[A.Length][];
for (int i = 0; i < A.Length; i++)
{
L[i] = new double[i + 1]; //L - треугольная матрица, поэтому в i-ой строке i+1 элементов
double temp;
//Сначала вычисляем значения элементов слева от диагонального элемента,
//так как эти значения используются при вычислении диагонального элемента.
for (int j = 0; j < i; j++)
{
temp = 0;
for (int k = 0; k < j; k++)
{
temp += L[i][k] * L[j][k];
}
L[i][j] = (A[i][j] - temp) / L[j][j];
}
//Находим значение диагонального элемента
temp = A[i][i];
for (int k = 0; k < i; k++)
{
temp -= L[i][k] * L[i][k];
}
L[i][i] = Math.Sqrt(temp);
}
return L;
}
public static double[][] InData(String file_name)
{
using (StreamReader f = new StreamReader(file_name))
{
int n = Convert.ToInt32(f.ReadLine());
double[][] matrix = new double[n][];
for (int i = 0; i < n; i++) {
matrix[i] = new double[n];
string[] s = f.ReadLine().Split(' ');
for (int j = 0; j < n; j++) {
matrix[i][j] = Convert.ToDouble(s[j]);
}
}
return matrix;
}
}
static void Main(string[] args)
{
double[][] result = CholeskyDecomposition(InData("test_a.txt"));
for (int i = 0; i < result.Length; i++)
{
for (int j = 0; j <= i; j++)
System.Console.Write("{0} ", result[i][j]);
System.Console.WriteLine();
}
double[] x, y, b;
int n = 0;
using (StreamReader f = new StreamReader("test_b.txt"))
{
n = Convert.ToInt32(f.ReadLine());
b = new double[n];
y = new double[n];
x = new double[n];
string[] s = f.ReadLine().Split(' ');
for (int j = 0; j < n; j++) {
b[j] = Convert.ToDouble(s[j]);
}
}
for (int i = 0; i < n; i++)
{
y[i] = b[i];
for (int j = 0; j < i; j++)
{
y[i] -= result[i][j] * y[j];
}
y[i] /= result[i][i];
}
for (int i = n-1; i >= 0; i--)
{
x[i] = y[i];
for (int j = n - 1; j > i; j--)
{
x[i] -= result[j][i] * x[j];
}
x[i] /= result[i][i];
}
for (int i = 0; i < n; i++)
{
System.Console.Write("{0} ", x[i]);
}
System.Console.ReadKey();
}
}
}