code: def psd_solve_with_chol(node):
      if node.op == solve:
          A, b = node.inputs  # result is solution Ax=b
          if is_psd(A):
              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 the L matrix during the
              #     stabilization
              Li_b = Solve("lower_triangular")(L, b)
              x = Solve("upper_triangular")(L.T, Li_b)
              return [x]
  
imports: from __future__ import absolute_import, print_function, division
  import logging
  from six import iteritems, integer_types
  from six.moves import xrange
  from theano.gof import Op, Apply
  from theano.tensor import DimShuffle, Dot
  from theano.tensor.blas import Dot22
  from theano import tensor
  import theano.tensor
  from theano.tensor.opt import register_stabilize, register_specialize, register_canonicalize
  from theano.gof import local_optimizer
  from theano.gof.opt import Optimizer
  from theano.tensor.nlinalg import MatrixInverse, matrix_inverse, extract_diag, trace, det
  from theano.tensor.slinalg import Cholesky, cholesky, Solve, solve, imported_scipy
label_desc: Use image node_id 12 for calling a global function with example usage: psd_solve_with_chol(node) and returns: list