code: def _hpsd_solve_fwd(matrix, rhs, cholesky_matrix):
      del matrix
      chol_linop = tf.linalg.LinearOperatorLowerTriangular(
          cholesky_matrix
      )
      solve_rhs = chol_linop.solve(chol_linop.solve(rhs), adjoint=True)
      return solve_rhs, (cholesky_matrix, rhs, solve_rhs)
  
imports: import numpy
  import tensorflow.compat.v1
  import tensorflow.compat.v2
  from tensorflow_probability.python.internal import assert_util
  from tensorflow_probability.python.internal import custom_gradient
  from tensorflow_probability.python.internal import distribution_util
  from tensorflow_probability.python.internal import dtype_util
  from tensorflow_probability.python.internal import prefer_static
  from tensorflow_probability.python.internal import tensorshape_util
  from tensorflow_probability.python.math import generic
label_desc: Use image node_id 25 for calling a global function with example usage: _hpsd_solve_fwd(matrix, rhs, cholesky_matrix) and returns: solve_rhs,  cholesky_matrix,  rhs,  solve_rhs