highest theshold rhodopsin distance matrix
wrapped test code into test()