Fix initial nonzero guess in TransportDG when solver type is 'preonly' (e.g. when direct solver is used).
Implemented P3 cubic finite elements.
Update in DG method (diffusive fluxes on dirichlet boundary).
Corrected mass balance.
Updated documentation and tests.