old_contents: # Standard imports
  import math
  import logging
  import numpy as np
  import emission.core.common as ec
  import emission.analysis.section_features as sf
  
  def calDistance(point1, point2):
      return ec.calDistance([point1.longitude, point1.latitude], [point2.longitude, point2.latitude])
  
  def calHeading(point1, point2):
      return sf.calHeading([point1.longitude, point1.latitude],
                      [point2.longitude, point2.latitude])
  
  def calHC(point1, point2, point3):
      return sf.calHC([point1.longitude, point1.latitude],
                      [point2.longitude, point2.latitude],
                      [point3.longitude, point3.latitude])
  
  def calSpeed(point1, point2):
      distanceDelta = calDistance(point1, point2)
      timeDelta = point2.mTime - point1.mTime
      # print "Distance delta = %s and time delta = %s" % (distanceDelta, timeDelta)
      # assert(timeDelta != 0)
      if (timeDelta == 0):
          logging.debug("timeDelta = 0, distanceDelta = %s, returning speed = 0")
          assert(distanceDelta < 0.01)
          return 0
  
      # TODO: Once we perform the conversions from ms to secs as part of the
      # usercache -> timeseries switch, we need to remove this division by 1000
      return distanceDelta/(float(timeDelta)/1000)
  
new_contents: # Standard imports
  import math
  import logging
  import numpy as np
  import emission.core.common as ec
  import emission.analysis.section_features as sf
  
  def calDistance(point1, point2):
      return ec.calDistance([point1.longitude, point1.latitude], [point2.longitude, point2.latitude])
  
  def calHeading(point1, point2):
      return sf.calHeading([point1.longitude, point1.latitude],
                      [point2.longitude, point2.latitude])
  
  def calHC(point1, point2, point3):
      return sf.calHC([point1.longitude, point1.latitude],
                      [point2.longitude, point2.latitude],
                      [point3.longitude, point3.latitude])
  
  def calSpeed(point1, point2):
      distanceDelta = calDistance(point1, point2)
      timeDelta = point2.ts - point1.ts
      # print "Distance delta = %s and time delta = %s" % (distanceDelta, timeDelta)
      # assert(timeDelta != 0)
      if (timeDelta == 0):
          logging.debug("timeDelta = 0, distanceDelta = %s, returning speed = 0")
          assert(distanceDelta < 0.01)
          return 0
  
      return distanceDelta/timeDelta
  
subject: Change the feature calculation to match the new unified format