Stack Overflow
Stack Overflow

Questions tagged [lapack]

LAPACK (Linear Algebra PACKage) is a software library package to solve linear algebra equations. LAPACK is written in Fortran 90.

LAPACK provides routines for solving systems of simultaneous linear equations, least-squares solutions of linear systems of equations, eigenvalue problems, and singular value problems. LAPACK is written in Fortran 90.

1025 questions
17
votes
1 answer

Is it possible to wrap a function from a shared library using F2PY?

I'm developing a package that requires Python bindings for the dgtsv subroutine from the LAPACK Fortran library. At the moment, I'm distributing the Fortran source file, dgtsv.f, alongside my Python code, and using numpy.distutils to automatically…
ali_m
  • 71,714
  • 23
  • 223
  • 298
16
votes
2 answers

c++11 std::array vs static array vs std::vector

First question, is it a good thing to start using c++11 if I will develop a code for the 3 following years? Then if it is, what is the "best" way to implement a matrix if I want to use it with Lapack? I mean, doing std::vector >…
PinkFloyd
  • 2,103
  • 4
  • 28
  • 47
16
votes
4 answers

What is a good free (open source) BLAS/LAPACK library for .net (C#)?

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…
Egil Hansen
  • 15,028
  • 8
  • 37
  • 54
14
votes
1 answer

Accessing submatrices using LAPACK

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?
user521968
13
votes
3 answers

Can I use Lapack for calculating the eigenvalues and eigenvectors of large sparse matrices?

If I had a square matrix that is 1,000 by 1,000 could Lapack calculate the eigenvectors and eigenvalues for this matrix? And if it can how long would it take? Also what about for a 10,000 by 10,000 matrix or even a 1,000,000 by 1,000,000 matrix?…
Spencer
  • 245
  • 2
  • 5
  • 9
12
votes
3 answers

How can I determine which matrix libraries my R install is using?

I am having a matrix error when using the computer cluster at my university that I cannot reproduce on my local machine. I think it might be due to a difference of matrix libraries (BLAS, LAPACK, ATLAS, etc.). I don't know much about these…
rcorty
  • 1,140
  • 1
  • 10
  • 28
12
votes
1 answer

Using Lapack with 128 bit precision

I am trying to use Lapack for a 128 bit precision calculation of a matrix singular value decomposition (SVD) and I found out that there is some black compiler magic to accomplish this. The Intel Fortran compiler (ifort) supports the option -r16…
Maxwell
  • 303
  • 1
  • 7
12
votes
2 answers

Parallel linear algebra for multicore system

I'm developing a program that needs to do heavy linear algebra calculations. Now I'm using LAPACK/BLAS routines, but I need to exploit my machine (24 core Xeon X5690). I've found projects like pblas and scalapack, but they all seem to focus on…
Patrik
  • 845
  • 2
  • 8
  • 20
11
votes
5 answers

Trouble installing scipy via pyCharm windows 8 - no lapack / blas resources found

I'm currently having trouble installing scipy via PyCharm's package manager. I have installed numpy successfully and do have the Microsoft Visual Studio C/C++ compiler in the System Variables. However, when it's time to install scipy in PyCharm, the…
MCKong
  • 113
  • 1
  • 1
  • 4
10
votes
1 answer

LAPACK: Are operations on packed storage matrices faster?

I want to tridiagonalize a real symmetric matrix using Fortran and LAPACK. LAPACK basically provides two routines, one operating on the full matrix, the other on the matrix in packed storage. While the latter surely uses less memory, I was wondering…
Michael Becker
  • 355
  • 1
  • 2
  • 10
10
votes
1 answer

LAPACK SVD (Singular Value Decomposition)

Do yo know any example to use LAPACK To calculate SVD?
edgarmtze
  • 24,683
  • 80
  • 235
  • 386
10
votes
1 answer

Pre-compiled Windows OMF BLAS/LAPACK?

Is there anywhere I can get pre-compiled BLAS and LAPACK binaries for Windows in OMF object format? I want to link some D language code to these. I'm aware of where I can get the relevant libs in COFF format. I'm also aware of the objconv tool,…
dsimcha
  • 67,514
  • 53
  • 213
  • 334
10
votes
1 answer

Will fortran's 'matmul' make use of MKL if I include the library?

I am writing some code right now and I have a placeholder with matmul that seems to be working pretty well, but I'd like to use a LAPACK dgemm implementation. I am only using gfortran right now and getting very good speeds with matmul, but I wonder…
drjrm3
  • 4,474
  • 10
  • 53
  • 91
10
votes
2 answers

Unable to import numpy: Error: /usr/lib/liblapack.so.3: undefined symbol: gotoblas

When I try to import numpy, I get the following error: /usr/local/lib/python2.7/dist-packages/numpy/linalg/__init__.py in () 49 from .info import __doc__ 50 ---> 51 from .linalg import * 52 53 from numpy.testing import…
vkmv
  • 1,345
  • 1
  • 14
  • 24
10
votes
4 answers

Shouldn't LAPACKs dsyevr function (for eigenvalues and eigenvectors) be thread-safe?

While trying to compute eigenvalues and eigenvectors of several matrices in parallel, I found that LAPACKs dsyevr function does not seem to be thread safe. Is this known to anyone? Is there something wrong with my code? (see minimal example…
dastrobu
  • 1,600
  • 1
  • 18
  • 32
1
2
3
…
68 69