Where can I find BLAS example code (in Fortran)?

Posted by Feynman on Stack Overflow See other posts from Stack Overflow or by Feynman
Published on 2011-01-09T05:47:28Z Indexed on 2011/01/09 5:53 UTC
Read the original article Hit count: 226

Filed under:

vector

|

matrix

|

fortran

|

linear-algebra

|

blas

I have been searching for decent documentation on blas, and I have found some 315 pages of dense material that ctrl-f does not work on. It provides all the information regarding what input arguments the routines take, but there are a LOT of input arguments and I could really use some example code. I am unable to locate any. I know there has to be some or no one would be able to use these libraries!

Specifically, I use ATLAS installed via macports on a mac osx 10.5.8 and I use gfortran from gcc 4.4 (also installed via macports). I am coding in Fortran 90. I am still quite new to Fortran, but I have a fair amount of experience with mathematica, matlab, perl, and shell scripting.

I would like to be able to initialize and multiply a dense complex vector by a dense symmetric (but not hermitian) complex matrix. The elements of the matrix are defined through a mathematical function of the indices--call it f(i,j).

Could anyone provide some code or a link to some code?

© Stack Overflow or respective owner

Related posts about vector

vector iterator not dereferencable at runtime on a vector<vector<vector<A*>*>*>

as seen on Stack Overflow - Search for 'Stack Overflow'
Hi, I have this destructor that create error at runtime "vector iterator not dereferencable". The gridMatrix is a std::vector<std::vector<std::vector<AtomsCell< Atom<T> * > * > * > * > I added the typename and also the typedef but I still have the error. I will move… >>> More
Adding Functions to an Implementation of Vector

as seen on Stack Overflow - Search for 'Stack Overflow'
I have this implementation of vector that I've been working on for a few days using examples from a textbook: #include <iostream> #include <string> #include <cassert> #include <algorithm> #include <cstring> // Vector.h using namespace std; template <class T> class… >>> More
Making only the outer vector in vector<vector<int>> fixed

as seen on Stack Overflow - Search for 'Stack Overflow'
I want to create a vector<vector<int>> where the outer vector is fixed (always containing the same vectors), but the inner vectors can be changed. For example: int n = 2; //decided at runtime assert(n>0); vector<vector<int>> outer(n); //outer vector contains n empty vectors outer… >>> More
(R) How can I get the complement of vector y in vector x

as seen on Stack Overflow - Search for 'Stack Overflow'
That's x \ y using mathematical notation. Suppose x <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,2,1,1,1,3) y <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1) How can I get a vector with ALL the values in x that are not in y. i.e the result should be: 2,1,1,3 There is a similar question here… >>> More
C++: Appending a vector to a vector

as seen on Stack Overflow - Search for 'Stack Overflow'
Assuming I have 2 STL vectors: vector<int> a; vector<int> b; Let's also say the both have around 30 elements. How do I add the vector b to the end of vector a? The dirty way would be iterating through b and adding each element via push_back, though I wouldn't like to do that! >>> More

Related posts about matrix

Find the "largest" dense sub matrix in a large sparse matrix

as seen on Stack Overflow - Search for 'Stack Overflow'
Given a large sparse matrix (say 10k+ by 1M+) I need to find a subset, not necessarily continuous, of the rows and columns that form a dense matrix (all non-zero elements). I want this sub matrix to be as large as possible (not the largest sum, but the largest number of elements) within some aspect… >>> More
How does a template class inherit another template class?

as seen on Stack Overflow - Search for 'Stack Overflow'
I have a "SquareMatrix" template class which inherits "Matrix" template class, like below: SquareMatrix.h: #ifndef SQUAREMATRIX_H #define SQUAREMATRIX_H #include "Matrix.h" template <class T> class SquareMatrix : public Matrix<T> { public: T GetDeterminant(); }; template… >>> More
iPhone OpenGL ES: How do I use gravity vector to correctly transform scene for augmented reality

as seen on Stack Overflow - Search for 'Stack Overflow'
I'm trying figure out how to get an OpenGL specified object to be displayed correctly according to the device orientation (ie. according to the gravity vector from the accelerometer, and heading from compass). The GLGravity sample project has an example which is almost like this (despite ignoring… >>> More
How do I use the gravity vector to correctly transform scene for augmented reality?

as seen on Stack Overflow - Search for 'Stack Overflow'
I'm trying figure out how to get an OpenGL specified object to be displayed correctly according to the device orientation (ie. according to the gravity vector from the accelerometer, and heading from compass). The GLGravity sample project has an example which is almost like this (despite ignoring… >>> More
C++ templates - matrix class

as seen on Stack Overflow - Search for 'Stack Overflow'
Hi I'm learning templates in C++ so I decied to write matrix class which would be a template class. In Matrix.h file I wrote #pragma once #include "stdafx.h" #include <vector> using namespace std; template<class T> class Matrix { public: Matrix(); ~Matrix(); GetDataVector(); … >>> More