Студопедия

КАТЕГОРИИ:


Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748)

Методы Крыловского типа





Выбор структуры разреженности

Поскольку предобусловливание является способом уменьшения количества арифметических операций при итерационном решении СЛАУ, за счет улучшения сходимости, следовательно, выбор структуры (портрета) разреженности матрицы предобусловливания является одной из важнейших практических задач. Одной из простейших стратегий является выбор заранее известной структуры (структурная предфильтрация), например, ленточной или блочно-диагональной т.е. элементы матрицы М совпадают с ленточной (блочно-диагональной) частью матрицы А, и равны нулю вне данной структуры. Данная предфильтрация характеризуется шириной ленты (k) (размером блока). Другой более часто используемой стратегией является динамическое определение структуры разреженности, основанное на некотором пороге (алгебраическая предфильтрация). Алгебраическая предфильтрация заключается в нахождении определенного порога, с помощью которого происходит отбрасывание (обнуление) малозначащих элементов путем сравнения модуля каждого элемента со значением полученного порога.

Упомянутый порог обнуления, в свою очередь, может быть получен многими способами. В данной версии системы реализована следующие стратегии

mij = 0, если |aij| < ei = ||ai*||2. t,  

которая основана на построчном нахождении евклидовой нормы строки далее умножением полученного значения на значение задаваемой точности обнуления (t) получим, собственно, значение порога (ei) для определения элементов матрицы М соответствующей строки. Далее следует перейти к следующей строке и осуществить вышеописанные действия, для текущей строки. Процесс закончится, когда все строки будут обработаны подобным образом. Таким образом, при таком подходе, порог обнуления является своим для каждой строки. Также реализовано еще два подхода. Первый способ аналогичен вышеописанному, только вместо строк обрабатываются столбцы матрицы

mij = 0, если |aij| < ej = ||a*j||2. t.  

Второй способ является более глобальным, поскольку порог обнуления, при таком подходе, является одним для всех элементов матрицы, т.е. первоначально вычисляется, евклидова норма матрицы далее, как и ранее умножением полученного значения на значение задаваемой точности получим искомый порог обнуления (e).

Ранее был построен класс проекционных методов, основой которых служит условие Петрова-Галеркина (2.13) и формула (2.16). Далее были приведены примеры построения нескольких методов для достаточно простого случая одномерных подпространств K и L. Однако выбор m=1, как правило, дает медленную сходимость, и соответствующие методы оказываются пригодными лишь для СЛАУ небольшой размерности. В настоящее время широкое распространение получили проекционные методы, которые в качестве K используют описанные ранее подпространства Крылова размерности больше единицы. Такие методы называются методами крыловского типа.





Дата добавления: 2014-12-27; Просмотров: 501; Нарушение авторских прав?


Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет



ПОИСК ПО САЙТУ:


Рекомендуемые страницы:

Читайте также:
studopedia.su - Студопедия (2013 - 2020) год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! Последнее добавление
Генерация страницы за: 0.007 сек.