Компактное хранение и сбор матрицы в МКЭ

[sergeik]

Здравствуйте!
При построении СЛАУ в МКЭ матрица получается симметричной и сильно разреженной. Есть несколько методов для компактного хранения разреженных матриц (CSR, Skyline и им подобные). Поделитесь, пожалуйста, опытом или ссылками на литературу кто и как собирает матрицу сразу в компактный вид, без первоначального сбора глобальной матрицы и ее дальнейшего сжатия в компактную? Как при этом внедряются условия типа Дирихле. Насколько я понимаю, основное время при сборе компактной матрицы без глобальной занимает поиск позиции соответствующего столбца, куда необходимо добавить вклад от другого узла. Можно сократить количество таких поисков, если сначала строить локальную матрицу,а потом ее элементы добавлять в глобальную.

P.S. Вариант с поузловым подходом не подходит, интересует поэлементный))