Параллельный решатель трёхдиагональных СЛУ

[Karoed]

Решаю задачку моделирования численными методами, которая сводится к набору трёхдиагональных СЛУ (с диагональным преобладанием всё нормально). Ну и как обычно упёрся в производительность.
Фактически мне нужен метод прогонки, но такой чтобы мог работать параллельно.
Я знаю, что есть много статей по распределению метода прогонки, а реализаций я что-то так сразу не нашел.
Может быть кто-то сталкивался с такой проблемой и может что-то подсказать?
Главное только, чтобы реализация была достаточно эффективная, язык и платформа сильно не важны.

Спасибо.

[wanderers]

Если у вас набор матриц, то запускайте решатель на каждом ядре отдельно. Другое дело, когда матрица одна, но невероятных размеров.

[Karoed]

Матрица конечно не одна, но у меня довольно подробные сетки (интересующий меня размер примерно 10 000x10 000).
В некоторых местах я действительно могу совершить прогонку по каждой из матриц независимо вдоль каждого измерения. В этом совсем нет сложности.
Но в данном случае вопрос в том как решать параллельно задачу прогонки или вернее чем решить.