How to run dgemms subroutine on mac computer

[missfangula]

Homework Statement

I wrote a program to run the dgemul subroutine on my mac to find the matrix product for matrix a and matrix b, outputted as matrix c, just like the reference manual told me to. Being a mac user, I added the accelerate framework to the file which contain all the lapack stuff, but it does not recognize 'dgemul' as a function. Can someone please enlighten me as to what I am doing wrong? Below is the completed program.

Homework Equations

The Attempt at a Solution

#import <Foundation/Foundation.h>
#include <Accelerate/Accelerate.h>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>

int main ()
{

double a[16] =
{
1.0, 2.0, 3.0, 4.0,
1.0, 2.0, 3.0, 4.0,
1.0, 2.0, 3.0, 4.0,
1.0, 2.0, 3.0, 4.0

};

double b[16] =
{
5.0, 4.0, 3.0, 2.0,
5.0, 4.0, 3.0, 2.0,
5.0, 4.0, 3.0, 2.0,
5.0, 4.0, 3.0, 2.0
};

double c[16] =
{
0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0
};

int lda = 4;
int ldb = 4;
char transa = 'N';
char transb = 'N';
int ldc = 4;
int l = 4;
int m = 4;
int n = 4;

dgemul(a, lda, transa, b, ldb, transb, c, ldc, l, m, n);
int i;
for (i=0; i<16; ++i) printf("%5.10f\n", c);


return 0;
}


Thank you!
-miss fangula

 

[missfangula]

sorry, i meant the dgemul subroutine in the title, not dgemms

 

[vela]

I took a quick look at lapack and didn't see a routine called dgemul included.

 

[missfangula]

where did you find that information? this is where I got the dgemul thing:

http://publib.boulder.ibm.com/infoc...luster.essl43.guideref.doc/am501_hsgemul.html

 

[missfangula]

what subroutine should i call to do a matrix-matrix product?

 

[vela]

where did you find that information? this is where I got the dgemul thing:

http://publib.boulder.ibm.com/infoc...luster.essl43.guideref.doc/am501_hsgemul.html

That page describes ESSL, not LAPACK.

You can check the LAPACK documentation for the routine you want, but as far as I can tell, it doesn't include a matrix-multiplication subroutine.