Triangular Arrays in julia?

You could just wrap a float in a new type and define the needed operations. This should dispatch to the generic versions of the linear algebra library.