Julia v1.0 interpolation indexing

[matthewware]

The microwave control logic maps waveform information to drive term amplitudes in

QSimulator.jl/src/controls.jl

Line 45 in 240fe35

return mc.sequence(t/mc.timeStep)*cos(2*pi*mc.freq*t + mc.phase)

Issues pop up when the mapping from scaled time to sequence index is close to zero (Julia uses 1-based indexing). This issue could have arisen from the Interpolations.jl changes (note OnGrid() is gone). More investigation is needed.

[matthewware]

Minimal working example:

function run_simulation(sys, rhoIn, measOp, times)
    populations = zeros(length(times)-1)
    U = eye(ComplexF64, dim(sys))
    dt = 0.01
    rho = rhoIn
    for ct = 1:length(times)-1
        U = unitary_propagator(sys, dt, times[ct], times[ct+1])
        populations[ct] = real(tr(U*rho*U'*measOp))
        rho = U*rho*U'
    end
    return populations
end

times = 0.0:200.0
populations = run_simulation(sys, rhoIn, measOp, times);

[matthewware]

Seems like forcing amplitue to 0.0 on under and over runs is doing the trick. I'll probably merge into master soon.