Variational formulations of time-dependent PDEs in space and time yield (d + 1)-dimensional problems to be solved numerically. This increases the number unknowns as well storage amount. On other hand, this approach enables adaptivity model reduction w.r.t. both type variables. In paper, we show that matrix oriented techniques can significantly reduce computational timings for solving arising li...