I forgot to mention in the previous email that we could add a "Get_vnl_vector" method to the Array class which would return a vnl_vector_ref of the appropriate type. This would let us use all the matrix math libraries from vnl with the Array and its derived Vector/Point classes without the overhead of the extra storage needed by vnl_vector. -Brad