  In this work we develop a methodology to approximate the covariance matrix
associated to the simulation of water diffusion inside the brain tissue. The
computation is based on an implementation of the Discontinuous Galerkin method
of the diffusion equation, in accord with the physical phenomenon. The
implementation in in parallel using GPUs in the CUDA language. Numerical
results are presented in 2D problems.
