Matrix product state algorithms for Gaussian fermionic states