def generate_XXZ_hamiltonian(num_qubits, delta):
#epsilon = 0.5
if num_qubits == 1:
raise(RuntimeError('Cannot generate Hamiltonian with 1 qubit'))
else:
hamiltonian = hcp.heisenberg_xyz_model(num_qubits, jx = 1, jy = 1, jz = delta)
return hamiltonian
import warnings
# warnings.filterwarnings("ignore", category=DeprecationWarning)
#Parameters
uptowhatK = 3
num_qubits = 3
endtime = 8
num_steps = 1001
optimizer = 'zvode'
p_invcond = 10**(-12)
#create initial state
initial_state = acp.Initialstate(num_qubits, "efficient_SU2", 123, 5)
#define Hamiltonian
hamiltonian = hcp.heisenberg_xyz_model(num_qubits)
#create Initial Ansatz for K = 0
ansatz = acp.initial_ansatz(num_qubits)
#finalresults
finalresults = []
finalresults.append(ansatz)
#Run QAS
for k in range(1, uptowhatK + 1):
print(k)
#Generate Ansatz for this round
ansatz = acp.gen_next_ansatz(ansatz, hamiltonian, num_qubits)