Matrices in Julia, unlike some other languages, aren't vectors of vectors - they are first class citizens of the multidimensional array system in Julia. Thus,
A = Array(Float64,1,1) creates a 1x1 matrix that holds Float64s. In most cases, it's even better to use A = zeros(Float64, 1, 1) to ensure that you don't have any garbage data (Array allocates, but doesn't wipe the memory). Whether that is possible in your case depends on the type parameter data. // T
