internal Register(ComputerModel compModel, Quantum.Register reg, RegisterModel regModel, int offsetToModel = 0)
{
_compModel = compModel;
_register = reg;
_model = regModel;
_offsetToModel = offsetToModel;
}
internal Register(ComputerModel compModel, Quantum.Register reg, RegisterModel regModel, int offsetToModel = 0)
{
_compModel = compModel;
_register = reg;
_model = regModel;
_offsetToModel = offsetToModel;
}
public Register this[int offset, int width]
{
get
{
ValidateOffsetWidth(offset, width);
Quantum.Register reg = _register[offset, width];
return(new Register(_compModel, reg, _model, _offsetToModel + offset));
}
}
public bool RunStep(IList <Gate> gates, bool runBackward = false, bool runComposite = false)
{
Quantum.Register root = _comp.GetSourceRoot();
bool somethingChanged = false;
int i = 0;
int maxI = (runComposite ? 1 : gates.Count);
while (i < maxI)
{
Gate gate = gates[i];
int? control = null;
if (gate.Control.HasValue)
{
control = gate.Control.Value.OffsetToRoot;
}
switch (gate.Name)
{
case GateName.Hadamard: // the same as inversion
root.Hadamard(gate.Target.OffsetToRoot, control);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.SigmaX: // the same as inversion
root.SigmaX(gate.Target.OffsetToRoot);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.SigmaY: // the same as inversion
root.SigmaY(gate.Target.OffsetToRoot, control);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.SigmaZ: // the same as inversion
root.SigmaZ(gate.Target.OffsetToRoot, control);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.SqrtX:
if (runBackward)
{
Complex[,] m = new Complex[2, 2];
Complex a = new Complex(0.5, -0.5);
Complex b = new Complex(0.5, 0.5);
m[0, 0] = b;
m[0, 1] = a;
m[1, 0] = a;
m[1, 1] = b;
root.Gate1(m, gate.Target.OffsetToRoot, control);
}
else
{
root.SqrtX(gate.Target.OffsetToRoot, control);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.RotateX:
RotateXGate rx = gate as RotateXGate;
if (runBackward)
{
root.RotateX(-rx.Gamma, rx.Target.OffsetToRoot, control);
}
else
{
root.RotateX(rx.Gamma, rx.Target.OffsetToRoot, control);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.RotateY:
RotateYGate ry = gate as RotateYGate;
if (runBackward)
{
root.RotateY(-ry.Gamma, ry.Target.OffsetToRoot, control);
}
else
{
root.RotateY(ry.Gamma, ry.Target.OffsetToRoot, control);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.RotateZ:
RotateZGate rz = gate as RotateZGate;
if (runBackward)
{
root.RotateZ(-rz.Gamma, rz.Target.OffsetToRoot, control);
}
else
{
root.RotateZ(rz.Gamma, rz.Target.OffsetToRoot, control);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.PhaseKick:
PhaseKickGate pk = gate as PhaseKickGate;
if (pk.Controls.Length > 0)
{
int[] controls1 = pk.Controls.Select <RegisterRefModel, int>(x => x.OffsetToRoot).ToArray <int>();
if (runBackward)
{
root.PhaseKick(-pk.Gamma, pk.Target.OffsetToRoot, controls1);
}
else
{
root.PhaseKick(pk.Gamma, pk.Target.OffsetToRoot, controls1);
}
}
else
{
if (runBackward)
{
root.PhaseKick(-pk.Gamma, pk.Target.OffsetToRoot);
}
else
{
root.PhaseKick(pk.Gamma, pk.Target.OffsetToRoot);
}
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.PhaseScale:
PhaseScaleGate ps = gate as PhaseScaleGate;
if (runBackward)
{
root.PhaseScale(-ps.Gamma, ps.Target.OffsetToRoot, control);
}
else
{
root.PhaseScale(ps.Gamma, ps.Target.OffsetToRoot, control);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.CPhaseShift:
CPhaseShiftGate cps = gate as CPhaseShiftGate;
if (cps.Controls.Length > 0)
{
int[] controls1 = cps.Controls.Select <RegisterRefModel, int>(x => x.OffsetToRoot).ToArray <int>();
if (runBackward)
{
root.InverseCPhaseShift(cps.Dist, cps.Target.OffsetToRoot, controls1);
}
else
{
root.CPhaseShift(cps.Dist, cps.Target.OffsetToRoot, controls1);
}
}
else
{
if (runBackward)
{
root.InverseCPhaseShift(cps.Dist, cps.Target.OffsetToRoot);
}
else
{
root.CPhaseShift(cps.Dist, cps.Target.OffsetToRoot);
}
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.InvCPhaseShift:
InvCPhaseShiftGate icps = gate as InvCPhaseShiftGate;
if (icps.Controls.Length > 0)
{
int[] controls1 = icps.Controls.Select <RegisterRefModel, int>(x => x.OffsetToRoot).ToArray <int>();
if (runBackward)
{
root.CPhaseShift(icps.Dist, icps.Target.OffsetToRoot, controls1);
}
else
{
root.InverseCPhaseShift(icps.Dist, icps.Target.OffsetToRoot, controls1);
}
}
else
{
if (runBackward)
{
root.CPhaseShift(icps.Dist, icps.Target.OffsetToRoot);
}
else
{
root.InverseCPhaseShift(icps.Dist, icps.Target.OffsetToRoot);
}
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.CNot:
CNotGate cn = gate as CNotGate;
root.CNot(cn.Target.OffsetToRoot, cn.Control.Value.OffsetToRoot);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.Toffoli:
ToffoliGate t = gate as ToffoliGate;
int[] controls = t.Controls.Select <RegisterRefModel, int>(x => x.OffsetToRoot).ToArray <int>();
root.Toffoli(t.Target.OffsetToRoot, controls);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.Unitary:
UnitaryGate u = gate as UnitaryGate;
if (runBackward)
{
Complex[,] m = new Complex[2, 2];
m[0, 0] = Complex.Conjugate(u.Matrix[0, 0]);
m[0, 1] = Complex.Conjugate(u.Matrix[1, 0]);
m[1, 0] = Complex.Conjugate(u.Matrix[0, 1]);
m[1, 1] = Complex.Conjugate(u.Matrix[1, 1]);
root.Gate1(m, gate.Target.OffsetToRoot, control);
}
else
{
root.Gate1(u.Matrix, u.Target.OffsetToRoot, control);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.Measure:
MeasureGate mg = gate as MeasureGate;
for (int j = mg.Begin; j <= mg.End; j++)
{
root.Measure(j);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.Parametric:
ParametricGate pg = gate as ParametricGate;
_comp.Evaluate(pg, runBackward);
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.Composite:
CompositeGate cg = gate as CompositeGate;
List <Gate> actual = _comp.GetActualGates(cg);
if (runBackward)
{
actual.Reverse();
}
foreach (Gate g in actual)
{
// TODO refactor RunStep() into RunGate()
List <Gate> toRun = new List <Gate>()
{
g
};
RunStep(toRun, runBackward, true);
}
somethingChanged = true;
i = gate.End + 1;
break;
case GateName.Empty:
default:
i++;
break;
}
}
return(somethingChanged);
}