From akhorshidi at live.com Sun Feb 2 09:49:56 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 2 Feb 2014 12:19:56 +0330 Subject: [Scilab-users] export graphics to EMF Message-ID: Hello; I wanna export Scilab graphics directly into emf format but xs2emf function returns unexpected error message! --> plot2d(); --> xs2emf(0,'plotex.emf') !--error 10000 xs2emf: Unable to execute pstoedit. at line 95 of function xs2emf called by : xs2emf(0,'plotex.emf') Before using this function, It's necessary to install the following software: ImageMagick-dynamic PStoEdit Ghostscript Moreover we need to edit path names to point to the correct files. So I modified the file "xs2emf.sci" as below: line# 52: if fileinfo("C:\Program Files\gs\gs9.07\bin\gsdll32.dll") == [] then line# 70: pstoeditPath = "C:\Program Files\pstoedit\pstoedit.exe"; Also there is an related bug report on bugzilla. Please see: http://bugzilla.scilab.org/show_bug.cgi?id=11853 So, how can I overcome this issue?! However, I can use pstoedit from other software such as GSview and it convert eps files to emf properly. Furthermore, it has been mentioned that PStoEdit has a DLL which can be used directly by GSview! (see GSview's homepage) Merci, Mehran _ -------------- next part -------------- An HTML attachment was scrubbed... URL: From akhorshidi at live.com Sun Feb 2 10:05:01 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 2 Feb 2014 12:35:01 +0330 Subject: [Scilab-users] clc() icon with additional functionality !!? Message-ID: Hi folk; The addition of the clc() icon was a good idea but clearing console by using clc() command is quicker than that icon. So I decided to change the usage of that icon on my Scilab so that the newer one can clear local and global variable as well as clearing console. Here is my changes: SCI>modules>gui>etc>main_toolbar.xml: After restarting Scilab, the new clc icon must can do all the above tasks but it cannot clear local and global variables! Actually, I think, the "instruction" attribute cannot perform some scilab commands sush as : clear(), clearglobal(),.. However there are many instructions which work fine: clc(), resethistory(), chdir(),.. So where is the problem? Thanks in advance. Be free and opensource! ;) Mehran _ -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Sun Feb 2 13:50:41 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 02 Feb 2014 13:50:41 +0100 Subject: [Scilab-users] clc() icon with additional functionality !!? In-Reply-To: References: Message-ID: <52EE3F21.2000701@free.fr> Hi Mehran, Le 02/02/2014 10:05, A Khorshidi a ?crit : > The addition of the clc() icon was a good idea but clearing console by > using clc() command is quicker than that icon. Even without the (), it is quicker for me to click the icon. And for newbies (as students, many students..) that do not know clc, the icon is explicit. May be are you using a padd rather than a mouse ;) > So I decided to change the usage of that icon on my Scilab so that the > newer one can clear local and global variable as well as clearing > console. > Are you coming from Matlab? Matlabers and former ones use very often /clear all/, /close all/ at the beginning of their scripts. I do not know where this "hygienism" is coming from. clear is much more brutal in Scilab than in Matlab. I never used them in the head of my scripts for 16 years, and i have never had problems for that. Some students or colleagues knowing a bit or a lot Matlab and this habit and coming to Scilab keeping the same have problems with that (but they do not know that this is the reason), since then libraries (sometimes local, loaded from ATOMS, etc) are also cleared, etc, and then many things get harder ;) Indeed, embedded Scilab libraries are automatically restored when first calling one of their functions, but it is not the case for other external libraries. So, definitively, i will keep the clear icon as is. Cheers Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Sun Feb 2 14:18:06 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 02 Feb 2014 14:18:06 +0100 Subject: [Scilab-users] clc() icon with additional functionality !!? In-Reply-To: References: Message-ID: <52EE458E.2050106@free.fr> Le 02/02/2014 10:05, A Khorshidi a ?crit : > > > After restarting Scilab, the new clc icon must can do all the above > tasks but it cannot clear local and global variables! > Actually, I think, the "instruction" attribute cannot perform some > scilab commands sush as : clear(), clearglobal(),.. > However there are many instructions which work fine: clc(), > resethistory(), chdir(),.. > > So where is the problem? My best guess would be that in the callback, you should drop the atomic bombs clear() (and may be clearglobal() if the callback instruction has survived to the clear() that it performed) as /last/ instructions. clear() also clears the /gcbo/ variable: So your proposal looks like a gcbo-callback suicide. And it is somewhat preferable to kiss one's mother /before/ suiciding oneself than /after/. AAaaa /clear all, close all/... Nice custom! :-) -------------- next part -------------- An HTML attachment was scrubbed... URL: From akhorshidi at live.com Sun Feb 2 14:42:05 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 2 Feb 2014 05:42:05 -0800 (PST) Subject: [Scilab-users] export graphics to EMF In-Reply-To: References: Message-ID: <1391348525659-4028496.post@n3.nabble.com> Thanks to the following article, now I can convert graphics format to each other in Scilab command line: http://how-to.wikia.com/wiki/How_to_convert_one_image_format_to_another_via_command_line by using "convert" command included in ImageMagick: unix_s('pstoedit -f emf plotex.eps plotex.emf') or by using "pstoedit" command from PStoEdit: unix_s('pstoedit -f emf plotex.eps plotex.emf') However, to execute the last command, it's necessary to change the working directory to PStoEdit installation directory. But above all, I create a patch which can fix that given bug. HTH Mehran _ -- View this message in context: http://mailinglists.scilab.org/Scilab-users-export-graphics-to-EMF-tp4028492p4028496.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From akhorshidi at live.com Sun Feb 2 14:43:20 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 2 Feb 2014 05:43:20 -0800 (PST) Subject: [Scilab-users] clc() icon with additional functionality !!? In-Reply-To: <52EE3F21.2000701@free.fr> References: <52EE3F21.2000701@free.fr> Message-ID: <1391348600814-4028497.post@n3.nabble.com> Samuel GOUGEON wrote > [..] > >> So I decided to change the usage of that icon on my Scilab so that the >> newer one can clear local and global variable as well as clearing >> console. >> > Are you coming from Matlab? Matlabers and former ones use very often > /clear all/, /close all/ at the beginning of their scripts. > I do not know where this "hygienism" is coming from. clear is much more > brutal in Scilab > than in Matlab. I never used them in the head of my scripts for 16 > years, and i have never had problems for that. > Some students or colleagues knowing a bit or a lot Matlab and this habit > and coming to Scilab keeping the same > have problems with that (but they do not know that this is the reason), > since then libraries (sometimes local, loaded > from ATOMS, etc) are also cleared, etc, and then many things get harder > ;) Indeed, embedded Scilab libraries are > automatically restored when first calling one of their functions, but it > is not the case for other external libraries. > > So, definitively, i will keep the clear icon as is. > [..] Samuel, I agree with you, and I do not intend to force my changes on any other user! ;) However, I wanna modify it only on my Scilab but there is a problem in calling some Scilab function from that XML file. -- View this message in context: http://mailinglists.scilab.org/Scilab-users-clc-icon-with-additional-functionality-tp4028493p4028497.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From JDHeinzmann at gmail.com Mon Feb 3 00:08:44 2014 From: JDHeinzmann at gmail.com (JD Heinzmann) Date: Sun, 2 Feb 2014 15:08:44 -0800 (PST) Subject: [Scilab-users] Trying to use old Scicos file with Xcos Message-ID: <1391382524436-4028498.post@n3.nabble.com> About 3 years ago I developed a Scicos simulation of a heat engine using ScicosLab-4.4b8 with Scicos 4.4. I have just downloaded Scilab 5.4.1 and would like to run this old simulation under Xcos. I can open the old .cos file in Xcos and see the blocks and navigate into the super-block hierarchy just fine. I can save this block diagram as .xcos and .zcos files just fine. I can load the context file into Scilab and run it to set up the workspace before running the simulation. But when I try to run the simulation under Xcos, Scilab and Xcos immediately quits silently (with no error messages). Are Scicos files not compatible with Xcos? If not, what do I need to do to adjust my model to work with Xcos? Thank you! JD -- View this message in context: http://mailinglists.scilab.org/Trying-to-use-old-Scicos-file-with-Xcos-tp4028498.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Mon Feb 3 09:17:53 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 3 Feb 2014 09:17:53 +0100 Subject: [Scilab-users] Intersection of wavy surfaces In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F326346C56@301EX00100.sidel.com> References: <3B5FFC67498DFF49AE7271A584867D16F3262F6C20@301EX00100.sidel.com> <52E95320.2090802@free.fr> <3B5FFC67498DFF49AE7271A584867D16F326346A0F@301EX00100.sidel.com> <3B5FFC67498DFF49AE7271A584867D16F326346C56@301EX00100.sidel.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3263E77D7@301EX00100.sidel.com> Hello, concerning the presumed "bug of SVG export" I mentioned, finally there's no bug, the export is just quite long because of the size of the file (5MB). Maybe I could suggest some kind of progress bar. Regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From Christophe.Dang at sidel.com Mon Feb 3 09:31:54 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 3 Feb 2014 09:31:54 +0100 Subject: [Scilab-users] clear() or not clear()? [was: clc() icon with additional functionality !!?] In-Reply-To: <52EE3F21.2000701@free.fr> References: <52EE3F21.2000701@free.fr> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3263E7844@301EX00100.sidel.com> Hello, > Matlabers and former ones use very often clear all, close all at the > beginning of their scripts. > I do not know where this "hygienism" is coming from. > clear is much more brutal in Scilab than in Matlab. > I never used them in the head of my scripts for 16 years Never been a Matlaber. Personally, I use clear() during the development of my scripts: I sometimes forget to initiallise some variables, and if it gets the value from a past execution, then it prevents me from seeing that the variable does not exist in a given case. After that, I don't to remove it in my final version... Probably I should. Maybe there's a beter way do spot such problems. -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From ian.h.bell at gmail.com Mon Feb 3 23:39:46 2014 From: ian.h.bell at gmail.com (Ian Bell) Date: Mon, 3 Feb 2014 23:39:46 +0100 Subject: [Scilab-users] Determine if 32-bit or 64-bit installation of scilab in code Message-ID: I need to determine in code if the active scilab version is 32-bit or 64-bit in order to load the right dynamically loaded library on windows. How can I determine whether scilab is 32-bit or 64-bit? Thanks, Ian -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Tue Feb 4 00:49:37 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 04 Feb 2014 00:49:37 +0100 Subject: [Scilab-users] Determine if 32-bit or 64-bit installation of scilab in code In-Reply-To: References: Message-ID: <52F02B11.2050303@free.fr> Le 03/02/2014 23:39, Ian Bell a ?crit : > I need to determine in code if the active scilab version is 32-bit or > 64-bit in order to load the right dynamically loaded library on windows. > > How can I determine whether scilab is 32-bit or 64-bit? maybe with [v, o] = getversion(), o(2)=="x64" Samuel From ian.h.bell at gmail.com Tue Feb 4 09:05:45 2014 From: ian.h.bell at gmail.com (Ian Bell) Date: Tue, 4 Feb 2014 09:05:45 +0100 Subject: [Scilab-users] Determine if 32-bit or 64-bit installation of scilab in code In-Reply-To: <52F02B11.2050303@free.fr> References: <52F02B11.2050303@free.fr> Message-ID: Brilliant, works great. Thanks. On Tue, Feb 4, 2014 at 12:49 AM, Samuel Gougeon wrote: > Le 03/02/2014 23:39, Ian Bell a ?crit : > > I need to determine in code if the active scilab version is 32-bit or >> 64-bit in order to load the right dynamically loaded library on windows. >> >> How can I determine whether scilab is 32-bit or 64-bit? >> > maybe with > [v, o] = getversion(), > o(2)=="x64" > > Samuel > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From ebmb.sci at gmail.com Mon Feb 3 13:27:59 2014 From: ebmb.sci at gmail.com (E Bmb) Date: Mon, 3 Feb 2014 10:27:59 -0200 Subject: [Scilab-users] Genectic algorithm question Message-ID: Hello All, I've been used the Genetic Algorithm in Scilab 5.4.1 (with the optim_ga macro). After a lot of tests, I would like to report a strange behaviour. That is, from the second time that I ran the algorithm, I always note a time increment to get the results. For example, to optimize a schedule problem with 40 tasks, at the first time the results could be given in 5 minutes, but from the second time it will be much more than 5 min. (This time is only a example, ok?!) I also noted that when I close the Scilab environment and start it again, I can get my results, for example, in the same 5 minutes. (But, I always will need close and start to get time and results (in terms of quality) corrects.) So, my questions are: 1. What is happening with my Scilab? 2. Is that behaviour correct? 3. Somebody who works with optimization already had problems like this? Thanks in advance for attention! Best regards, Eduardo. -------------- next part -------------- An HTML attachment was scrubbed... URL: From ebmb.sci at gmail.com Mon Feb 3 13:33:29 2014 From: ebmb.sci at gmail.com (ebmb.sci) Date: Mon, 3 Feb 2014 04:33:29 -0800 (PST) Subject: [Scilab-users] Genetic Algorithm question Message-ID: <1391430808990-4028501.post@n3.nabble.com> Hello All, I've been used the Genetic Algorithm in Scilab 5.4.1 (with the optim_ga macro). After a lot of tests, I would like to report a strange behaviour. That is, from the second time that I ran the algorithm, I always note a time increment to get the results. For example, to optimize a schedule problem with 40 tasks, at the first time the results could be given in 5 minutes, but from the second time it will be much more than 5 min. (This time is only a example, ok?!) I also noted that when I close the Scilab environment and start it again, I can get my results, for example, in the same 5 minutes. (But, I always will need close and start to get time and results (in terms of quality) corrects.) So, my questions are: 1. What is happening with my Scilab? 2. Is that behaviour correct? 3. Somebody who works with optimization already had problems like this? Thanks in advance for attention! Best regards, Eduardo. -- View this message in context: http://mailinglists.scilab.org/Genetic-Algorithm-question-tp4028501.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From clement.david at scilab-enterprises.com Tue Feb 4 10:43:13 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Tue, 04 Feb 2014 10:43:13 +0100 Subject: [Scilab-users] Trying to use old Scicos file with Xcos In-Reply-To: <1391382524436-4028498.post@n3.nabble.com> References: <1391382524436-4028498.post@n3.nabble.com> Message-ID: <1391506993.5710.5.camel@paros> Hello, It should work fine but it is model related. Can you provide your model ? -- Cl?ment Le dimanche 02 f?vrier 2014 ? 15:08 -0800, JD Heinzmann a ?crit : > About 3 years ago I developed a Scicos simulation of a heat engine using > ScicosLab-4.4b8 with Scicos 4.4. I have just downloaded Scilab 5.4.1 and > would like to run this old simulation under Xcos. > > I can open the old .cos file in Xcos and see the blocks and navigate into > the super-block hierarchy just fine. I can save this block diagram as .xcos > and .zcos files just fine. I can load the context file into Scilab and run > it to set up the workspace before running the simulation. But when I try to > run the simulation under Xcos, Scilab and Xcos immediately quits silently > (with no error messages). > > Are Scicos files not compatible with Xcos? If not, what do I need to do to > adjust my model to work with Xcos? > > Thank you! > > JD > > > > -- > View this message in context: http://mailinglists.scilab.org/Trying-to-use-old-Scicos-file-with-Xcos-tp4028498.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From paul.bignier at scilab-enterprises.com Tue Feb 4 11:39:39 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Tue, 04 Feb 2014 11:39:39 +0100 Subject: [Scilab-users] Determine if 32-bit or 64-bit installation of scilab in code In-Reply-To: References: Message-ID: <52F0C36B.5020903@scilab-enterprises.com> Hello, [a b]=getversion(); b(2) should give you the answer :) Regards, Paul On 02/03/2014 11:39 PM, Ian Bell wrote: > I need to determine in code if the active scilab version is 32-bit or > 64-bit in order to load the right dynamically loaded library on windows. > > How can I determine whether scilab is 32-bit or 64-bit? > > Thanks, > Ian > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Scilab Engineer & Xcos Developer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.bignier at scilab-enterprises.com Tue Feb 4 12:02:57 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Tue, 04 Feb 2014 12:02:57 +0100 Subject: [Scilab-users] Genectic algorithm question In-Reply-To: References: Message-ID: <52F0C8E1.8010903@scilab-enterprises.com> Hi Eduardo, What version of Scilab are you running, and under which OS please? Could you please provide a test case so we can reproduce the bug? As far as I know, this behavior is abnormal. Don't hesitate to report it on our bugzilla . Regards, Paul On 02/03/2014 01:27 PM, E Bmb wrote: > Hello All, > I've been used the Genetic Algorithm in Scilab 5.4.1 (with the > optim_ga macro). After a lot of tests, I would like to report a > strange behaviour. That is, from the second time that I ran the > algorithm, I always note a time increment to get the results. For > example, to optimize a schedule problem with 40 tasks, at the first > time the results could be given in 5 minutes, but from the second time > it will be much more than 5 min. (This time is only a example, ok?!) > > I also noted that when I close the Scilab environment and start it > again, I can get my results, for example, in the same 5 minutes. (But, > I always will need close and start to get time and results (in terms > of quality) corrects.) > > So, my questions are: > > 1. What is happening with my Scilab? > 2. Is that behaviour correct? > 3. Somebody who works with optimization already had problems like this? > > Thanks in advance for attention! > > Best regards, > Eduardo. > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Scilab Engineer & Xcos Developer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From contact at pierre-vuillemin.fr Wed Feb 5 09:56:01 2014 From: contact at pierre-vuillemin.fr (Pierre Vuillemin) Date: Wed, 05 Feb 2014 09:56:01 +0100 Subject: [Scilab-users] Interfacing fortran code with scilab Message-ID: <1391590561.4264.4.camel@Pelisse.lan> Hello all, I am currently learning fortran in order to speed up some part o my code and I would like to interface it with Scilab. To make some tests, I have used this fortran code : subroutine matmul(C,A,B,m,n,p) integer ::m,n,p,i,j,k real(kind=8),dimension(m,n),intent(in)::A real(kind=8),dimension(n,p),intent(in)::B real(kind=8),dimension(m,p),intent(out)::C real(kind=8):: temp C = 0 do j =1,p do k=1,n temp = B(k,j) do i=1,m C(i,j)=C(i,j) + A(i,k)*temp end do end do end do end subroutine matmul which performs a matrix multiplication. Then I have used 'ilib_for_link' to create the link. It works but I have some questions : - I did the same thing in Python with f2py and the code is seemingly faster in python, I was wondering why ? (the precision is the same) - If I change the 'kind = 8' to 'kind = 4' in the fortran code, the function returns a result which completely wrong. I expected it to be wrong, but not that much : for instance, with A = 10 and B = 20, the call to the fortran function gives me '6.36D-314' What is the reason? (My comprehension of how fortran types work with the kind parameter is still limited, I know that kind=8 and kind=4 is not a very portable expression though) - If I replace real(kind=8),dimension(m,p),intent(out)::C by real(kind=8),dimension(:,:),intent(out)::C I get a segfault in Scilab or it just closes. I was wondering why ? (I thought it was something valid since gfortran still compiles) - Is it possible to keep the upper case in the name of fortran functions when linked to scilab ? (At the beginning, my function was called 'matMul' and it took me some time to figure out what the error message 'matMul is not an entry point' meant) Best regards, Pierre Vuillemin From paul.carrico at esterline.com Wed Feb 5 13:38:40 2014 From: paul.carrico at esterline.com (Carrico, Paul) Date: Wed, 5 Feb 2014 12:38:40 +0000 Subject: [Scilab-users] Interfacing fortran code with scilab In-Reply-To: <1391590561.4264.4.camel@Pelisse.lan> References: <1391590561.4264.4.camel@Pelisse.lan> Message-ID: <3A6B7233274DB449A2A0053A47684F95042502FC@BGS-EX01.auxitrol.ad> All before using fortan code, did you have a look to the "vectorization" in order to avoid loops (that drastically decrease speed) and to increase speedup consequently ? Just an advice Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Pierre Vuillemin Envoy??: mercredi 5 f?vrier 2014 09:56 ??: users at lists.scilab.org Objet?: [Scilab-users] Interfacing fortran code with scilab Hello all, I am currently learning fortran in order to speed up some part o my code and I would like to interface it with Scilab. To make some tests, I have used this fortran code : subroutine matmul(C,A,B,m,n,p) integer ::m,n,p,i,j,k real(kind=8),dimension(m,n),intent(in)::A real(kind=8),dimension(n,p),intent(in)::B real(kind=8),dimension(m,p),intent(out)::C real(kind=8):: temp C = 0 do j =1,p do k=1,n temp = B(k,j) do i=1,m C(i,j)=C(i,j) + A(i,k)*temp end do end do end do end subroutine matmul which performs a matrix multiplication. Then I have used 'ilib_for_link' to create the link. It works but I have some questions : - I did the same thing in Python with f2py and the code is seemingly faster in python, I was wondering why ? (the precision is the same) - If I change the 'kind = 8' to 'kind = 4' in the fortran code, the function returns a result which completely wrong. I expected it to be wrong, but not that much : for instance, with A = 10 and B = 20, the call to the fortran function gives me '6.36D-314' What is the reason? (My comprehension of how fortran types work with the kind parameter is still limited, I know that kind=8 and kind=4 is not a very portable expression though) - If I replace real(kind=8),dimension(m,p),intent(out)::C by real(kind=8),dimension(:,:),intent(out)::C I get a segfault in Scilab or it just closes. I was wondering why ? (I thought it was something valid since gfortran still compiles) - Is it possible to keep the upper case in the name of fortran functions when linked to scilab ? (At the beginning, my function was called 'matMul' and it took me some time to figure out what the error message 'matMul is not an entry point' meant) Best regards, Pierre Vuillemin _______________________________________________ users mailing list users at lists.scilab.org https://urldefense.proofpoint.com/v1/url?u=http://lists.scilab.org/mailman/listinfo/users&k=b2vlTQszY8VIpYRvaG%2By2A%3D%3D%0A&r=SzPQe21KEY%2BPdUhJge5w%2FdtbtLgIXw5YTsnbzx%2F7JYE%3D%0A&m=ULVia0lWcxnisV%2BYAJ894mfD5PBI7UEZDeqeSqE9PQw%3D%0A&s=4501ea13ac011e234cb1dd1cdf1810c4388736b832cf43301a84c41f80b767f5 -------------------------------------------------------------------------------- Le pr?sent mail et ses pi?ces jointes sont confidentiels et destin?s ? la personne ou aux personnes vis?e(s) ci-dessus. Si vous avez re?u cet e-mail par erreur, veuillez contacter imm?diatement l'exp?diteur et effacer le message de votre syst?me. Toute divulgation, copie ou distribution de cet e-mail est strictement interdite. This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error, please contact the sender and delete the email from your system. If you are not the named addressee you should not disseminate, distribute or copy this email. From contact at pierre-vuillemin.fr Wed Feb 5 14:08:27 2014 From: contact at pierre-vuillemin.fr (Pierre Vuillemin) Date: Wed, 05 Feb 2014 14:08:27 +0100 Subject: [Scilab-users] Interfacing fortran code with scilab In-Reply-To: <3A6B7233274DB449A2A0053A47684F95042502FC@BGS-EX01.auxitrol.ad> References: <1391590561.4264.4.camel@Pelisse.lan> <3A6B7233274DB449A2A0053A47684F95042502FC@BGS-EX01.auxitrol.ad> Message-ID: <1391605707.4264.25.camel@Pelisse.lan> I am quite used to work with Matlab, hence my code is fully vectorized, however : - some parts of my code cannot be vectorized (optimization loops, etc) - some other parts are really not natural when written in a vectorized way. Thus writting, reading, debugging and maintaining them is a pain. But basically, I just would like to see the speed gain I can achieve with a lower level language on some parts of my code. Then depending on the performances of the forthcoming JIT compiler of Scilab, I may consider to translate (or not) my code in fortran. Pierre Vuillemin Le mercredi 05 f?vrier 2014 ? 12:38 +0000, Carrico, Paul a ?crit : > All > > before using fortan code, did you have a look to the "vectorization" in order to avoid loops (that drastically decrease speed) and to increase speedup consequently ? > > Just an advice > > Paul > > -----Message d'origine----- > De : users [mailto:users-bounces at lists.scilab.org] De la part de Pierre Vuillemin > Envoy? : mercredi 5 f?vrier 2014 09:56 > ? : users at lists.scilab.org > Objet : [Scilab-users] Interfacing fortran code with scilab > > Hello all, > > I am currently learning fortran in order to speed up some part o my code and I would like to interface it with Scilab. > > To make some tests, I have used this fortran code : > > subroutine matmul(C,A,B,m,n,p) > integer ::m,n,p,i,j,k > real(kind=8),dimension(m,n),intent(in)::A > real(kind=8),dimension(n,p),intent(in)::B > real(kind=8),dimension(m,p),intent(out)::C > real(kind=8):: temp > C = 0 > do j =1,p > do k=1,n > temp = B(k,j) > do i=1,m > C(i,j)=C(i,j) + A(i,k)*temp > end do > end do > end do > end subroutine matmul > > which performs a matrix multiplication. Then I have used 'ilib_for_link' > to create the link. > It works but I have some questions : > > - I did the same thing in Python with f2py and the code is seemingly faster in python, I was wondering why ? (the precision is the same) > > - If I change the 'kind = 8' to 'kind = 4' in the fortran code, the function returns a result which completely wrong. > I expected it to be wrong, but not that much : for instance, with A = 10 and B = 20, the call to the fortran function gives me '6.36D-314' > What is the reason? (My comprehension of how fortran types work with the kind parameter is still limited, I know that kind=8 and kind=4 is not a very portable expression though) > > - If I replace > real(kind=8),dimension(m,p),intent(out)::C > by > real(kind=8),dimension(:,:),intent(out)::C > I get a segfault in Scilab or it just closes. I was wondering why ? (I thought it was something valid since gfortran still compiles) > > - Is it possible to keep the upper case in the name of fortran functions when linked to scilab ? (At the beginning, my function was called 'matMul' and it took me some time to figure out what the error message 'matMul is not an entry point' meant) > > Best regards, > > Pierre Vuillemin > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v1/url?u=http://lists.scilab.org/mailman/listinfo/users&k=b2vlTQszY8VIpYRvaG%2By2A%3D%3D%0A&r=SzPQe21KEY%2BPdUhJge5w%2FdtbtLgIXw5YTsnbzx%2F7JYE%3D%0A&m=ULVia0lWcxnisV%2BYAJ894mfD5PBI7UEZDeqeSqE9PQw%3D%0A&s=4501ea13ac011e234cb1dd1cdf1810c4388736b832cf43301a84c41f80b767f5 > -------------------------------------------------------------------------------- > > > Le pr?sent mail et ses pi?ces jointes sont confidentiels et destin?s ? la personne ou aux personnes vis?e(s) ci-dessus. Si vous avez re?u cet e-mail par erreur, veuillez contacter imm?diatement l'exp?diteur et effacer le message de votre syst?me. Toute divulgation, copie ou distribution de cet e-mail est strictement interdite. > > This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error, please contact the sender and delete the email from your system. If you are not the named addressee you should not disseminate, distribute or copy this email. > From jdheinzmann at gmail.com Wed Feb 5 17:52:55 2014 From: jdheinzmann at gmail.com (John David Heinzmann) Date: Wed, 5 Feb 2014 11:52:55 -0500 Subject: [Scilab-users] Trying to use old Scicos file with Xcos In-Reply-To: <1391506993.5710.5.camel@paros> References: <1391382524436-4028498.post@n3.nabble.com> <1391506993.5710.5.camel@paros> Message-ID: Hi, Cl?ment, I am new here, so would you tell me the best way to provide my model? Simply as an email attachment? .zcos format to save space? Thanks. JD On Tue, Feb 4, 2014 at 4:43 AM, Cl?ment David < clement.david at scilab-enterprises.com> wrote: > Hello, > > It should work fine but it is model related. Can you provide your > model ? > > -- > Cl?ment > > Le dimanche 02 f?vrier 2014 ? 15:08 -0800, JD Heinzmann a ?crit : > > About 3 years ago I developed a Scicos simulation of a heat engine using > > ScicosLab-4.4b8 with Scicos 4.4. I have just downloaded Scilab 5.4.1 and > > would like to run this old simulation under Xcos. > > > > I can open the old .cos file in Xcos and see the blocks and navigate into > > the super-block hierarchy just fine. I can save this block diagram as > .xcos > > and .zcos files just fine. I can load the context file into Scilab and > run > > it to set up the workspace before running the simulation. But when I > try to > > run the simulation under Xcos, Scilab and Xcos immediately quits silently > > (with no error messages). > > > > Are Scicos files not compatible with Xcos? If not, what do I need to do > to > > adjust my model to work with Xcos? > > > > Thank you! > > > > JD > > > > > > > > -- > > View this message in context: > http://mailinglists.scilab.org/Trying-to-use-old-Scicos-file-with-Xcos-tp4028498.html > > Sent from the Scilab users - Mailing Lists Archives mailing list archive > at Nabble.com. > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > http://lists.scilab.org/mailman/listinfo/users > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From kvinith7 at gmail.com Wed Feb 5 09:23:30 2014 From: kvinith7 at gmail.com (Vinith Kumar) Date: Wed, 5 Feb 2014 13:53:30 +0530 Subject: [Scilab-users] regarding the need of help manual Message-ID: Iam doing my final year under graduate project in interfacing Arduino with scilab The help manual is available in French language. I want the soft copy of help manual for interfacing scilab with Arduino in English language. So that i and my project mates can understand the Arduino Atom interfacing with scilab and it will be useful for my project also. with regards, [1]. J.HARIHARASUDHAN [2]. K.DINESH KUMAR [3]. B.VINITH KUMAR Department of Electronics and Communication Engineering, Institute of Road and Transport Technology, Erode Tamilnadu India -------------- next part -------------- An HTML attachment was scrubbed... URL: From ebmb.sci at gmail.com Wed Feb 5 11:24:41 2014 From: ebmb.sci at gmail.com (ebmb.sci) Date: Wed, 5 Feb 2014 02:24:41 -0800 (PST) Subject: [Scilab-users] Genectic algorithm question In-Reply-To: <52F0C8E1.8010903@scilab-enterprises.com> References: <52F0C8E1.8010903@scilab-enterprises.com> Message-ID: <1391595881063-4028511.post@n3.nabble.com> Dear Paul, first of all, thanks for your answer... The version of my Scilab is 5.4.1 1364571296, under MS Windows 8 6.2 64bits. I will contact you by email to send a test case, ok?! Thanks a lot for your attention! Cheers, Eduardo. -- View this message in context: http://mailinglists.scilab.org/Scilab-users-Genectic-algorithm-question-tp4028507p4028511.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Thu Feb 6 08:57:59 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Thu, 6 Feb 2014 08:57:59 +0100 Subject: [Scilab-users] regarding the need of help manual In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3264C90AF@301EX00100.sidel.com> Hello, > De la part de Vinith Kumar > Envoy? : mercredi 5 f?vrier 2014 09:24 > > I want the soft copy of help > manual for interfacing scilab with Arduino in English language. What do you mean by "soft copy"? A PDF? You normally have the help in english provided with the software. You can start the help navigator just by typing help() in the console, or using the dedicated button (?). The help is available at http://help.scilab.org/ you also have a wiki in english at http://wiki.scilab.org/ Finally, you probably have the help files intalled on your computer. In my case (Windows XP), you have the .chm files at ~/modules/helptools where ~ is the installation folder. Hope this helps. -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From ddu at bluewin.ch Thu Feb 6 09:58:18 2014 From: ddu at bluewin.ch (Skydive) Date: Thu, 6 Feb 2014 00:58:18 -0800 (PST) Subject: [Scilab-users] uicontrol not printed out in xs2pdf (and all) Message-ID: <1391677098235-4028520.post@n3.nabble.com> Hello I have some issues with the uicontrol. When i want to export a figure with axes and uicontrol's i get only the axes out of it. here an example: ///------------------------------------------- clc clear close(gcf()) f=gcf(); plot2d(); t1=uicontrol(f,'style',"text", 'string', "hello", 'position', [10 100 100 50]); filename='test'; xs2pdf(f, filename) ///------------------------------------------- I'm using scilab5.4.1 on an xp 32bit machine I found nowhere some answer... Thanks for your help in advance greets domi -- View this message in context: http://mailinglists.scilab.org/Scilab-users-uicontrol-not-printed-out-in-xs2pdf-and-all-tp4028520.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From boljeleel at yahoo.com Thu Feb 6 12:41:06 2014 From: boljeleel at yahoo.com (Debola Abduljeleel) Date: Thu, 6 Feb 2014 03:41:06 -0800 (PST) Subject: [Scilab-users] Help on how to plot double y-axes on thesame x-axis using Scilab! Message-ID: <1391686866.34021.YahooMailNeo@web160304.mail.bf1.yahoo.com> Hi everyone, Please how can I plot a graph with 2 y-axes i.e. Left & right with just one single x-axis. Regards. ? ABDULJELEEL A. OSUNKUNLE SCIENTIFIC OFFICER PROJECTS DEVELOPMENT INSTITUTE, P.M.B. 01609, EMENE INDUSTRIAL LAYOUT, ENUGU, NIGERIA. +2347058408440. -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Thu Feb 6 13:05:25 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Thu, 6 Feb 2014 13:05:25 +0100 Subject: [Scilab-users] Help on how to plot double y-axes on thesame x-axis using Scilab! In-Reply-To: <1391686866.34021.YahooMailNeo@web160304.mail.bf1.yahoo.com> References: <1391686866.34021.YahooMailNeo@web160304.mail.bf1.yahoo.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3264C9765@301EX00100.sidel.com> Hello, > De la part de Debola Abduljeleel > Envoy? : jeudi 6 f?vrier 2014 12:41 > > Please how can I plot a graph with 2 y-axes i.e. Left & right with > just one single x-axis. We already answered you in last september, please see http://mailinglists.scilab.org/Scilab-users-Double-y-axis-on-a-single-plot-td4027480.html I personally invited you to have a look at https://commons.wikimedia.org/wiki/File:Trace_ln_sqrt_1_2_deux_echelles_scilab.svg?uselang=en Best regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From mariajovera at icloud.com Fri Feb 7 23:39:33 2014 From: mariajovera at icloud.com (samaelkreutz) Date: Fri, 7 Feb 2014 14:39:33 -0800 (PST) Subject: [Scilab-users] Loading multiple text files Message-ID: <1391812772237-4028528.post@n3.nabble.com> Hello, I'm newbie in scilab... I have one folder /Users/konstantina/Desktop/CHECK/export_5-20/export_5-20HZ_C*.txt where * takes the values= 1:1:5 I need open each text file and make a plot. So, I was thinking in something like for i=1:1:4// DM(i)=fscanfMat('export_5-20HZ_C(i).txt') figure(i) plot(DM(i)(:,7),DM(i)(:,4),'ro*',DM(i)(:,7),DM(i)(:,5),'o') end but is obvious that it doesn't work, hehehe. Please, please any suggestions!!! Ill be very thankful :) Best wishes!!! -- View this message in context: http://mailinglists.scilab.org/Loading-multiple-text-files-tp4028528.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Sat Feb 8 01:29:06 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sat, 08 Feb 2014 01:29:06 +0100 Subject: [Scilab-users] Loading multiple text files In-Reply-To: <1391812772237-4028528.post@n3.nabble.com> References: <1391812772237-4028528.post@n3.nabble.com> Message-ID: <52F57A52.3040701@free.fr> Le 07/02/2014 23:39, samaelkreutz a ?crit : > for i=1:1:4// > DM(i)=fscanfMat('export_5-20HZ_C(i).txt') > figure(i) > plot(DM(i)(:,7),DM(i)(:,4),'ro*',DM(i)(:,7),DM(i)(:,5),'o') > end for i=1:1:4// DM=fscanfMat('export_5-20HZ_C(i).txt') figure(i) plot(DM(:,7),DM(:,4),'ro*',DM(:,7),DM(:,5),'o') end may work From mariajovera at icloud.com Sat Feb 8 02:21:40 2014 From: mariajovera at icloud.com (samaelkreutz) Date: Fri, 7 Feb 2014 17:21:40 -0800 (PST) Subject: [Scilab-users] Loading multiple text files In-Reply-To: <52F57A52.3040701@free.fr> References: <1391812772237-4028528.post@n3.nabble.com> <52F57A52.3040701@free.fr> Message-ID: <1391822500814-4028534.post@n3.nabble.com> :( nop, it doesn't work seems the program doesn't recognize the variable (i) -- View this message in context: http://mailinglists.scilab.org/Loading-multiple-text-files-tp4028528p4028534.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.carrico at free.fr Sat Feb 8 08:52:36 2014 From: paul.carrico at free.fr (Paul CARRICO) Date: Sat, 8 Feb 2014 08:52:36 +0100 Subject: [Scilab-users] Loading multiple text files In-Reply-To: <1391812772237-4028528.post@n3.nabble.com> References: <1391812772237-4028528.post@n3.nabble.com> Message-ID: <003001cf24a2$bc578d30$3506a790$@carrico@free.fr> I haven't tested specifically your example, but I would use something like : for i=1:1:4 str1 = "... DM("+ string(i) +")=fscanfMat(''export_5-20HZ_C("+ string(i) +").txt''),... figure("+ string(i) +"),... plot(DM("+ string(i) +")(:,7),DM("+ string(i) +")(:,4),''ro*'',DM("+ string(i) +")(:,7),DM("+ string(i) +")(:,5),''o''),... "; execstr(str1) ; end Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de samaelkreutz Envoy??: vendredi 7 f?vrier 2014 23:40 ??: users at lists.scilab.org Objet?: [Scilab-users] Loading multiple text files Hello, I'm newbie in scilab... I have one folder /Users/konstantina/Desktop/CHECK/export_5-20/export_5-20HZ_C*.txt where * takes the values= 1:1:5 I need open each text file and make a plot. So, I was thinking in something like for i=1:1:4// DM(i)=fscanfMat('export_5-20HZ_C(i).txt') figure(i) plot(DM(i)(:,7),DM(i)(:,4),'ro*',DM(i)(:,7),DM(i)(:,5),'o') end but is obvious that it doesn't work, hehehe. Please, please any suggestions!!! Ill be very thankful :) Best wishes!!! -- View this message in context: http://mailinglists.scilab.org/Loading-multiple-text-files-tp4028528.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com From paul.carrico at free.fr Sat Feb 8 09:06:23 2014 From: paul.carrico at free.fr (Paul CARRICO) Date: Sat, 8 Feb 2014 09:06:23 +0100 Subject: [Scilab-users] Loading multiple text files In-Reply-To: <003001cf24a2$bc578d30$3506a790$@carrico@free.fr> References: <1391812772237-4028528.post@n3.nabble.com> <003001cf24a2$bc578d30$3506a790$@carrico@free.fr> Message-ID: <003401cf24a4$a8edffc0$fac9ff40$@carrico@free.fr> Note the "plot" line has been cut in 2 ones by outlook !!!! -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Paul CARRICO Envoy??: samedi 8 f?vrier 2014 08:53 ??: 'International users mailing list for Scilab.' Objet?: Re: [Scilab-users] Loading multiple text files I haven't tested specifically your example, but I would use something like : for i=1:1:4 str1 = "... DM("+ string(i) +")=fscanfMat(''export_5-20HZ_C("+ string(i) +").txt''),... figure("+ string(i) +"),... plot(DM("+ string(i) +")(:,7),DM("+ string(i) +")(:,4),''ro*'',DM("+ string(i) +")(:,7),DM("+ string(i) +")(:,5),''o''),... "; execstr(str1) ; end Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de samaelkreutz Envoy??: vendredi 7 f?vrier 2014 23:40 ??: users at lists.scilab.org Objet?: [Scilab-users] Loading multiple text files Hello, I'm newbie in scilab... I have one folder /Users/konstantina/Desktop/CHECK/export_5-20/export_5-20HZ_C*.txt where * takes the values= 1:1:5 I need open each text file and make a plot. So, I was thinking in something like for i=1:1:4// DM(i)=fscanfMat('export_5-20HZ_C(i).txt') figure(i) plot(DM(i)(:,7),DM(i)(:,4),'ro*',DM(i)(:,7),DM(i)(:,5),'o') end but is obvious that it doesn't work, hehehe. Please, please any suggestions!!! Ill be very thankful :) Best wishes!!! -- View this message in context: http://mailinglists.scilab.org/Loading-multiple-text-files-tp4028528.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com From sgougeon at free.fr Sat Feb 8 13:14:12 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sat, 08 Feb 2014 13:14:12 +0100 Subject: [Scilab-users] Loading multiple text files In-Reply-To: <1391822500814-4028534.post@n3.nabble.com> References: <1391812772237-4028528.post@n3.nabble.com> <52F57A52.3040701@free.fr> <1391822500814-4028534.post@n3.nabble.com> Message-ID: <52F61F94.5020101@free.fr> Le 08/02/2014 02:21, samaelkreutz a ?crit : > :( nop, it doesn't work > > Yes, i did not pay attention to the management of files names in your script. How are numbered the names of your files? Assuming that they are: export_5-20HZ_C1.txt, export_5-20HZ_C2.txt, etc, then for i=1:1:4// DM = fscanfMat(msprintf("export_5-20HZ_C%d.txt",i)); figure(i); plot(DM(:,7), DM(:,4), 'ro*', DM(:,7), DM(:,5), 'o') end should work, provided that files are set in the working directory (run /pwd, dir/ to see it and them) and are well formatted to be accepted by fscanfMat(). Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From mariajovera at icloud.com Sat Feb 8 16:23:10 2014 From: mariajovera at icloud.com (samaelkreutz) Date: Sat, 8 Feb 2014 07:23:10 -0800 (PST) Subject: [Scilab-users] Loading multiple text files In-Reply-To: <52F61F94.5020101@free.fr> References: <1391812772237-4028528.post@n3.nabble.com> <52F57A52.3040701@free.fr> <1391822500814-4028534.post@n3.nabble.com> <52F61F94.5020101@free.fr> Message-ID: <1391872990998-4028546.post@n3.nabble.com> *THANKS ALOT!!! NOW ITWORKS :) :) * -- View this message in context: http://mailinglists.scilab.org/Loading-multiple-text-files-tp4028528p4028546.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From limipl2 at orange.fr Fri Feb 7 10:01:11 2014 From: limipl2 at orange.fr (limipl) Date: Fri, 7 Feb 2014 01:01:11 -0800 (PST) Subject: [Scilab-users] graph format Message-ID: <1391763671323-4028525.post@n3.nabble.com> Hi, I save graphs from scilab on my desktop but the file has no extension and i'm now unable to open it. Is there a solution ? Thanks -- View this message in context: http://mailinglists.scilab.org/graph-format-tp4028525.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From noelec57 at yahoo.fr Sun Feb 9 16:14:09 2014 From: noelec57 at yahoo.fr (noelec57) Date: Sun, 9 Feb 2014 07:14:09 -0800 (PST) Subject: [Scilab-users] Scilab vs Xcos Message-ID: <1391958848959-4028555.post@n3.nabble.com> Bonjour Je rejoins en tant qu'enseignant les utilisateurs de Scilab et Xcos. Pour le moment j'essaie de l'utiliser pour la mod?lisation de syst?mes asservis continus avec PID. Je suis parvenu ? obtenir la r?ponse de mon syst?me en utilisant Xcos (mode graphique) et scilab (ligne de code). Malheureusement le m?me syst?me ?tudi? avec la m?me fonction de transfert ne fournit pas exactement la m?me r?ponse. Par exemple, j'ai cherch? le point d'instabilit? avec la m?thode Ziegler Nichols. j'obtiens deux valeurs diff?rentes de gain critique. Idem, l'application d'un PID avec les m?mes valeurs ne fournit pas la m?e r?ponse( variation d'environ 10% !!). Je ne parviens pas ? trouver la raison de cette diff?rence. A titre d'exemple, avec la fonction de transfert en BO suivante : 60 / (1+8*s) * (1+2*s) * (1+s) le syst?me en BF ? retour unitaire provoque des oscillations stables dans un cas avec un gain critique de 0,281 sous scilab et 0,287 sous xcos. Peut ?tre que quelqu'un peut me dire, lequel des deux r?sultats je dois prendre et surtout pourquoi cet ?cart ? Merci. -- View this message in context: http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From akhorshidi at live.com Mon Feb 10 17:00:05 2014 From: akhorshidi at live.com (A Khorshidi) Date: Mon, 10 Feb 2014 19:30:05 +0330 Subject: [Scilab-users] Modeling a damped driven pendulum by using Coselica blocks Message-ID: Hi; Is it possible to model a damped-nonlinear pendulum under external excitation by using Coselica blocks? Suppose we have a simple pendulum with a bob of mass "m" and a weightless rod of length "L". The following 2nd-order ODE describes the motion of the pendulum: "$ \ddot\theta + \frac{c}{mL^{2}}\dot{\theta}(t) + \frac{g}{L} \sin {\theta(t)} = \frac{\tau (t)}{mL^{2}} $" g: the gravitational acceleration (9.81 m/s^2) ct: the (rotational) damping constant (due to air resistance and friction of the rod at its pivot point) tau(t): the external torque applied to the pendulum The attached files show how we can use Scilab ode function and Xcos blocks to solve the pendulum's ODE. So my question is, How can I solve this problem by using Modelica "standard"blocks? I'm debating whether or not I should write my own Modelica blocks. Merci, Mehran _ -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: equation.gif Type: image/gif Size: 1260 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Pendulum_xcos.xcos Type: application/octet-stream Size: 85653 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Pendulum_ode.sci Type: application/octet-stream Size: 551 bytes Desc: not available URL: From Christophe.Dang at sidel.com Tue Feb 11 09:20:30 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Tue, 11 Feb 2014 09:20:30 +0100 Subject: [Scilab-users] graph format In-Reply-To: <1391763671323-4028525.post@n3.nabble.com> References: <1391763671323-4028525.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3265AD6BA@301EX00100.sidel.com> Hello, > De la part de limipl > Envoy? : vendredi 7 f?vrier 2014 10:01 > > I save graphs from scilab on my desktop but the file has no extension > and i'm now unable to open it. > Is there a solution ? Mosts OS use the extension to determine the type of file. Maybe you could just add the correct extension. Now the question is: how do you save your graph? -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From limipl2 at orange.fr Tue Feb 11 16:30:49 2014 From: limipl2 at orange.fr (limipl) Date: Tue, 11 Feb 2014 07:30:49 -0800 (PST) Subject: [Scilab-users] graph format In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F3265AD6BA@301EX00100.sidel.com> References: <1391763671323-4028525.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F3265AD6BA@301EX00100.sidel.com> Message-ID: <1392132649563-4028580.post@n3.nabble.com> Hi, That was the question... but I've successfully open the file with xload command. Thanks Dang -- View this message in context: http://mailinglists.scilab.org/graph-format-tp4028525p4028580.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Wed Feb 12 11:05:46 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 12 Feb 2014 11:05:46 +0100 Subject: [Scilab-users] graph format In-Reply-To: <1392132649563-4028580.post@n3.nabble.com> References: <1391763671323-4028525.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F3265AD6BA@301EX00100.sidel.com> <1392132649563-4028580.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3265FACE4@301EX00100.sidel.com> Hello, > De la part de limipl > Envoy? : mardi 11 f?vrier 2014 16:31 > > but I've successfully open the file with xload command. OK, so you probably saved the graphics with the xsave() command. Unfortunately, the format is not described in the help pages. Probably a specific format. If you want to process your graphics with another application, you should use cammands like xs2png() (raster) or xs2svg() (vector graphics). see e.g. http://help.scilab.org/docs/5.4.1/en_US/xs2svg.html Hope this helps. Regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From julian.estevez at ehu.es Wed Feb 12 12:49:56 2014 From: julian.estevez at ehu.es (Jeibros) Date: Wed, 12 Feb 2014 03:49:56 -0800 (PST) Subject: [Scilab-users] 3d curve plot Message-ID: <1392205796285-4028615.post@n3.nabble.com> Hi all, I'm totally new to Scilab (I used to work with Matlab), and I have a rather simple question. I have two arrays (x,y). Next, I make a conversion so that to get a new variable u out of x and y. And finally, I have a z,u relationship. *I would like to make a x,y,z plot*, and I don't need a surface, but a line. Which is the solution please? imagine that the code I have is this: x=0:0.1:2; y=0:0.2:4; u=x+y; z=3*u + 8; plot3(x,y,z) //this is obviously wrong Thanks a lot and sorry for the simple doubt -- View this message in context: http://mailinglists.scilab.org/3d-curve-plot-tp4028615.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Serge.Steer at inria.fr Thu Feb 13 10:51:25 2014 From: Serge.Steer at inria.fr (Serge Steer) Date: Thu, 13 Feb 2014 10:51:25 +0100 Subject: [Scilab-users] 3d curve plot In-Reply-To: <1392205796285-4028615.post@n3.nabble.com> References: <1392205796285-4028615.post@n3.nabble.com> Message-ID: <52FC959D.60309@inria.fr> Le 12/02/2014 12:49, Jeibros a ?crit : > Hi all, > > I'm totally new to Scilab (I used to work with Matlab), and I have a rather > simple question. I have two arrays (x,y). Next, I make a conversion so that > to get a new variable u out of x and y. And finally, I have a z,u > relationship. > > *I would like to make a x,y,z plot*, and I don't need a surface, but a line. > Which is the solution please? > > imagine that the code I have is this: > > x=0:0.1:2; > y=0:0.2:4; > u=x+y; > z=3*u + 8; > plot3(x,y,z) //this is obviously wrong Why? > > Thanks a lot and sorry for the simple doubt > > > > > -- > View this message in context: http://mailinglists.scilab.org/3d-curve-plot-tp4028615.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From stephane.mottelet at utc.fr Thu Feb 13 10:56:49 2014 From: stephane.mottelet at utc.fr (=?ISO-8859-1?Q?St=E9phane_Mottelet?=) Date: Thu, 13 Feb 2014 10:56:49 +0100 Subject: [Scilab-users] 3d curve plot In-Reply-To: <52FC959D.60309@inria.fr> References: <1392205796285-4028615.post@n3.nabble.com> <52FC959D.60309@inria.fr> Message-ID: <52FC96E1.6010105@utc.fr> Le 13/02/2014 10:51, Serge Steer a ?crit : > Le 12/02/2014 12:49, Jeibros a ?crit : >> Hi all, >> >> I'm totally new to Scilab (I used to work with Matlab), and I have a >> rather >> simple question. I have two arrays (x,y). Next, I make a conversion >> so that >> to get a new variable u out of x and y. And finally, I have a z,u >> relationship. >> >> *I would like to make a x,y,z plot*, and I don't need a surface, but >> a line. >> Which is the solution please? >> >> imagine that the code I have is this: >> >> x=0:0.1:2; >> y=0:0.2:4; >> u=x+y; >> z=3*u + 8; >> plot3(x,y,z) //this is obviously wrong > Why? For a parametric curve in 3D, the scilab macro is "param3d". S. >> >> Thanks a lot and sorry for the simple doubt >> >> >> >> >> -- >> View this message in context: >> http://mailinglists.scilab.org/3d-curve-plot-tp4028615.html >> Sent from the Scilab users - Mailing Lists Archives mailing list >> archive at Nabble.com. >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From martin at null-a.de Thu Feb 13 10:52:53 2014 From: martin at null-a.de (Martin Helm) Date: Thu, 13 Feb 2014 10:52:53 +0100 Subject: [Scilab-users] 3d curve plot In-Reply-To: <1392205796285-4028615.post@n3.nabble.com> References: <1392205796285-4028615.post@n3.nabble.com> Message-ID: <52FC95F5.60607@null-a.de> Am 12.02.2014 12:49, schrieb Jeibros: > plot3(x,y,z) //this is obviously wrong use param3d(x,y,z) From paul.bignier at scilab-enterprises.com Thu Feb 13 13:44:19 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Thu, 13 Feb 2014 13:44:19 +0100 Subject: [Scilab-users] Scilab vs Xcos In-Reply-To: <1391958848959-4028555.post@n3.nabble.com> References: <1391958848959-4028555.post@n3.nabble.com> Message-ID: <52FCBE23.8060708@scilab-enterprises.com> Bonjour, Pourriez-vous svp nous envoyer en pi?ce jointe vos diagramme Xcos et script Scilab, pour que nous puissions reproduire l'erreur ? Merci, Bonne journ?e, Paul On 02/09/2014 04:14 PM, noelec57 wrote: > Bonjour > > Je rejoins en tant qu'enseignant les utilisateurs de Scilab et Xcos. Pour le > moment j'essaie de l'utiliser pour la mod?lisation de syst?mes asservis > continus avec PID. Je suis parvenu ? obtenir la r?ponse de mon syst?me en > utilisant Xcos (mode graphique) et scilab (ligne de code). Malheureusement > le m?me syst?me ?tudi? avec la m?me fonction de transfert ne fournit pas > exactement la m?me r?ponse. Par exemple, j'ai cherch? le point d'instabilit? > avec la m?thode Ziegler Nichols. j'obtiens deux valeurs diff?rentes de gain > critique. Idem, l'application d'un PID avec les m?mes valeurs ne fournit pas > la m?e r?ponse( variation d'environ 10% !!). > > Je ne parviens pas ? trouver la raison de cette diff?rence. > > A titre d'exemple, avec la fonction de transfert en BO suivante : 60 / > (1+8*s) * (1+2*s) * (1+s) > > le syst?me en BF ? retour unitaire provoque des oscillations stables dans > un cas avec un gain critique de 0,281 sous scilab et 0,287 sous xcos. > > Peut ?tre que quelqu'un peut me dire, lequel des deux r?sultats je dois > prendre et surtout pourquoi cet ?cart ? > > Merci. > > > > -- > View this message in context: http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Scilab Engineer & Xcos Developer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com From antoine.monmayrant at laas.fr Thu Feb 13 14:07:49 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Thu, 13 Feb 2014 14:07:49 +0100 Subject: [Scilab-users] Scilab vs Xcos In-Reply-To: <52FCBE23.8060708@scilab-enterprises.com> References: <1391958848959-4028555.post@n3.nabble.com> <52FCBE23.8060708@scilab-enterprises.com> Message-ID: <52FCC3A5.7000805@laas.fr> Bonjour Paul et noelec57, Juste une remarque en passant, cette liste de mail est destin?e ? tous les utilisateurs de Scilab (utilisateurs 'internationaux') et l'anglais est normalement de rigueur. Vous aurez beaucoup plus de chance d'avoir des r?ponses ? vos questions si vous pouvez les poser en anglais. Il existe une liste francophone ici: users-fr at lists.scilab.org . Je ne sais pas si elle est aussi active que la liste internationale... Cordialement, Antoine ----- Hi Paul and noelec57 (???), Just a remark, this user list is supposed to be in English. If you post your questions here, you might get way more answers if you do it in English... There is a French user list here: users-fr at lists.scilab.org . I don't know whether the French list is as active as the English one... Cheers, Antoine ----- On 02/13/2014 01:44 PM, Paul Bignier wrote: > > Bonjour, > > Pourriez-vous svp nous envoyer en pi?ce jointe vos diagramme Xcos et > script Scilab, pour que nous puissions reproduire l'erreur ? > > Merci, > Bonne journ?e, > Paul > > > On 02/09/2014 04:14 PM, noelec57 wrote: >> Bonjour >> >> Je rejoins en tant qu'enseignant les utilisateurs de Scilab et Xcos. >> Pour le >> moment j'essaie de l'utiliser pour la mod?lisation de syst?mes asservis >> continus avec PID. Je suis parvenu ? obtenir la r?ponse de mon >> syst?me en >> utilisant Xcos (mode graphique) et scilab (ligne de code). >> Malheureusement >> le m?me syst?me ?tudi? avec la m?me fonction de transfert ne fournit pas >> exactement la m?me r?ponse. Par exemple, j'ai cherch? le point >> d'instabilit? >> avec la m?thode Ziegler Nichols. j'obtiens deux valeurs diff?rentes >> de gain >> critique. Idem, l'application d'un PID avec les m?mes valeurs ne >> fournit pas >> la m?e r?ponse( variation d'environ 10% !!). >> >> Je ne parviens pas ? trouver la raison de cette diff?rence. >> >> A titre d'exemple, avec la fonction de transfert en BO suivante : 60 / >> (1+8*s) * (1+2*s) * (1+s) >> >> le syst?me en BF ? retour unitaire provoque des oscillations stables >> dans >> un cas avec un gain critique de 0,281 sous scilab et 0,287 sous xcos. >> >> Peut ?tre que quelqu'un peut me dire, lequel des deux r?sultats je dois >> prendre et surtout pourquoi cet ?cart ? >> >> Merci. >> >> >> >> -- >> View this message in context: >> http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555.html >> Sent from the Scilab users - Mailing Lists Archives mailing list >> archive at Nabble.com. >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > -- +++++++++++++++++++++++++++++++++++++++++++++++++++++++ Antoine Monmayrant LAAS - CNRS 7 avenue du Colonel Roche BP 54200 31031 TOULOUSE Cedex 4 FRANCE Tel:+33 5 61 33 64 59 email : antoine.monmayrant at laas.fr permanent email : antoine.monmayrant at polytechnique.org +++++++++++++++++++++++++++++++++++++++++++++++++++++++ -------------- next part -------------- An HTML attachment was scrubbed... URL: From Serge.Steer at inria.fr Fri Feb 14 14:38:22 2014 From: Serge.Steer at inria.fr (Serge Steer) Date: Fri, 14 Feb 2014 14:38:22 +0100 Subject: [Scilab-users] Scilab vs Xcos In-Reply-To: <1391958848959-4028555.post@n3.nabble.com> References: <1391958848959-4028555.post@n3.nabble.com> Message-ID: <52FE1C4E.9060607@inria.fr> Le 09/02/2014 16:14, noelec57 a ?crit : > Bonjour > > Je rejoins en tant qu'enseignant les utilisateurs de Scilab et Xcos. Pour le > moment j'essaie de l'utiliser pour la mod?lisation de syst?mes asservis > continus avec PID. Je suis parvenu ? obtenir la r?ponse de mon syst?me en > utilisant Xcos (mode graphique) et scilab (ligne de code). Malheureusement > le m?me syst?me ?tudi? avec la m?me fonction de transfert ne fournit pas > exactement la m?me r?ponse. Par exemple, j'ai cherch? le point d'instabilit? > avec la m?thode Ziegler Nichols. j'obtiens deux valeurs diff?rentes de gain > critique. Idem, l'application d'un PID avec les m?mes valeurs ne fournit pas > la m?e r?ponse( variation d'environ 10% !!). > > Je ne parviens pas ? trouver la raison de cette diff?rence. > > A titre d'exemple, avec la fonction de transfert en BO suivante : 60 / > (1+8*s) * (1+2*s) * (1+s) > > le syst?me en BF ? retour unitaire provoque des oscillations stables dans > un cas avec un gain critique de 0,281 sous scilab et 0,287 sous xcos. Vous ne dites pas comment vous avez obtenu ces deux valeurs. pour calculer la valeur du gain critique vous pouvez utiliser la fonction kpure de Scilab h=60 /((1+8*%s) * (1+2*%s) * (1+%s)); Ki=kpure(h) //->0.28125 //V?rification H=(Ki*h) /.(1.0); roots(H.den) Si sous Xcos vous avez obtenu ce gain par simulation il faut savoir que l'integration num?rique repose sur des tol?rances. Dans xcos ces tol?rances sont par d?faut "tol?rantes". pour am?liorer la pr?cision des r?sultats vous devez sp?cifier des tol?rance absolue et relatives plus strictes (menu: Simulation/Configurer) > > Peut ?tre que quelqu'un peut me dire, lequel des deux r?sultats je dois > prendre et surtout pourquoi cet ?cart ? > > Merci. > > > > -- > View this message in context: http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From antoine.monmayrant at laas.fr Fri Feb 14 17:08:31 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Fri, 14 Feb 2014 17:08:31 +0100 Subject: [Scilab-users] SaveAndExecuteFileAction -> ChangeDirectoryAndSaveAndExecuteFileAction Message-ID: <52FE3F7F.3080509@laas.fr> Hi everyone, I would like to add a new button / new action to Scinotes that behaves in a slightly different way than the 'F5 Save and execute'. Basically, I would like to first change to the directory of the file before executing it: 1) Change to the directory that contains the file active in scinotes (the from which I launch the action). 2) Save the file. 3) Execute the file. Is that possible without dissecting org.scilab.modules.scinotes.jar ? As an alternative, if someone knows a way to guess from within a script file the path to the said script, that would also do the trick. (and no, 'get_absolute_file_path' does not work as the script file is not necessarily opened in scilab) Antoine From rm.engineer84 at gmail.com Sun Feb 16 17:54:10 2014 From: rm.engineer84 at gmail.com (R M) Date: Sun, 16 Feb 2014 11:54:10 -0500 Subject: [Scilab-users] how to install modules scilab 5.1.1 Message-ID: Hi, I have scilab 5.1.1 installed. I want to install some new tool boxes. How would I go about doing the same. Thanks and Regards, RM -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Mon Feb 17 08:56:21 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 17 Feb 2014 08:56:21 +0100 Subject: [Scilab-users] how to install modules scilab 5.1.1 In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32670F02A@301EX00100.sidel.com> Hello, > De la part de R M > Envoy? : dimanche 16 f?vrier 2014 17:54 > > I want to install some new tool boxes. > > How would I go about doing the same. Do you mean Atoms modules? Maybe you could have a look at http://wiki.scilab.org/Modules http://atoms.scilab.org/ Hope this helps regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From communication at scilab-enterprises.com Mon Feb 17 09:16:10 2014 From: communication at scilab-enterprises.com (Scilab Communications) Date: Mon, 17 Feb 2014 09:16:10 +0100 Subject: [Scilab-users] ScilabTEC 2014 program is now online! Message-ID: <5301C54A.1080408@scilab-enterprises.com> Is this email not displaying correctly? View it in your browser scilabtec ScilabTEC 2014 program is now online Two-day conference with: 2 keynote speakers: - Alan Edelman, Professor of Applied Mathematics, MIT? - Philippe Marchal, Deputy Director On Board/Soil Systems and Space Dynamics, CNES and also high quality presentations and innovative applications in various fields with the interventions of INRA, AREVA, Noesis Solutions, Sysfera, Structures, Esterline, Upper Austria University of Applied Sciences, CNES, University of Luxembourg, Enginsoft, CEREMA, Cospas-Sarsat, TBSS.c3e, Inria and GEOLAB-CNRS. Discover the complete program and register on http://www.scilabtec.com/ The first evening, participants may meet over dinner with the organizers during a moment of conviviality and exchange between all. REGISTER NOW (limited seating) scilabtec Communication Department, Scilab Enterprises | communication at scilab-enterprises.com 143bis rue Yves Le Coz - 78000 Versailles | www.scilab-enterprises.com - www.scilab.org Scilabtec2014 -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/png Size: 22295 bytes Desc: http://www.scilabtec.com/ URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/png Size: 16786 bytes Desc: http://scilabtec.com/ URL: From noelec57 at yahoo.fr Sun Feb 16 11:54:23 2014 From: noelec57 at yahoo.fr (noelec57) Date: Sun, 16 Feb 2014 02:54:23 -0800 (PST) Subject: [Scilab-users] Scilab vs Xcos In-Reply-To: <52FE1C4E.9060607@inria.fr> References: <1391958848959-4028555.post@n3.nabble.com> <52FE1C4E.9060607@inria.fr> Message-ID: <1392548017.37558.YahooMailNeo@web172506.mail.ir2.yahoo.com> Bonjour et Merci J'ai pris note de votre message mais je n'ai pas eu encore assez de temps pour m'y pencher s?rieusement.? N?anmoins je pense que je dois commettre un autre type d'erreur due ? ma m?connaissance du logiciel. Car en plus des petits ?carts, j'obtiens parfois des r?sultats tr?s surprenants lorsque j'essaie d'utiliser un correcteur PID. Les r?sultats sous Xcos et Scilab sont trop diff?rents. Je me permets de vous joindre le fruit de ma r?flexion actuelle sous la forme de fichiers xcos et Scilab, ainsi qu'une petite note explicative. Je n'ai nullement l'intention de vous faire travailler ? ma place, mais peut ?tre que vous verrez instantan?ment ce qui ne va pas. Le point de d?part ?tait de faire ?tudier des fonctions de transfert avec correcteurs PID mais je voulais ?pargner aux ?tudiants le code scilab et donc utiliser la version "graphique" en utilisant Xcos et ses sch?mas blocs. Mais actuellement je suis un peu bloqu?. Cordialement Noel DEMISSY Enseignant. ________________________________ De?: Serge Steer [via Scilab / Xcos - Mailing Lists Archives] ??: noelec57 Envoy? le : Vendredi 14 f?vrier 2014 14h39 Objet?: Re: Scilab vs Xcos Le 09/02/2014 16:14, noelec57 a ?crit : > Bonjour > > Je rejoins en tant qu'enseignant les utilisateurs de Scilab et Xcos. Pour le > moment j'essaie de l'utiliser pour la mod?lisation de syst?mes asservis > continus avec PID. Je suis parvenu ? obtenir la r?ponse de mon syst?me en > utilisant Xcos (mode graphique) et scilab (ligne de code). Malheureusement > le m?me syst?me ?tudi? avec la m?me fonction de transfert ne fournit pas > exactement la m?me r?ponse. Par exemple, j'ai cherch? le point d'instabilit? > avec la m?thode Ziegler Nichols. j'obtiens deux valeurs diff?rentes de gain > critique. Idem, l'application d'un PID avec les m?mes valeurs ne fournit pas > la m?e r?ponse( variation d'environ 10% !!). > > Je ne parviens pas ? trouver la raison de cette diff?rence. > > A titre d'exemple, avec la fonction de transfert en BO suivante : ?60 / > (1+8*s) * (1+2*s) * (1+s) > > le syst?me en BF ? retour unitaire ?provoque des oscillations stables dans > un cas avec un gain critique de 0,281 sous scilab et 0,287 sous xcos. Vous ne dites pas comment vous avez obtenu ces deux valeurs. pour calculer la valeur du gain critique vous pouvez utiliser la fonction kpure de Scilab h=60 /((1+8*%s) * (1+2*%s) * (1+%s)); Ki=kpure(h) ?//->0.28125 //V?rification H=(Ki*h) /.(1.0); roots(H.den) Si sous Xcos vous avez obtenu ce gain par simulation il faut savoir que l'integration num?rique repose sur des tol?rances. Dans xcos ces tol?rances sont par d?faut "tol?rantes". pour am?liorer la pr?cision des r?sultats vous devez sp?cifier des tol?rance absolue et relatives plus strictes (menu: Simulation/Configurer) > ? > Peut ?tre que quelqu'un peut me dire, lequel des deux r?sultats je dois > prendre et surtout pourquoi cet ?cart ? > > Merci. > > > > -- > View this message in context: http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > [hidden email] > http://lists.scilab.org/mailman/listinfo/users _______________________________________________ users mailing list [hidden email] http://lists.scilab.org/mailman/listinfo/users ________________________________ If you reply to this email, your message will be added to the discussion below:http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555p4028640.html To unsubscribe from Scilab vs Xcos, click here. NAML Scilab vs Xcos.doc (374K) systeme_C_ZN.xcos (90K) systeme_PID.xcos (187K) -- View this message in context: http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555p4028652.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.bingham at tiscali.co.uk Mon Feb 17 21:51:15 2014 From: john.bingham at tiscali.co.uk (john) Date: Mon, 17 Feb 2014 12:51:15 -0800 (PST) Subject: [Scilab-users] Xcos PID controller problem Message-ID: <1392670275497-4028669.post@n3.nabble.com> Hello, I'm new to Scilab and Xcos so may be making a simple mistake, but I have the following problem. I created an Xcos model of a simple mass-spring-damper with PID control using the data from the Matlab example shown here . The Xcos results are the same as the Matlab example except for the case when all three gains are used as follows Kp=350, Ki=300 and Kd=50 . My result looks like this (red line) and is clearly giving a different response to the Matlab example for the same gain values. I have tried different solvers and time steps, but this doesn't seem to be the problem. I have also modelled the problem in Mathcad 7, Mathcad Prime 3 and using csim in Scilab and all those results agree broadly with the Matlab example. Any advice would be greatly appreciated and I have attached the Xcos file for information. spring-damper-PID.zcos -- View this message in context: http://mailinglists.scilab.org/Xcos-PID-controller-problem-tp4028669.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Serge.Steer at inria.fr Mon Feb 17 17:11:53 2014 From: Serge.Steer at inria.fr (Serge Steer) Date: Mon, 17 Feb 2014 17:11:53 +0100 Subject: [Scilab-users] Fwd: Re: Scilab vs Xcos In-Reply-To: <5301BD62.5050105@scilab-enterprises.com> References: <1392547493.36083.YahooMailNeo@web172502.mail.ir2.yahoo.com> <5301BD62.5050105@scilab-enterprises.com> Message-ID: <530234C9.5070606@inria.fr> Le 17/02/2014 08:42, Paul Bignier a ?crit : > > -------- Original Message -------- > Subject: Re: [Scilab-users] Scilab vs Xcos > Date: Sun, 16 Feb 2014 10:44:53 +0000 (GMT) > From: noel demissy > Reply-To: noel demissy > To: Paul Bignier > > > > Bonjour > > Merci de votre r?ponse. Je vous joins 3 fichiers : scilab, xcos et un > "word" explicatif. > > Comme je l'ai mentionn?, je d?couvre ces logiciels. Initialement j'ai > pr?vu de faire travailler des ?tudiants sur des mod?les de syst?mes > asservis mais je voulais porter le travail fait initialement avec > Scilab sous Xcos car je trouvais plus simple de leur faire manipuler > des blocs sans rentrer dans le code. > Je confirme ce que j'avais dit dans un precedent mail. Le probleme avac la simulation Xcos est les r?glage des tol?rances de l'integrateur avec une tol absolue de 1e-8 et une relative de 1e-9 la gain de 0.281 assure une quasi stabilit?. Je vous attache une version de votre sch?ma parametr?e avec le gain calcul? par kpure Serge Steer > Mes erreurs proviennent peut ?tre sans aucun doute de la > m?connaissance du logiciel et en partie des param?tres d'?tude. > > En vous remerciant d'avance. > > Noel DEMISSY > > ------------------------------------------------------------------------ > *De :* Paul Bignier > *? :* International users mailing list for Scilab. > ; noelec57 at yahoo.fr > *Envoy? le :* Jeudi 13 f?vrier 2014 13h44 > *Objet :* Re: [Scilab-users] Scilab vs Xcos > > > Bonjour, > > Pourriez-vous svp nous envoyer en pi?ce jointe vos diagramme Xcos et > script Scilab, pour que nous puissions reproduire l'erreur ? > > Merci, > Bonne journ?e, > Paul > > > On 02/09/2014 04:14 PM, noelec57 wrote: > > Bonjour > > > > Je rejoins en tant qu'enseignant les utilisateurs de Scilab et Xcos. > Pour le > > moment j'essaie de l'utiliser pour la mod?lisation de syst?mes asservis > > continus avec PID. Je suis parvenu ? obtenir la r?ponse de mon > syst?me en > > utilisant Xcos (mode graphique) et scilab (ligne de code). > Malheureusement > > le m?me syst?me ?tudi? avec la m?me fonction de transfert ne fournit pas > > exactement la m?me r?ponse. Par exemple, j'ai cherch? le point > d'instabilit? > > avec la m?thode Ziegler Nichols. j'obtiens deux valeurs diff?rentes > de gain > > critique. Idem, l'application d'un PID avec les m?mes valeurs ne > fournit pas > > la m?e r?ponse( variation d'environ 10% !!). > > > > Je ne parviens pas ? trouver la raison de cette diff?rence. > > > > A titre d'exemple, avec la fonction de transfert en BO suivante : 60 / > > (1+8*s) * (1+2*s) * (1+s) > > > > le syst?me en BF ? retour unitaire provoque des oscillations > stables dans > > un cas avec un gain critique de 0,281 sous scilab et 0,287 sous xcos. > > > > Peut ?tre que quelqu'un peut me dire, lequel des deux r?sultats je dois > > prendre et surtout pourquoi cet ?cart ? > > > > Merci. > > > > > > > > -- > > View this message in context: > http://mailinglists.scilab.org/Scilab-vs-Xcos-tp4028555.html > > Sent from the Scilab users - Mailing Lists Archives mailing list > archive at Nabble.com. > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > http://lists.scilab.org/mailman/listinfo/users > > -- > Paul BIGNIER > Scilab Engineer & Xcos Developer > ----------------------------------------------------------- > Scilab Enterprises > 143bis rue Yves Le Coz - 78000 Versailles, France > Phone: +33.1.80.77.04.69 > http://www.scilab-enterprises.com > > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: systeme_C_ZN.xcos Type: text/xml Size: 67909 bytes Desc: not available URL: From mobash2007 at gmail.com Tue Feb 18 00:33:35 2014 From: mobash2007 at gmail.com (Mohammed Bashir) Date: Mon, 17 Feb 2014 15:33:35 -0800 Subject: [Scilab-users] inquiry In-Reply-To: References: Message-ID: is scilab equivalent to matlab On Mon, Nov 25, 2013 at 3:06 PM, Mohammed Bashir wrote: > > is scilab equivalent to matlab > -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.bignier at scilab-enterprises.com Tue Feb 18 09:06:40 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Tue, 18 Feb 2014 09:06:40 +0100 Subject: [Scilab-users] inquiry In-Reply-To: References: Message-ID: <53031490.5030609@scilab-enterprises.com> Hi Mohammed, Scilab and Matlab are meant for the same purposes and share many similar features, but their semantic, function names and conventions are not exactly the same. So no, they are not directly equivalent, but a Scilab job is often translatable in Matlab, and vice versa. Check out the Matlab-Scilab equivalents help page for more information. Regards, Paul On 02/18/2014 12:33 AM, Mohammed Bashir wrote: > is scilab equivalent to matlab > > > > On Mon, Nov 25, 2013 at 3:06 PM, Mohammed Bashir > wrote: > > > is scilab equivalent to matlab > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Scilab Engineer & Xcos Developer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From Serge.Steer at inria.fr Tue Feb 18 12:04:49 2014 From: Serge.Steer at inria.fr (Serge Steer) Date: Tue, 18 Feb 2014 12:04:49 +0100 Subject: [Scilab-users] Xcos PID controller problem In-Reply-To: <1392670275497-4028669.post@n3.nabble.com> References: <1392670275497-4028669.post@n3.nabble.com> Message-ID: <53033E51.2070802@inria.fr> The PID controller is an improper transfer function (degree of num>degree of den). In Xcos the derivative term of the PID bloc is computed by an approximation which seems to fail in your context. A solution should be to replace the PID bloc with a transfert function block made proper adding a very small term in s^2 for the denominator (sse the attched diagram) Please report the problem in http://bugzilla.scilab.org/ Serge Steer Le 17/02/2014 21:51, john a ?crit : > Hello, > I'm new to Scilab and Xcos so may be making a simple mistake, but I have the > following problem. I created an Xcos model of a simple mass-spring-damper > with PID control using the data from the Matlab example shown here > > . > > The Xcos results are the same as the Matlab example except for the case when > all three gains are used as follows Kp=350, Ki=300 and Kd=50 . > > My result looks like this (red line) and is clearly giving a different > response to the Matlab example for the same gain values. > > I have tried different solvers and time steps, but this doesn't seem to be > the problem. I have also modelled the problem in Mathcad 7, Mathcad Prime 3 > and using csim in Scilab and all those results agree broadly with the Matlab > example. > > Any advice would be greatly appreciated and I have attached the Xcos file > for information. spring-damper-PID.zcos > > > > > > -- > View this message in context: http://mailinglists.scilab.org/Xcos-PID-controller-problem-tp4028669.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- A non-text attachment was scrubbed... Name: spring-damper-PID.zcos Type: application/octet-stream Size: 8189 bytes Desc: not available URL: From matiasb at gmail.com Tue Feb 18 20:59:10 2014 From: matiasb at gmail.com (Matiasb) Date: Tue, 18 Feb 2014 11:59:10 -0800 (PST) Subject: [Scilab-users] How to get a reguar sampled vector in Scilab Message-ID: <1392753550766-4028712.post@n3.nabble.com> I'm trying to program function (or even better it it already exists) in scilab that calculates a regular timed samples of values. IE: I have a vector 'values' which defines contains the value of a signal at different times. This times are in the vector 'times'. So at time times(N), the signal has value values(N). At the moment the times are not regular, so the variable 'times' and 'values' can look like: times = [0, 2, 6, 8, 14] values= [5, 9, 10, 1, 6] This represents that the signal had value 5 from second 0 to second 2. Value 9 from second 2 to second 6, etc. Therefore, if I want to calculate the signal average value I can not just calculate the average of vector 'values'. This is because for example the signal can be for a long time with the same value, but there will be only one value in the vector. One option is to take the deltaT to calculate the media, but I will also need to perform other calculations:average, etc. Other option is to create a function that given a deltaT, samples the time and values vectors to produce an equally spaced timeVector and corresponding values. For example, with deltaT=2 and the previous vectors, [sampledTime, sampledValues] = regularSample(times, values, 2) sampledTime = [0, 2, 4, 6, 8, 10, 12, 14] sampledValues = [5, 9, 9, 10, 1, 1, 1, 6] This is easy if deltaT is small enough to fit exactly with all the times. If the deltaT is bigger, then the average of values or some aproximation must be done... Is there anything already done in Scilab? How can this function be programmed? Thanks a lot! PS: I don't know if this is the correct forum to post scilab questions, so any pointer would also be useful. -- View this message in context: http://mailinglists.scilab.org/How-to-get-a-reguar-sampled-vector-in-Scilab-tp4028712.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From serge.steer at inria.fr Tue Feb 18 22:26:03 2014 From: serge.steer at inria.fr (Serge Steer) Date: Tue, 18 Feb 2014 22:26:03 +0100 (CET) Subject: [Scilab-users] How to get a reguar sampled vector in Scilab In-Reply-To: <1392753550766-4028712.post@n3.nabble.com> References: <1392753550766-4028712.post@n3.nabble.com> Message-ID: <1213788196.680525.1392758762995.JavaMail.zimbra@inria.fr> I think there is a little mistake in the sampledValues you are expecting: times = [0, 2, 6, 8, 14] values= [5, 9, 10, 1, 6] sampledTime = [0, 2, 4, 6, 8, 10, 12, 14] sampledValues = [5, 9, 9, 10, 1, 1, 1, 6] According to the rule used for the previous values, the last value of sampledValues must be 1 and not 6 Here is a code doing what you expect (but with a different last point) times = [0, 2, 6, 8, 14]; values= [5, 9, 10, 1, 6]; sampletimes=0:14; V=zeros(sampletimes); for k=1:size(times,'*')-1 V(sampletimes>=times(k)&sampletimestimes($))=timesvalues($); clf;plot(sampletimes,V,'-+') //You can also use linear interpolation V1=interpln([times;values],sampletimes); plot(sampletimes,V1,'r-+') //or Splin interpolation V2=interp(sampletimes,times,values,splin(times,values)) plot(sampletimes,V2,'m-*') Serge Steer ----- Mail original ----- > De: "Matiasb" > ?: users at lists.scilab.org > Envoy?: Mardi 18 F?vrier 2014 20:59:10 > Objet: [Scilab-users] How to get a reguar sampled vector in Scilab > > I'm trying to program function (or even better it it already exists) in > scilab that calculates a regular timed samples of values. IE: I have a > vector 'values' which defines contains the value of a signal at different > times. This times are in the vector 'times'. So at time times(N), the signal > has value values(N). At the moment the times are not regular, so the > variable 'times' and 'values' can look like: > > times = [0, 2, 6, 8, 14] > values= [5, 9, 10, 1, 6] > This represents that the signal had value 5 from second 0 to second 2. Value > 9 from second 2 to second 6, etc. Therefore, if I want to calculate the > signal average value I can not just calculate the average of vector > 'values'. This is because for example the signal can be for a long time with > the same value, but there will be only one value in the vector. One option > is to take the deltaT to calculate the media, but I will also need to > perform other calculations:average, etc. > > Other option is to create a function that given a deltaT, samples the time > and values vectors to produce an equally spaced timeVector and corresponding > values. For example, with deltaT=2 and the previous vectors, > > [sampledTime, sampledValues] = regularSample(times, values, 2) > sampledTime = [0, 2, 4, 6, 8, 10, 12, 14] > sampledValues = [5, 9, 9, 10, 1, 1, 1, 6] > > This is easy if deltaT is small enough to fit exactly with all the times. If > the deltaT is bigger, then the average of values or some aproximation must > be done... > > Is there anything already done in Scilab? How can this function be > programmed? > > Thanks a lot! PS: I don't know if this is the correct forum to post scilab > questions, so any pointer would also be useful. > > > > -- > View this message in context: > http://mailinglists.scilab.org/How-to-get-a-reguar-sampled-vector-in-Scilab-tp4028712.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at > Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From Christophe.Dang at sidel.com Wed Feb 19 09:20:05 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 19 Feb 2014 09:20:05 +0100 Subject: [Scilab-users] inquiry In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3267A40F6@301EX00100.sidel.com> Hello, > De la part de Mohammed Bashir > Envoy? : mardi 18 f?vrier 2014 00:34 > > is scilab equivalent to matlab You might have a look at Alan Edelman, Open Source and Traditional Technical Computing, MIT (2010) https://www.scilab.org/content/download/395/2850/file/ScilabTec_Keynote.pdf Best regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From yohann.morille at univ-angers.fr Tue Feb 18 15:06:48 2014 From: yohann.morille at univ-angers.fr (Yohann) Date: Tue, 18 Feb 2014 06:06:48 -0800 (PST) Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq Message-ID: <1392732408401-4028696.post@n3.nabble.com> Hi all !! First of all, thank you to all contributors of Scilab. However, I have a problem to calculate the estimated error on each parameter calculated with leastsq. For example if I have the simple following function to fit: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! function y = fun2fit(x, c) y = c(1)*x*x + c(2)*x + c(3); endfunction function e = myerror(c, x, y) e = fun2fit(x, c) - y; endfunction // Data x = [ 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6]'; y = [5.02 6.08 3.33 -0.93 -0.22 7.83 16.52 15.55 2.67 -11.42 -11.78 5.09 25.25]'; // First guess c0 = [1 1 1]'; y0 = fun2fit(x, c0); // least square [f, copt] = leastsq(list(myerror, x, y), c0) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! I'd like to have an estimated error on each of the 3 element of copt. Any idea ? Thanks ! Y -- View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From ronan_2 at hotmail.com Tue Feb 18 20:12:10 2014 From: ronan_2 at hotmail.com (Ronan_) Date: Tue, 18 Feb 2014 11:12:10 -0800 (PST) Subject: [Scilab-users] Installation of SIVP fo Scilab Message-ID: <1392750730616-4028711.post@n3.nabble.com> I can't execute the builder file of SIVP because "this module need a compileur" -->exec('builder.sce') -->mode(-1); Cr?ation des macros... -- Cr?ation de [SIVPlib] (Macros) -- !--error 10000 tbx_builder_gateway : Ce module requiert un compilateur C mais celui-ci n'a pas pu ?tre trouv?. at line 27 of function tbx_builder_gateway called by : tbx_builder_gateway(toolbox_dir); at line 35 of exec file called by : exec('builder.sce') * on the read file it is written that we need a version of scilab older than the 2.0 version but there are not avalaible anymore, maybe this could pose a problem ? * i have installed Eclipse (with mingw) for the occasion, maybe i can create a link between Eclipse and Scilabe to solve the problem ? * I also installed OPEN CV but i don't know what i have to do with it ? it put itself automatically in the PATH but that's all for SIVP use ? Thank you for your answer, Sorry for my english, Ronan SY -- View this message in context: http://mailinglists.scilab.org/Installation-of-SIVP-fo-Scilab-tp4028711.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From julian.estevez at ehu.es Wed Feb 19 08:20:12 2014 From: julian.estevez at ehu.es (Jeibros) Date: Tue, 18 Feb 2014 23:20:12 -0800 (PST) Subject: [Scilab-users] 3d curve plot In-Reply-To: <52FC96E1.6010105@utc.fr> References: <1392205796285-4028615.post@n3.nabble.com> <52FC959D.60309@inria.fr> <52FC96E1.6010105@utc.fr> Message-ID: <1392794412387-4028731.post@n3.nabble.com> Hi, I tested it and worked perfectly. Thanks a lot, J -- View this message in context: http://mailinglists.scilab.org/3d-curve-plot-tp4028615p4028731.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From antoine.monmayrant at laas.fr Wed Feb 19 10:32:20 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Wed, 19 Feb 2014 10:32:20 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392732408401-4028696.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> Message-ID: <53047A24.90004@laas.fr> On 02/18/2014 03:06 PM, Yohann wrote: > Hi all !! > > First of all, thank you to all contributors of Scilab. > However, I have a problem to calculate the estimated error on each parameter > calculated with leastsq. > > For example if I have the simple following function to fit: > !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! > function y = fun2fit(x, c) > y = c(1)*x*x + c(2)*x + c(3); > endfunction > > function e = myerror(c, x, y) > e = fun2fit(x, c) - y; > endfunction > > // Data > x = [ 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6]'; > y = [5.02 6.08 3.33 -0.93 -0.22 7.83 16.52 15.55 2.67 -11.42 -11.78 > 5.09 25.25]'; > > // First guess > c0 = [1 1 1]'; > y0 = fun2fit(x, c0); > > // least square > [f, copt] = leastsq(list(myerror, x, y), c0) > !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! > > I'd like to have an estimated error on each of the 3 element of copt. > Any idea ? What do you mean? Do you want a confidence interval or any kind of "error bar" for each parameter? Or any way to see how your fit is dependent on each parameter? If it's the case, you might be interested in bootstraping or jackknifing: http://en.wikipedia.org/wiki/Resampling_(statistics) These methods are quite ressource hungry but work quite well to measure the quality and robustness of your fit. Cheers, Antoine > > Thanks ! > Y > > > > > > -- > View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- +++++++++++++++++++++++++++++++++++++++++++++++++++++++ Antoine Monmayrant LAAS - CNRS 7 avenue du Colonel Roche BP 54200 31031 TOULOUSE Cedex 4 FRANCE Tel:+33 5 61 33 64 59 email : antoine.monmayrant at laas.fr permanent email : antoine.monmayrant at polytechnique.org +++++++++++++++++++++++++++++++++++++++++++++++++++++++ From arctica1963 at gmail.com Wed Feb 19 15:15:51 2014 From: arctica1963 at gmail.com (Lester Anderson) Date: Wed, 19 Feb 2014 14:15:51 +0000 Subject: [Scilab-users] Rounding error or precision setting? Message-ID: Hello, Just a quick query. Not been using Scilab for a while, so maybe I have missed something, but is there a built-in precision setting for working with real and imaginary functions? Tried the Euler formula (special case) in Scilab 5.4.1: -->x=cos(%pi)+(%i*sin(%pi)) x = - 1. + 1.225D-16i -->y=%e^(%i*%pi) y = - 1. + 1.225D-16i The real part is fine, but the complex has a measured value (approaching zero). Is there a way to tell Scilab if this gets beyond a limit to round off at zero say? Sure this is an easy one! Cheers Lester -------------- next part -------------- An HTML attachment was scrubbed... URL: From vogt at centre-cired.fr Wed Feb 19 15:24:33 2014 From: vogt at centre-cired.fr (Adrien Vogt-Schilb) Date: Wed, 19 Feb 2014 15:24:33 +0100 Subject: [Scilab-users] Rounding error or precision setting? In-Reply-To: References: Message-ID: <5304BEA1.20007@centre-cired.fr> Hi The function *clean *does that: clean(cos(%pi)+%i*sin(%pi)) Note that the limited precision of floating point numbers is a general issue in computer science. http://en.wikipedia.org/wiki/Floating_point#Representable_numbers.2C_conversion_and_rounding You should be aware of that for instance when testing numbers, as in: if cos(%pi)+%i*sin(%pi) == -1 disp("i am a mathematician") else disp("i am a computer") end if abs(cos(%pi)+%i*sin(%pi) +1)<1d-8 disp("yep") else disp("nope") end or when you calculate indefinite integrals. Also note that some modules on ATOMS and scilabfilexchange allow to make computation with arbitrary (but not infenitesimal) precision. On 19/02/2014 15:15, Lester Anderson wrote: > Hello, > Just a quick query. Not been using Scilab for a while, so maybe I have > missed something, but is there a built-in precision setting for > working with real and imaginary functions? > Tried the Euler formula (special case) in Scilab 5.4.1: > -->x=cos(%pi)+(%i*sin(%pi)) > x = > > - 1. + 1.225D-16i > > -->y=%e^(%i*%pi) > y = > > - 1. + 1.225D-16i > The real part is fine, but the complex has a measured value > (approaching zero). Is there a way to tell Scilab if this gets beyond > a limit to round off at zero say? > Sure this is an easy one! > Cheers > > Lester > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Adrien Vogt-Schilb PhD Student (Cired) -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Wed Feb 19 15:25:41 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 19 Feb 2014 15:25:41 +0100 Subject: [Scilab-users] Rounding error or precision setting? In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3267A4903@301EX00100.sidel.com> Hello, > De la part de Lester Anderson > Envoy? : mercredi 19 f?vrier 2014 15:16 > > Is there a way to tell Scilab if this gets beyond a limit to round off > at zero say? Maybe format() or clean() ? http://help.scilab.org/docs/5.4.1/en_US/format.html http://help.scilab.org/docs/5.4.1/en_US/clean.html Best regards -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From arctica1963 at gmail.com Wed Feb 19 15:30:25 2014 From: arctica1963 at gmail.com (Lester Anderson) Date: Wed, 19 Feb 2014 14:30:25 +0000 Subject: [Scilab-users] Rounding error or precision setting? In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F3267A4903@301EX00100.sidel.com> References: <3B5FFC67498DFF49AE7271A584867D16F3267A4903@301EX00100.sidel.com> Message-ID: Thanks for the pointer - knew it was a simple one! Cheers On 19 February 2014 14:25, Dang, Christophe wrote: > Hello, > > > De la part de Lester Anderson > > Envoy? : mercredi 19 f?vrier 2014 15:16 > > > > Is there a way to tell Scilab if this gets beyond a limit to round off > > at zero say? > > Maybe format() or clean() ? > > http://help.scilab.org/docs/5.4.1/en_US/format.html > http://help.scilab.org/docs/5.4.1/en_US/clean.html > > Best regards > > -- > Christophe Dang Ngoc Chan > Mechanical calculation engineer > > ______________________________________________________________________ > > This e-mail may contain confidential and/or privileged information. If you > are not the intended recipient (or have received this e-mail in error), > please notify the sender immediately and destroy this e-mail. Any > unauthorized copying, disclosure or distribution of the material in this > e-mail is strictly forbidden. > ______________________________________________________________________ > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From denis.crete at thalesgroup.com Wed Feb 19 16:05:34 2014 From: denis.crete at thalesgroup.com (CRETE Denis) Date: Wed, 19 Feb 2014 16:05:34 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392732408401-4028696.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> Message-ID: <22361_1392822336_5304C840_22361_9512_1_8F1D40232A0E68409E3FC23A30C3266201690C8E195C@THSONEA01CMS04P.one.grp> Dear Yohann The (x,y) dataset does not look like a parabola... (it resembles more a type of x*sin(x) function). Trying to fit this dataset with a parabola, will not only give unreliable result, but also will put a high challenge to error estimation algorithms... Denis -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Yohann Envoy??: mardi 18 f?vrier 2014 15:07 ??: users at lists.scilab.org Objet?: [Scilab-users] evaluate error on each parameter calculated with leastsq Hi all !! First of all, thank you to all contributors of Scilab. However, I have a problem to calculate the estimated error on each parameter calculated with leastsq. For example if I have the simple following function to fit: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! function y = fun2fit(x, c) y = c(1)*x*x + c(2)*x + c(3); endfunction function e = myerror(c, x, y) e = fun2fit(x, c) - y; endfunction // Data x = [ 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6]'; y = [5.02 6.08 3.33 -0.93 -0.22 7.83 16.52 15.55 2.67 -11.42 -11.78 5.09 25.25]'; // First guess c0 = [1 1 1]'; y0 = fun2fit(x, c0); // least square [f, copt] = leastsq(list(myerror, x, y), c0) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! I'd like to have an estimated error on each of the 3 element of copt. Any idea ? Thanks ! Y -- View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users From cfuttrup at gmail.com Wed Feb 19 17:32:07 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Wed, 19 Feb 2014 17:32:07 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392732408401-4028696.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> Message-ID: <5304DC87.4020507@gmail.com> Hi Yohan I see that I'm not alone wondering a bit about what you'd like to do. A mathematically solid and readily available technique to corner mathematical rounding errors and other errors of operation is called Interval Analysis. For matlab there's INTLAB. See: http://www.ti3.tu-harburg.de/~rump/intlab/ I'm sure it could be converted to Scilab without big trouble. Maybe I even have an old one laying around (??) somewhere ... maybe. Another approach would be in a specific situation to apply the techniques of interval analysis, only centered around your concern ... and leave the rest of your script normal. Best regards, Claus On 2/18/2014 15:06, Yohann wrote: > Hi all !! > > First of all, thank you to all contributors of Scilab. > However, I have a problem to calculate the estimated error on each parameter > calculated with leastsq. > > For example if I have the simple following function to fit: > !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! > function y = fun2fit(x, c) > y = c(1)*x*x + c(2)*x + c(3); > endfunction > > function e = myerror(c, x, y) > e = fun2fit(x, c) - y; > endfunction > > // Data > x = [ 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6]'; > y = [5.02 6.08 3.33 -0.93 -0.22 7.83 16.52 15.55 2.67 -11.42 -11.78 > 5.09 25.25]'; > > // First guess > c0 = [1 1 1]'; > y0 = fun2fit(x, c0); > > // least square > [f, copt] = leastsq(list(myerror, x, y), c0) > !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! > > I'd like to have an estimated error on each of the 3 element of copt. > Any idea ? > > Thanks ! > Y > > > > > > -- > View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > . > --- This email is free from viruses and malware because avast! Antivirus protection is active. http://www.avast.com From barbanera at crbnet.it Wed Feb 19 11:31:47 2014 From: barbanera at crbnet.it (bomber1978) Date: Wed, 19 Feb 2014 02:31:47 -0800 (PST) Subject: [Scilab-users] To find maximum value in a random curve Message-ID: <1392805906969-4028741.post@n3.nabble.com> Hi all, I have a random 2D-curve which is defined by a set of points (I don't know the equation of the curve), obtained by interpolation. How can I find the maximum, minimum, & inflexion points on this curve. Thank you very much! -- View this message in context: http://mailinglists.scilab.org/To-find-maximum-value-in-a-random-curve-tp4028741.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From ml-node+s994242n4028741h12 at n3.nabble.com Wed Feb 19 11:31:47 2014 From: ml-node+s994242n4028741h12 at n3.nabble.com (bomber1978 [via Scilab / Xcos - Mailing Lists Archives]) Date: Wed, 19 Feb 2014 02:31:47 -0800 (PST) Subject: [Scilab-users] To find maximum value in a random curve Message-ID: <1392805906969-4028741.post@n3.nabble.com> Hi all, I have a random 2D-curve which is defined by a set of points (I don't know the equation of the curve), obtained by interpolation. How can I find the maximum, minimum, & inflexion points on this curve. Thank you very much! ______________________________________ If you reply to this email, your message will be added to the discussion below: http://mailinglists.scilab.org/To-find-maximum-value-in-a-random-curve-tp4028741.html This email was sent by bomber1978 (via Nabble) To receive all replies by email, subscribe to this discussion: http://mailinglists.scilab.org/template/NamlServlet.jtp?macro=subscribe_by_code&node=4028741&code=dXNlcnNAbGlzdHMuc2NpbGFiLm9yZ3w0MDI4NzQxfC0xODEwOTc0OTQz -------------- next part -------------- An HTML attachment was scrubbed... URL: From romain.desbats at ifsttar.fr Wed Feb 19 11:34:40 2014 From: romain.desbats at ifsttar.fr (Romain Desbats) Date: Wed, 19 Feb 2014 11:34:40 +0100 (CET) Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <365589264.1926236.1392804055518.JavaMail.root@ifsttar.fr> Message-ID: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> Dear Scilab users, I am running a simple electrical model in xcos (see attached). This model reads a .csv input file and outputs a variable to the workspace. The model works fine when I compile it and run it by using the menus in xcos. What I would like is to run it in batch mode from Scilab in order optimise the model parameters. My plan is to use xcos_simulate to compile the diagram, to modify the cpr variable generated by the compilation (to change the model parameters) and than to simulate the model with scicosim . The problem is that calling the cpr variable after the compilation results in kind of "freezing" the workspace. - The variable browser shows "N/A" for all the variables (from image 1 to image 2 attached) - Calling cpr does not give any result (image 3 attached) Does anyone has any idea of what is happening? I am using Scilab 5.4.1. My code is as follows: // parameters data=csvRead('input.csv',';'); input_current=evstr(data(:,1)); size1=length(input_current); C=1/10000; R=12900; T=0.0001; // timestep time=[0:T:(size1-1)*T]'; // in s end_time=time($)+T; // end time in s // loading Xcos functions loadScicos; loadXcosLibs(); exec(loadpallibs,-1); // loading model importXcosDiagram('graphic_model.xcos') // compilation: scs_m is the result of the import, xcos_simulate(scs_m,4) // '4' is to perform compilation tf=end_time; atol=1.e-8; rtol=1.e-9; ttol=1.e-10; deltat=tf; scale=0; hmax=0; tol=[atol,rtol,ttol,deltat,scale,hmax]; // PROBLEM COMES WHEN EXECUTING NEXT LINE state=cpr.state; // cpr variable came from the compilation sim=cpr.sim; cpr.sim.rpar(5)=T; // timestep cpr.sim.rpar(5)=0; // initialisation cpr.sim.rpar(7)=R; // resistance cpr.sim.rpar(8)=C; // capacitor // initialisation [state,t]=scicosim(state,0,tf,sim,'start',tol); // run [state,t]=scicosim(state,t,tf,sim,'run',tol) //end [state,t]=scicosim(state,t,tf,sim,'finish',tol); graph_voltage=graph_voltage.values; // output variable Thanks a lot for your help. Best regards, Romain -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 1.png Type: image/png Size: 121580 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 2.png Type: image/png Size: 121447 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 3.png Type: image/png Size: 120644 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: input.csv Type: text/csv Size: 2222 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: graphic_model.xcos Type: application/xml Size: 74159 bytes Desc: not available URL: From yohann.morille at univ-angers.fr Wed Feb 19 14:31:17 2014 From: yohann.morille at univ-angers.fr (Yohann) Date: Wed, 19 Feb 2014 05:31:17 -0800 (PST) Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <53047A24.90004@laas.fr> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> Message-ID: <1392816677490-4028742.post@n3.nabble.com> Hi Antoine, thank you for your answer but what I need is a confidence interval on each parameter ! Cheers Yohann -- View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028742.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From yohann.morille at univ-angers.fr Wed Feb 19 17:10:57 2014 From: yohann.morille at univ-angers.fr (Yohann) Date: Wed, 19 Feb 2014 08:10:57 -0800 (PST) Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <22361_1392822336_5304C840_22361_9512_1_8F1D40232A0E68409E3FC23A30C3266201690C8E195C@THSONEA01CMS04P.one.grp> References: <1392732408401-4028696.post@n3.nabble.com> <22361_1392822336_5304C840_22361_9512_1_8F1D40232A0E68409E3FC23A30C3266201690C8E195C@THSONEA01CMS04P.one.grp> Message-ID: <1392826257532-4028748.post@n3.nabble.com> yes I know Denis, It was just an example to illustrate my question. My real dataset and function to fit are completely different and more complex. Thank you -- View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028748.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.bignier at scilab-enterprises.com Thu Feb 20 10:35:55 2014 From: paul.bignier at scilab-enterprises.com (paul.bignier at scilab-enterprises.com) Date: Thu, 20 Feb 2014 10:35:55 +0100 Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> References: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> Message-ID: <6a0f175ed0d6c0c5408403abfefaa077@scilab-enterprises.com> Hello Romain, Running your code in 5.4.1 just tells me that 'cpr' variable is unknown when you call it at line "state=cpr.state;". So I replaced "xcos_simulate(scs_m,4)" by "cpr = xcos_simulate(scs_m,4)" and it runs smooth ("-->cpr" yields an output) until the end. I noted that you modifed cpr.sim.rpar(5) twice, but you probably meant cpr.sim.rpar(6) the second time? cpr.sim.rpar(5)=T; // timestep cpr.sim.rpar(5)=0; // initialisation Otherwise, I didn't manage to reproduce your bug. Regards, Paul PS: you may want to use the Context to modify a diagram from batch, see xcos_simulate help page [3]. It is safer to fully recompile Modelica diagrams between simulations. On 02/19/2014 11:34 AM, Romain Desbats wrote: > Dear Scilab users, > > I am running a simple electrical model in xcos (see attached). This model reads a .csv input file and outputs a variable to the workspace. The model works fine when I compile it and run it by using the menus in xcos. > > What I would like is to run it in batch mode from Scilab in order optimise the model parameters. > > My plan is to use _xcos_simulate_ to compile the diagram, to modify the _cpr_ variable generated by the compilation (to change the model parameters) and than to simulate the model with _scicosim_. > > The problem is that calling the _cpr_ variable after the compilation results in kind of "freezing" the workspace. > > - The variable browser shows "N/A" for all the variables (from image 1 to image 2 attached) > - Calling _cpr_ does not give any result (image 3 attached) > > Does anyone has any idea of what is happening? > > I am using Scilab 5.4.1. > > My code is as follows: > > // parameters > data=csvRead('input.csv',';'); > input_current=evstr(data(:,1)); > size1=length(input_current); > > C=1/10000; > R=12900; > > T=0.0001; // timestep > > time=[0:T:(size1-1)*T]'; // in s > > end_time=time($)+T; // end time in s > > // loading Xcos functions > loadScicos; > loadXcosLibs(); > exec(loadpallibs,-1); > > // loading model > importXcosDiagram('graphic_model.xcos') > > // compilation: scs_m is the result of the import, > xcos_simulate(scs_m,4) // '4' is to perform compilation > > tf=end_time; > > atol=1.e-8; > rtol=1.e-9; > ttol=1.e-10; > deltat=tf; > scale=0; > hmax=0; > > tol=[atol,rtol,ttol,deltat,scale,hmax]; > > // PROBLEM COMES WHEN EXECUTING NEXT LINE > > state=cpr.state; // cpr variable came from the compilation > sim=cpr.sim; > > cpr.sim.rpar(5)=T; // timestep > cpr.sim.rpar(5)=0; // initialisation > cpr.sim.rpar(7)=R; // resistance > cpr.sim.rpar(8)=C; // capacitor > > // initialisation > [state,t]=scicosim(state,0,tf,sim,'start',tol); > // run > [state,t]=scicosim(state,t,tf,sim,'run',tol) > //end > [state,t]=scicosim(state,t,tf,sim,'finish',tol); > > graph_voltage=graph_voltage.values; // output variable > > Thanks a lot for your help. > > Best regards, > > Romain > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users [1] -- Paul BIGNIER Development engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com [2] Links: ------ [1] http://lists.scilab.org/mailman/listinfo/users [2] http://www.scilab-enterprises.com [3] http://help.scilab.org/docs/5.4.1/en_US/xcos_simulate.html -------------- next part -------------- An HTML attachment was scrubbed... URL: From Serge.Steer at inria.fr Thu Feb 20 12:07:46 2014 From: Serge.Steer at inria.fr (Serge Steer) Date: Thu, 20 Feb 2014 12:07:46 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392816677490-4028742.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> <1392816677490-4028742.post@n3.nabble.com> Message-ID: <5305E202.60906@inria.fr> Le 19/02/2014 14:31, Yohann a ?crit : > Hi Antoine, > thank you for your answer but > what I need is a confidence interval on each parameter ! The confidence interval on the computed parameters depends on the confidence interval on your x and y data. In your case where you are looking for polynomial fitting then perturbations on y produces perturbation on the computed parameters roughly proportionnal to the conditionning of the matrix A=[x*x x ones(x)] if c is the mean square solution of the probleme for y and c1 is the mean square solution of the problem of y1 then norm(c-c1)/norm(c)<=cond(A)*norm(y-y1)/norm(y) x = [ 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6]'; y = [5.02 6.08 3.33 -0.93 -0.22 7.83 16.52 15.55 2.67 -11.42 -11.78 5.09 25.25]'; A=[x.^2 x ones(x)]; y1=y;y1=y1.*(1+0.000001*rand(y1));ey=norm(y-y1)/norm(y) c1=A\y1;ec=norm(c-c1)/norm(c) ec<=cond(A)*ey the dependance with respect to perturbation on x is more difficult. > Cheers > Yohann > > > > > -- > View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028742.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From antoine.monmayrant at laas.fr Thu Feb 20 09:04:58 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Thu, 20 Feb 2014 09:04:58 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392816677490-4028742.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> <1392816677490-4028742.post@n3.nabble.com> Message-ID: <5305B72A.3090500@laas.fr> On 02/19/2014 02:31 PM, Yohann wrote: > Hi Antoine, > thank you for your answer but > what I need is a confidence interval on each parameter ! Well, so resampling can do the trick! Antoine > Cheers > Yohann > > > > > -- > View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028742.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- +++++++++++++++++++++++++++++++++++++++++++++++++++++++ Antoine Monmayrant LAAS - CNRS 7 avenue du Colonel Roche BP 54200 31031 TOULOUSE Cedex 4 FRANCE Tel:+33 5 61 33 64 59 email : antoine.monmayrant at laas.fr permanent email : antoine.monmayrant at polytechnique.org +++++++++++++++++++++++++++++++++++++++++++++++++++++++ From stephane.mottelet at gmail.com Thu Feb 20 10:00:38 2014 From: stephane.mottelet at gmail.com (=?ISO-8859-1?Q?St=E9phane_Mottelet?=) Date: Thu, 20 Feb 2014 10:00:38 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: References: <1392732408401-4028696.post@n3.nabble.com> <22361_1392822336_5304C840_22361_9512_1_8F1D40232A0E68409E3FC23A30C3266201690C8E195C@THSONEA01CMS04P.one.grp> <1392826257532-4028748.post@n3.nabble.com> Message-ID: Here is a simple example where the derivative is approximated with the "derivative" macro of Scilab. S. 2014-02-20 9:43 GMT+01:00 St?phane Mottelet : > Hello, > > If you can make the hypothesis that your data is corrupted by gaussian > noise, then you can approximate the covariance matrix of your estimated > parameters. Let p be the vector of parameters and r(p) the residual vector > given by > > r(p)=sigma^(-1)*(y-Y(p)) > > where y is your measuement vector, Y(p) the "simulated" measurement, sigma > a diagonal matrix with the std error for each measurement. If we denote by > > drdp(p) the derivative (or jacobian matrix) of r with respect to p then > the covariance matrice C of parameters can be estimated by > > > C=F^(-1) > > where > > F=drdp(p)' * drdp(p) > > is the Fisher information matrix. The diagonal terms of V give you the > variance of the parameters. Of course, even in the gaussian case, this is a > crude approximation.... > > S. > > > 2014-02-19 17:10 GMT+01:00 Yohann : > > yes I know Denis, It was just an example to illustrate my question. >> My real dataset and function to fit are completely different and more >> complex. >> Thank you >> >> >> >> -- >> View this message in context: >> http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028748.html >> Sent from the Scilab users - Mailing Lists Archives mailing list archive >> at Nabble.com. >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> > > -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: premier_ordre_stats.sce Type: application/octet-stream Size: 3039 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: premier_ordre_data.dat Type: application/octet-stream Size: 544 bytes Desc: not available URL: From aweeks at hidglobal.com Fri Feb 21 11:45:30 2014 From: aweeks at hidglobal.com (aweeks at hidglobal.com) Date: Fri, 21 Feb 2014 10:45:30 +0000 Subject: [Scilab-users] Switching windows 'focus' in SciLab. Message-ID: Hello All, I have SciLab set up with a docked SciNotes and console. Calling 'plot', in the SciLab script, creates a separate graphics window which has the 'windows focus' on it. Is there a way, within the running SciLab script, to switch the focus back onto the console (specifically for keyboard input) ? (SciLab 5.4.0 running in Windows XP SP3 32-bit). Many thanks in anticipation. Adrian Weeks Development Engineer HID Global, 3 Cae Gwyrdd, Green Meadow Springs Business Park, Cardiff CF15 7AB, United Kingdom. +44 (0)29 20528523 (Office) aweeks at hidglobal.com www.hidglobal.com -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 13947181.gif Type: image/gif Size: 4559 bytes Desc: not available URL: From aweeks at hidglobal.com Fri Feb 21 12:00:25 2014 From: aweeks at hidglobal.com (aweeks at hidglobal.com) Date: Fri, 21 Feb 2014 11:00:25 +0000 Subject: [Scilab-users] Linking DLL from MS Visual C++ to SciLab. Message-ID: Hello All, I can create a DLL using Microsoft Visual C++ and I can link to the exported functions in SciLab using 'link' and then call them using 'call'. If I understand correctly, the exported functions have to return void and pass values through the argument list pointers. I have a third-party DLL that I would like to use in SciLab and where the functions return types other than void. The more elegant method would seem to be to use the API_Scilab gateway to interface to the exported functions. Can I / how do I link a DLL (compiled in Visual C++) to SciLab at the 'ilib_build' stage ? (SciLab 5.4.0, Windows XP SP3 32-bit, Visual C++ 2008 Express). Many thanks in anticipation. Adrian Weeks Development Engineer HID Global, 3 Cae Gwyrdd, Green Meadow Springs Business Park, Cardiff CF15 7AB, United Kingdom. +44 (0)29 20528523 (Office) aweeks at hidglobal.com www.hidglobal.com -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 12748956.gif Type: image/gif Size: 4559 bytes Desc: not available URL: From sgougeon at free.fr Wed Feb 19 20:55:21 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Wed, 19 Feb 2014 20:55:21 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <5304DC87.4020507@gmail.com> References: <1392732408401-4028696.post@n3.nabble.com> <5304DC87.4020507@gmail.com> Message-ID: <53050C29.1000308@free.fr> Le 19/02/2014 17:32, Claus Futtrup a ?crit : > Hi Yohan > > I see that I'm not alone wondering a bit about what you'd like to do. > > A mathematically solid and readily available technique to corner > mathematical rounding errors and other errors of operation is called > Interval Analysis. > > For matlab there's INTLAB. See: > http://www.ti3.tu-harburg.de/~rump/intlab/ > > I'm sure it could be converted to Scilab without big trouble. Maybe I > even have an old one laying around (??) somewhere ... maybe. Yes indeed, here : http://www-sop.inria.fr/coprin/logiciels/Int4Sci/ Samuel From jasper at amsterchem.com Thu Feb 20 09:31:49 2014 From: jasper at amsterchem.com (jasper van baten) Date: Thu, 20 Feb 2014 09:31:49 +0100 Subject: [Scilab-users] To find maximum value in a random curve In-Reply-To: <1392805906969-4028741.post@n3.nabble.com> References: <1392805906969-4028741.post@n3.nabble.com> Message-ID: <5305BD75.8030801@amsterchem.com> You could determine the coefficients of a cubic spline, then determine analytically the extrema and inflection points on the cubic equation that defines each interval (there where the roots of the parabolic that is the first order derivative of the cubic are within the interval). See here for setting up the spline: http://en.wikipedia.org/wiki/Spline_interpolation Best wishes, Jasper. On 2/19/2014 11:31, bomber1978 [via Scilab / Xcos - Mailing Lists Archives] wrote: > Hi all, > > I have a random 2D-curve which is defined by a set of points (I don't > know the equation of the curve), obtained by interpolation. > How can I find the maximum, minimum, & inflexion points on this curve. > > Thank you very much! > > ------------------------------------------------------------------------ > If you reply to this email, your message will be added to the > discussion below: > http://mailinglists.scilab.org/To-find-maximum-value-in-a-random-curve-tp4028741.html > > This email was sent by bomber1978 > > (via Nabble) > To receive all replies by email, subscribe to this discussion > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From grocer.toolbox at gmail.com Thu Feb 20 09:08:42 2014 From: grocer.toolbox at gmail.com (Eric Dubois) Date: Thu, 20 Feb 2014 09:08:42 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392816677490-4028742.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> <1392816677490-4028742.post@n3.nabble.com> Message-ID: Hello. You could have ? look at the nls function in my toolbox Grocer, avaolable undet Atoms ( see also my web page dubois.ensae.net/grocer.html). ?ric Le 20 f?vr. 2014 08:32, "Yohann" a ?crit : > Hi Antoine, > thank you for your answer but > what I need is a confidence interval on each parameter ! > Cheers > Yohann > > > > > -- > View this message in context: > http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028742.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive > at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at gmail.com Thu Feb 20 09:43:46 2014 From: stephane.mottelet at gmail.com (=?ISO-8859-1?Q?St=E9phane_Mottelet?=) Date: Thu, 20 Feb 2014 09:43:46 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392826257532-4028748.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> <22361_1392822336_5304C840_22361_9512_1_8F1D40232A0E68409E3FC23A30C3266201690C8E195C@THSONEA01CMS04P.one.grp> <1392826257532-4028748.post@n3.nabble.com> Message-ID: Hello, If you can make the hypothesis that your data is corrupted by gaussian noise, then you can approximate the covariance matrix of your estimated parameters. Let p be the vector of parameters and r(p) the residual vector given by r(p)=sigma^(-1)*(y-Y(p)) where y is your measuement vector, Y(p) the "simulated" measurement, sigma a diagonal matrix with the std error for each measurement. If we denote by drdp(p) the derivative (or jacobian matrix) of r with respect to p then the covariance matrice C of parameters can be estimated by C=F^(-1) where F=drdp(p)' * drdp(p) is the Fisher information matrix. The diagonal terms of V give you the variance of the parameters. Of course, even in the gaussian case, this is a crude approximation.... S. 2014-02-19 17:10 GMT+01:00 Yohann : > yes I know Denis, It was just an example to illustrate my question. > My real dataset and function to fit are completely different and more > complex. > Thank you > > > > -- > View this message in context: > http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028748.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive > at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From antoine.monmayrant at laas.fr Fri Feb 21 16:06:23 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Fri, 21 Feb 2014 16:06:23 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <53050C29.1000308@free.fr> References: <1392732408401-4028696.post@n3.nabble.com> <5304DC87.4020507@gmail.com> <53050C29.1000308@free.fr> Message-ID: <53076B6F.1070306@laas.fr> On 02/19/2014 08:55 PM, Samuel Gougeon wrote: > Le 19/02/2014 17:32, Claus Futtrup a ?crit : >> Hi Yohan >> >> I see that I'm not alone wondering a bit about what you'd like to do. >> >> A mathematically solid and readily available technique to corner >> mathematical rounding errors and other errors of operation is called >> Interval Analysis. >> >> For matlab there's INTLAB. See: >> http://www.ti3.tu-harburg.de/~rump/intlab/ >> >> I'm sure it could be converted to Scilab without big trouble. Maybe I >> even have an old one laying around (??) somewhere ... maybe. > Yes indeed, here : http://www-sop.inria.fr/coprin/logiciels/Int4Sci/ The thing is I don't see how to use these tools to do what Yoahn wants to do. I think I have the same kind of issues than Yoahn (I've done a fit, how good is it as compare to another one, and what kind of "confidence" or "error bar" should I associate to each parameter). Could you tell us how we can use interval analysis to adress this (just curious, it's a bit far from my field of expertise)? Cheers, Antoine From reitertamas at gmail.com Thu Feb 20 09:57:42 2014 From: reitertamas at gmail.com (=?ISO-8859-1?Q?Tam=E1s_Reiter?=) Date: Thu, 20 Feb 2014 08:57:42 +0000 Subject: [Scilab-users] help, please Message-ID: Hi all, I would appreciate if anyone could suggest a solution for extracting the content of a column vector to a pre-defined matrix size (identical row numbers) within a for loop. I keep getting an error message:Submatrix incorrecly defined... my loop is as follows: for j=1:numpcs do t=X(:,1);//pick the score vector (first column of X) t_pre=t-t; while norm(t-t_pre) > threshold do p=(X'*t)/norm(X'*t);//calculate and normalise the loading vector t_pre=t;//saving the original vector into t_pre t=X*p;//calculating the new score vector (product of eigenvalue and eigenvector) end //extracting components explored X=X-t*p';//calculating the residual matrix *T(:,j)=t(:,:);--THIS IS WHERE I AM GETTING THIS ERRORP(:,j)=p(:,:);--THIS IS WHERE I AM GETTING THIS ERROR*end //T(:,j+1:numpcs)=[]; //P(:,j+1:numpcs)=[]; S=diag(T'*T); Very much appriciate any help, Kind regards, Tamas -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at gmail.com Thu Feb 20 11:15:30 2014 From: stephane.mottelet at gmail.com (=?ISO-8859-1?Q?St=E9phane_Mottelet?=) Date: Thu, 20 Feb 2014 11:15:30 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392816677490-4028742.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> <1392816677490-4028742.post@n3.nabble.com> Message-ID: Hello, If you can make the hypothesis that your data is corrupted by gaussian noise, then you can approximate the covariance matrix of your estimated parameters. Let p be the vector of parameters and r(p) the residual vector given by r(p)=sigma^(-1)*(y-Y(p)) where y is your measuement vector, Y(p) the "simulated" measurement, sigma a diagonal matrix with the std error for each measurement. If we denote by drdp(p) the derivative (or jacobian matrix) of r with respect to p then the covariance matrice C of parameters can be estimated by C=F^(-1) where F=drdp(p)' * drdp(p) is the Fisher information matrix. The diagonal terms of V give you the variance of the parameters. Of course, even in the gaussian case, this is a crude approximation.... Here is (attached) a simple example where the derivative is approximated with the "derivative" macro of Scilab. S. 2014-02-19 14:31 GMT+01:00 Yohann : > Hi Antoine, > thank you for your answer but > what I need is a confidence interval on each parameter ! > Cheers > Yohann > > > > > -- > View this message in context: > http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028742.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive > at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: premier_ordre_stats.sce Type: application/octet-stream Size: 3039 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: premier_ordre_data.dat Type: application/octet-stream Size: 544 bytes Desc: not available URL: From mariajovera at icloud.com Fri Feb 21 06:18:14 2014 From: mariajovera at icloud.com (samaelkreutz) Date: Thu, 20 Feb 2014 21:18:14 -0800 (PST) Subject: [Scilab-users] finding similar values with scilab Message-ID: <1392959894551-4028792.post@n3.nabble.com> Hello!!! I have a big question ! I have a biiiiiiiiiiig data archive with text files with the form: 18.87 2.6 0.00545558 19.98 2.6 0.00225349 18.87 2.6 0.00405905 13.32 2.6 0.01338288 19.98 2.6 0.01537532 18.87 2.6 0.00481375 19.98 2.6 0.00936207 12.21 2.6 0.00558517 I need eliminate extreme values, for example: If I plot a random text file I obtain a image like this: As you can see theres points extremes. If i make a zoom DN=find(DM(:,3)<=0.0001 & DM(:,3)>=0.00002 ) dn=DM(DN,:) figure(2) plot(dn(:,1),dn(:,3),'O') I have this picture My question is: Has scilab a command that can help me to find the "dense data" , I mean find all the values that are similar. I ask this cause I have about 20,000 text files and is imposible for me make "artesal zoms". Please any suggestion, answer going to help me. -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From romain.desbats at ifsttar.fr Thu Feb 20 11:21:52 2014 From: romain.desbats at ifsttar.fr (R0M1) Date: Thu, 20 Feb 2014 02:21:52 -0800 (PST) Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <6a0f175ed0d6c0c5408403abfefaa077@scilab-enterprises.com> References: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> <6a0f175ed0d6c0c5408403abfefaa077@scilab-enterprises.com> Message-ID: <1392891712865-4028769.post@n3.nabble.com> Hi Paul, Thanks a lot for your answer!! Actually, changing "xcos_simulate(scs_m,4)" to "cpr=xcos_simulate(scs_m,4)" gave me the error "Probl?me avec la pile" ("Problem with the stack" in English I assume). But this was maybe due to the initial parameters of timestep / en time in the diagram. But anyway, using the context to change the parameters is much easier and it works perfectly!! Exactly as in the example of the help page you mentioned. Thanks again for this. Regards, Romain -- View this message in context: http://mailinglists.scilab.org/Scilab-users-Problem-running-xcos-in-batch-mode-from-Scilab-tp4028765p4028769.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From matiashs at hotmail.com Fri Feb 21 10:04:35 2014 From: matiashs at hotmail.com (matacosta) Date: Fri, 21 Feb 2014 01:04:35 -0800 (PST) Subject: [Scilab-users] Complex text file problem getting matrices from it Message-ID: <1392973475489-4028793.post@n3.nabble.com> Dear All, I need to obtain two matrices from a text file. Normally this should not be a problem with fscanfMat but because I have headings in the middle of the file this function does not work. I though up to now in two solutions. One is using mgetl and obtain the part where the first numbers are and then the second (I have a problem with this since I am not sure if mgetl is able to start getting the data from a certain line on or only beginning). Therefore, I got properly the first half (which normally the fscanfMat analysis properly) but the fscanfMat does not recognize the 3 columns, only 1 when I save with write. Is something of formatting but I do not get why (related to tab and space, I checked manually). The other to read the whole part of interested as 1 string with mgetstr. Then separate this string with strsplit which I am having problems with the matched_separators vector. Then I would transform that matrix of strings into double. This would allow me to get both matrices but again i would have the same problem of writing function. Any help? better way of doing it? A simplified version of the data is down: "BW: 4" "SWEEP AVERAGING: OFF" "POINT AVERAGING: OFF" "TRACE: A" "FORMAT: LINEAR" "Frequency" "Data Trace Real" "Data Trace Imag" 1.00000000000e+03 1.883005e-11 0.000000e+00 1.00000000000e+03 1.876519e-11 0.000000e+00 1.00000000000e+03 1.880004e-11 0.000000e+00 "TRACE: B" "FORMAT: LINEAR" "Frequency" "Data Trace Real" "Data Trace Imag" 1.00000000000e+03 1.906406e-02 0.000000e+00 1.00000000000e+03 1.790656e-02 0.000000e+00 1.00000000000e+03 1.664777e-02 0.000000e+00 1.00000000000e+03 1.774269e-02 0.000000e+00 Thanks a lot!! -- View this message in context: http://mailinglists.scilab.org/Complex-text-file-problem-getting-matrices-from-it-tp4028793.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Mike at Page-One.Waitrose.com Fri Feb 21 18:28:39 2014 From: Mike at Page-One.Waitrose.com (Mike Page) Date: Fri, 21 Feb 2014 17:28:39 -0000 Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <1392959894551-4028792.post@n3.nabble.com> Message-ID: Hi, If you want to remove outliers from a data set, you can use the modified Thompson Tau method. Basically tau provides a limit based on mean and standard deviation so that values outside (mean +/- n * std dev) are rejected. The method is unfortunately iterative, but it works well. You can remove the located outliers using something like x = x(find(delta < tau * sd)), where delta is the absolute difference from the mean. Hope that helps, Mike. -----Original Message----- From: users [mailto:users-bounces at lists.scilab.org]On Behalf Of samaelkreutz Sent: 21 February 2014 05:18 To: users at lists.scilab.org Subject: [Scilab-users] finding similar values with scilab Hello!!! I have a big question ! I have a biiiiiiiiiiig data archive with text files with the form: 18.87 2.6 0.00545558 19.98 2.6 0.00225349 18.87 2.6 0.00405905 13.32 2.6 0.01338288 19.98 2.6 0.01537532 18.87 2.6 0.00481375 19.98 2.6 0.00936207 12.21 2.6 0.00558517 I need eliminate extreme values, for example: If I plot a random text file I obtain a image like this: As you can see theres points extremes. If i make a zoom DN=find(DM(:,3)<=0.0001 & DM(:,3)>=0.00002 ) dn=DM(DN,:) figure(2) plot(dn(:,1),dn(:,3),'O') I have this picture My question is: Has scilab a command that can help me to find the "dense data" , I mean find all the values that are similar. I ask this cause I have about 20,000 text files and is imposible for me make "artesal zoms". Please any suggestion, answer going to help me. -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792. html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users ----- No virus found in this message. Checked by AVG - www.avg.com Version: 2014.0.4335 / Virus Database: 3705/7111 - Release Date: 02/20/14 From stephane.mottelet at gmail.com Fri Feb 21 19:22:45 2014 From: stephane.mottelet at gmail.com (=?ISO-8859-1?Q?St=E9phane_Mottelet?=) Date: Fri, 21 Feb 2014 19:22:45 +0100 Subject: [Scilab-users] Test Message-ID: <53079975.3080307@gmail.com> This is just a test, sorry for the disturbance. S. From chuox.bece at gmail.com Fri Feb 21 19:31:36 2014 From: chuox.bece at gmail.com (Chuox) Date: Fri, 21 Feb 2014 10:31:36 -0800 (PST) Subject: [Scilab-users] Linking DLL from MS Visual C++ to SciLab. In-Reply-To: References: Message-ID: <1393007496941-4028813.post@n3.nabble.com> Using a third-party DLL in SciLab might be tricky. I had a problem linking EPANET toolkit library, what i did was create a dll to call the epanet2.dll (the library form EPANET) and return the values in SciLab. Maybe is not the best solution and there is a better way, but i didn't find it. This is and example of how i did it: /////The C code///// #include #include "epanet2.h" int SetLinkValue(int *index, int *param, float *Val, int *errcode ) { *errcode=ENsetlinkvalue(*index, *param, *Val); return (*errcode); } /////SciLab Code///// //Link the dll used to cal the external dll unSetLink=link(Ubicacion+"enlace.dll","SetLinkValue","c"); //This return and error [err]=call("SetLinkValue",TubIndex,1,"i",0,2,"i",diametro,3,"r","out",[1,1],4,"i"); -- View this message in context: http://mailinglists.scilab.org/Scilab-users-Linking-DLL-from-MS-Visual-C-to-SciLab-tp4028800p4028813.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From stephane.mottelet at gmail.com Fri Feb 21 19:34:04 2014 From: stephane.mottelet at gmail.com (=?ISO-8859-1?Q?St=E9phane_Mottelet?=) Date: Fri, 21 Feb 2014 19:34:04 +0100 Subject: [Scilab-users] Test Message-ID: This a test, sorry. S. -------------- next part -------------- An HTML attachment was scrubbed... URL: From ian.h.bell at gmail.com Fri Feb 21 19:39:24 2014 From: ian.h.bell at gmail.com (Ian Bell) Date: Fri, 21 Feb 2014 19:39:24 +0100 Subject: [Scilab-users] Linking DLL from MS Visual C++ to SciLab. In-Reply-To: <1393007496941-4028813.post@n3.nabble.com> References: <1393007496941-4028813.post@n3.nabble.com> Message-ID: Yes, using third-party DLL with scilab does seem to be hard. Two crucial things are not supported: - All inputs to the function are passed by *reference*, not value, but often times DLL are build with inputs being passed by value - Functions can not return a value. IMHO these are very serious problems with the DLL integration. I couldn't use my DLL in scilab because of these issues. And the documentation on this topic is quite weak. Especially the always passing by reference part. That took me a few hours of debugging to realize that. Ian On Fri, Feb 21, 2014 at 7:31 PM, Chuox wrote: > Using a third-party DLL in SciLab might be tricky. > I had a problem linking EPANET toolkit library, what i did was create a dll > to call the epanet2.dll (the library form EPANET) and return the values in > SciLab. Maybe is not the best solution and there is a better way, but i > didn't find it. > This is and example of how i did it: > > /////The C code///// > #include > #include "epanet2.h" > int SetLinkValue(int *index, int *param, float *Val, int *errcode ) > { > *errcode=ENsetlinkvalue(*index, *param, *Val); > return (*errcode); > } > > /////SciLab Code///// > //Link the dll used to cal the external dll > unSetLink=link(Ubicacion+"enlace.dll","SetLinkValue","c"); > > //This return and error > > [err]=call("SetLinkValue",TubIndex,1,"i",0,2,"i",diametro,3,"r","out",[1,1],4,"i"); > > > > -- > View this message in context: > http://mailinglists.scilab.org/Scilab-users-Linking-DLL-from-MS-Visual-C-to-SciLab-tp4028800p4028813.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive > at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Fri Feb 21 20:50:30 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Fri, 21 Feb 2014 20:50:30 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <53076B6F.1070306@laas.fr> References: <1392732408401-4028696.post@n3.nabble.com> <5304DC87.4020507@gmail.com> <53050C29.1000308@free.fr> <53076B6F.1070306@laas.fr> Message-ID: <5307AE06.5090002@gmail.com> Hi Antoine Initial message from Yohan was not sufficient to determine whether Interval Analysis could do it or not. As I said in my previous answer - IntLab is for determining mathematical rounding errors and other errors of operation ... for example if you make almost division-by-zero, then you have a big gap of potentially correct solutions. This is what IntLab is for, robust mathematics, not curve fitting or analyzing measurement data. Best regards, Claus On 2/21/2014 16:06, Antoine Monmayrant wrote: > On 02/19/2014 08:55 PM, Samuel Gougeon wrote: >> Le 19/02/2014 17:32, Claus Futtrup a ?crit : >>> Hi Yohan >>> >>> I see that I'm not alone wondering a bit about what you'd like to do. >>> >>> A mathematically solid and readily available technique to corner >>> mathematical rounding errors and other errors of operation is called >>> Interval Analysis. >>> >>> For matlab there's INTLAB. See: >>> http://www.ti3.tu-harburg.de/~rump/intlab/ >>> >>> I'm sure it could be converted to Scilab without big trouble. Maybe >>> I even have an old one laying around (??) somewhere ... maybe. >> Yes indeed, here : http://www-sop.inria.fr/coprin/logiciels/Int4Sci/ > > The thing is I don't see how to use these tools to do what Yoahn wants > to do. > I think I have the same kind of issues than Yoahn (I've done a fit, > how good is it as compare to another one, and what kind of > "confidence" or "error bar" should I associate to each parameter). > Could you tell us how we can use interval analysis to adress this > (just curious, it's a bit far from my field of expertise)? > > Cheers, > > Antoine > > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > . > --- This email is free from viruses and malware because avast! Antivirus protection is active. http://www.avast.com From sgougeon at free.fr Fri Feb 21 21:28:23 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 21 Feb 2014 21:28:23 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <53076B6F.1070306@laas.fr> References: <1392732408401-4028696.post@n3.nabble.com> <5304DC87.4020507@gmail.com> <53050C29.1000308@free.fr> <53076B6F.1070306@laas.fr> Message-ID: <5307B6E7.1010103@free.fr> Le 21/02/2014 16:06, Antoine Monmayrant a ?crit : > On 02/19/2014 08:55 PM, Samuel Gougeon wrote: >> Le 19/02/2014 17:32, Claus Futtrup a ?crit : >>> Hi Yohan >>> >>> I see that I'm not alone wondering a bit about what you'd like to do. >>> >>> A mathematically solid and readily available technique to corner >>> mathematical rounding errors and other errors of operation is called >>> Interval Analysis. >>> >>> For matlab there's INTLAB. See: >>> http://www.ti3.tu-harburg.de/~rump/intlab/ >>> >>> I'm sure it could be converted to Scilab without big trouble. Maybe >>> I even have an old one laying around (??) somewhere ... maybe. >> Yes indeed, here : http://www-sop.inria.fr/coprin/logiciels/Int4Sci/ > > The thing is I don't see how to use these tools to do what Yoahn wants > to do. Neither do i. I just got the opportunity to say that it would be nice to have this external toolbox in ATOMS... Samuel From sgougeon at free.fr Fri Feb 21 21:46:07 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 21 Feb 2014 21:46:07 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <1392816677490-4028742.post@n3.nabble.com> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> <1392816677490-4028742.post@n3.nabble.com> Message-ID: <5307BB0F.3030100@free.fr> Le 19/02/2014 14:31, Yohann a ?crit : > Hi Antoine, > thank you for your answer but > what I need is a confidence interval on each parameter ! A raw empirical approach: Let R be the RMS residue for the best set P of fitting parameters. * change the value of P(1) by -- say -- +1%, for the parameter #1. * calculate the new residue R1p, and the change dR1p = abs(R-R1p) you may symmetrize the process : change P(1) => P(1)*0.99 (-1%) calculate R1m with this biased fitting model, and get dR1=dR1m+dR1p * These dR somewhat represent sensitivities of the fit w.r.t. to each parameter. * Assume that the confidence interval for parameter # is proportional to 1/dR# For "stiff" parameters, dR# will be big (easily divergent), and 1/dR# small. DISCLAIMER : You may use this approach under your own responsability ;) Its main (may be unique) merit is to be simple. Samuel From sgougeon at free.fr Fri Feb 21 21:51:25 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 21 Feb 2014 21:51:25 +0100 Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <5307BB0F.3030100@free.fr> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> <1392816677490-4028742.post@n3.nabble.com> <5307BB0F.3030100@free.fr> Message-ID: <5307BC4D.9030201@free.fr> This makes me also thinking about /Design Of Experiment/. You may get useful ressources about DOE through the module http://atoms.scilab.org/toolboxes/scidoe -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.carrico at free.fr Fri Feb 21 22:36:54 2014 From: paul.carrico at free.fr (Paul CARRICO) Date: Fri, 21 Feb 2014 22:36:54 +0100 Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <1392959894551-4028792.post@n3.nabble.com> References: <1392959894551-4028792.post@n3.nabble.com> Message-ID: <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> Did you have a look to "awk" and "sed" tools that can be used in conjunction with Scilab ? -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de samaelkreutz Envoy??: vendredi 21 f?vrier 2014 06:18 ??: users at lists.scilab.org Objet?: [Scilab-users] finding similar values with scilab Hello!!! I have a big question ! I have a biiiiiiiiiiig data archive with text files with the form: 18.87 2.6 0.00545558 19.98 2.6 0.00225349 18.87 2.6 0.00405905 13.32 2.6 0.01338288 19.98 2.6 0.01537532 18.87 2.6 0.00481375 19.98 2.6 0.00936207 12.21 2.6 0.00558517 I need eliminate extreme values, for example: If I plot a random text file I obtain a image like this: As you can see theres points extremes. If i make a zoom DN=find(DM(:,3)<=0.0001 & DM(:,3)>=0.00002 ) dn=DM(DN,:) figure(2) plot(dn(:,1),dn(:,3),'O') I have this picture My question is: Has scilab a command that can help me to find the "dense data" , I mean find all the values that are similar. I ask this cause I have about 20,000 text files and is imposible for me make "artesal zoms". Please any suggestion, answer going to help me. -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792. html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com From Christophe.Dang at sidel.com Thu Feb 20 09:47:48 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Thu, 20 Feb 2014 09:47:48 +0100 Subject: [Scilab-users] To find maximum value in a random curve In-Reply-To: <1392805906969-4028741.post@n3.nabble.com> References: <1392805906969-4028741.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F3267F1160@301EX00100.sidel.com> Hello, > De la part de bomber1978 > Envoy? : mercredi 19 f?vrier 2014 11:32 > > I have a random 2D-curve which is defined by a set of points > (I don't know the equation of the curve) [...] > How can I find the maximum, minimum, & inflexion points on this curve. This is more a methodology problem than a Scilab problem. You don't have a model to fit your data, so you have two solutions: perform a non-parametric analysis, e.g. the min() and max() functions give you the extrema of the values you have, but you don't have things like inflection point. smooth the data, extract the second derivative and detect when it goes through 0. You can look at the Savitzky-Golay algorithm, which gives the derivatives. You have an example of its implementation here when the points are evenly spaced : http://commons.wikimedia.org/wiki/File:Savitzky-golay_pic_gaussien_bruite.svg If the point sare not evenly spaced, you then have to perform a polynomial regression on a sliding interval "manually". Hope this helps. -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From mariajovera at icloud.com Sat Feb 22 00:25:59 2014 From: mariajovera at icloud.com (samaelkreutz) Date: Fri, 21 Feb 2014 15:25:59 -0800 (PST) Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> References: <1392959894551-4028792.post@n3.nabble.com> <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> Message-ID: <1393025159051-4028865.post@n3.nabble.com> I need compute an iteration , but (obviously I don't know how to do it =/ ) Read thousand of text files... then what I do is check the mean, stand. deviation, and variation coefficient ( std/mean ). If the variation coefficient is <= 0.4, stop the process. If not, eliminate the big value and recalculate all the process again. But i'm lost!!! because Im reading multiple text files and is confusing!!!!! Heeelp!!!! clc clear all z=[]; // genero mi matriz vacia que se ira llenando con el for for i=103 //35:2:37// for j=2.6:0.1:3.9// DM = fscanfMat(msprintf("new-C1_%d_%3.1f.txt",i,j)); //disp(i,j) med=mean(DM(:,3)); std=st_deviation(DM(:,3)); CV=std/abs(med); Z=[med std CV CV*100]; z=[z;Z]; // Matriz con informaci?n requerida //figure(i); //plot(DM(:,1), DM(:,3), 'ro*') end end z -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792p4028865.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.carrico at free.fr Sat Feb 22 08:48:28 2014 From: paul.carrico at free.fr (Paul CARRICO) Date: Sat, 22 Feb 2014 08:48:28 +0100 Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <1393025159051-4028865.post@n3.nabble.com> References: <1392959894551-4028792.post@n3.nabble.com> <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> <1393025159051-4028865.post@n3.nabble.com> Message-ID: <000301cf2fa2$7a713280$6f539780$@carrico@free.fr> Do you have to do this on a single file and to repeat it X times (where X is the thousands of file you're speaking) or do you have to calculate the deviation and so on) on the X files at the same time ? Obviously the first case is much easier than the second one ..... -----Message d'origine----- De : users [mailto:users-bounces at lists.scilab.org] De la part de samaelkreutz Envoy? : samedi 22 f?vrier 2014 00:26 ? : users at lists.scilab.org Objet : Re: [Scilab-users] finding similar values with scilab I need compute an iteration , but (obviously I don't know how to do it =/ ) Read thousand of text files... then what I do is check the mean, stand. deviation, and variation coefficient ( std/mean ). If the variation coefficient is <= 0.4, stop the process. If not, eliminate the big value and recalculate all the process again. But i'm lost!!! because Im reading multiple text files and is confusing!!!!! Heeelp!!!! clc clear all z=[]; // genero mi matriz vacia que se ira llenando con el for for i=103 //35:2:37// for j=2.6:0.1:3.9// DM = fscanfMat(msprintf("new-C1_%d_%3.1f.txt",i,j)); //disp(i,j) med=mean(DM(:,3)); std=st_deviation(DM(:,3)); CV=std/abs(med); Z=[med std CV CV*100]; z=[z;Z]; // Matriz con informaci?n requerida //figure(i); //plot(DM(:,1), DM(:,3), 'ro*') end end z -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792p4028865.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com From paul.carrico at free.fr Sat Feb 22 09:04:28 2014 From: paul.carrico at free.fr (Paul CARRICO) Date: Sat, 22 Feb 2014 09:04:28 +0100 Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <000301cf2fa2$7a713280$6f539780$@carrico@free.fr> References: <1392959894551-4028792.post@n3.nabble.com> <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> <1393025159051-4028865.post@n3.nabble.com> <000301cf2fa2$7a713280$6f539780$@carrico@free.fr> Message-ID: <000001cf2fa4$b690aaf0$23b200d0$@carrico@free.fr> I've ever suggested to use "awk" / "sed" and other linux tool (even under Windows) in conjunction to Scilab ; 1rst step : concatenate all the files if necessary (with Scilab for example) Average : http://azimuth.biz/2009/03/30/calculate-average-using-awk/ standard deviation : http://azimuth.biz/2010/02/25/calculate-standard-deviation-using-awk/ Min-Max: http://azimuth.biz/2010/10/20/find-minimum-and-maximum-using-awk/ supress values greater than : (ditto for min value) http://www.unix.com/shell-programming-scripting/178904-awk-show-lines-where-nth-column-greater-than-some-number.html http://oreilly.com/catalog/unixnut3/chapter/ch11.html and so on I'm not a awk/sed specialist, but I've a look to such tools when my input files have million's of lines ... and it works fine and fast !!! Just an idea NB: tutos http://www.grymoire.com/Unix/Awk.html http://www.cs.unibo.it/~renzo/doc/awk/nawkA4.pdf etc. ... -----Message d'origine----- De : users [mailto:users-bounces at lists.scilab.org] De la part de Paul CARRICO Envoy? : samedi 22 f?vrier 2014 08:48 ? : 'International users mailing list for Scilab.' Objet : Re: [Scilab-users] finding similar values with scilab Do you have to do this on a single file and to repeat it X times (where X is the thousands of file you're speaking) or do you have to calculate the deviation and so on) on the X files at the same time ? Obviously the first case is much easier than the second one ..... -----Message d'origine----- De : users [mailto:users-bounces at lists.scilab.org] De la part de samaelkreutz Envoy? : samedi 22 f?vrier 2014 00:26 ? : users at lists.scilab.org Objet : Re: [Scilab-users] finding similar values with scilab I need compute an iteration , but (obviously I don't know how to do it =/ ) Read thousand of text files... then what I do is check the mean, stand. deviation, and variation coefficient ( std/mean ). If the variation coefficient is <= 0.4, stop the process. If not, eliminate the big value and recalculate all the process again. But i'm lost!!! because Im reading multiple text files and is confusing!!!!! Heeelp!!!! clc clear all z=[]; // genero mi matriz vacia que se ira llenando con el for for i=103 //35:2:37// for j=2.6:0.1:3.9// DM = fscanfMat(msprintf("new-C1_%d_%3.1f.txt",i,j)); //disp(i,j) med=mean(DM(:,3)); std=st_deviation(DM(:,3)); CV=std/abs(med); Z=[med std CV CV*100]; z=[z;Z]; // Matriz con informaci?n requerida //figure(i); //plot(DM(:,1), DM(:,3), 'ro*') end end z -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792p4028865.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users --- Ce courrier ?lectronique ne contient aucun virus ou logiciel malveillant parce que la protection avast! Antivirus est active. http://www.avast.com From mariajovera at icloud.com Sat Feb 22 15:01:07 2014 From: mariajovera at icloud.com (samaelkreutz) Date: Sat, 22 Feb 2014 06:01:07 -0800 (PST) Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <000001cf2fa4$b690aaf0$23b200d0$@carrico@free.fr> References: <1392959894551-4028792.post@n3.nabble.com> <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> <1393025159051-4028865.post@n3.nabble.com> <000301cf2fa2$7a713280$6f539780$@carrico@free.fr> <000001cf2fa4$b690aaf0$23b200d0$@carrico@free.fr> Message-ID: <1393077667016-4028875.post@n3.nabble.com> Is the second option... with a for cycle im calling the files... I have about 2,000 files each one with hundred of lines. And yep... I manipulate the files with awk, at the beginning I was programming with bash the process becomes slow and two... (it takes about half hour, ok is not so much... but for me it is). And second... I need eliminate the outliers, because if a take a random file and compute the mean and std, the mean is so sensitive to extreme values and sometimes those values aren representative. I have to eliminate them and recalculate the process again without that outlier and test again... is coefficient std/mean less that 0.4 ?? if not, find again other outlier, eliminate and repeat the process again and so on until the coefficient is 0.4 or less. Save the information required (mean, std, coefficient,) into matrix and then take the following text file repeat the process before and save the results in the matrix. At the end of the process I pretend obtain a representative matrix with the information of all my text files. As a user posted here, I think the modified Thompson Tau method sounds great... I think I'm going to try this... seriously I never though that statistics were useful until my thesis obviously... =/ -- View this message in context: http://mailinglists.scilab.org/finding-similar-values-with-scilab-tp4028792p4028875.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From yohann.morille at univ-angers.fr Fri Feb 21 10:35:30 2014 From: yohann.morille at univ-angers.fr (Yohann) Date: Fri, 21 Feb 2014 01:35:30 -0800 (PST) Subject: [Scilab-users] evaluate error on each parameter calculated with leastsq In-Reply-To: <53047A24.90004@laas.fr> References: <1392732408401-4028696.post@n3.nabble.com> <53047A24.90004@laas.fr> Message-ID: <1392975330818-4028795.post@n3.nabble.com> Contrary to what I understood at the first reading of Antoine's respons and what can do Bootstrapping, it seems that this method is well suited to do the job. http://en.wikipedia.org/wiki/Bootstrapping_%28statistics%29#Resampling_residuals If I have well understood this link, after applying this algorithm you obtain a large number of occurrences (Nboot) of each one of your parameter (C1 C2 and C3 in my case) and then you can determine a value and a confidence interval: C1 = mean(C1boot)+/- std(C1boot)*1.95996 C2 = mean(C2boot)+/- std(C2boot)*1.95996 C3 = mean(C3boot)+/- std(C3boot)*1.95996 I think 1.95996 is for a confidence interval of 97.5% (from Student law), but not sure... If it is good, it looks easy to implement and understand. -- View this message in context: http://mailinglists.scilab.org/evaluate-error-on-each-parameter-calculated-with-leastsq-tp4028696p4028795.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From aweeks at hidglobal.com Mon Feb 24 08:58:23 2014 From: aweeks at hidglobal.com (aweeks at hidglobal.com) Date: Mon, 24 Feb 2014 07:58:23 +0000 Subject: [Scilab-users] Linking DLL from MS Visual C++ to SciLab. In-Reply-To: <1393007496941-4028813.post@n3.nabble.com> References: <1393007496941-4028813.post@n3.nabble.com> Message-ID: Thanks Chuox, I'll try that and see what happens. Adrian Weeks Development Engineer HID Global, 3 Cae Gwyrdd, Green Meadow Springs Business Park, Cardiff CF15 7AB, United Kingdom. +44 (0)29 20528523 (Office) aweeks at hidglobal.com www.hidglobal.com -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 13483180.gif Type: image/gif Size: 4559 bytes Desc: not available URL: From Christophe.Dang at sidel.com Mon Feb 24 10:36:16 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 24 Feb 2014 10:36:16 +0100 Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <1393077667016-4028875.post@n3.nabble.com> References: <1392959894551-4028792.post@n3.nabble.com> <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> <1393025159051-4028865.post@n3.nabble.com> <000301cf2fa2$7a713280$6f539780$@carrico@free.fr> <000001cf2fa4$b690aaf0$23b200d0$@carrico@free.fr> <1393077667016-4028875.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32689F5CB@301EX00100.sidel.com> Hello, > De la part de samaelkreutz > Envoy? : samedi 22 f?vrier 2014 15:01 > > I need eliminate the outliers, because if a take a random file and > compute the mean and std, the mean is so sensitive to extreme values > and sometimes those values aren representative. Another approach is to use another position criterion that is less sensitive to outliers, e.g. the median. Best regards. -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From Christophe.Dang at sidel.com Mon Feb 24 10:46:01 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 24 Feb 2014 10:46:01 +0100 Subject: [Scilab-users] help, please In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32689F60D@301EX00100.sidel.com> Hello, > De la part de Tam?s Reiter > Envoy? : jeudi 20 f?vrier 2014 09:58 > > I keep getting an error message:Submatrix incorrecly defined... > [...] > T(:,j)=t(:,:);--THIS IS WHERE I AM GETTING THIS ERROR > P(:,j)=p(:,:);--THIS IS WHERE I AM GETTING THIS ERROR I'm not sure about what you try to do (I couldn't run your code because we don't have all the values), but the overall syntax is correct: A = [1 1 ; 1 1] b = |2 ; 2] A(:, 1) = b(:, :) works well. Mind that t(:, :) is the same as t, so you can simplify your syntax. Now, t must be a column vector. Can you check it is before the assignment? Hope this helps -- Christophe Dang Ngoc Chan Mechanical calculation engineer ______________________________________________________________________ This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ______________________________________________________________________ From cfuttrup at gmail.com Mon Feb 24 11:52:14 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Mon, 24 Feb 2014 11:52:14 +0100 Subject: [Scilab-users] finding similar values with scilab In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F32689F5CB@301EX00100.sidel.com> References: <1392959894551-4028792.post@n3.nabble.com> <001801cf2f4d$0b2eafb0$218c0f10$@carrico@free.fr> <1393025159051-4028865.post@n3.nabble.com> <000301cf2fa2$7a713280$6f539780$@carrico@free.fr> <000001cf2fa4$b690aaf0$23b200d0$@carrico@free.fr> <1393077667016-4028875.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32689F5CB@301EX00100.sidel.com> Message-ID: <530B245E.20401@gmail.com> On 2/24/2014 10:36, Dang, Christophe wrote: > Hello, > >> De la part de samaelkreutz >> Envoy? : samedi 22 f?vrier 2014 15:01 >> >> I need eliminate the outliers, because if a take a random file and >> compute the mean and std, the mean is so sensitive to extreme values >> and sometimes those values aren representative. > Another approach is to use another position criterion that is less > sensitive to outliers, e.g. the median. > > Best regards. > I agree with Christophe ... it sounds like you should take a look at median filtering. /Claus --- This email is free from viruses and malware because avast! Antivirus protection is active. http://www.avast.com From W.Schulz at ove.at Tue Feb 25 12:55:51 2014 From: W.Schulz at ove.at (SCHULZ Wolfgang) Date: Tue, 25 Feb 2014 11:55:51 +0000 Subject: [Scilab-users] Plotting density values Message-ID: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> Hello, I'm asking for an intelligent method to solve the following problem without a for loop. I have 2 vectors x an y with e.g. 3000 values. The range of x values is between 0 and 360 and the range of y values is between -15 and 15. I would like to plot the density of individual grid cells for e.g. dx=10, dy=1 and therefore I would need a matrix with 36x30 values where each value corresponds to the number of data points in the cell. Any idea who I can do that without a for loop? Thanks a lot Wolfgang From paul.bignier at scilab-enterprises.com Tue Feb 25 13:54:17 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Tue, 25 Feb 2014 13:54:17 +0100 Subject: [Scilab-users] Plotting density values In-Reply-To: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> References: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> Message-ID: <530C9279.2010206@scilab-enterprises.com> Hi, It sounds to me like you are trying to use histplot (y,x)? Regards, Paul On 02/25/2014 12:55 PM, SCHULZ Wolfgang wrote: > Hello, > I'm asking for an intelligent method to solve the following problem without a for loop. > > I have 2 vectors x an y with e.g. 3000 values. The range of x values is between 0 and 360 and the range of y values is between -15 and 15. > I would like to plot the density of individual grid cells for e.g. dx=10, dy=1 and therefore I would need a matrix with 36x30 values where each value corresponds to the number of data points in the cell. > > Any idea who I can do that without a for loop? > > Thanks a lot > Wolfgang > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Development engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From Serge.Steer at inria.fr Tue Feb 25 15:57:20 2014 From: Serge.Steer at inria.fr (Serge Steer) Date: Tue, 25 Feb 2014 15:57:20 +0100 Subject: [Scilab-users] Plotting density values In-Reply-To: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> References: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> Message-ID: <530CAF50.50001@inria.fr> The following function does the job but with one loop. It is based on the 1D classification function dsearch Serg Steer function occ=dsearch2D(xy,dx,dy,opt) //x,y coordinates of the points of the cloud //dx, discretization of the x axis //dy , discretization of the y axis //occ, table such as occ(i,j) contains the number of points in the //"pixel" [dx(i) dx(i+1)) x [dy(j) dy(j+1)) if argn(2)<4 then opt='c';end if and(opt<>["c","d"]) then error(msprintf(_("%: unknown char specifier (must be ''c'' or ''d'')\n"),"dsearch2D")) end if size(xy,2)<>2 then error("Wrong dimension for first argument") end [indx,occx]=dsearch(xy(:,1),dx,opt); occ=[]; dy=matrix(dy,1,-1); for k=1:length(dx)-1 i=find(indx==k); [indy,occy]=dsearch(xy(i,2),dy,opt); occ=[occ; occy]; end endfunction Le 25/02/2014 12:55, SCHULZ Wolfgang a ?crit : > Hello, > I'm asking for an intelligent method to solve the following problem without a for loop. > > I have 2 vectors x an y with e.g. 3000 values. The range of x values is between 0 and 360 and the range of y values is between -15 and 15. > I would like to plot the density of individual grid cells for e.g. dx=10, dy=1 and therefore I would need a matrix with 36x30 values where each value corresponds to the number of data points in the cell. > > Any idea who I can do that without a for loop? > > Thanks a lot > Wolfgang > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From romain.desbats at ifsttar.fr Tue Feb 25 16:15:09 2014 From: romain.desbats at ifsttar.fr (Romain Desbats) Date: Tue, 25 Feb 2014 16:15:09 +0100 (CET) Subject: [Scilab-users] Increment of parameters in optim function In-Reply-To: <1984960779.3708825.1393340406433.JavaMail.root@ifsttar.fr> Message-ID: <879307472.3717216.1393341309894.JavaMail.root@ifsttar.fr> Dear Scilab Community, I am using the optim function to optimise a two-variable problem (called x1 and x2 ). During the optimisation I monitor the values of those variables as well as the value of the cost function to minimise. I find that after 36 iterations that the values of the parameters "did not change much": iteration 1 : x1=765 x2=12900 f=1.402D+08 iteration 36: x1=764.82264 x2=12900 f=1.400D+08 Is there a way to force the increment value of the parameters that are being optimised? For instance x1 could only be equal to 765, 764, 763... My question is similar to this one but I did not find further discussion: http://lists.scilab.org/pipermail/bugzilla/2010-February/000638.html Thanks a lot for your help. Romain -------------- next part -------------- An HTML attachment was scrubbed... URL: From denis.crete at thalesgroup.com Tue Feb 25 16:24:19 2014 From: denis.crete at thalesgroup.com (CRETE Denis) Date: Tue, 25 Feb 2014 16:24:19 +0100 Subject: [Scilab-users] Plotting density values In-Reply-To: <530C9279.2010206@scilab-enterprises.com> References: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> <530C9279.2010206@scilab-enterprises.com> Message-ID: <2438_1393341860_530CB5A4_2438_16838_1_8F1D40232A0E68409E3FC23A30C3266201690CBE5E41@THSONEA01CMS04P.one.grp> It looks like you wanted to suggest hist3d Denis De : users [mailto:users-bounces at lists.scilab.org] De la part de Paul Bignier Envoy? : mardi 25 f?vrier 2014 13:54 ? : International users mailing list for Scilab.; W.Schulz at ove.at Objet : Re: [Scilab-users] Plotting density values Hi, It sounds to me like you are trying to use histplot(y,x)? Regards, Paul On 02/25/2014 12:55 PM, SCHULZ Wolfgang wrote: Hello, I'm asking for an intelligent method to solve the following problem without a for loop. I have 2 vectors x an y with e.g. 3000 values. The range of x values is between 0 and 360 and the range of y values is between -15 and 15. I would like to plot the density of individual grid cells for e.g. dx=10, dy=1 and therefore I would need a matrix with 36x30 values where each value corresponds to the number of data points in the cell. Any idea who I can do that without a for loop? Thanks a lot Wolfgang _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Development engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From romain.desbats at ifsttar.fr Tue Feb 25 17:52:39 2014 From: romain.desbats at ifsttar.fr (R0M1) Date: Tue, 25 Feb 2014 08:52:39 -0800 (PST) Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> References: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> Message-ID: <1393347159977-4028965.post@n3.nabble.com> Dear Users, Trying different compilation methods for this model, I found some discrepancies in the results that I do not explain so far. With the same .xcos model and .csv input file as in the first post, I compile using /xcos_simulate/ and /scicos_simulate/: The plots shows the difference between the results from the two methods. Any idea on what is causing this ? Thanks a lot. Here is the code: clear // parameters data=csvRead('input.csv',';'); input_current=evstr(data(:,1)); size1=length(input_current); C=1/12900; R=765; T=0.0001; // timestep time=[0:T:(size1-1)*T]'; // in s end_time=time($)+T; // end time in s // loading Xcos functions loadScicos; loadXcosLibs(); exec(loadpallibs,-1); // loading model: creates scs_m importXcosDiagram('graphic_model.xcos') // first compilation method scs_m.props.context=strcat(["C=",string(C),"; R=",string(R),"; T=",string(T)"; end_time=",string(end_time)]); xcos_simulate(scs_m,4) // '4' is to perform compilation graph_voltage_1=graph_voltage.values; // output variable // second compilation method %scicos_context.C=C; %scicos_context.R=R; %scicos_context.T=T; %scicos_context.end_time=end_time; Info=scicos_simulate(scs_m,list(),%scicos_context,flag='nw'); graph_voltage_2=graph_voltage.values; // output variable set("current_figure",10) plot(time,graph_voltage_1,'-b') plot(time,graph_voltage_2,'-r') xtitle('graph title','time (s)','voltage (V)'); hl=legend(['xcos_simulate';'scicos_simulate'],"in_upper_right"); -- View this message in context: http://mailinglists.scilab.org/Scilab-users-Problem-running-xcos-in-batch-mode-from-Scilab-tp4028765p4028965.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From contact at pierre-vuillemin.fr Tue Feb 25 19:22:50 2014 From: contact at pierre-vuillemin.fr (Pierre Vuillemin) Date: Tue, 25 Feb 2014 19:22:50 +0100 Subject: [Scilab-users] Increment of parameters in optim function In-Reply-To: <879307472.3717216.1393341309894.JavaMail.root@ifsttar.fr> References: <879307472.3717216.1393341309894.JavaMail.root@ifsttar.fr> Message-ID: <1393352570.3923.32.camel@Pelisse.lan> Dear Romain, Just some hints concerning one of the underlying algorithm of optim, the "quasi newton" one with a Wolfe type line search : this algorithm is a descent algorithm, i.e given an initial point, it will find a local minimun of your objective function "near" (in some sense) to the initial point. The main ideas are the following : - at each iteration, the gradient of the objective function at the current position x is computed (and the approximation of the Hessian which explains the name of "quasi Newton" method). The descent direction, ie the direction which decreases the function, is given by the opposite of the gradient, hence the algorithm has to go this way. But how far? The length of the step in that direction is given by a line search algorithm - in this case, it is a Wolfe type line search which means that the step length is choosen so that the next point ensures (i) a sufficient decrease in the objective function AND (ii) leads to a sufficient decrease in the slope of the function (some variations exist between the Wolfe and strong Wolfe conditions). The algorithm stops when a local minimum is reached. Concerning your question, you cannot force the algorithm to move further between two iterations. This means that in your case, your initial point may be close to a local minimum. Some hints : - since you only have 2 optimization parameters, you can plot your objective function to see where are the local minima, or to give you some indications on what the initial point should be. - try different initial points to (hopefully) find different local minima and choose the best one. - you may want to try evolutionary algorithms (genetic algorithms or particle swarm optimization) which do not get stuck in local minima. I think that there is a built-in function or an atom module containing such algorithm. (- there exist plenty of other optimization algorithms which may (or may not) be better suited for your problem depending on its structure) Hope this helps, Best regards, Pierre Le mardi 25 f?vrier 2014 ? 16:15 +0100, Romain Desbats a ?crit : > Dear Scilab Community, > > I am using the optim function to optimise a two-variable problem > (called x1 and x2). During the optimisation I monitor the values of > those variables as well as the value of the cost function to minimise. > > I find that after 36 iterations that the values of the parameters "did > not change much": > > iteration 1 : x1=765 x2=12900 f=1.402D+08 > iteration 36: x1=764.82264 x2=12900 f=1.400D+08 > > Is there a way to force the increment value of the parameters that are > being optimised? For instance x1 could only be equal to 765, 764, > 763... > > My question is similar to this one but I did not find further > discussion: > http://lists.scilab.org/pipermail/bugzilla/2010-February/000638.html > > Thanks a lot for your help. > > Romain > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From n_hizume at nifty.com Wed Feb 26 07:08:59 2014 From: n_hizume at nifty.com (Nobuo HIZUME) Date: Wed, 26 Feb 2014 15:08:59 +0900 Subject: [Scilab-users] Algebraic loop error in Xcos step response Message-ID: Dear Xcos users who are interested in step response: I have recently tried step response of dc motor circuit in Xcos. But to get the attached result, slight modification was necessary because original circuit can?t be simulated showing ?Algebraic loop? error . I needed to change one of the Gains(value=22) into 1st order delay(gain=22, time constant=0.1ms). Regarding detail refer to the attached ZIP file, which contains Scilab step response too, and the results are the same. I?d like to know what ?Algebraic loop? error means, and is there another method to avoid the error? Regards, Nobuo HIZUME -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dc_motor_circuit.ZIP Type: application/x-zip-compressed Size: 24366 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dcmotor_simulation_by xcos.pdf Type: application/pdf Size: 13759 bytes Desc: not available URL: From paul.bignier at scilab-enterprises.com Wed Feb 26 09:57:13 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Wed, 26 Feb 2014 09:57:13 +0100 Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <1393347159977-4028965.post@n3.nabble.com> References: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> <1393347159977-4028965.post@n3.nabble.com> Message-ID: <530DAC69.20908@scilab-enterprises.com> Hello Romain, The difference is due to the solver type. Your diagram currently selects the ODE solver Sundials/CVode - BDF - Newton. "xcos_simulate(scs_m, 4);" detects two continuous states and switches to a DAE solver (Sundials/IDA). "scicos_simulate(scs_m);" only detects one continuous state and keeps the ODE solver. This must be due to a Modelica compilation difference between xcos and scicos. I don't know the equations of your model so I cannot assess which method is correct. Also, I noticed that the context variables (C, R, T, end_time) are not present in the diagram, so what are you trying to change when you modify the context? You might want ot use the END block to set the end time. Hope this helps, Regards, Paul On 02/25/2014 05:52 PM, R0M1 wrote: > Dear Users, > > Trying different compilation methods for this model, I found some > discrepancies in the results that I do not explain so far. > > With the same .xcos model and .csv input file as in the first post, I > compile using /xcos_simulate/ and /scicos_simulate/: > > The plots shows the difference between the results from the two methods. > > Any idea on what is causing this ? > > Thanks a lot. > > Here is the code: > > clear > > // parameters > data=csvRead('input.csv',';'); > input_current=evstr(data(:,1)); > size1=length(input_current); > > C=1/12900; > R=765; > > T=0.0001; // timestep > time=[0:T:(size1-1)*T]'; // in s > end_time=time($)+T; // end time in s > > // loading Xcos functions > loadScicos; > loadXcosLibs(); > exec(loadpallibs,-1); > > // loading model: creates scs_m > importXcosDiagram('graphic_model.xcos') > > // first compilation method > scs_m.props.context=strcat(["C=",string(C),"; R=",string(R),"; > T=",string(T)"; end_time=",string(end_time)]); > xcos_simulate(scs_m,4) // '4' is to perform compilation > > graph_voltage_1=graph_voltage.values; // output variable > > // second compilation method > %scicos_context.C=C; > %scicos_context.R=R; > %scicos_context.T=T; > %scicos_context.end_time=end_time; > Info=scicos_simulate(scs_m,list(),%scicos_context,flag='nw'); > > graph_voltage_2=graph_voltage.values; // output variable > > set("current_figure",10) > plot(time,graph_voltage_1,'-b') > plot(time,graph_voltage_2,'-r') > xtitle('graph title','time (s)','voltage (V)'); > hl=legend(['xcos_simulate';'scicos_simulate'],"in_upper_right"); > > > > > -- > View this message in context: http://mailinglists.scilab.org/Scilab-users-Problem-running-xcos-in-batch-mode-from-Scilab-tp4028765p4028965.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Development engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From romain.desbats at ifsttar.fr Wed Feb 26 10:21:30 2014 From: romain.desbats at ifsttar.fr (R0M1) Date: Wed, 26 Feb 2014 01:21:30 -0800 (PST) Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <530DAC69.20908@scilab-enterprises.com> References: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> <1393347159977-4028965.post@n3.nabble.com> <530DAC69.20908@scilab-enterprises.com> Message-ID: <1393406490356-4028993.post@n3.nabble.com> Hi Paul, Thank you very much for your answer. I have changed the model after your first answer but I forgot to attach it again. I did use the END block. Thanks for highlighting the difference of solver type. It seems that Sundials/IDA gives the result I expect. Can you tell me where to check the solver that have been used for the simulation? (The only variable generated by running the simulation is %scicos_prob but I can not open it). -- View this message in context: http://mailinglists.scilab.org/Scilab-users-Problem-running-xcos-in-batch-mode-from-Scilab-tp4028765p4028993.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.bignier at scilab-enterprises.com Wed Feb 26 10:40:09 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Wed, 26 Feb 2014 10:40:09 +0100 Subject: [Scilab-users] Problem running xcos in batch mode from Scilab In-Reply-To: <1393406490356-4028993.post@n3.nabble.com> References: <1313506737.1940942.1392806080120.JavaMail.root@ifsttar.fr> <1393347159977-4028965.post@n3.nabble.com> <530DAC69.20908@scilab-enterprises.com> <1393406490356-4028993.post@n3.nabble.com> Message-ID: <530DB679.4090404@scilab-enterprises.com> During the compilation, a WARNING is displayed to notify the solver switch. xcos_simulate(): WARNING: Diagram has been compiled for implicit solver switching to implicit Solver scicos_simualte(): no warning, so no switching. If the diagram contains continuous states, xcos and scicos are supposed to detect the problem type (ODE or DAE) and switch to the default one (Sundials/BDF-Newton for ODE and Sundials/IDA for DAE) if the Setup solver is not adapted. On 02/26/2014 10:21 AM, R0M1 wrote: > Hi Paul, > > Thank you very much for your answer. I have changed the model after your > first answer but I forgot to attach it again. I did use the END block. > > Thanks for highlighting the difference of solver type. It seems that > Sundials/IDA gives the result I expect. Can you tell me where to check the > solver that have been used for the simulation? (The only variable generated > by running the simulation is %scicos_prob but I can not open it). > > > > -- > View this message in context: http://mailinglists.scilab.org/Scilab-users-Problem-running-xcos-in-batch-mode-from-Scilab-tp4028765p4028993.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Development engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com From stephane.mottelet at gmail.com Wed Feb 26 11:39:28 2014 From: stephane.mottelet at gmail.com (=?UTF-8?B?U3TDqXBoYW5lIE1vdHRlbGV0?=) Date: Wed, 26 Feb 2014 11:39:28 +0100 Subject: [Scilab-users] Test In-Reply-To: <1393352570.3923.32.camel@Pelisse.lan> References: <879307472.3717216.1393341309894.JavaMail.root@ifsttar.fr> <1393352570.3923.32.camel@Pelisse.lan> Message-ID: <530DC460.5050800@gmail.com> Sorry, this is a test... S. From W.Schulz at ove.at Wed Feb 26 15:44:01 2014 From: W.Schulz at ove.at (SCHULZ Wolfgang) Date: Wed, 26 Feb 2014 14:44:01 +0000 Subject: [Scilab-users] Plotting density values In-Reply-To: <530CAF50.50001@inria.fr> References: <07A8AF089844A549B0F0CCC727AF31AA56F95F7D@OVEKBG.ove.at> <530CAF50.50001@inria.fr> Message-ID: <07A8AF089844A549B0F0CCC727AF31AA56F968D6@OVEKBG.ove.at> Wow - this was exactly what I was searching for. Thanks a lot Wolfgang > -----Urspr?ngliche Nachricht----- > Von: users [mailto:users-bounces at lists.scilab.org] Im Auftrag von Serge Steer > Gesendet: Dienstag, 25. Februar 2014 15:57 > An: users at lists.scilab.org > Betreff: Re: [Scilab-users] Plotting density values > > The following function does the job but with one loop. It is based on the 1D > classification function dsearch Serg Steer > > > function occ=dsearch2D(xy,dx,dy,opt) > //x,y coordinates of the points of the cloud > //dx, discretization of the x axis > //dy , discretization of the y axis > //occ, table such as occ(i,j) contains the number of points in the > //"pixel" [dx(i) dx(i+1)) x [dy(j) dy(j+1)) > if argn(2)<4 then opt='c';end > if and(opt<>["c","d"]) then > error(msprintf(_("%: unknown char specifier (must be ''c'' or > ''d'')\n"),"dsearch2D")) > end > if size(xy,2)<>2 then > error("Wrong dimension for first argument") > end > [indx,occx]=dsearch(xy(:,1),dx,opt); > occ=[]; > dy=matrix(dy,1,-1); > for k=1:length(dx)-1 > i=find(indx==k); > [indy,occy]=dsearch(xy(i,2),dy,opt); > occ=[occ; occy]; > end > endfunction > > > > > Le 25/02/2014 12:55, SCHULZ Wolfgang a ?crit : > > Hello, > > I'm asking for an intelligent method to solve the following problem without a > for loop. > > > > I have 2 vectors x an y with e.g. 3000 values. The range of x values is between > 0 and 360 and the range of y values is between -15 and 15. > > I would like to plot the density of individual grid cells for e.g. dx=10, dy=1 and > therefore I would need a matrix with 36x30 values where each value > corresponds to the number of data points in the cell. > > > > Any idea who I can do that without a for loop? > > > > Thanks a lot > > Wolfgang > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From W.Schulz at ove.at Wed Feb 26 15:56:25 2014 From: W.Schulz at ove.at (SCHULZ Wolfgang) Date: Wed, 26 Feb 2014 14:56:25 +0000 Subject: [Scilab-users] Bug? Message-ID: <07A8AF089844A549B0F0CCC727AF31AA56F968ED@OVEKBG.ove.at> Hello, I'm not sure this is a bug - therefore I want to ask her first: I'm using Scilab 5.4.1 64 bit under Win7: Depending on the size of the figure the axis are disappearing or not. Is it a bug or can I avoid/circumvent this feature somehow? Thanks a lot Wolfgang -------------- next part -------------- A non-text attachment was scrubbed... Name: test.png Type: image/png Size: 10274 bytes Desc: test.png URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: test1.png Type: image/png Size: 10917 bytes Desc: test1.png URL: From calixte.denizet at scilab-enterprises.com Wed Feb 26 16:07:52 2014 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Wed, 26 Feb 2014 16:07:52 +0100 Subject: [Scilab-users] Bug? In-Reply-To: <07A8AF089844A549B0F0CCC727AF31AA56F968ED@OVEKBG.ove.at> References: <07A8AF089844A549B0F0CCC727AF31AA56F968ED@OVEKBG.ove.at> Message-ID: <530E0348.5050008@scilab-enterprises.com> Hi Wolfgang, It is probably a bug. Could you report it on bugzilla.scilab.org with a little test case please ? Thanks Best regards Calixte On 26/02/2014 15:56, SCHULZ Wolfgang wrote: > Hello, > I'm not sure this is a bug - therefore I want to ask her first: > I'm using Scilab 5.4.1 64 bit under Win7: > Depending on the size of the figure the axis are disappearing or not. > > Is it a bug or can I avoid/circumvent this feature somehow? > Thanks a lot > Wolfgang > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Calixte Denizet Software Development Engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From antoine.monmayrant at laas.fr Wed Feb 26 16:28:41 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Wed, 26 Feb 2014 16:28:41 +0100 Subject: [Scilab-users] Bug? In-Reply-To: <530E0348.5050008@scilab-enterprises.com> References: <07A8AF089844A549B0F0CCC727AF31AA56F968ED@OVEKBG.ove.at> <530E0348.5050008@scilab-enterprises.com> Message-ID: <530E0829.5060604@laas.fr> On 02/26/2014 04:07 PM, Calixte Denizet wrote: > Hi Wolfgang, > > It is probably a bug. > Could you report it on bugzilla.scilab.org with a little test case > please ? And if you post the test case here, I'll try to reproduce it on linux... Antoine > > Thanks > > Best regards > > Calixte > > On 26/02/2014 15:56, SCHULZ Wolfgang wrote: >> Hello, >> I'm not sure this is a bug - therefore I want to ask her first: >> I'm using Scilab 5.4.1 64 bit under Win7: >> Depending on the size of the figure the axis are disappearing or not. >> >> Is it a bug or can I avoid/circumvent this feature somehow? >> Thanks a lot >> Wolfgang >> >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > > -- > Calixte Denizet > Software Development Engineer > ----------------------------------------------------------- > Scilab Enterprises > 143bis rue Yves Le Coz - 78000 Versailles, France > http://www.scilab-enterprises.com > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From athanasiou at dica.unict.it Wed Feb 26 16:13:19 2014 From: athanasiou at dica.unict.it (Anastasia Athanasiou) Date: Wed, 26 Feb 2014 16:13:19 +0100 Subject: [Scilab-users] fmincon won't work Message-ID: Hello everyone! I am a new Scilab 5.4.1 user (in Ubuntu 12.04.4). I installed fmincon 0.3.1-1 using atoms, but when I tried it on the first help example I obtained a rather weird output: fval = 0. x = 10^(-318) * Inf Inf Could anyone help me please? Thank you in advance, Natasha From paul.bignier at scilab-enterprises.com Wed Feb 26 16:46:24 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Wed, 26 Feb 2014 16:46:24 +0100 Subject: [Scilab-users] fmincon won't work In-Reply-To: References: Message-ID: <530E0C50.4050809@scilab-enterprises.com> Hello, Did you check that your problem matched the cases described on fmincon home page ? Regards, Paul On 02/26/2014 04:13 PM, Anastasia Athanasiou wrote: > Hello everyone! I am a new Scilab 5.4.1 user (in Ubuntu 12.04.4). > I installed fmincon 0.3.1-1 using atoms, but when I tried it on the first > help example I obtained a rather weird output: > > fval = 0. > > x = 10^(-318) * Inf Inf > > Could anyone help me please? > > Thank you in advance, Natasha > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Paul BIGNIER Development engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.69 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Wed Feb 26 20:25:38 2014 From: stephane.mottelet at utc.fr (=?ISO-8859-1?Q?=22St=E9phane_Mottelet_=28UTC=29=22?=) Date: Wed, 26 Feb 2014 20:25:38 +0100 Subject: [Scilab-users] fmincon won't work In-Reply-To: References: <530E0C50.4050809@scilab-enterprises.com> Message-ID: <530E3FB2.5050301@utc.fr> Hello, Currently, sci_ipopt seems to be broken on Linux and MacOSX when using the atoms package for Scilab 5.4.1. Since fmincon relies on sci_ipot, fmicon is equally broken... BTW, users needing an optimization toolbox with general non-linear constraints are not lucky with Scilab. As far as I am concerned, I still use fsqp, which cannot be easily released as an Atoms module because of proprietary code. But Ipopt is said to really outperform fsqp, so I am rather impatient to see it coming. As it is the only module of this class/kind (optimization + complementary module), the Scilab team should take care of the functionnality of it when new versions of Scilab are released. St?phane Le 26/02/2014 20:13, St?phane Mottelet a ?crit : > > > ---------- Forwarded message ---------- > From: *Paul Bignier* > > Date: 2014-02-26 16:46 GMT+01:00 > Subject: Re: [Scilab-users] fmincon won't work > To: athanasiou at dica.unict.it , > "International users mailing list for Scilab." > > > > > Hello, > > Did you check that your problem matched the cases described on fmincon > home page ? > > Regards, > Paul > > > > On 02/26/2014 04:13 PM, Anastasia Athanasiou wrote: >> Hello everyone! I am a new Scilab 5.4.1 user (in Ubuntu 12.04.4). >> I installed fmincon 0.3.1-1 using atoms, but when I tried it on the first >> help example I obtained a rather weird output: >> >> fval = 0. >> >> x = 10^(-318) * Inf Inf >> >> Could anyone help me please? >> >> Thank you in advance, Natasha >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > -- > Paul BIGNIER > Development engineer > ----------------------------------------------------------- > Scilab Enterprises > 143bis rue Yves Le Coz - 78000 Versailles, France > Phone:+33.1.80.77.04.69 > http://www.scilab-enterprises.com > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From W.Schulz at ove.at Wed Feb 26 20:38:30 2014 From: W.Schulz at ove.at (SCHULZ Wolfgang) Date: Wed, 26 Feb 2014 19:38:30 +0000 Subject: [Scilab-users] Bug? In-Reply-To: <530E0348.5050008@scilab-enterprises.com> References: <07A8AF089844A549B0F0CCC727AF31AA56F968ED@OVEKBG.ove.at> <530E0348.5050008@scilab-enterprises.com> Message-ID: <07A8AF089844A549B0F0CCC727AF31AA56F96A8F@OVEKBG.ove.at> I reported it - bug ID 13255. Thanks Wolfgang Von: users [mailto:users-bounces at lists.scilab.org] Im Auftrag von Calixte Denizet Gesendet: Mittwoch, 26. Februar 2014 16:08 An: users at lists.scilab.org Betreff: Re: [Scilab-users] Bug? Hi Wolfgang, It is probably a bug. Could you report it on bugzilla.scilab.org with a little test case please ? Thanks Best regards Calixte On 26/02/2014 15:56, SCHULZ Wolfgang wrote: Hello, I'm not sure this is a bug - therefore I want to ask her first: I'm using Scilab 5.4.1 64 bit under Win7: Depending on the size of the figure the axis are disappearing or not. Is it a bug or can I avoid/circumvent this feature somehow? Thanks a lot Wolfgang _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -- Calixte Denizet Software Development Engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From W.Schulz at ove.at Wed Feb 26 20:39:22 2014 From: W.Schulz at ove.at (SCHULZ Wolfgang) Date: Wed, 26 Feb 2014 19:39:22 +0000 Subject: [Scilab-users] Bug? In-Reply-To: <530E0829.5060604@laas.fr> References: <07A8AF089844A549B0F0CCC727AF31AA56F968ED@OVEKBG.ove.at> <530E0348.5050008@scilab-enterprises.com> <530E0829.5060604@laas.fr> Message-ID: <07A8AF089844A549B0F0CCC727AF31AA56F96A9C@OVEKBG.ove.at> Hello, the test case is attached. Wolfgang Von: users [mailto:users-bounces at lists.scilab.org] Im Auftrag von Antoine Monmayrant Gesendet: Mittwoch, 26. Februar 2014 16:29 An: International users mailing list for Scilab. Betreff: Re: [Scilab-users] Bug? On 02/26/2014 04:07 PM, Calixte Denizet wrote: Hi Wolfgang, It is probably a bug. Could you report it on bugzilla.scilab.org with a little test case please ? And if you post the test case here, I'll try to reproduce it on linux... Antoine Thanks Best regards Calixte On 26/02/2014 15:56, SCHULZ Wolfgang wrote: Hello, I'm not sure this is a bug - therefore I want to ask her first: I'm using Scilab 5.4.1 64 bit under Win7: Depending on the size of the figure the axis are disappearing or not. Is it a bug or can I avoid/circumvent this feature somehow? Thanks a lot Wolfgang _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -- Calixte Denizet Software Development Engineer ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France http://www.scilab-enterprises.com _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Scilab_bugreport.zip Type: application/x-zip-compressed Size: 1793 bytes Desc: Scilab_bugreport.zip URL: From athanasiou at dica.unict.it Thu Feb 27 09:58:11 2014 From: athanasiou at dica.unict.it (Anastasia Athanasiou) Date: Thu, 27 Feb 2014 09:58:11 +0100 Subject: [Scilab-users] fmincon won't work In-Reply-To: <530E3FB2.5050301@utc.fr> References: <530E0C50.4050809@scilab-enterprises.com> <530E3FB2.5050301@utc.fr> Message-ID: Thank you St?phane, actually I tried fmincon in the previous Scilab version (5.3.3) and it works there! Have a nice day! Regards, Anastasia > Hello, > > Currently, sci_ipopt seems to be broken on Linux and MacOSX when using > the atoms package for Scilab 5.4.1. Since fmincon relies on sci_ipot, > fmicon is equally broken... BTW, users needing an optimization toolbox > with general non-linear constraints are not lucky with Scilab. As far as > I am concerned, I still use fsqp, which cannot be easily released as an > Atoms module because of proprietary code. But Ipopt is said to really > outperform fsqp, so I am rather impatient to see it coming. > > As it is the only module of this class/kind (optimization + > complementary module), the Scilab team should take care of the > functionnality of it when new versions of Scilab are released. > > St?phane > > > Le 26/02/2014 20:13, St?phane Mottelet a ?crit : >> >> >> ---------- Forwarded message ---------- >> From: *Paul Bignier* > > >> Date: 2014-02-26 16:46 GMT+01:00 >> Subject: Re: [Scilab-users] fmincon won't work >> To: athanasiou at dica.unict.it , >> "International users mailing list for Scilab." > > >> >> >> >> Hello, >> >> Did you check that your problem matched the cases described on fmincon >> home page ? >> >> Regards, >> Paul >> >> >> >> On 02/26/2014 04:13 PM, Anastasia Athanasiou wrote: >>> Hello everyone! I am a new Scilab 5.4.1 user (in Ubuntu 12.04.4). >>> I installed fmincon 0.3.1-1 using atoms, but when I tried it on the >>> first >>> help example I obtained a rather weird output: >>> >>> fval = 0. >>> >>> x = 10^(-318) * Inf Inf >>> >>> Could anyone help me please? >>> >>> Thank you in advance, Natasha >>> >>> _______________________________________________ >>> users mailing list >>> users at lists.scilab.org >>> http://lists.scilab.org/mailman/listinfo/users >> >> -- >> Paul BIGNIER >> Development engineer >> ----------------------------------------------------------- >> Scilab Enterprises >> 143bis rue Yves Le Coz - 78000 Versailles, France >> Phone:+33.1.80.77.04.69 >> http://www.scilab-enterprises.com >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> >> > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From clement.david at scilab-enterprises.com Thu Feb 27 15:45:59 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Thu, 27 Feb 2014 15:45:59 +0100 Subject: [Scilab-users] Trying to use old Scicos file with Xcos In-Reply-To: References: <1391382524436-4028498.post@n3.nabble.com> <1391506993.5710.5.camel@paros> Message-ID: <1393512359.2204.3.camel@paros> Hello, Sorry for the lag, please provide it as a zcos email attachment. -- Cl?ment Le mercredi 05 f?vrier 2014 ? 11:52 -0500, John David Heinzmann a ?crit : > Hi, Cl?ment, > > > I am new here, so would you tell me the best way to provide my model? > Simply as an email attachment? .zcos format to save space? > > > Thanks. > > > JD > > > On Tue, Feb 4, 2014 at 4:43 AM, Cl?ment David > wrote: > Hello, > > It should work fine but it is model related. Can you provide > your > model ? > > -- > Cl?ment > > Le dimanche 02 f?vrier 2014 ? 15:08 -0800, JD Heinzmann a > ?crit : > > About 3 years ago I developed a Scicos simulation of a heat > engine using > > ScicosLab-4.4b8 with Scicos 4.4. I have just downloaded > Scilab 5.4.1 and > > would like to run this old simulation under Xcos. > > > > I can open the old .cos file in Xcos and see the blocks and > navigate into > > the super-block hierarchy just fine. I can save this block > diagram as .xcos > > and .zcos files just fine. I can load the context file into > Scilab and run > > it to set up the workspace before running the simulation. > But when I try to > > run the simulation under Xcos, Scilab and Xcos immediately > quits silently > > (with no error messages). > > > > Are Scicos files not compatible with Xcos? If not, what do > I need to do to > > adjust my model to work with Xcos? > > > > Thank you! > > > > JD > > > > > > > > -- > > View this message in context: > http://mailinglists.scilab.org/Trying-to-use-old-Scicos-file-with-Xcos-tp4028498.html > > Sent from the Scilab users - Mailing Lists Archives mailing > list archive at Nabble.com. > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > http://lists.scilab.org/mailman/listinfo/users > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From athanasiou at dica.unict.it Thu Feb 27 19:57:54 2014 From: athanasiou at dica.unict.it (natasha) Date: Thu, 27 Feb 2014 10:57:54 -0800 (PST) Subject: [Scilab-users] CMA-ES boundary handling in scilab Message-ID: <1393527474008-4029013.post@n3.nabble.com> Hello everyone, I am currently working on translating an optimization problem (solved with CMA-ES) from Matlab to Scilab. If I have understood properly cma_optim does not contain any Matlab equivalent option of optsLB, opts.UB . In the source page there is a comment on boundary handling https://www.lri.fr/~hansen/cmaes_inmatlab.html . Has anyone ever tried establishing a search space for the CMA-ES in Scilab? Thank you very much for your attention, Anastasia -- View this message in context: http://mailinglists.scilab.org/CMA-ES-boundary-handling-in-scilab-tp4029013.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From cfuttrup at gmail.com Thu Feb 27 20:07:29 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Thu, 27 Feb 2014 20:07:29 +0100 Subject: [Scilab-users] CMA-ES boundary handling in scilab In-Reply-To: <1393527474008-4029013.post@n3.nabble.com> References: <1393527474008-4029013.post@n3.nabble.com> Message-ID: <530F8CF1.5080302@gmail.com> Hi Natasha I have several years of experience with CMA-ES. I am also the author of the Scilab ATOMS module you can find here: https://atoms.scilab.org/toolboxes/CMA-ES I am well aware the code is not working well with Scilab (i.e. if you install the ATOMS module it does something bad to your Scilab and the Scilab GUI doesn't work). I don't know what to do about it (need help from Scilab experts). You could say, Scilab needs to document their ATOMS properly, so that independent people (like me) can make fully working ATOMS modules. Regarding search space, I have successfully implemented much of this. You can write me in private ( cfuttrup -at- gmail.com ), not to unnecessarily disturb this Scilab forum. Best regards, Claus On 2/27/2014 19:57, natasha wrote: > Hello everyone, > > I am currently working on translating an optimization problem (solved with > CMA-ES) from Matlab to Scilab. > > If I have understood properly cma_optim does not contain any Matlab > equivalent option of optsLB, opts.UB . > > In the source page there is a comment on boundary handling > https://www.lri.fr/~hansen/cmaes_inmatlab.html > . > > Has anyone ever tried establishing a search space for the CMA-ES in Scilab? > > Thank you very much for your attention, Anastasia > > > > > > -- > View this message in context: http://mailinglists.scilab.org/CMA-ES-boundary-handling-in-scilab-tp4029013.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > --- This email is free from viruses and malware because avast! Antivirus protection is active. http://www.avast.com From rasleendeol at yahoo.com Fri Feb 28 08:07:28 2014 From: rasleendeol at yahoo.com (rasleendeol) Date: Thu, 27 Feb 2014 23:07:28 -0800 (PST) Subject: [Scilab-users] operations on a DAG Message-ID: <1393571248887-4029021.post@n3.nabble.com> i have implemented a directed acyclic graph in Matlab n would like to know how to perform operations such as to find critical path, bottom-level, top-level...etc. also, clustering needs to be done based upon the critical paths. -- View this message in context: http://mailinglists.scilab.org/operations-on-a-DAG-tp4029021.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Fri Feb 28 11:18:37 2014 From: sgougeon at free.fr (sgougeon at free.fr) Date: Fri, 28 Feb 2014 11:18:37 +0100 (CET) Subject: [Scilab-users] operations on a DAG In-Reply-To: <1393571248887-4029021.post@n3.nabble.com> Message-ID: <1220328945.1039149172.1393582717804.JavaMail.root@zimbra75-e12.priv.proxad.net> Hello, Metanet might be your friend: http://atoms.scilab.org/toolboxes/metanet It should be available through the ATOMS module manager (open box icon in the console toolbar). Regards Samuel ----- Mail original ----- De: "rasleendeol" ?: users at lists.scilab.org Envoy?: Vendredi 28 F?vrier 2014 08:07:28 Objet: [Scilab-users] operations on a DAG i have implemented a directed acyclic graph in Matlab n would like to know how to perform operations such as to find critical path, bottom-level, top-level...etc. also, clustering needs to be done based upon the critical paths.