| 14 | using namespace af; |
| 15 | |
| 16 | int main(int argc, char* argv[]) { |
| 17 | try { |
| 18 | // Select a device and display arrayfire info |
| 19 | int device = argc > 1 ? atoi(argv[1]) : 0; |
| 20 | af::setDevice(device); |
| 21 | af::info(); |
| 22 | |
| 23 | printf("Running QR InPlace\n"); |
| 24 | array in = randu(5, 8); |
| 25 | af_print(in); |
| 26 | |
| 27 | array qin = in.copy(); |
| 28 | |
| 29 | array tau; |
| 30 | qrInPlace(tau, qin); |
| 31 | |
| 32 | af_print(qin); |
| 33 | af_print(tau); |
| 34 | |
| 35 | printf("Running QR with Q and R factorization\n"); |
| 36 | array q, r; |
| 37 | qr(q, r, tau, in); |
| 38 | |
| 39 | af_print(q); |
| 40 | af_print(r); |
| 41 | af_print(tau); |
| 42 | |
| 43 | } catch (af::exception& e) { |
| 44 | fprintf(stderr, "%s\n", e.what()); |
| 45 | throw; |
| 46 | } |
| 47 | |
| 48 | return 0; |
| 49 | } |