Решить полосчатую матричную систему уравнений

Мне нужно решить двумерное уравнение Пуассона, то есть систему уравнений в for AX=B, где A-матрица n-на-n, А B-вектор n-на-1. Будучи матрицей дискретизации для двумерной задачи Пуассона, я знаю, что только 5 диагоналей не будут нулевыми. Lapack не предоставляет функций для решения этой конкретной задачи, но у него есть функции для решения связанных матричных систем уравнений, а именно DGBTRF (для факторизации LU) и DGBTRS. Теперь 5 диагоналей: главная диагональ, первая диагональ. диагонали выше и ниже главной и двум диагоналям выше и ниже по диагонали м относительно главной диагонали. Прочитав документацию lapack о band storage, я узнал, что мне нужно создать матрицу (3*m+1)-by-n для хранения A в формате band storage, назовем эту матрицу AB. Теперь вопросы:

1) в чем разница между dgbtrs и dgbtrs_? Intel MKL предоставляет оба, но я не могу понять, почему

2) dgbtrf требует, чтобы Матрица хранения каналов была массивом. Должен Ли Я линеаризовать AB по строкам или по столбцам?

3) является ли это правильным способом вызова двух функций?

int n, m;
double *AB;
/*... fill n, m, AB, with appropriate numbers */

int *pivots;
int nrows = 3 * m + 1, info, rhs = 1;
dgbtrf_(&n, &n, &m, &m, AB, &nrows, pivots, &info);
char trans = 'N';
dgbtrs_(&trans, &n, &m, &m, &rhs, AB, &nrows, pivots, B, &n, &info);