"> файл зі значеннями вектора швидкості в кожній вузловій точці.
Вихідні дані :
файл зі значеннями концентрації в кожній вузловій точці.
Основні реалізовані функції:
Звернення матриці загального вигляду
bool inverse (double ** a, int n);
Вхідні параметри :
A - матриця. масив з нумерацією елементів [1 .. N, 1 .. N]
N - розмірність матриці A
Вихідні параметри :
A - матриця, обернена до вихідної. Масив з нумерацією елементів [1 .. N, 1 .. N]
Результат :
True, якщо вихідна матриця невироджена.
False, якщо вихідна матриця вироджена.
LU-розкладання матриці загального вигляду розміру M x N
void ludecomposition (double ** a, int m, int n, int * pivots);
Функція обчислює LU-розкладання прямокутної матриці загального вигляду з частковим вибором ведучого елементу (з перестановками рядків).
Вхідні параметри :
A - матриця A. Нумерація елементів: [1 .. M, 1 .. N]
M - число рядків в матриці A
N - число стовпців в матриці A
Вихідні параметри :
A - матриці L і U в компактній формі (див. нижче).
Нумерація елементів: [1 .. M, 1 .. N]
Pivots - матриця пере становок в компактній формі (див. нижче).
Нумерація елементів: [1 .. Min (M, N)]
Матриця A представляється, як A = P * L * U, де P - матриця перестановок, матриця L - ніжнетреугольная (або ніжнетрапецоідальная, якщо M> N) матриця, U - верхнетреугольная (або верхнетрапецоідальная, якщо M
Обчислення визначника матриці, заданої LU-азложеніем.
double determinantlu (double ** a, int * pivots, int n);
Вхідні параметри :
A - LU-розкладання матриці (результат роботи підпрограми LUDecomposition).
Pivots - таблиця перестановок, вироблених в ході LU-розкладання. (Результат роботи підпрограми LUDecomposition). p align="justify"> N - розмірність матриці
Результат:
визначник матриці
LDL ^ T розкладання симетричної матриці
<...
