| 161 | |
| 162 | template<typename T> |
| 163 | Array<T> solveLU(const Array<T> &A, const Array<int> &pivot, const Array<T> &b, |
| 164 | const af_mat_prop options) { |
| 165 | UNUSED(options); |
| 166 | int N = A.dims()[0]; |
| 167 | int NRHS = b.dims()[1]; |
| 168 | |
| 169 | Array<T> B = copyArray<T>(b); |
| 170 | |
| 171 | mapped_ptr<T> aPtr = A.getMappedPtr(); |
| 172 | mapped_ptr<T> bPtr = B.getMappedPtr(); |
| 173 | mapped_ptr<int> pPtr = pivot.getMappedPtr(); |
| 174 | |
| 175 | getrs_func<T>()(AF_LAPACK_COL_MAJOR, 'N', N, NRHS, aPtr.get(), |
| 176 | A.strides()[1], pPtr.get(), bPtr.get(), B.strides()[1]); |
| 177 | |
| 178 | return B; |
| 179 | } |
| 180 | |
| 181 | template<typename T> |
| 182 | Array<T> triangleSolve(const Array<T> &A, const Array<T> &b, |
nothing calls this directly
no test coverage detected