12

Working on a project using Tensorflow. However, I can't seem to reproduce my results.

I have tried setting the graph level seed, numpy random seed and even operation level seeds. However, it still not reproducible.

On searching Google, most people point to the reduce_sum function as the culprit as the reduce_sum function has a non-deterministic property on gpu even after setting the seeds. However, since I am working on a project for a paper, I need to reproduce the results. Is there any other efficient function that can work around this?

Another suggestion was to use CPU. However, I'm working on bug data and such CPU is not an option. How do people working on complex projects using Tensorflow work around this issue? Or it is acceptable to reviewers to load the saved model checkpoint file for result verification?

Martin Thoma
  • 124,992
  • 159
  • 614
  • 958
uchman21
  • 678
  • 3
  • 6
  • 22
  • 4
    I think it's because of sum operation for floats is not really accociative, i.e. `a+(b+c)` not always equals to `(a+b)+c`. So in case of any parallel computing when an order of operations is not guaranteed, you can get different results each time. It can begin with small differences, but over time they grow. – stop-cran Aug 24 '17 at 15:36
  • 2
    If you are not sure about what may or may not be acceptable for a reviewer you can ask more senior members of your department, or your supervisor if you're a student; [here is a related question in Academia SE](https://academia.stackexchange.com/questions/88576/reproducible-and-exact-calculations-vs-approximate-and-slightly-random-results). What I can tell is that TensorFlow has been used in many accepted academic publications. Probably it depends on the variability of the results, the size of the data and your particular field/problem, among other things. – jdehesa Aug 24 '17 at 15:47
  • @stop-cran I noticed that... I noticed the more I increase the dataset size, the wider the gap between the results. – uchman21 Aug 24 '17 at 16:18
  • Can you try in nightly version? A recent commit (d93a55b8) is supposed to make reduce_sum deterministic on GPU – Yaroslav Bulatov Aug 24 '17 at 17:37
  • @YaroslavBulatov, Never heard of that. Thanks will look at it to see how it works. – uchman21 Aug 24 '17 at 20:20
  • The nightly build didn't help... maybe just reduced the difference a bit. – uchman21 Aug 25 '17 at 11:36
  • Beside `tf.set_random_seed(1234)`, have you tried to set the number of threads to 1 (as the order in parallelism in `reduce_sum` is not guaranteed)? : `session_conf = tf.ConfigProto(intra_op_parallelism_threads=1, inter_op_parallelism_threads=1)` ; `sess = tf.Session(graph=tf.get_default_graph(), config=session_conf)` – TNg Nov 21 '18 at 16:28
  • @thanhtang, that would not be utilizing the computing power of GPU – uchman21 Jan 09 '19 at 12:22
  • @uchman21: Right, I know. I was just trying to see if anybody else could reproduce the result with those settings (because I couldn't). – TNg Jan 10 '19 at 02:31

1 Answers1

6

Cool, that you want to make your results reproducible! However, there are many things to note here:

I call a paper reproducible if one can obtain exactly the same numbers as found in the paper by executing exactly the same steps. This means if one had access to the same environment, the same software, hardware and data, one would be able to get the same results. In contrast, a paper is called replicatable if one can achieve the same results if one only follows the textual description in the paper. Hence replicability is harder to achieve, but also a more powerful indicator of the quality of the paper

You want to achieve that the training results on a bit-wise identical model. The holy grail would be to write your paper in a way that if people ONLY have the paper, they can still confirm your results.

Please also note that in many important papers results are practically impossible to reproduce:

  • Datasets are often not available: JFT-300M
  • Massive usage of computational power: For one of the AutoML/Architecture Search papers by Google I asked the author how many GPU-hours they spent on one of the experiments. At the time, if I wanted that many GPU-hours it would have costed me around 250,000 USD.

If that is a problem, depends very much on the context. As a comparison, think of CERN / LHC: It is impossible to have completely identical experiments. Only very few institutions on earth have the instruments to check the results. Still it is not a problem. So ask your advisor / people who have already published in that journal / conference.

Achieving Replicatability

This is super hard. I think the following is helpful:

  • Make sure the quality metrics you mention don't have too many digits
  • As the training likely depends on random initialization, you might also want to give rather an interval than a single number
  • Try minor variations
  • Re-implement things from scratch (maybe with another library?)
  • Ask colleagues to read your paper and then explain back to you what they think you did.

Getting Bit-Wise identical Model

It seems to me that you already do the important things:

  • Setting all seeds: numpy, tensorflow, random, ...
  • Making sure the Training-Test split is consistent
  • Making sure the training data is loaded in the same order

Please note that there might be factors out of your control:

  • Bitflips: B. Schroeder, E. Pinheiro, and W.-D. Weber, “Dram errors in the wild: a large-scale field study”
  • Inherent Hardware/Software reproducibility problems: Floating point multiplication is not associative and different cores on a GPU might finish computations at different times. Thus each single run could lead to different results. (I'd be happy if somebody could give an authorative reference here)
Martin Thoma
  • 124,992
  • 159
  • 614
  • 958