Stefan-Endres / DWPM-Mixture-Model

Phase seperation calculation using the DWPM mixture rule.
7 stars 3 forks source link

Error in math in tgo.py #14

Closed alchemyst closed 8 years ago

alchemyst commented 8 years ago

In tgo.py:

    k_c = numpy.floor((-(ep - 1) + numpy.sqrt((ep - 1.0)--2 + 80.0 * ep))
                      / 2.0)

There is a --2 which I am pretty sure is not supposed to be there. Is it perhaps supposed to be **2 instead?

Stefan-Endres commented 8 years ago

Ah yes thank you, it should be **2 (this wouldn't have affected the results in nComp.py since a forced k = 5 option is used in the routines). Fixed with 9f3027a.