Симметричное аппроксимированное минимальное сочетание степени
p = symamd(S)
p = symamd(S,knobs)
[p,stats] = symamd(...)
p = symamd(S)
для симметричного положительного определенного матричного S
, возвращает вектор сочетания p
таким образом, что S(p,p)
имеет тенденцию иметь более разреженный Фактор Холецкого, чем S
. Найти упорядоченное расположение для S
symamd
создает матричный M
таким образом, что spones(M'*M) = spones (S)
, и затем вычисляет p = colamd(M)
. symamd
функция может также работать хорошо на симметричные неопределенные матрицы.
S
должно быть квадратным; только на строго нижнюю треугольную часть ссылаются.
p = symamd(S,knobs)
где knobs
скаляр. Если S
n
- n
, строки и столбцы с больше, чем knobs*n
записи удалены до упорядоченного расположения и упорядочены в последний раз в выходном сочетании p
. Если knobs
параметр не присутствует, затем knobs = spparms('wh_frac')
.
[p,stats] = symamd(...)
дает дополнительный вектор stats
это обеспечивает данные об упорядоченном расположении и валидности матричного S
.
stats(1) | Количество плотных или пустых строк проигнорировано |
stats(2) | Количество плотных или пустых столбцов проигнорировано |
stats(3) | Количество сборок мусора, выполняемых на внутренней структуре данных, используется |
stats(4) | 0 |
stats(5) | Индекс крайнего правого столбца, который не отсортирован или содержит дублирующиеся записи или |
stats(6) | В последний раз замеченная дублирующаяся или неисправная строка индексирует в индексе столбца, данном |
stats(7) | Количество дублирующихся и неисправных индексов строки |
Несмотря на то, что, MATLAB® встроенные функции генерируют допустимые разреженные матрицы, пользователь может создать недопустимую разреженную матрицу с помощью MATLAB C или API Фортрана и передать его symamd
. Поэтому symamd
проверяет тот S
isvalid:
Если индекс строки появляется два или больше раза в том же столбце, symamd
игнорирует дублирующиеся записи, продолжает обрабатывать и предоставляет информацию о дублирующихся записях в stats(4:7)
.
Если индексы строки в столбце не работают, symamd
виды каждый столбец его внутренней копии матричного S
(но не восстанавливает входную матрицу S
), продолжает обрабатывать и предоставляет информацию о неисправных записях в stats(4:7)
.
Если S
недопустимо любым другим способом, symamd
не может продолжиться. Это распечатывает сообщение об ошибке и не возвращает выходных аргументов (p
или stats
).
Упорядоченное расположение сопровождается симметричным поступорядоченным расположением дерева устранения.
Авторы кода для symamd
Стефан Ай. Лэримор и Тимоти А. Дэвис (davis@cise.ufl.edu
), Флоридский университет. Алгоритм был разработан в сотрудничестве с Джоном Гильбертом, Xerox PARC, и Эсмондом Ыном, Окриджской национальной лабораторией. Исследование Алгоритмов Разреженной матрицы во Флоридском университете: https://www.cise.ufl.edu/research/sparse/