code: def find_full_domain(self, obj):
          r
          if 'Subdomain' not in obj._mro():
              # Network, Phase, Alg
              return obj
          else:
              if obj._isa() == 'geometry':
                  # Geom
                  return self.network
              else:
                  # Phys
                  return self.find_phase(obj)
docs: r"""
          Find the full domain object associated with a given object.
          For geometry the network is found, for physics the phase is found and
          for all other objects which are defined for for the full domain,
          themselves are found.
  
          Parameters
          ----------
          obj : OpenPNM Object
              Can be any object
  
          Returns
          -------
          An OpenPNM object