Fill the domain with colors of the current palette corresponding to the current component of the solution. Coloring is done by linear interpolation between the grid points.