Diego,

There aren't many linear solvers that are bit-consistent, where the answer is the same no matter how many cores or processes you use. Intel's version of Pardiso is bit-consistent and I think MUMPS 5.0 might be, but that's all. You should assume your answer will not be exactly the same as you change the number of cores or processes, although you should reach the same overall error tolerance in approximately the same number of iterations.

Damien

On 2015-10-28 3:51 PM, Diego Avesani wrote:
dear Andreas, dear all,
The code is quite long. It is a conjugate gradient algorithm to solve a complex system.

I have noticed that when a do cycle is small, let's say
do i=1,3

enddo

the results are identical. If the cycle is big, let's say do i=1,20, the results are different and the difference increase with the number of iterations.

What do you think?



Diego


On 28 October 2015 at 22:32, Andreas Schäfer <gent...@gmx.de <mailto:gent...@gmx.de>> wrote:

    On 22:03 Wed 28 Oct     , Diego Avesani wrote:
    > When I use a single CPU a get a results, when I use 4 CPU I get
    another
    > one. I do not think that very is a bug.

    Sounds like a bug to me, most likely in your code.

    > Do you think that these small differences are normal?

    It depends on what small means. Floating point operations in a
    computer are generally not commutative, so parallelization may in deed
    lead to different results.

    > Is there any way to get the same results? is some align problem?

    Impossible to say without knowing your code.

    Cheers
    -Andreas


    --
    ==========================================================
    Andreas Schäfer
    HPC and Grid Computing
    Department of Computer Science 3
    Friedrich-Alexander-Universität Erlangen-Nürnberg, Germany
    +49 9131 85-27910 <tel:%2B49%209131%2085-27910>
    PGP/GPG key via keyserver
    http://www.libgeodecomp.org
    ==========================================================

    (\___/)
    (+'.'+)
    (")_(")
    This is Bunny. Copy and paste Bunny into your
    signature to help him gain world domination!

    _______________________________________________
    users mailing list
    us...@open-mpi.org <mailto:us...@open-mpi.org>
    Subscription: http://www.open-mpi.org/mailman/listinfo.cgi/users
    Link to this post:
    http://www.open-mpi.org/community/lists/users/2015/10/27933.php




_______________________________________________
users mailing list
us...@open-mpi.org
Subscription: http://www.open-mpi.org/mailman/listinfo.cgi/users
Link to this post: 
http://www.open-mpi.org/community/lists/users/2015/10/27934.php

Reply via email to