Then<br>Scilab 5.3.2. spends about 0.4sec<br>Scilab 5.3.1. spends about 5 sec<br>Scilab 4.1.2. spends about 0.15sec<br><br>However, the code result is different from the previous one:<br>
<br>n=100000;<br>sigeps=1;<br>mi=grand(n,1,'nor',0,sigeps);<br>for i=2:n;<br>    mi(i)=mi(i)+mi(i-1);<br>end<br><br>produces a random walk, but<br>
<br>n=100000;<br>sigeps=1;<br>mi=grand(n,1,'nor',0,sigeps);<br>mi(2:n)=mi(2:n)+mi(1:n-1);<br><br>does not.<br>
<br>
Gin.<br><br><br><br><div class="gmail_quote">On Wed, Jul 20, 2011 at 11:50 AM, Allan CORNET <span dir="ltr"><<a href="mailto:allan.cornet@scilab.org">allan.cornet@scilab.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div link="blue" vlink="purple" lang="FR"><div><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">Hi,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">You know that you can remove loop in  your script ?<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">Try:<u></u><u></u></span></p><div class="im"><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">tic()<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">sigeta=2; sigeps=1;<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">mi=[];<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">c=[];<u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11.0pt;color:#1F497D">y=[];<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D">skait=[];<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D">sauc2=[];<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D">rho=[];<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D">f1=1.2;<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">f2=-0.5;<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">n=102000;<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">k=2000;<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">c=grand(n,1,'nor',0,sigeta);<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">c(2)=f1*c(1)+c(2);<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">mi=grand(n,1,'nor',0,sigeps);<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p>
</div><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">c(3:n) = c(3:n)+f1*c(2:n-1)+f2*c(1:n-2);<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">mi(2:n)=mi(2:n)+mi(1:n-1);<u></u><u></u></span></p><div class="im"><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">y=c+mi;<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">y(1:k)=[];<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">c(1:k)=[];<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">i=1:length(y);<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">skait=sum((y(i)-mean(y)).*(c(i)-mean(c)));<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">sauc2=sum((y(i)-mean(y))^2)*sum((c(i)-mean(c))^2);<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">rho=skait/sqrt(sauc2)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">toc()<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p></div><p class="MsoNormal"><span style="font-size:11.0pt;color:#1F497D" lang="EN-US">Allan<u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11.0pt;color:#1F497D" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><b><span style="font-size:10.0pt">De :</span></b><span style="font-size:10.0pt"> Ginters Bušs [mailto:<a href="mailto:ginters.buss@gmail.com" target="_blank">ginters.buss@gmail.com</a>] <br>
<b>Envoyé :</b> mercredi 20 juillet 2011 10:35<br><b>À :</b> <a href="mailto:users@lists.scilab.org" target="_blank">users@lists.scilab.org</a><br><b>Objet :</b> [scilab-Users] Re: speed of Scilab 5.3.2 vs Scilab 5.3.1<u></u><u></u></span></p>
<div><div></div><div class="h5"><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal" style="margin-bottom:12.0pt">And about 1.65sec. on Scilab 4.1.2.<br><br>So, still Scilab 4 is faster, but only marginally.<br><br>
<u></u><u></u></p><div><p class="MsoNormal">On Wed, Jul 20, 2011 at 11:23 AM, Ginters Bušs <<a href="mailto:ginters.buss@gmail.com" target="_blank">ginters.buss@gmail.com</a>> wrote:<u></u><u></u></p><p class="MsoNormal">
Hi,<br><br>Just to let you know:<br>I just found out that version 5.3.2. runs faster than 5.3.1. <br><br>Say, this code:<br><br>tic()<br>sigeta=2; sigeps=1;<br>mi=[]; c=[]; y=[]; skait=[]; sauc2=[]; rho=[];<br>f1=1.2; f2=-0.5;<br>
n=102000; k=2000;<br>c=grand(n,1,'nor',0,sigeta);<br>c(2)=f1*c(1)+c(2);<br>mi=grand(n,1,'nor',0,sigeps);<br>for i=3:n;<br>    c(i)=c(i)+f1*c(i-1)+f2*c(i-2);<br>end<br>for i=2:n;<br>    mi(i)=mi(i)+mi(i-1);<br>
end<br>y=c+mi;<br>y(1:k)=[];<br>c(1:k)=[];<br>i=1:length(y);<br>skait=sum((y(i)-mean(y)).*(c(i)-mean(c)));<br>sauc2=sum((y(i)-mean(y))^2)*sum((c(i)-mean(c))^2);<br>rho=skait/sqrt(sauc2)<br>toc()<br><br>takes approx. 1.75sec. (on my PC) on Scilab 5.3.2, while approx. 6 sec. on Scilab 5.3.1. <br>
<br>Gin.<u></u><u></u></p></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div></blockquote></div><br>