def METHOD_NAME(self):
    """Test when using a user defined initial point."""
    initial_point = np.asarray([-7.35250290e-05, -9.73079292e-02, -5.43346282e-05])
    estimator = Estimator()
    optimizer = COBYLA(maxiter=1)
    initial_state = VSCF(self.vibrational_problem.num_modals, self.qubit_mapper)
    ansatz = UVCCSD(
        self.vibrational_problem.num_modals, self.qubit_mapper, initial_state=initial_state
    )
    solver = VQE(estimator, ansatz, optimizer, initial_point=initial_point)
    gsc = GroundStateEigensolver(self.qubit_mapper, solver)
    esc = QEOM(gsc, estimator, "sd")
    results = esc.solve(self.vibrational_problem)
    np.testing.assert_array_almost_equal(
        results.raw_result.ground_state_raw_result.optimal_point, initial_point
    )