Understanding LAPACK calls in C++ with a simple example
You need to factor the matrix (by calling dgetrf) before you can solve the system using dgetrs. Alternatively, you can use the dgesv routine, which does both steps for you. By the way, you don’t need to declare the interfaces yourself, they are in the Accelerate headers: // LAPACK test code #include <iostream> #include <vector> … Read more