# Copyright 2021 Huawei Technologies Co., Ltd
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
# ============================================================================
"""Circuit for time evolution."""
from ..circuit.utils import decompose_single_term_time_evolution
from ..parameterresolver import ParameterResolver
from .qubit_operator import QubitOperator
[文档]class TimeEvolution: # pylint: disable=too-few-public-methods
r"""
The time evolution operator that can generate a corresponded circuit.
The time evolution operator will do the following evolution:
.. math::
\left|\varphi(t)\right>=e^{-iHt}\left|\varphi(0)\right>
Note:
The hamiltonian should be a parameterized or non parameterized
QubitOperator. If the QubitOperator has multiple terms, the
first order trotter decomposition will be used.
Args:
ops (QubitOperator): The qubit operator hamiltonian, could be parameterized or
non parameterized.
time (Union[numbers.Number, dict, ParameterResolver]): The evolution time,
could be a number or a parameter resolver. If ``None``, the time will be set to 1.
Default: ``None``.
Examples:
>>> from mindquantum.core.operators import TimeEvolution, QubitOperator
>>> q1 = QubitOperator('Z0 Y1', 'a')
>>> q2 = QubitOperator('X0 Z1', 'b')
>>> ops1 = q1 + q2
>>> ops2 = q2 + q1
>>> TimeEvolution(ops1).circuit
┏━━━┓ ┏━━━┓
q0: ────────────────■─────────────────■───┨ H ┠──────────■─────────────────■───┨ H ┠───
┃ ┃ ┗━━━┛ ┃ ┃ ┗━━━┛
┏━━━━━━━━━┓ ┏━┻━┓ ┏━━━━━━━━━┓ ┏━┻━┓ ┏━━━━━━━━━━┓ ┏━┻━┓ ┏━━━━━━━━━┓ ┏━┻━┓
q1: ──┨ RX(π/2) ┠─┨╺╋╸┠─┨ RZ(2*a) ┠─┨╺╋╸┠─┨ RX(7π/2) ┠─┨╺╋╸┠─┨ RZ(2*b) ┠─┨╺╋╸┠─────────
┗━━━━━━━━━┛ ┗━━━┛ ┗━━━━━━━━━┛ ┗━━━┛ ┗━━━━━━━━━━┛ ┗━━━┛ ┗━━━━━━━━━┛ ┗━━━┛
>>> TimeEvolution(ops2).circuit
┏━━━┓ ┏━━━┓
q0: ──┨ H ┠───■─────────────────■───┨ H ┠─────────■─────────────────■──────────────────
┗━━━┛ ┃ ┃ ┗━━━┛ ┃ ┃
┏━┻━┓ ┏━━━━━━━━━┓ ┏━┻━┓ ┏━━━━━━━━━┓ ┏━┻━┓ ┏━━━━━━━━━┓ ┏━┻━┓ ┏━━━━━━━━━━┓
q1: ────────┨╺╋╸┠─┨ RZ(2*b) ┠─┨╺╋╸┠─┨ RX(π/2) ┠─┨╺╋╸┠─┨ RZ(2*a) ┠─┨╺╋╸┠─┨ RX(7π/2) ┠───
┗━━━┛ ┗━━━━━━━━━┛ ┗━━━┛ ┗━━━━━━━━━┛ ┗━━━┛ ┗━━━━━━━━━┛ ┗━━━┛ ┗━━━━━━━━━━┛
"""
def __init__(self, ops: QubitOperator, time=None):
"""Initialize a TimeEvolution object."""
# pylint: disable=import-outside-toplevel
from mindquantum.utils.type_value_check import _num_type
if time is None:
time = 1
if not isinstance(time, _num_type):
if not isinstance(time, (dict, ParameterResolver)):
raise TypeError(f"time requires a number or a dict or a ParameterResolver, but get {type(time)}")
time = ParameterResolver(time)
self.time = time
if isinstance(time, dict):
self.time = ParameterResolver(time)
self.ops = ops
@property
def circuit(self):
"""Get the first order trotter decomposition circuit of this time evolution operator."""
from ..circuit import Circuit # pylint: disable=import-outside-toplevel
from ..gates import GlobalPhase # pylint: disable=import-outside-toplevel
circ = Circuit()
for k, v in self.ops.terms.items():
pr_tmp = self.time * v
if k:
tmp_circ = decompose_single_term_time_evolution(k, pr_tmp)
else:
tmp_circ = GlobalPhase(pr_tmp).on(0)
circ += tmp_circ
return circ