def as_late_as_possible(circuit: QuantumCircuit, schedule_config: ScheduleConfig) -> Schedule: """ Return the pulse Schedule which implements the input circuit using an "as late as possible" (alap) scheduling policy. Circuit instructions are first each mapped to equivalent pulse Schedules according to the command definition given by the schedule_config. Then, this circuit instruction-equivalent Schedule is appended at the latest time that it can be without allowing unnecessary time between instructions or allowing instructions with common qubits to overlap. This method should improves the outcome fidelity over ASAP scheduling, because we may maximize the time that the qubit remains in the ground state. Args: circuit: The quantum circuit to translate. schedule_config: Backend specific parameters used for building the Schedule. Returns: A schedule corresponding to the input ``circuit`` with pulses occurring as late as possible. """ qubit_time_available = defaultdict(int) def update_times(inst_qubits: List[int], time: int = 0) -> None: """Update the time tracker for all inst_qubits to the given time.""" for q in inst_qubits: qubit_time_available[q] = time rev_stop_times = [] circ_pulse_defs = lower_gates(circuit, schedule_config) for circ_pulse_def in reversed(circ_pulse_defs): start_time = max(qubit_time_available[q] for q in circ_pulse_def.qubits) stop_time = start_time if not isinstance(circ_pulse_def.schedule, Barrier): stop_time += circ_pulse_def.schedule.duration rev_stop_times.append(stop_time) update_times(circ_pulse_def.qubits, stop_time) last_stop = max( t for t in qubit_time_available.values()) if qubit_time_available else 0 start_times = [last_stop - t for t in reversed(rev_stop_times)] timed_schedules = [(time, cpd.schedule) for time, cpd in zip(start_times, circ_pulse_defs) if not isinstance(cpd.schedule, Barrier)] schedule = Schedule.initialize_from(circuit) for time, inst in timed_schedules: schedule.insert(time, inst, inplace=True) return schedule
def as_soon_as_possible(circuit: QuantumCircuit, schedule_config: ScheduleConfig) -> Schedule: """ Return the pulse Schedule which implements the input circuit using an "as soon as possible" (asap) scheduling policy. Circuit instructions are first each mapped to equivalent pulse Schedules according to the command definition given by the schedule_config. Then, this circuit instruction-equivalent Schedule is appended at the earliest time at which all qubits involved in the instruction are available. Args: circuit: The quantum circuit to translate. schedule_config: Backend specific parameters used for building the Schedule. Returns: A schedule corresponding to the input ``circuit`` with pulses occurring as early as possible. """ qubit_time_available = defaultdict(int) def update_times(inst_qubits: List[int], time: int = 0) -> None: """Update the time tracker for all inst_qubits to the given time.""" for q in inst_qubits: qubit_time_available[q] = time start_times = [] circ_pulse_defs = lower_gates(circuit, schedule_config) for circ_pulse_def in circ_pulse_defs: start_time = max(qubit_time_available[q] for q in circ_pulse_def.qubits) stop_time = start_time if not isinstance(circ_pulse_def.schedule, Barrier): stop_time += circ_pulse_def.schedule.duration start_times.append(start_time) update_times(circ_pulse_def.qubits, stop_time) timed_schedules = [(time, cpd.schedule) for time, cpd in zip(start_times, circ_pulse_defs) if not isinstance(cpd.schedule, Barrier)] schedule = Schedule.initialize_from(circuit) for time, inst in timed_schedules: schedule.insert(time, inst, inplace=True) return schedule
def sequence(scheduled_circuit: QuantumCircuit, schedule_config: ScheduleConfig) -> Schedule: """ Return the pulse Schedule which implements the input scheduled circuit. Assume all measurements are done at once at the last of the circuit. Schedules according to the command definition given by the schedule_config. Args: scheduled_circuit: The scheduled quantum circuit to translate. schedule_config: Backend specific parameters used for building the Schedule. Returns: A schedule corresponding to the input ``circuit``. Raises: QiskitError: If invalid scheduled circuit is supplied. """ circ_pulse_defs = lower_gates(scheduled_circuit, schedule_config) # find the measurement start time (assume measurement once) def _meas_start_time(): _qubit_time_available = defaultdict(int) for inst, qubits, _ in scheduled_circuit.data: if isinstance(inst, Measure): return _qubit_time_available[qubits[0]] for q in qubits: _qubit_time_available[q] += inst.duration return None meas_time = _meas_start_time() # restore start times qubit_time_available = {} start_times = [] out_circ_pulse_defs = [] for circ_pulse_def in circ_pulse_defs: active_qubits = [ q for q in circ_pulse_def.qubits if q in qubit_time_available ] start_time = max((qubit_time_available[q] for q in active_qubits), default=0) for q in active_qubits: if qubit_time_available[q] != start_time: # print(q, ":", qubit_time_available[q], "!=", start_time) raise QiskitError("Invalid scheduled circuit.") stop_time = start_time if not isinstance(circ_pulse_def.schedule, Barrier): stop_time += circ_pulse_def.schedule.duration delay_overlaps_meas = False for q in circ_pulse_def.qubits: qubit_time_available[q] = stop_time if (meas_time is not None and circ_pulse_def.schedule.name == "delay" and stop_time > meas_time): qubit_time_available[q] = meas_time delay_overlaps_meas = True # skip to delays overlapping measures and barriers if not delay_overlaps_meas and not isinstance(circ_pulse_def.schedule, Barrier): start_times.append(start_time) out_circ_pulse_defs.append(circ_pulse_def) timed_schedules = [(time, cpd.schedule) for time, cpd in zip(start_times, out_circ_pulse_defs)] sched = Schedule(*timed_schedules, name=scheduled_circuit.name) return pad(sched)