cirqやqiskitを使えば快適に計算できると思いますが、blueqatSDKでも実装が可能で、まだチューニングしてませんが、比較的高速にいろいろ計算ができそうです。
blueqatSDKからcuStateVecをつないでNVIDIA A100 メモリ80Gのマシンから実行してみました。
単体のマシンからは29量子ビットがメモリの関係で最大でした。
マルチGPUですともっといけますので、持ってる方はぜひチャレンジしてみてください。その際にはcirq+qsim+cuStateVecになるとおもいます。
from blueqat import Circuit
import time
N = 29
ite = 100
start = time.time()
circ = Circuit(N)
for j in range(ite):
circ.x[:].y[:].z[:].h[:]
for i in range(N):
circ.cx[i,i+1]
result = circ.run(backend="cusv")
t = time.time() - start
print(result)
[ 9.9995494e-01+0.j -1.2907941e-07+0.j 9.6412592e-08+0.j ...
1.9410919e-08+0.j -8.7635840e-08+0.j -8.3126771e-08+0.j]
440.07840061187744
7分ほどかかりましたが複雑なもつれ計算も十分できました。弊社では要望があればチューニングをして高速化をします。以上です。