Hermitean Clifford analysis encompasses the theory of holomorphic functions of several complex variables; the theory relies on the introduction of a complex structure in Euclidean space and is developed via an order two circulant matrix approach. Quaternionic Hermitean Clifford analysis involves several complex structures, refining the system. This project aims at developing the corresponding function theory by a higher order circulant matrix approach.