Qiskitで任意のユニタリー行列から量子ゲート(量子回路)を作る

ここ数ヶ月何度か同じコードが必要になる(聞かれる)ケースがあったので、コードサンプルを載せておきます。
何かしようとすると意外と使いたい機能な気がします。

コードは2020/3/13現在のものです。
もし間違いがあったら教えてください。

from qiskit import QuantumCircuit
from qiskit.quantum_info.operators.operator import Operator
import numpy as np

# 任意の数の量子ビットと古典ビットを持つ量子回路を用意(今回は適当に2qubit)
qc = QuantumCircuit(2,2)

# ユニタリー行列を用意(今回は適当に単位行列)
unitary = np.array([[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1]])

# ユニタリー行列からオペレーターを作成
operator = Operator(unitary)

# 量子回路(qubit0番と1番)にオペレーターを作用(ユニタリーでなかった場合はエラーが出る)
qc.append(operator, [0,1])

# 実行する部分は省略