code: def psd_solve_with_chol(fgraph, node):
      """
      This utilizes a boolean `psd` tag on matrices.
      """
      if (
          isinstance(node.op.core_op, Solve)
          and node.op.core_op.b_ndim == 2
      ):
          A, b = node.inputs  # result is solution Ax=b
          if getattr(A.tag, "psd", None) is True:
              L = cholesky(A)
              # N.B. this can be further reduced to a yet-unwritten cho_solve Op
              #     __if__ no other Op makes use of the L matrix during the
              #     stabilization
              Li_b = solve_triangular(L, b, lower=True, b_ndim=2)
              x = solve_triangular(_T(L), Li_b, lower=False, b_ndim=2)
              return [x]
  
imports: import logging
  from typing import cast
  from pytensor.graph.rewriting.basic import copy_stack_trace, node_rewriter
  from pytensor.tensor.basic import TensorVariable, diagonal, swapaxes
  from pytensor.tensor.blas import Dot22
  from pytensor.tensor.blockwise import Blockwise
  from pytensor.tensor.elemwise import DimShuffle
  from pytensor.tensor.math import Dot, Prod, _matrix_matrix_matmul, log, prod
  from pytensor.tensor.nlinalg import MatrixInverse, det
  from pytensor.tensor.rewriting.basic import register_canonicalize, register_specialize, register_stabilize
  from pytensor.tensor.slinalg import Cholesky, Solve, SolveBase, cholesky, solve, solve_triangular
label_desc: Use image node_id 8 for calling a global function with example usage: psd_solve_with_chol(fgraph, node) and returns: list