Search:

Type: Posts; User: collymitch

Search: Search took 0.00 seconds.

  1. Replies
    4
    Views
    3,987

    It was working...

    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...
  2. Replies
    4
    Views
    3,987

    Lookin good now thanks.

    Lookin good now thanks.
  3. Replies
    4
    Views
    3,987

    comparing two matrices or double values

    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 )
    {...
  4. OpenMP parallel threads for matrix multiplier

    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...
  5. Like I said they can be linked fine. It's just...

    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.

    ...
  6. External Functions and Compiler linking issues

    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...
  7. compare double value to 0 - unreliable comparisons warning

    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...
  8. Replies
    13
    Views
    2,579

    Excellent. That seems to be working now. I could...

    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...
  9. Replies
    13
    Views
    2,579

    Mmmm.... I modified it to take in the a[m][n]...

    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...
  10. Replies
    13
    Views
    2,579

    Yes that's right - i'm comparing my result with...

    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...
  11. Replies
    13
    Views
    2,579

    Ok Dave, Thanks, that seems to work fine. I...

    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...
  12. Replies
    13
    Views
    2,579

    Ok - i got a bit confused by that so had to ask...

    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...
  13. Replies
    13
    Views
    2,579

    matrix mult - double ptr in col major

    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...
  14. Hey I've got that set up and working fine now by...

    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.
  15. Ok.... Yeah the 'void main' bit was a typing...

    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....
  16. Hey thanks, Ok that helped a bit. I read that...

    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...
  17. Matrix multipliers - using a double pointer instead of an array

    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. ...
Results 1 to 17 of 17