
Linear algebraic expressions are the essence of many computationally
intensive problems, including scientific simulations and machine learning
applications. However, translating highlevel formulations of these expressions
to efficient machinelevel representations is far from trivial: developers
should be assisted by automatic optimization tools so that they can focus their
attention on highlevel problems, rather than lowlevel details. The
tractability of these optimizations is highly dependent on the choice of the
primitive constructs in terms of which the computations are to be expressed. In
this work we propose to describe operations on multidimensional arrays using a
selection of higherorder functions, inspired by functional programming, and we
present rewrite rules for these such that they can be automatically optimized
for modern hierarchical and heterogeneous architectures. Using this formalism
we systematically construct and analyse different subdivisions and permutations
of the dense matrix multiplication problem.

In this paper the emergence of the Chiral Magnetic Effect (CME) and the
related anomalous current is investigated using the real time
DiracHeisenbergWigner formalism. This method is widely used for describing
strong field physics and QED vacuum tunneling phenomena as well as pair
production in heavyion collisions. We extend earlier investigations of the CME
in constant flux tube configuration by considering time dependent fields. In
this model we can follow the formation of axial charge separation, formation of
axial current and then the emergence of the anomalous electric current.
Qualitative results are shown for special field configurations that help to
interpret the predictions of CME related effects in heavyion collisions at the
RHIC Beam Energy Scan program.

We investigate the process of Abelian pair production in the presence of
strong inhomogeneous and timedependent external electric fields. The spatial
dependence of the external field is motivated by a nonAbelian color flux tube
in heavyion collisions. We show that the inhomogeneity significantly increase
the particle yield compared to that in the commonly used models with a constant
and homogeneous field. Moreover our results indicate that in contrast to the
latter, most of the particles are produced at the interface of the field
profile in accordance with Heisenberg's prediction.

In ultrarelativistic heavy ion collisions the baryonmeson ratio is not fully
understood in the intermediatetohigh transverse momentum region. Although
quark coalescence methods combined with jet fragmentation yield results close
to the experimental data, some tendencies are not reproduced and require
further investigation. We propose a new channel, namely extra quarkantiquark
and diquarkantidiquark pair production from coherent gluon field created in
the early stage of heavy ion collisions. This process is described
quantitatively by the timedependent Schwingermechanism. The extra
diquarkquark coalescence yields are able to explain the measured anomalous
baryon production.

In this article we review a theoretical framework for pair production from
strong external electromagnetic fields. We propose a numerical method to solve
the resulting equations of motion and present results for both cases of
spatially homogeneous and inhomogeneous electric fields.