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
6
votes
1 answer

errors using lapack C header in C++ with visual studio 2010

Please help me! It takes me hours to look up in the internet and I haven't found a solution.... I am trying to use the call lapack function from C++ functions but I failed at the very beginning. Here is my code: #include "stdafx.h" #include…
Jason
  • 1,200
  • 1
  • 10
  • 25
6
votes
2 answers

Wrapping a LAPACKE function using Cython

I'm trying to wrap the LAPACK function dgtsv (a solver for tridiagonal systems of equations) using Cython. I came across this previous answer, but since dgtsv is not one of the LAPACK functions that are wrapped in scipy.linalg I don't think I can…
ali_m
  • 71,714
  • 23
  • 223
  • 298
6
votes
2 answers

Use Emscripten with Fortran: LAPACK binding

My goal is to use LAPACK with Emscripten. My question is: how to port LAPACK to JS? The are two ways I can think of: CLAPACK to JS where my question is: does anybody know an unofficial version that is later than 3.2.1? And the other way to think of…
japedo
  • 368
  • 3
  • 15
6
votes
1 answer

gfortran LAPACK "undefined reference" error

I installed LAPACK on Ubuntu by following the instruction, sudo apt-get install liblapack-dev thus I can find /usr/lib/libblas/libblas.a and /usr/lib/lapack/liblapack.a, and then tested it in gfortran with the randomsys1 example, gfortran…
Randel
  • 320
  • 1
  • 6
  • 17
6
votes
2 answers

Stack overflow when solving large system using Julia

I am trying to solve a random linear system with a large square system matrix using Octave and Julia. Because the syntax of Octave and Julia are quite similar I run the following code in both a Octave shell and a Julia shell: N = 5000; A = rand(N,…
Aeronaelius
  • 1,291
  • 3
  • 12
  • 31
6
votes
3 answers

Undefined Reference to LAPACK wrapper of Armadillo

I have a question regarding the use of Armadillo. I'm using Ubuntu 12.10, and the gcc compiler in Code::Blocks. I installed LAPACK and BLAS using the synaptic package manager. I also installed Armadillo once using the synaptic package manager, once…
user2477647
  • 61
  • 1
  • 1
  • 3
6
votes
4 answers

(Pseudo)-Inverse of N by N matrix with zero determinant

I would like to take the inverse of a nxn matrix to use in my GraphSlam. The issues that I encountered: .inverse() Eigen-library (3.1.2) doesn't allow zero values, returns NaN The LAPACK (3.4.2) library doesn't allow to use a zero determinant, but…
Toon
  • 19
  • 1
  • 1
  • 3
5
votes
2 answers

BLAS equivalent of a LAPACK function for GPUs

In LAPACK there is this function for diagonalization SUBROUTINE DSPGVX( ITYPE, JOBZ, RANGE, UPLO, N, AP, BP, VL, VU, $ IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, $ IFAIL, INFO ) * I am looking for its GPU…
Open the way
  • 26,225
  • 51
  • 142
  • 196
5
votes
1 answer

Using BLAS or LAPACK with Xcode

I'm looking for a good guide on how to incorporate BLAS or LAPACK functions into my Objective C Program developed through Xcode. The only sources I can find online of programs in BLAS/LAPACK are written in Fortran - this is confusing. How does one…
Shaddy
  • 155
  • 1
  • 3
  • 6
5
votes
1 answer

how to use lapack under windows

I want to use lapack and make C++ matrix wrapper for it, but lapack is written in Fortran, there are some clapack but I want to use it from source. compile firstly *.f and *.cpp files to object files then link it into a application.. the following…
nurmurat
  • 151
  • 1
  • 11
5
votes
2 answers

numpy matrix mult does not work when it is parallized on HPC

I have two dense matrices with the sizes (2500, 208) and (208, 2500). I want to calculate their product. It works fine and fast when it is a single process but when it is in a multiprocessing block, the processes stuck in there for hours. I do…
rando
  • 365
  • 3
  • 12
5
votes
0 answers

Eigensystem of parametrized, hermitian matrix in python

Suppose we are interested in the eigenvalues and eigenvectors of a hermitian matrix h(t) that depends on a parameter t. My matrix is large and sparse and hence needs to be treated numerically. A naive approach is to evaluate the matrix h(t_k) at…
tree frog
  • 53
  • 4
5
votes
1 answer

Best performance method when solving 7000x7000 linear system with python

I’m in the need of an efficient method to inverse a 7000x7000 aerodynamic influence coefficient (dense) matrix in python. I had started before a FORTRAN routine to handle the problem using the LU decompostition routines from LAPACK, which I had seen…
Pedro Secchi
  • 109
  • 5
5
votes
1 answer

Alpine Linux can't install lapack-dev on python:3.5-alpine3.4

I'm using docker python:3.5-alpine3.4 image and trying to install lapack-dev but it keeps failing. It is complaining that it can't find libgfortran.so.5. However, I've tried installing libgfortran and that does not seem to fix the problem. (1/1)…
James
  • 2,488
  • 2
  • 28
  • 45
5
votes
2 answers

is R creating too many threads on startup

Every invocation of R is creating 63 sub processes Rscript --vanilla -e 'Sys.sleep(5)' & pstree -p $! | grep -c '{R}' # 63 where pstree looks something like this R(2562809)─┬─{R}(2562818) …
Will
  • 1,206
  • 9
  • 22
1 2 3
…
68 69