Description 
Since quaternions are not regular vectors, but rather representations of orientation, an average quaternion cannot just be obtained by taking a weighted mean. This function implements the work done by F. Landis Merkley to calculate the average quaternion. The algorithm is explained by F. Landis Markley at:
http://www.acsu.buffalo.edu/~johnc/ave_quat07.pdf
For this particular implementation, I would also like to reference Mandar Harshe:
http://wwwsop.inria.fr/members/Mandar.Harshe/kneejoint/html/index.html
While being basic and straightforward, this algorithm is compared with rotqrmean from VoiceBox and found to produce quite similar results, yet it is more elegant, much simpler to implement and follow. (Though, there might be difference in signs)
