/**
* @param {function(!CircuitEvalContext)} updateAction
* @param {!Matrix} matrix
* @param {!number|undefined=} forcedTime
*/
function assertThatCircuitMutationActsLikeMatrix_single(updateAction, matrix, forcedTime=undefined) {
let qubitSpan = Math.round(Math.log2(matrix.height()));
let extraWires = Math.floor(Math.random()*5);
let time = Math.random();
let qubitIndex = Math.floor(Math.random() * extraWires);
if (USE_SIMPLE_VALUES) {
extraWires = 0;
time = 0;
qubitIndex = 0;
}
if (forcedTime !== undefined) {
time = forcedTime;
}
let wireCount = qubitSpan + extraWires;
let controls = Controls.NONE;
for (let i = 0; i < extraWires; i++) {
if (Math.random() < 0.5) {
controls = controls.and(Controls.bit(i + (i < qubitIndex ? 0 : qubitSpan), Math.random() < 0.5));
}
}
let ampCount = 1 << wireCount;
let inVec = Matrix.generate(1, ampCount, () => USE_SIMPLE_VALUES ?
(Math.random() < 0.5 ? 1 : 0) :
new Complex(Math.random()*10 - 5, Math.random()*10 - 5));
let tex = Shaders.vec2Data(inVec.rawBuffer()).toVec2Texture(wireCount);
let trader = new WglTextureTrader(tex);
let controlsTexture = CircuitShaders.controlMask(controls).toBoolTexture(wireCount);
let ctx = new CircuitEvalContext(
time,
qubitIndex,
wireCount,
controls,
controlsTexture,
trader,
new Map());
updateAction(ctx);
controlsTexture.deallocByDepositingInPool();
let outData = KetTextureUtil.tradeTextureForVec2Output(trader);
let outVec = new Matrix(1, ampCount, outData);
let expectedOutVec = matrix.applyToStateVectorAtQubitWithControls(inVec, qubitIndex, controls);
assertThat(outVec).withInfo({matrix, inVec, ctx}).isApproximatelyEqualTo(expectedOutVec, 0.005);
}
/**
* @param {!int} wireCount The number of wires in the circuit.
* @param {!function(!CircuitEvalContext)} updateAction The actual update action.
* @param {!function(!int) : !int} permutation The expected permutation.
* @param {*} permuteInfo Debug info included when the assertion fails.
*/
function assertThatCircuitUpdateActsLikePermutation(wireCount, updateAction, permutation, permuteInfo=undefined) {
let time = Math.random();
let ampCount = 1 << wireCount;
let inVec = Matrix.generate(1, ampCount, r => new Complex(r + Math.random(), Math.random()*1000));
let tex = Shaders.vec2Data(inVec.rawBuffer()).toVec2Texture(wireCount);
let trader = new WglTextureTrader(tex);
let controlsTexture = CircuitShaders.controlMask(Controls.NONE).toBoolTexture(wireCount);
let ctx = new CircuitEvalContext(
time,
0,
wireCount,
Controls.NONE,
controlsTexture,
trader,
new Map());
updateAction(ctx);
controlsTexture.deallocByDepositingInPool();
let outData = KetTextureUtil.tradeTextureForVec2Output(trader);
let outVec = new Matrix(1, ampCount, outData);
for (let i = 0; i < ampCount; i++) {
let j = permutation(i);
let inVal = inVec.cell(0, i);
let outVal = outVec.cell(0, j);
if (!outVal.isApproximatelyEqualTo(inVal, 0.001)) {
let actualIn = Math.floor(outVec.cell(0, j).real);
let actualOut = Seq.range(ampCount).filter(k => Math.floor(outVec.cell(0, k).real) === i).first('[NONE]');
assertThat(outVal).
withInfo({i, j, actualIn, actualOut, permuteInfo}).
isApproximatelyEqualTo(inVal, 0.01);
}
}
// Increment assertion count.
assertTrue(true);
}