Assigning return value of a function to multidimensional arrays in VHDL fails with Vivado

Question

I have the following package which defines a 2D array type and some random function which returns an initialized 2D array.

-- 
library IEEE;
use