PETSc FDTD implementation for the solution of Maxwell's electrodynamic equations in anistropic, dispersive and time-varying materials