Часть (и) почти блочно-диагональной матрицы
[nb,rows,ncols,last,blocks] = bkbrk(blokmat)
bkbrk(blokmat)
[nb,rows,ncols,last,blocks] = bkbrk(blokmat)
возвращает детали почти блочно-диагональной матрицы, содержащейся в blokmat
, с rows
и last
nb
-векторы и blocks
матрица размера [sum(rows),ncols]
.
Эта служебная программа вряд ли будет интересна случайному пользователю. Он используется в slvblk
для декодирования информации, предоставляемой spcol
, о сплайн-коллокации матрицы в почти блоке диагонали форме, особенно подходящей для сплайнов. Но bkbrk
может также декодировать почти блочно-диагональную форму, используемую в [1].
bkbrk(blokmat)
ничего не возвращает, но детали распечатываются. Это полезно при попытке понять, что пошло не так с такой матрицей.
[1] C. de Boor and R. Weiss. «SOLVEBLOK: Пакет для решения почти блок-диагональных линейных систем». ACM Trans. Mathem. Программное обеспечение 6 (1980), 80-87.