В поиске метода решения больших разреженных матриц (МКЭ) столкнулся с методом сопряженных градиентов с ILU-предобуславливанием. Разобрался как сделать ILU-разложение для исходной матрицы, разобрался в методе сопряженных градиентов. Но мне не понятно, как использовать предобуславливание (ILU-разложение) для итерационных методов... Ведь предобуславливание используется для улучшения и ускорения сходимости (для итерационных методов), но тогда получается две матрицы L и U, как здесь применять метод сопряженных градиентов? Или его надо использовать сначала для прямого хода, а затем для обратного, как в прямых методах? Вопрос для многих элементарный, я уверен, но уже неделю не могу разобраться с этим)
|