Skip to content

Instantly share code, notes, and snippets.

Avatar

David Smith davidssmith

  • Vanderbilt University
  • Nashville, USA
View GitHub Profile
@davidssmith
davidssmith / tv_fista_demo.jl
Last active Dec 29, 2015
This is 10x slower in Julia than in Python/Anaconda. Why?
View tv_fista_demo.jl
function tv_fista(X::Array{Complex{Float64},2}; mu=0.1, niter=10)
m, n = size(X);
P = Any[zeros(Complex{Float64}, (m-1,n)), zeros(Complex{Float64}, (m,n-1))];
R = Any[zeros(Complex{Float64}, (m-1,n)), zeros(Complex{Float64}, (m,n-1))];
tk = 1;
tkp1 = 1;
D = zeros(Complex{Float64}, (m,n));
fval = Inf;
obj = zeros(niter);
for t = 1:niter