Stack Overflow
Stack Overflow

Questions tagged [scalapack]

ScaLAPACK is a library of high-performance linear algebra routines for parallel distributed memory machines. ScaLAPACK solves dense and banded linear systems, least squares problems, eigenvalue problems, and singular value problems.

The goals of the ScaLAPACK project are the same than the one’s of LAPACK, namely:

efficiency (to run as fast as possible),

scalability (as the problem size and number of processors grow),

reliability (including error bounds),

portability (across all important parallel machines),

flexibility (so users can construct new routines from well-designed parts),

and ease of use (by making the interface to LAPACK and ScaLAPACK look as similar as possible).

62 questions
0
votes
1 answer

Are ScaLAPACK subroutines 'blocking' in MPI lingo?

I feel like this is a dumb question, but I really don't know how to google this, since every combination of the words 'ScaLAPACK' and 'block' just bombards you with information about the block-cyclic matrix distribution ScaLAPACK uses. However, what…
Jakob Filser
  • 131
  • 9
0
votes
0 answers

Writing a block-cyclic distributed matrix to a binary file

I am using Scalapack for diagonalizing matrices (25k × 25k) and need the eigenvectors. I am using pzlawrite to write out the distributed matrices to a file. But the speed of writing is quite slow in ASCII for such huge files. I suppose writing the…
Unpotating Potato
  • 21
  • 2
0
votes
1 answer

Compile the first example of ScaLapack

I am trying to compile the first example of ScaLapack. For that I am using a makefile to compile and I am having these error's message: I already did all the procedures to install ScaLapack but I think that I have to do some modification in my…
Mateus
  • 149
  • 7
0
votes
0 answers

How to customize makefile file to compile BLACS library

I have a makefile file to compile my codes in Fortran. This makefile already compile Lapack library as well. In view of I already have BLACS library install in my system how to modify this makefile in order to compile BLACS library as well. My…
Mateus
  • 149
  • 7
0
votes
1 answer

"ld: library not found for -lblacsF77init"

I am trying to compile my code using mpi on my macbook but I get the following error: "ld: library not found for -lblacsF77init". I don't understand where I can find this library and how to create a correct path. This is the makefile: # serial…
David
  • 1
0
votes
0 answers

Kronecker tensor product with identity matrix in fortran using mpi

I need to take Kronecker product of two matrices one of them is an identity matrix, my system size happens to be of the order of 256*256 for both the matrices. Although I've written a code which works fine but the problem is with the system size. It…
Puneet Sharma
  • 1
  • 1
0
votes
1 answer

python ctypes link multiple shared library with example gsl gslcblas

I want to use some functions from a shared library in python. From the python doc, I know ctypes is a good choice. However such library have some undefined symbols and I should link it to another shared library to get the symbols. In g++, it is…
ken
  • 318
  • 1
  • 11
0
votes
1 answer

Blocks of different sizes in ScaLAPACK?

I am performing a Cholesky factorization with Intel-MKL, which uses ScaLAPACK. I distributed the matrix, based on this example, where the matrix is distributed in blocks, which are of equal size (i.e. Nb x Mb). I tried to make it so that every block…
gsamaras
  • 71,951
  • 46
  • 188
  • 305
0
votes
1 answer

On entry to NIT parameter number 9 had an illegal value

I go this ex1.c file from Intel 11. However, when I execute it, it fails: gsamaras@pythagoras:~/konstantis$ ../mpich-install/bin/mpicc -o test ex1.c -I../intel/mkl/include ../intel/mkl/lib/intel64/libmkl_scalapack_ilp64.a -Wl,--start-group…
gsamaras
  • 71,951
  • 46
  • 188
  • 305
0
votes
1 answer

Inplace restoration of upper/lower part of distributed lower/upper triangular matrix

I wonder if there exists any Scalapack routine that may convert a triangular matrix back to a general one, i.e. copies the lower/upper part to the upper/lower part - I couldn't find it yet. I know there exists an easy workaraound using p?geadd with…
zp3
  • 5
  • 2
0
votes
1 answer

Adapt dgemm example code to use sgemm (scalapack)

I need to make the following program (from http://www.netlib.org/scalapack/examples/pblas.tgz) work with SGEMM. What do I need to change to make it work? My knowledge of Fortran is quite limited, I'm pretty much treating this as a black-box and…
pldimitrov
  • 1,597
  • 2
  • 16
  • 21
0
votes
1 answer

Writing a configurable scalapack linear system solver that prints execution time

I'm trying to adapt the following example program to use as a coarse-grained parallel benchmark in my experiments. I added the following lines to the code: START_TIME = MPI_WTIME() * <- added this CALL PDGESV( N, NRHS, MEM( IPA ), 1, 1,…
pldimitrov
  • 1,597
  • 2
  • 16
  • 21
0
votes
1 answer

Problems compiling example scalapack application

When I do: mpif77 example1.f -L scalapack/scalapack-1.8.0/ -lscalapack -L scalapack/blacs/BLACS/LIB -l:blacsF77init_MPI-LINUX-0.a -l:blacs_MPI-LINUX-0.a -l:blacsF77init_MPI-LINUX-0.a -L scalapack/blas/BLAS/ -l:blas_LINUX.a -L scalapack/lapack/…
pldimitrov
  • 1,597
  • 2
  • 16
  • 21
0
votes
1 answer

results of linear system using scalapack changing with processors number

I have used Jonathan's code and tried to solve A*X=B, but my results are changing based on the number of processors. Can anyone please help me in this regard. ! Use MPI-IO to read a diagonal matrix distributed block-cycliCALLy, ! use Scalapack to…
ziv
  • 189
  • 1
  • 9
0
votes
2 answers

Calling BLACS with more processes than used

I want to create a parallel program, which makes heavy use of SCALAPACK. The basis of SCALAPACK is BLACS, which itself relies on MPI for interprocess communication. I want to start the program with a defined number of processes (e.g. the number of…
Stefan
  • 2,460
  • 1
  • 17
  • 33
1 2 3
4
5