docs: def _build(self, inputs, prev_state, is_training=None, test_local_stats=True):
      """Connects the LSTM module into the graph.
  
      If this is not the first time the module has been connected to the graph,
      the Tensors provided as inputs and state must have the same final
      dimension, in order for the existing variables to be the correct size for
      their corresponding multiplications. The batch size may differ for each
      connection.
  
      Args:
        inputs: Tensor of size `[batch_size, input_size]`.
        prev_state: Tuple (prev_hidden, prev_cell), or if batch norm is enabled
          and `max_unique_stats > 1`, then (prev_hidden, prev_cell, time_step).
          Here, prev_hidden and prev_cell are tensors of size
          `[batch_size, hidden_size]`, and time_step is used to indicate the
          current RNN step.
        is_training: Boolean indicating whether we are in training mode (as
          opposed to testing mode), passed to the batch norm
          modules. Note to use this you must wrap the cell via the
          `with_batch_norm_control` function.
        test_local_stats: Boolean indicating whether to use local batch statistics
          in test mode. See the `BatchNorm` documentation for more on this.
  
      Returns:
        A tuple (output, next_state) where 'output' is a Tensor of size
        `[batch_size, hidden_size]` and 'next_state' is a tuple
        (next_hidden, next_cell) or (next_hidden, next_cell, time_step + 1),
        where next_hidden and next_cell have size `[batch_size, hidden_size]`.
  
      Raises:
        ValueError: If connecting the module into the graph any time after the
          first time, and the inferred size of the inputs does not match previous
          invocations.
      """
      if is_training is None:
        raise ValueError("Boolean is_training flag must be explicitly specified "
                         "when using batch normalization.")
  
      if self._max_unique_stats == 1:
        prev_hidden, prev_cell = prev_state
        time_step = None
      else:
        prev_hidden, prev_cell, time_step = prev_state
  
      # pylint: disable=invalid-unary-operand-type
      if self._hidden_clip_value is not None:
        prev_hidden = tf.clip_by_value(
            prev_hidden, -self._hidden_clip_value, self._hidden_clip_value)
      if self._cell_clip_value is not None:
        prev_cell = tf.clip_by_value(
            prev_cell, -self._cell_clip_value, self._cell_clip_value)
      # pylint: enable=invalid-unary-operand-type
  
      self._create_gate_variables(inputs.get_shape(), inputs.dtype)
      self._create_batch_norm_variables(inputs.dtype)
  
      # pylint false positive: calling module of same file;
      # pylint: disable=not-callable
  
      if self._use_batch_norm_h or self._use_batch_norm_x:
        gates_h = tf.matmul(prev_hidden, self._w_h)
        gates_x = tf.matmul(inputs, self._w_x)
        if self._use_batch_norm_h:
          gates_h = self._gamma_h * self._batch_norm_h(gates_h,
                                                       time_step,
                                                       is_training,
                                                       test_local_stats)
        if self._use_batch_norm_x:
          gates_x = self._gamma_x * self._batch_norm_x(gates_x,
                                                       time_step,
                                                       is_training,
                                                       test_local_stats)
        gates = gates_h + gates_x
      else:
        # Parameters of gates are concatenated into one multiply for efficiency.
        inputs_and_hidden = tf.concat([inputs, prev_hidden], 1)
        gates = tf.matmul(inputs_and_hidden, self._w_xh)
  
      gates += self._b
  
      # i = input_gate, j = next_input, f = forget_gate, o = output_gate
      i, j, f, o = tf.split(value=gates, num_or_size_splits=4, axis=1)
  
      if self._use_peepholes:  # diagonal connections
        self._create_peephole_variables(inputs.dtype)
        f += self._w_f_diag * prev_cell
        i += self._w_i_diag * prev_cell
  
      forget_mask = tf.sigmoid(f + self._forget_bias)
      next_cell = forget_mask * prev_cell + tf.sigmoid(i) * tf.tanh(j)
      cell_output = next_cell
      if self._use_batch_norm_c:
        cell_output = (self._beta_c
                       + self._gamma_c * self._batch_norm_c(cell_output,
                                                            time_step,
                                                            is_training,
                                                            test_local_stats))
      if self._use_peepholes:
        cell_output += self._w_o_diag * cell_output
      next_hidden = tf.tanh(cell_output) * tf.sigmoid(o)
  
      if self._max_unique_stats == 1:
        return next_hidden, (next_hidden, next_cell)
      else:
        return next_hidden, (next_hidden, next_cell, time_step + 1)