| 49 | |
| 50 | template<typename T> |
| 51 | Array<T> inverse(const Array<T> &in) { |
| 52 | int M = in.dims()[0]; |
| 53 | int N = in.dims()[1]; |
| 54 | |
| 55 | if (M != N) { |
| 56 | Array<T> I = identity<T>(in.dims()); |
| 57 | return solve(in, I); |
| 58 | } |
| 59 | |
| 60 | Array<T> A = copyArray<T>(in); |
| 61 | Array<int> pivot = lu_inplace<T>(A, false); |
| 62 | |
| 63 | auto func = [=](Param<T> A, Param<int> pivot, int M) { |
| 64 | getri_func<T>()(AF_LAPACK_COL_MAJOR, M, A.get(), A.strides(1), |
| 65 | pivot.get()); |
| 66 | }; |
| 67 | getQueue().enqueue(func, A, pivot, M); |
| 68 | |
| 69 | return A; |
| 70 | } |
| 71 | |
| 72 | #define INSTANTIATE(T) template Array<T> inverse<T>(const Array<T> &in); |
| 73 | |