import numpy as np
import matplotlib
import matplotlib.pyplot as plt

#S3  1276 http://browser.geekbench.com/geekbench3/8570680
#S4: 1772 https://browser.geekbench.com/android_devices/108
#S5: 2414 https://browser.geekbench.com/android_devices/163
#S6: 4049 https://browser.geekbench.com/android_devices/209
#S7: 5348 https://browser.geekbench.com/android_devices/220
#S8: 6428 https://browser.geekbench.com/android_devices/376
# S9: 8945 https://www.tek.no/artikler/test-samsung-galaxy-s9/431900/2
samsung_scores = np.array([1276, 1772, 2414, 4049, 5348, 6428, 8945])

# Release dates of Galaxy S (Rounded down to month)
# (source: google "Samsung Galaxy S<version> release date")
samsung_versions = np.array([2012+5./12, 2013+4/12., 2014+4/12., 2015+4/12.,
                             2016+3./12, 2017+4/12., 2018+3/12.])

p_samsung =  plt.plot(samsung_versions, samsung_scores, '*',
         label='Samsung Galaxy S <Year>')

# Now we fit the data, except the latest, and see how good a prediction it is
samsung_versions_shifted = samsung_versions - samsung_versions[0]*np.ones_like(samsung_versions)
samsung_poly_coeff = np.polyfit(samsung_versions_shifted[:-2],
    np.log(samsung_scores[:-2]), 1)


plt.plot(samsung_versions,
         np.exp(samsung_poly_coeff[1])*np.exp(samsung_versions_shifted*samsung_poly_coeff[0]),
         '--', color=p_samsung[0].get_color(),
         label='$%.2f (%.2f)^{\\mathrm{Year} - %.2f}$'
         % (np.exp(samsung_poly_coeff[1]),
          np.exp(samsung_poly_coeff[0]),
          samsung_versions[0]))
# iphone 4s 491 https://browser.geekbench.com/ios_devices/6
# iphone 5 1207 https://browser.geekbench.com/ios_devices/20
# iphone 5s 2144 https://browser.geekbench.com/ios_devices/28
# iphone 6 2307 https://browser.geekbench.com/ios_devices/33
# iphone 6s 3813 https://browser.geekbench.com/ios_devices/38
# iphone 7 5694 https://browser.geekbench.com/ios_devices/44
# iphone 8 10129 https://browser.geekbench.com/ios_devices/50

iphone_scores = np.array([491, 1207, 2144, 2307, 3813, 5694, 10129])

# Release dates of iPhone (Rounded down to month)
# (source: google "Iphone <version> release date")
iphone_versions = np.array(
    [2011. + 10/12., 2012 + 9/12., 2013+ 9/12., 2014+9/12.,
                   2015+9/12., 2016+9/12., 2017+9/12.])


p_apple = plt.plot(iphone_versions, iphone_scores,
                   'o', label='Apple iPhone <Year>')

iphone_versions_shifted = iphone_versions - iphone_versions[0]*np.ones_like(iphone_versions)
apple_poly_coeff = np.polyfit(iphone_versions_shifted[:-2],
    np.log(iphone_scores[:-2]), 1)


plt.plot(iphone_versions,
         np.exp(apple_poly_coeff[1])*np.exp(iphone_versions_shifted*apple_poly_coeff[0]),
         '--', color=p_apple[0].get_color(),
         label='$%.2f (%.2f)^{\\mathrm{Year} - %.2f}$'
         % (np.exp(apple_poly_coeff[1]),
          np.exp(apple_poly_coeff[0]),
          iphone_versions[0]))


plt.xlabel('Year')
plt.ylabel('3D Mark Ice Storm')
plt.legend()
plt.grid("on")
plt.show()