def METHOD_NAME(self, solver):
    mat = get_base_matrix(use_tril=False)
    zero_mat = mat.copy()
    zero_mat.data.fill(0)
    stat = solver.do_symbolic_factorization(zero_mat)
    self.assertEqual(stat.status, LinearSolverStatus.successful)
    stat = solver.do_numeric_factorization(mat)
    self.assertEqual(stat.status, LinearSolverStatus.successful)
    x_true = np.array([1, 2, 3], dtype=np.double)
    rhs = mat * x_true
    x, res = solver.do_back_solve(rhs)
    self.assertTrue(np.allclose(x, x_true))
    x_true = np.array([4, 2, 3], dtype=np.double)
    rhs = mat * x_true
    x, res = solver.do_back_solve(rhs)
    self.assertTrue(np.allclose(x, x_true))