It was working fine like this:
boolean compare_matrices(double *matrixC, double *matrixD, int rows, int cols)^M
{
const double EPSILON = 0.0000000001;
int row, col;
boolean flag...
Type: Posts; User: collymitch
It was working fine like this:
boolean compare_matrices(double *matrixC, double *matrixD, int rows, int cols)^M
{
const double EPSILON = 0.0000000001;
int row, col;
boolean flag...
Lookin good now thanks.
I'm trying to compare two matrices to see if they are equal. They are stored in the same manner, in that I use a double pointer.
This is my code
for ( row = 0; row < rows; ++row )
{...
Hey,
Background:
I've built a matrix multiplier in C based on the fortran DGEMM function.
It works great the only remaining modification I need to make is to add OpenMP pragmas so that it can...
Like I said they can be linked fine. It's just with the intel compiler that I can't get it to compile. f77 and g77 work fine
>f77 -o mm mm.c mmult.c dgemm.o xerbla.o lsame.o
compiles just file.
...
Hey,
Ok here goes...
I needed to separate a function (mmult) I had made from the driver program so I put it in its own .c file. The C driver though calls mmult.c and other Fortran functions so i...
Hey - Quick one for ya...
I get the warning message below at compile time because i'm trying to compare a double value to 0.
What should I use to compare it? I'm working on a matrix multiplier...
Excellent. That seems to be working now. I could have sworn that I'd tried that! maybe i didn't save it right or something. thanks Dave you've been a great help.
That makes it easier to read and...
Mmmm....
I modified it to take in the a[m][n] x b[n][p] = c[m][p] sizes.
I've not got it just right, it's close but close is wrong so...
void mmult(size_t m_rows, size_t n_common, size_t...
Yes that's right - i'm comparing my result with the Blas implementation which uses the DGEMM routine as described Here.
Ah-ha. That makes sense to me now.
So i'm gonna have
void...
Ok Dave, Thanks, that seems to work fine.
I was unsure how you got this part:
I went thorugh it step by step and i see what it's doing but i'm not completely sure why that works. Thing is...
Ok - i got a bit confused by that so had to ask my supervisor again...
Swoopy you were on the right track - I need to have a single pointer instead of a double so
would work.
This means that...
Hey,
I'm creating a matrix-matrix multiplier in C that will use the same format as the GOTO BLAS and ATLAS implementations. I need to allocate memory and store the matrices in COL MAJOR format not...
Hey I've got that set up and working fine now by the way. It was a great help it got my program kicked off and have made a good deal of progress so thanks for your help quzah!
Colm.
Ok....
Yeah the 'void main' bit was a typing mistake - i had called the function something else but changed it to main but forgot to change the rest. I was just trying to show a quick example....
Hey thanks,
Ok that helped a bit. I read that an array name is just a pointer to the beginning of the allocated memory space. Does this mean I treat the matrix just as an array?
I know that If...
Hey,
I'm working on an openMP enabled matrix-matrix multiplier that will multiply two matrices, A and B, storing the result in C. I have been told to write it in C, of which I have no experience. ...