In my test, I change the code and define Xb::Vector, o::Vector,... etc but the result doesn't change. Is there some special way to define them?
Regards Em quarta-feira, 26 de novembro de 2014 17h50min04s UTC-2, Stefan Karpinski escreveu: > > It turns out that doesn't actually matter but there are other unrelated > ill-typed variables in the code: > > Xb, o, y, r, e, Nj, ind, q, t > > > I haven't determined why this is happening, but something is ill-typed > here. > > On Wed, Nov 26, 2014 at 1:20 PM, Tim Holy <tim....@gmail.com <javascript:> > > wrote: > >> Stefan's point is more relevant than you realize; if the type of the >> variable >> is uncertain, the loop will be very slow. You can make that resize >> operation a >> separate function, from which you call your main routine. >> >> --Tim >> >> On Wednesday, November 26, 2014 09:26:35 AM Emerson Vitor Castelani wrote: >> > Humm, It's useful but I guess that the problem is something in loop >> while... >> > >> > Em quarta-feira, 26 de novembro de 2014 15h05min03s UTC-2, Stefan >> Karpinski >> > >> > escreveu: >> > > I'm not sure if this is the problem, but changing the type of a >> variable >> > > in a function body causes problems for type inference. For that >> reason, >> > > the >> > > first two lines of this code may cause performance issue if you call >> this >> > > with b and c as matrices. A more Julian idiom for this is to do >> something >> > > like this: >> > > >> > > function simplex(A::Matrix, b::Vector, c::Vector) >> > > >> > > (m,n) = size(A); >> > > ... >> > > >> > > end >> > > simplex(A::Matrix, b::Matrix, c::Matrix) = simplex(A, vec(b), vec(c)) >> > > >> > > >> > > This also avoids copying the data of b and c needlessly. >> > > >> > > On Wednesday, November 26, 2014, Emerson Vitor Castelani < >> > > >> > > emerso...@gmail.com <javascript:>> wrote: >> > >> Have you considered the example scsd8.mat? >> > >> >> > >> Em quarta-feira, 26 de novembro de 2014 14h39min14s UTC-2, Pileas >> > >> >> > >> escreveu: >> > >>> Result with tic() toc() at the very beginning and at the very end: >> > >>> >> > >>> elapsed time: 0.025719973 seconds >> > >>> >> > >>> >> > >>> Τη Τετάρτη, 26 Νοεμβρίου 2014 11:06:26 π.μ. UTC-5, ο χρήστης Emerson >> > >>> >> > >>> Vitor Castelani έγραψε: >> > >>>> Well, I am tried to implement a simple version of simplex in Julia >> and >> > >>>> I have had some troubles. In Julia, my algorithm spends about 30 >> sec >> > >>>> and in >> > >>>> matlab/octave 3 sec for the same problem. So, I saw some tips in >> order >> > >>>> to >> > >>>> get a better performance but the best that I got in Julia was 17-20 >> > >>>> sec. >> > >>>> The codes are in attachment. I am new in Julia and the algorithms >> are >> > >>>> little roughly implemented but they are very similar. >> > >>>> >> > >>>> Thanks >> >> >