code: def _hpsd_quadratic_form_solve_jvp(primals, tangents):
      """JVP for hpsd_quadratic_form_solve."""
      matrix, rhs, cholesky_matrix = primals
      gmatrix, grhs, gcholesky_matrix = tangents
      del gcholesky_matrix
      output, (
          cholesky_matrix,
          _,
          solve_rhs,
      ) = _hpsd_quadratic_form_solve_fwd(matrix, rhs, cholesky_matrix)
      chol_linop = tf.linalg.LinearOperatorLowerTriangular(
          cholesky_matrix
      )
      full_solve = chol_linop.solve(solve_rhs, adjoint=True)
      gmatrix_linop = tf.linalg.LinearOperatorFullMatrix(gmatrix)
      jvp = tf.linalg.matmul(
          full_solve,
          2.0 * grhs - gmatrix_linop.matmul(full_solve),
          transpose_a=True,
      )
      return output, jvp
  
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 33 for calling a global function with example usage: _hpsd_quadratic_form_solve_jvp(primals, tangents) and returns: output,  jvp