Implement partial trace of fermionic operators

There is an old implementation here. We need to update it so that it fits we the package.