Hessian of vector-valued functions

Could someone explain how the Hessian of a vector-valued function f = (f1,f2,...,fn) is calculated or point me towards a reference where the information is available? I understand the array of second partial derivatives is not a matrix, but a tensor of rank 3....
.