Форум программистов, компьютерный форум CyberForum.ru

Поиск универсального алгоритма LU разложения - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Исправить код, реализующий алгоритм сортировки http://www.cyberforum.ru/cpp-beginners/thread830317.html
Доброе утро. Сделал попытку реализовать функцию сортировки простым двухпутевым слиянием, но не вышло. При запуске происходит ошибка, имправить которую никак не удаётся. Вот код с комментариями: //Метод простого двухпутевого слияния void SimpleMerge (int * data, int size) { int * p_write = new int ; //вспомогательный массив для записи int * p_read(data);
C++ Работа с матрицами. Переписать вторую матрицу, так , чтобы элемент, находящийся на пересечении 1 строки и 2 столбца равнялся 1... у меня есть пустая матрица a. И есть матрица вида 1 2 1 4 1 7 2 4 Как мне переписать вторую матрицу, так , чтобы элемент, находящийся на пересечении 1 строки и 2 столбца равнялся 1, на пересечении 1 строки и 4 столбца равнялся 1 и т.д., до конца второй матрицы. А все оставшиеся эементы равнялись 0? http://www.cyberforum.ru/cpp-beginners/thread830313.html
C++ Проверить корректность ввода элементов матрицы
Ввод элементов матрицы А(m, n) осуществляется в произвольном порядке тройками чисел < i, j, Aij>. Признаком конца ввода служат три нуля: <0, 0, 0>. Проверить корректность такого ввода: все ли элементы введены, нет ли попытки повторного ввода или указания несуществующих координат i и j. Указание. Разрешается выделение дополнительного (рабочего) массива такой же размерности, что и исходный...
C++ Передача неопределенного числа параметров в функцию
В общем вот программка реализовано так,чтобы узнать когда конец переданных параметров, последним передается 0, и проверяется в ф-ии на 0. #include<iostream.h> #include<conio.h> float avg(float a...) { float* p; p=&a; float sum=0,count=0;
C++ Не пойму как из цикла выйти http://www.cyberforum.ru/cpp-beginners/thread830301.html
в учебнике такой пример. #include <stdio.h> main() { long nc; nc = 0; while ( getchar() != EOF)
C++ Класс // Класс "Армия", содержащий все типы боевых единиц class Army { public: vector<Infantryman> vi; //Infantryman, Archer, Horseman, Catapult, Elephant это классы с методом info, который //просто выводит названия войск vector<Archer> va; //это вектор типа класс Archer vector<Horseman> vh; vector<Catapult> vc; vector<Elephant> ve; подробнее

Показать сообщение отдельно
meraxujiep
 Аватар для meraxujiep
40 / 40 / 4
Регистрация: 31.01.2013
Сообщений: 133
07.04.2013, 16:08  [ТС]     Поиск универсального алгоритма LU разложения
Цитата Сообщение от palva Посмотреть сообщение
Обычно ищут максимальный по модулю элемент. Вы же находите первый не нулевой. Если он окажется очень близок к нулю, то теория предсказывает неприятности.
Согласен, промах.

Цитата Сообщение от palva Посмотреть сообщение
Информацию о перестановке строк вы не запоминаете? Тогда какая ценность этого алгоритма?
Задача стоит сделать LU-разложение матрицы. По итогу получаю **PointU на верхнюю треугольную матрицу и **PointL на нижнюю треугольную матрицу с единичками на диагонали. Как мне кажется этот кусок можно использовать для нахождения определителя например, или для решения СЛУ(это собираюсь попробовать реализовать).


Цитата Сообщение от palva Посмотреть сообщение
Надеюсь, вы понимаете что делаете, когда объявляете параметр double **PointU и в курсе, куда на самом деле обратится программа, когда увидит PointU[i][i] То, как вы размещаете матрицы в памяти довольно странно.
Не понял что вы имели ввиду в этой части.
 
Текущее время: 17:24. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru