code: def strip_hydrogens(
      coords: np.ndarray, mol: Union[RDKitMol, MolecularFragment]
  ) -> Tuple[np.ndarray, MolecularFragment]:
      """Strip the hydrogens from input molecule
  
      Parameters
      ----------
      coords: np.ndarray
          The coords must be of shape (N, 3) and correspond to coordinates of mol.
      mol: rdkit.Chem.rdchem.Mol or MolecularFragment
          The molecule to strip
  
      Returns
      -------
      Tuple[np.ndarray, MolecularFragment]
          A tuple of `(coords, mol_frag)` where `coords` is a numpy array of
          coordinates with hydrogen coordinates. `mol_frag` is a `MolecularFragment`.
  
      Notes
      -----
      This function requires RDKit to be installed.
      """
      mol_atoms = mol.GetAtoms()
      atomic_numbers = [atom.GetAtomicNum() for atom in mol_atoms]
      atom_indices_to_keep = [
          ind
          for (ind, atomic_number) in enumerate(atomic_numbers)
          if (atomic_number != 1)
      ]
      return get_mol_subset(coords, mol, atom_indices_to_keep)
  
imports: import itertools
  import numpy
  from typing import List, Optional, Sequence, Set, Tuple, Union
  import logging
  from deepchem.utils.typing import RDKitAtom, RDKitMol
  from deepchem.utils.geometry_utils import compute_pairwise_distances
label_desc: Use image node_id 4 for calling a global function with example usage: strip_hydrogens(coords, mol) and returns: get_mol_subset