lapack

I am writing a Python module that includes Cython extensions and uses LAPACK (and BLAS ). I am open to using either clapack or lapacke , or some kind of f2c or f2py solution if necessary. What is important is that I am able to call lapack and blas routines from Cython in tight loops without Python call overhead. I've found one example here . However, that example depends on SAGE. I want my module to be installable without installing SAGE, since my users are not likely to want or need SAGE for anything else. My users are likely to have packages like numpy, scipy, pandas, and scikit learn

I have a project written in C# where I need to do various linear algebraic operations on matrices (like LU-factorization). Since the program is mainly a prototype created to confirm a theory, a C# implementation will suffice (compared to a possibly speedier C++ one), but I would still like a good BLAS or LAPACK library available to save me some coding. Long story short, can anybody recommend a free/open source BLAS or LAPACK library for use with .net? Best regards Egil. Update: Found Math.NET Numerics today, looks interesting, anybody have any experience with that? AMD's ACML is a free

Is there a function in LAPACK, which will give me the elements of a particular submatrix? If so how what is the syntax in C++? Or do I need to code it up? Stephen Canon There is no function for accessing a submatrix. However, because of the way matrix data is stored in LAPACK routines, you don't need one. This saves a lot of copying, and the data layout was (partially) chosen for this reason: Recall that a dense (i.e., not banded, triangular, hermitian, etc) matrix in LAPACK is defined by four values: a pointer to the top left corner of the matrix the number of rows in the matrix the number of

I don't understand how BLAS, LAPACK and ATLAS are related and how I should use them together! I have been looking through all of their manuals and I have a general idea of BLAS and LAPACK and how to use them with the very few examples I find, but I can't find any actual examples using ATLAS to see how it is related with these two. I am trying to do some low level work on matrixes and my primary language is C. First I wanted to use GSL, but it says that if you want the best performance you should use BLAS and ATLAS. Is there any good webpage giving some nice examples of how to use these (in C)

I installed atlas, blas and lapack x86_64 packages via yum install atlas.x86_64 blas.x86_64 lapack.x86_64 on a Redhat 6.6 (ii) distro which installs a shared library but am having problems compiling and linking. For example, if I try to compile the minimal working example: program main print *, 'hello world' end program main using gfortran -L. main.f90 -llapack -lblas -o main the compiler doesn't find the libraries and I get the error message: /usr/bin/ld: cannot find -llapack collect2: ld returned 1 exit status I'm relatively new to fortran and linux so I'm probably missing something obvious.

I'm new to numerical linear algebra, and I have just started using LAPACK and BLAS. Is there a routine that can copy/convert a symmetric matrix between packed and full storage? I have found dtrttp , which I can use to convert a double-precision full symmetric matrix to packed storage. However, these routines are meant for triangular matrices, so the corresponding dtpttr only fills a triangle of the full matrix. How can I fill the other half? The obvious solution is to symmetrize the matrix by a "home-made/diy" code, the risk being to re-invent the wheel. It is pretty easy to write the for