Search:

Type: Posts; User: KBriggs

Page 1 of 20 1 2 3 4

Search: Search took 0.01 seconds.

  1. Replies
    30
    Views
    4,246

    As you can tell, the documentation is purely...

    As you can tell, the documentation is purely theoretical at this point ^_^. But I added a README with an acknowledgement :).

    Your debugging suggestion is a good one, I will implement something...
  2. Replies
    30
    Views
    4,246

    The other compiler flags did not really do...

    The other compiler flags did not really do anything either, but using a more liberal definition of zero seems to make the difference.

    I want to thank you for taking the time to help be understand...
  3. Replies
    30
    Views
    4,246

    So rerunning with the previous settings resulted...

    So rerunning with the previous settings resulted in -ffloat-store not doing anything after all. What changed before was that I changed the value of EPS (defined in utils.h) to 1e-10 from 1e-30. From...
  4. Replies
    30
    Views
    4,246

    I thought that fixed it but later realized I had...

    I thought that fixed it but later realized I had changed an important setting so the test was not valid. I am rerunning all my comparisons again, I will update shortly.
  5. Replies
    30
    Views
    4,246

    Further investigation found a whole bunch of...

    Further investigation found a whole bunch of unsafe floating point operations in lmmin_int64.c. Stuff like


    for (j = 0; j < n; j++) {
    temp3 = fjac[j*m+j];
    if...
  6. Replies
    30
    Views
    4,246

    nice find! I also found a whole bunch of...

    nice find!

    I also found a whole bunch of unsafe floating point comparisons (== and !=) in that file, which I fixed in a commit a couple of days ago. I'm testing at my end after removing the ! you...
  7. Replies
    30
    Views
    4,246

    I think the fact that I get different results...

    I think the fact that I get different results with the serial program using gcc vs mingw64 is probably telling. I am starting to suspect that the differences with different thread counts are a red...
  8. Replies
    30
    Views
    4,246

    OK, I will test that as well. Out of...

    OK, I will test that as well.

    Out of curiosity, where is the problem with the function as implemented before? (num>EPS) and (num<-EPS) are both just either 1 or 0, so I'm not sure where the...
  9. Replies
    30
    Views
    4,246

    what do you mean by long hand? Something like: ...

    what do you mean by long hand? Something like:


    int signum(double num)
    {
    if (num > EPS)
    {
    return 1;
    }
    else if (num < -EPS)
  10. Replies
    30
    Views
    4,246

    OK, an update on results with a larger data set....

    OK, an update on results with a larger data set. This makes less and less sense.

    on Windows 7, using x86_64-w64-mingw32-gcc (mingw64) to compile without any optimizations, my results DO depend on...
  11. Replies
    30
    Views
    4,246

    Yes, that's what I am seeing as well. I am...

    Yes, that's what I am seeing as well. I am rerunning the threaded comparisons (on a much bigger data set, to amplify errors if they are present) with the most recent code to see if the number of...
  12. Replies
    30
    Views
    4,246

    So, basically, out of my control. Good to know. ...

    So, basically, out of my control. Good to know.

    Thanks for your help fixing the memory issues. There's still the main bug to hunt down, but it's in a much better place to be able to diagnose than...
  13. Replies
    30
    Views
    4,246

    A few more memory leaks, this time picked up...

    A few more memory leaks, this time picked up using a 4-threaded run. I've been unable to figure out where the problems actually are, since the lines in main.c which valgrind indicates are problematic...
  14. Replies
    30
    Views
    4,246

    I have fixed and pushed all of those memory...

    I have fixed and pushed all of those memory leaks, thanks. That last one is odd - if you look at the master branch in that spot you will see that cpos is free'd correctly. Must have been some...
  15. Replies
    30
    Views
    4,246

    Ah. That does seem like it might be overkill at...

    Ah. That does seem like it might be overkill at this point. I'll look into whether or not that problems applies to later versions of gcc - they mention 4.2 and 4.3 but I think I am using 5.x at the...
  16. Replies
    30
    Views
    4,246

    My valgrind run on the serial version revealed...

    My valgrind run on the serial version revealed one memory leak, which I fixed (and certainly doesn't matter to the outcome). I am running it again to make sure I fixed it. My compilation and valgrind...
  17. Replies
    30
    Views
    4,246

    Apologies for double post, but I think it is...

    Apologies for double post, but I think it is warranted here.

    I tried to find a small data set that would demonstrate the problem, but it seems to be a cumulative issue. Using a long data set...
  18. Replies
    30
    Views
    4,246

    Thank you for taking the time to put that...

    Thank you for taking the time to put that together!

    I can certainly provide test data and config files if you would like (I honestly didn't expect anyone to try to compile and run it, I was just...
  19. Replies
    30
    Views
    4,246

    debugging multi-threaded code with OpenMP

    I have a piece of signal processing code that I've been working on for the past 2 years or so whenever I have some downtime. It works pretty well, and I've been slowly optimizing it, more as a way to...
  20. Replies
    12
    Views
    783

    You could do the same thing with ...

    You could do the same thing with

    scanf("%5d",x)
  21. Replies
    7
    Views
    776

    You were correct. Malloc in the head node case...

    You were correct. Malloc in the head node case was causing the reference to the actual head node to be lost. Thanks!

    My error was assuming that if I declare a pointer and then malloc it, that...
  22. Replies
    7
    Views
    776

    I don't think I understand - the first branch of...

    I don't think I understand - the first branch of the if happens if there is already a node in existence, and it adds a new one and returns a pointer to the new one. The head node should be stored...
  23. Replies
    7
    Views
    776

    Yes, the format specifiers are correct: PRIu64 is...

    Yes, the format specifiers are correct: PRIu64 is a macro for printing uint64_t types. This code compiles without error or warnings with -Wall and -Wextra flags enabled so I think it would complain...
  24. Replies
    7
    Views
    776

    Help understanding linked list behavior

    Hi all,

    I have the following two structs which form linked lists:


    struct Cusumlevel
    {
    double current;
    uint64_t length;
    struct Cusumlevel *next;
  25. Interesting, I had misunderstood something basic...

    Interesting, I had misunderstood something basic about pointer assignment apparently. Thanks for the help!
Results 1 to 25 of 487
Page 1 of 20 1 2 3 4