From rfabbri at gmail.com Sun Jun 1 00:28:21 2014 From: rfabbri at gmail.com (Ricardo Fabbri) Date: Sat, 31 May 2014 19:28:21 -0300 Subject: [Scilab-users] Image processing for color positioning and distance from a specific point. In-Reply-To: References: Message-ID: Hi, take a look at siptoolbox.sf.net. Download the Git version, and there is a function called color_classify in macros/color_classify.sci It takes an input image and classifies each pixel into 6 broad classes: Red, Green, Blue, Yellow, Black, or 'others'. I haven't documented it properly yet, but you can play with it. The Red color was tuned to be sensitive, but you can change the parameters inside the function to suit your detection needs. SIP is designed for Linux and Mac OSX but you should be able to use this specific function in another OS. Best, -- Dr Ricardo Fabbri Professor of Computer Engineering GNU/Linux registered user #175401 www.lems.brown.edu/~rfabbri pt.wikipedia.org/wiki/IPRJ labmacambira.sf.net On Wed, May 28, 2014 at 5:35 AM, Varun Singh wrote: > Hello > > I am having problems with a program regarding the color detection and > location of that particular color in the picture. I browsed through many > tutorials but was unable to put it into effect for the desired output of > mine. Basically i want to locate the red points in the attached image and > also find the distance between the red and green parts of the image using a > scilab program. > > So can u please send me the commands and instructions regarding how i can > proceed through this as i m new to scilab, asap and oblige. > > Regards > Varun Singh > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Sun Jun 1 00:32:49 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 01 Jun 2014 00:32:49 +0200 Subject: [Scilab-users] How to get SIVP working with Scilab on the Mac? In-Reply-To: <1401538032926-4030634.post@n3.nabble.com> References: <1401538032926-4030634.post@n3.nabble.com> Message-ID: <538A5891.8020702@free.fr> Hello, Le 31/05/2014 14:07, mstoth a ?crit : > I may have put this on the wrong list originally, perhaps someone here can > help? > > atoms has SIVP on Windows but not on Mac for some reason. AFAIK, SIVP uses TcL/Tk windows rather than Java ones. Now, TcL/Tk is not supported on Mac (for some reason that could be found in the archives). I am afraid that there is really no chance to run SIVP on Mac. Regards Samuel From antoine.elias at scilab-enterprises.com Sun Jun 1 11:41:24 2014 From: antoine.elias at scilab-enterprises.com (Antoine ELIAS) Date: Sun, 01 Jun 2014 11:41:24 +0200 Subject: [Scilab-users] Help write simple GUI In-Reply-To: <538A34FD.1050505@gmail.com> References: <538A34FD.1050505@gmail.com> Message-ID: <538AF544.9080806@scilab-enterprises.com> Hi Claus, For small dialogs to ask inputs, use *x_mdialog* function ( http://help.scilab.org/docs/5.5.0/fr_FR/x_mdialog.html ) function result = compute_voltage() txt = ['Power (Watt)';'Impedance (ohm)']; values = *x_mdialog*('System values',txt,['100';'8']) if values <> [] then w = values(1); o = values(2); w = eval(w); o = eval(o); if typeof(w) <> "constant" | int(w) <> w | typeof(o) <> "constant" | int(o) <> o then error("bad value, integer values expected"); end result = sqrt(w * o); else result = []; end endfunction voltage = compute_voltage() Antoine Le 31/05/2014 22:01, Claus Futtrup a ?crit : > Hi there > > I hope someone can help me make a simple GUI front end to my script, preferably without using GUI Builder or guimaker, because I think I'll understand better without using additional tools. > > The GUI - a window containing two variable input: > 1. text + input box + text > Again the input is a number, typically an integer (as in "system power : 100 Watt"). 100 being the user variable input. > 2. text + input box + text > A number, typically an intereger (as in "impedance : 8 ohm"). 8 being the user variable input. > > The script will then based on 100 Watt into 8 ohm calculate the voltage = sqrt(100*8). I will round the number to one decimal, convert it to a string and create an answer, like "Input level : 28.3 > Volt". > > The GUI should contain two buttons, e.g. an OK + a QUIT button. The OK button runs the calculations. The quit button quits the script (may as well quit Scilab). > > I've studied the GUI controls in Scilab already and I understand the the old way of making a new window is: > > f = figure(); > > but it seems there's a new command in Scilab 5.5.0: > > f = createWindow(); > > The new method probably has its advantages, but the Scilab manual doesn't really explain... can anybody explain to me? > > Can you use the same uicontrol commands with both methods? (can I use examples found on the internet using "figure" and just replace figure with createWindow?). It seems to me that the examples in > the createWindow gives uicontrol some inputs named "units" - "normalized" - etc. - which are not present when looking at old figure() examples. > > My next question will be, how do I create a uicontrol which accepts a string input (as I need 2 of those). The uicontrol documentation says the "Edit" style will do that. Funny, but also a bit > strange, that you can find all sorts of examples for using uicontrol with checboxes, listboxes, radiobuttons, pushbuttons, tables etc. - but no examples are provided with this most useful editable > string input ... (???) > > ... hereafter I expect to evaluate the two strings (or if they're returned as numbers, that's fine of course, I just don't need the script to crash if something is not right). Then execute my code. > > Best regards, > Claus > > --- > This email is free from viruses and malware because avast! Antivirus protection is active. > http://www.avast.com > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Antoine ELIAS Software developer ----------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles Phone: 01.80.77.04.70 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Sun Jun 1 12:31:22 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sun, 01 Jun 2014 12:31:22 +0200 Subject: [Scilab-users] Help write simple GUI In-Reply-To: <538AF544.9080806@scilab-enterprises.com> References: <538A34FD.1050505@gmail.com> <538AF544.9080806@scilab-enterprises.com> Message-ID: <538B00FA.8020609@gmail.com> Hi Antoine Thank you for your suggestion, which almost does what I need. The users are not necessarily Scilab users, so I'd like a user interface where the result is also displayed. I realize I could show a messagebox with the result. What I'd like to do requires a single window/figure where input and output is shown + a QUIT button so that the user always only sees one single frame / user interface, as an application / for simplicity (no Scilab). Later I will add a noise spectrum simulation + a filter - e.g. 2. order high-pass butterworth at 2500 Hz. The filtered response will reduce the voltage, but how much? The user should be informed about before + after filter (3 inputs, 2 outputs) and dialog boxes is not as nice when it's more complicated like this. Best regards, Claus On 01-Jun-14 11:41, Antoine ELIAS wrote: > Hi Claus, > > For small dialogs to ask inputs, use *x_mdialog* function ( > http://help.scilab.org/docs/5.5.0/fr_FR/x_mdialog.html ) > > function result = compute_voltage() > txt = ['Power (Watt)';'Impedance (ohm)']; > values = *x_mdialog*('System values',txt,['100';'8']) > if values <> [] then > w = values(1); > o = values(2); > w = eval(w); > o = eval(o); > > if typeof(w) <> "constant" | int(w) <> w | typeof(o) <> > "constant" | int(o) <> o then > error("bad value, integer values expected"); > end > > result = sqrt(w * o); > else > result = []; > end > endfunction > > voltage = compute_voltage() > > Antoine > > > Le 31/05/2014 22:01, Claus Futtrup a ?crit : >> Hi there >> >> I hope someone can help me make a simple GUI front end to my script, >> preferably without using GUI Builder or guimaker, because I think >> I'll understand better without using additional tools. >> >> The GUI - a window containing two variable input: >> 1. text + input box + text >> Again the input is a number, typically an integer (as in "system >> power : 100 Watt"). 100 being the user variable input. >> 2. text + input box + text >> A number, typically an intereger (as in "impedance : 8 ohm"). 8 >> being the user variable input. >> >> The script will then based on 100 Watt into 8 ohm calculate the >> voltage = sqrt(100*8). I will round the number to one decimal, >> convert it to a string and create an answer, like "Input level : 28.3 >> Volt". >> >> The GUI should contain two buttons, e.g. an OK + a QUIT button. The >> OK button runs the calculations. The quit button quits the script >> (may as well quit Scilab). >> >> I've studied the GUI controls in Scilab already and I understand the >> the old way of making a new window is: >> >> f = figure(); >> >> but it seems there's a new command in Scilab 5.5.0: >> >> f = createWindow(); >> >> The new method probably has its advantages, but the Scilab manual >> doesn't really explain... can anybody explain to me? >> >> Can you use the same uicontrol commands with both methods? (can I use >> examples found on the internet using "figure" and just replace figure >> with createWindow?). It seems to me that the examples in the >> createWindow gives uicontrol some inputs named "units" - "normalized" >> - etc. - which are not present when looking at old figure() examples. >> >> My next question will be, how do I create a uicontrol which accepts a >> string input (as I need 2 of those). The uicontrol documentation says >> the "Edit" style will do that. Funny, but also a bit strange, that >> you can find all sorts of examples for using uicontrol with >> checboxes, listboxes, radiobuttons, pushbuttons, tables etc. - but no >> examples are provided with this most useful editable string input ... >> (???) >> >> ... hereafter I expect to evaluate the two strings (or if they're >> returned as numbers, that's fine of course, I just don't need the >> script to crash if something is not right). Then execute my code. >> >> Best regards, >> Claus >> >> --- >> This email is free from viruses and malware because avast! Antivirus >> protection is active. >> http://www.avast.com >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > > -- > Antoine ELIAS > Software developer > ----------------------- > Scilab Enterprises > 143bis rue Yves Le Coz - 78000 Versailles > Phone: 01.80.77.04.70 > http://www.scilab-enterprises.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 -------------- next part -------------- An HTML attachment was scrubbed... URL: From sterg at teikav.edu.gr Sun Jun 1 15:29:03 2014 From: sterg at teikav.edu.gr (Stergios Papadimitriou) Date: Sun, 01 Jun 2014 16:29:03 +0300 Subject: [Scilab-users] ScalaLab and GroovyLab support for caling SciLab Message-ID: <538B2A9F.1030108@teikav.edu.gr> Hi all, I like to inform the group (mostly concerns the Java funs of it) that I interfaced my projects ScalaLab and GroovyLab to cooperate with SciLab scripts: https://code.google.com/p/scalalab/wiki/SciLabInScalaLab https://code.google.com/p/jlabgroovy/wiki/SciLabInGroovyLab Best Regards Stergios From michael at virtualpianist.com Mon Jun 2 10:51:13 2014 From: michael at virtualpianist.com (mstoth) Date: Mon, 2 Jun 2014 01:51:13 -0700 (PDT) Subject: [Scilab-users] How to get SIVP working with Scilab on the Mac? In-Reply-To: <538A5891.8020702@free.fr> References: <1401538032926-4030634.post@n3.nabble.com> <538A5891.8020702@free.fr> Message-ID: <1401699073288-4030646.post@n3.nabble.com> Thank you for your response. I guess I have to wonder what people use to do image processing on the mac with Scilab. -- View this message in context: http://mailinglists.scilab.org/How-to-get-SIVP-working-with-Scilab-on-the-Mac-tp4030634p4030646.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Mon Jun 2 13:24:25 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 02 Jun 2014 13:24:25 +0200 Subject: [Scilab-users] How to get SIVP working with Scilab on the Mac? In-Reply-To: <1401699073288-4030646.post@n3.nabble.com> References: <1401538032926-4030634.post@n3.nabble.com> <538A5891.8020702@free.fr> <1401699073288-4030646.post@n3.nabble.com> Message-ID: <538C5EE9.4020504@free.fr> Le 02/06/2014 10:51, mstoth a ?crit : > Thank you for your response. > > I guess I have to wonder what people use to do image processing on the mac with Scilab. Right... In a recent message, Riccardo Fabri told that SIP can be used with OSX: http://mailinglists.scilab.org/Scilab-users-Image-processing-for-color-positioning-and-distance-from-a-specific-point-tp4030626.html Since image display is now possible and easy through Matlplot() with Scilab 5.5.0, this might invite authors of image processing packages to upgrade their display function, and give up with the TCl/Tk technology for doing this. Regards -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Mon Jun 2 13:41:40 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 02 Jun 2014 13:41:40 +0200 Subject: [Scilab-users] Figure In-Reply-To: <538639E7.7080102@durietz.se> References: <53838A7A.2080405@durietz.se> <53839BB2.3020005@free.fr> <538639E7.7080102@durietz.se> Message-ID: <538C62F4.8020109@free.fr> Le 28/05/2014 21:32, Stefan Du Rietz a ?crit : > On 2014-05-26 21:53, Samuel Gougeon wrote: >> Le 26/05/2014 20:39, Stefan Du Rietz a ?crit : >>> Hello, >>> I am trying to get my figures from 5.4.1 in 5.5.0. >>> >>> Why is everything red? (no help found and no answer to an equal >>> question on the Internet) >> Like here?: >> http://bugzilla.scilab.org/attachment.cgi?id=3561 >> http://bugzilla.scilab.org/show_bug.cgi?id=12788 >> > Thank you, Samuel. I have made a comment to your bug report. > > Another question: > Is it possible to get rid of the extra title bar in a figure? It just > repeats the text in the main title bar and makes the axes space lower, > to no avail whatsoever. It is used as handle for the docking system. When a window is set as .dockable="off", this bar could be masked (since then it becomes useless). Scilab 5.5.0 allows to create a window without this bar, through the new function createWindow(). HTH Samuel From flavio.alves at gmail.com Mon Jun 2 14:38:53 2014 From: flavio.alves at gmail.com (Flavio Castro Alves Filho) Date: Mon, 2 Jun 2014 09:38:53 -0300 Subject: [Scilab-users] Running Xcos on Beaglebone Black Message-ID: Hello, I'm trying to execute Xcos on a Beaglebone Black (ARM based embedded computer) and it is not working. I posted an issue on the Bugzilla website: http://bugzilla.scilab.org/show_bug.cgi?id=13451 Does anybody faced similar problem? How should I solve this? Best regards, Flavio -- Flavio de Castro Alves Filho flavio.alves at gmail.com www.linuxembarcado.com Twitter: http://twitter.com/#!/fraviofii LinkedIn profile: www.linkedin.com/in/flaviocastroalves From sdr at durietz.se Mon Jun 2 19:02:50 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Mon, 02 Jun 2014 19:02:50 +0200 Subject: [Scilab-users] Figure In-Reply-To: <538C62F4.8020109@free.fr> References: <53838A7A.2080405@durietz.se> <53839BB2.3020005@free.fr> <538639E7.7080102@durietz.se> <538C62F4.8020109@free.fr> Message-ID: <538CAE3A.5090904@durietz.se> On 2014-06-02 13:41, Samuel Gougeon wrote: -------------------- > Le 28/05/2014 21:32, Stefan Du Rietz a ?crit : >> On 2014-05-26 21:53, Samuel Gougeon wrote: >>> Le 26/05/2014 20:39, Stefan Du Rietz a ?crit : >>>> Hello, >>>> I am trying to get my figures from 5.4.1 in 5.5.0. >>>> >>>> Why is everything red? (no help found and no answer to an equal >>>> question on the Internet) >>> Like here?: >>> http://bugzilla.scilab.org/attachment.cgi?id=3561 >>> http://bugzilla.scilab.org/show_bug.cgi?id=12788 >>> >> Thank you, Samuel. I have made a comment to your bug report. >> >> Another question: >> Is it possible to get rid of the extra title bar in a figure? It >> just repeats the text in the main title bar and makes the axes space >> lower, to no avail whatsoever. > It is used as handle for the docking system. When a window is set as > .dockable="off", this bar could be masked (since then it becomes > useless). > Scilab 5.5.0 allows to create a window without this bar, through the > new function createWindow(). > Excellent! Thank you very much! Is there any help on docking? Stefan From cfuttrup at gmail.com Mon Jun 2 19:16:28 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Mon, 02 Jun 2014 19:16:28 +0200 Subject: [Scilab-users] createWindow and uicontrol Message-ID: <538CB16C.9020305@gmail.com> Hi there My prevoius email didn't get a really good answer, so I try again - to get partially a good answer: f = createWindow(); The new method probably has its advantages, but the Scilab manual doesn't really explain... can anybody explain to me? Can you use the same uicontrol commands with both methods? (can I use examples found on the internet using "figure" and just replace figure with createWindow?). It seems to me that the examples in the createWindow gives uicontrol some inputs named "units" - "normalized" - etc. - which are not present when looking at old figure() examples. My next question will be, how do I create a uicontrol which accepts a string input (as I need 2 of those). The uicontrol documentation says the "Edit" style will do that, but how? /Claus --- This email is free from viruses and malware because avast! Antivirus protection is active. http://www.avast.com From sdr at durietz.se Mon Jun 2 20:05:43 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Mon, 02 Jun 2014 20:05:43 +0200 Subject: [Scilab-users] Figure In-Reply-To: <538CAE3A.5090904@durietz.se> References: <53838A7A.2080405@durietz.se> <53839BB2.3020005@free.fr> <538639E7.7080102@durietz.se> <538C62F4.8020109@free.fr> <538CAE3A.5090904@durietz.se> Message-ID: <538CBCF7.4010601@durietz.se> On 2014-06-02 19:02, Stefan Du Rietz wrote: -------------------- > On 2014-06-02 13:41, Samuel Gougeon wrote: > -------------------- >> Le 28/05/2014 21:32, Stefan Du Rietz a ?crit : >>> On 2014-05-26 21:53, Samuel Gougeon wrote: >>>> Le 26/05/2014 20:39, Stefan Du Rietz a ?crit : >>>>> Hello, >>>>> I am trying to get my figures from 5.4.1 in 5.5.0. >>>>> >>>>> Why is everything red? (no help found and no answer to an equal >>>>> question on the Internet) >>>> Like here?: >>>> http://bugzilla.scilab.org/attachment.cgi?id=3561 >>>> http://bugzilla.scilab.org/show_bug.cgi?id=12788 >>>> >>> Thank you, Samuel. I have made a comment to your bug report. >>> >>> Another question: >>> Is it possible to get rid of the extra title bar in a figure? It >>> just repeats the text in the main title bar and makes the axes space >>> lower, to no avail whatsoever. >> It is used as handle for the docking system. When a window is set as >> .dockable="off", this bar could be masked (since then it becomes >> useless). >> Scilab 5.5.0 allows to create a window without this bar, through the >> new function createWindow(). >> > Excellent! Thank you very much! > > Is there any help on docking? > And how can I plot anything in this figure window? gcf() returns a handle to this figure, but newaxes() and plot() creates a new window ... However, I can put GUI:s in it. Stefan From antoine.monmayrant at laas.fr Tue Jun 3 08:53:10 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Tue, 03 Jun 2014 08:53:10 +0200 Subject: [Scilab-users] createWindow and uicontrol In-Reply-To: <538CB16C.9020305@gmail.com> References: <538CB16C.9020305@gmail.com> Message-ID: <538D70D6.1020105@laas.fr> On 06/02/2014 07:16 PM, Claus Futtrup wrote: > Hi there > > My prevoius email didn't get a really good answer, so I try again - to > get partially a good answer: > > f = createWindow(); > > The new method probably has its advantages, but the Scilab manual > doesn't really explain... can anybody explain to me? > > Can you use the same uicontrol commands with both methods? Yes, both methods are almost identical. You can see that by looking at the source code of createWindow. Just type 'edit createWindow': function f = createWindow() f = figure( ... "dockable", "off", ... "infobar_visible", "off", ... "toolbar_visible", "off", ... "menubar_visible", "off", ... "default_axes", "off"); endfunction So createWindow is just a wrapper around figure, with different default values. > (can I use examples found on the internet using "figure" and just > replace figure with createWindow?). > > It seems to me that the examples in the createWindow gives uicontrol > some inputs named "units" - "normalized" - etc. - which are not > present when looking at old figure() examples. > > My next question will be, how do I create a uicontrol which accepts a > string input (as I need 2 of those). The uicontrol documentation says > the "Edit" style will do that, but how? Have a look at the Demos->Gui->Uicontrols1 or Uicontrols2 : there are text input uicontrols. Click on the 'view code' menu to see what's under the hood. Hope it helps, Antoine > > /Claus > > --- > This email is free from viruses and malware because avast! Antivirus > protection is active. > http://www.avast.com > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From TANTHIAMHUAT at e.ntu.edu.sg Tue Jun 3 07:33:24 2014 From: TANTHIAMHUAT at e.ntu.edu.sg (#TAN THIAM HUAT#) Date: Tue, 3 Jun 2014 05:33:24 +0000 Subject: [Scilab-users] [Scilab] Fwd: Scilab - Symbolic toolbox integration In-Reply-To: References: <1401686695094.87998@e.ntu.edu.sg> <1401767806620.24733@e.ntu.edu.sg>, Message-ID: <1401773604683.7956@e.ntu.edu.sg> ?Hi All, I would appreciate if anyone could guide me in installing Symbolic Toolbox in Scilab 5.5.0 under Windows 7. thanks. Sincerely, Thiam Huat ________________________________ From: Jayaram R Pai Sent: Tuesday, June 3, 2014 1:04 PM To: #TAN THIAM HUAT# Cc: Scilab.in Team; scilab group Subject: Re: [Scilab] Fwd: Scilab - Symbolic toolbox integration Dear Tan, I am not an Scilab Expert, I don't use Windows for programming. So it would be better if you send your query to the Official Scilab mailing lists: https://www.scilab.org/development/ml The above steps that I gave you worked for me in GNU/Linux. Please refer this Scimax README too: http://forge.scilab.org/index.php/p/scimax/source/tree/master/README Thanks, Jayaram R Pai On Tue, Jun 3, 2014 at 9:26 AM, #TAN THIAM HUAT# > wrote: Hi Jayaram, thanks for your kind reply. I am using Scilab 5.5.0 on Windows 7. I could not find the Overload Toolbox in ATOMS. Is it possible to have the Symbolic Toolbox installed? The link you provided seems to say it is only able to install in Scilab V 5.3.3, and in Linux system. I would appreciate your prompt response. Sincerely, Thiam Huat? ________________________________ From: Jayaram R Pai > Sent: Monday, June 2, 2014 5:49 PM To: #TAN THIAM HUAT# Cc: scilab group; Scilab.in Team Subject: Re: [Scilab] Fwd: Scilab - Symbolic toolbox integration Dear TAN THIAM HUAT, Please have a look at this wiki page: http://es.wikibooks.org/wiki/Manual_de_Scilab/Xcos/C%C3%A1lculo_simb%C3%B3lico_con_SciMax The above page is in spanish. Please translate it using google translate. It it one of the best material for Scimax integration. Note: Scimax can currently be build only using Scilab v5.3.3. But once its build, using a workaround you can easily make it work with Scilab v5.4.1. In case you get an error in step: Loading toolbox of the above link. Comment the below line of the file: macros/maxfun/personnalizedOverload.sce inside the scimax folder. Here is the diff. ------------------------ < overload('printf',["c"]) --- > //overload('printf',["c"]) For further knowledge, please check this repository: https://github.com/sengupta/scilab-maxima Please let me know if you are facing anymore issues. Thanks, Jayaram R Pai On Mon, Jun 2, 2014 at 11:00 AM, Scilab.in Team > wrote: Can someone address the query mentioned below? Shamika ---------- Forwarded message ---------- From: #TAN THIAM HUAT# > Date: Mon, Jun 2, 2014 at 10:54 AM Subject: Scilab - Symbolic toolbox integration To: "contact at scilab.in" > Hi, I saw something related to Symbolic toolbox from http://www.scilab.in/maxima http://www.scilab.in/node/135 http://www.scilab.in/node/706 http://www.scilab.in/node/76 Can I know how can I integrate Symbolic Toolbox into Scilab such that it is possible to have a screen like in http://forge.scilab.org/index.php/p/scimax/ Sincerely, Thiam Huat _______________________________________________ Scilab mailing list Scilab at fossee.in http://fossee.in/cgi-bin/mailman/listinfo/scilab -------------- next part -------------- An HTML attachment was scrubbed... URL: From clement.david at scilab-enterprises.com Tue Jun 3 11:06:43 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Tue, 03 Jun 2014 11:06:43 +0200 Subject: [Scilab-users] Running Xcos on Beaglebone Black In-Reply-To: References: Message-ID: <1401786403.2431.6.camel@paros> Hello Flavio, I don't know the current status of the Debian ARM JVM but last time I checked to used Scilab on Fedora ARM the JVM crashed on some JNI functions. Are you able to call others Java implemented functions ? Does JIMS work ? -- Cl?ment Le lundi 02 juin 2014 ? 09:38 -0300, Flavio Castro Alves Filho a ?crit : > Hello, > > I'm trying to execute Xcos on a Beaglebone Black (ARM based embedded > computer) and it is not working. > > I posted an issue on the Bugzilla website: > > http://bugzilla.scilab.org/show_bug.cgi?id=13451 > > Does anybody faced similar problem? How should I solve this? > > Best regards, > > Flavio > From clement.david at scilab-enterprises.com Tue Jun 3 11:34:43 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Tue, 03 Jun 2014 11:34:43 +0200 Subject: [Scilab-users] ScalaLab and GroovyLab support for caling SciLab In-Reply-To: <538B2A9F.1030108@teikav.edu.gr> References: <538B2A9F.1030108@teikav.edu.gr> Message-ID: <1401788083.2431.13.camel@paros> Hi Stergios, Great job ! I did not known scalalab nor jlabgroovy before. Do you know that Scilab can call any Java library using JIMS ? Do you think that it may be easily possible to use groovy are scala instead ? -- Cl?ment Le dimanche 01 juin 2014 ? 16:29 +0300, Stergios Papadimitriou a ?crit : > Hi all, > > I like to inform the group (mostly concerns the Java funs of it) > that I interfaced my projects ScalaLab and GroovyLab to > cooperate with SciLab scripts: > > https://code.google.com/p/scalalab/wiki/SciLabInScalaLab > > https://code.google.com/p/jlabgroovy/wiki/SciLabInGroovyLab > > > Best Regards > Stergios > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From sterg at teikav.edu.gr Tue Jun 3 11:46:46 2014 From: sterg at teikav.edu.gr (Stergios Papadimitriou) Date: Tue, 03 Jun 2014 12:46:46 +0300 Subject: [Scilab-users] ScalaLab and GroovyLab support for caling SciLab In-Reply-To: <1401788083.2431.13.camel@paros> References: <538B2A9F.1030108@teikav.edu.gr> <1401788083.2431.13.camel@paros> Message-ID: <538D9986.1080909@teikav.edu.gr> Hi Cl?ment, I need to study the Scilab-Java interface in order to improve also the interface of Scilab - ScalaLab/GroovyLab. (... this week I'm very pressed, perhaps I can do more next week!) Do you know (or somebody else) where to find good staff on SciLab-Java interfacing? I do not know JIMS, a possible link for it? .. perhaps, it is good to have both the compactness of SciLab and the general purpose scripting of Scala/Groovy, this is why I like to improve on that interfacing work ... Best Regards Stergios ???? 3/6/2014 12:34 ??, ?/? Cl?ment David ??????: > Hi Stergios, > > Great job ! I did not known scalalab nor jlabgroovy before. > > Do you know that Scilab can call any Java library using JIMS ? Do you > think that it may be easily possible to use groovy are scala instead ? > > -- > Cl?ment > > Le dimanche 01 juin 2014 ? 16:29 +0300, Stergios Papadimitriou a ?crit : >> Hi all, >> >> I like to inform the group (mostly concerns the Java funs of it) >> that I interfaced my projects ScalaLab and GroovyLab to >> cooperate with SciLab scripts: >> >> https://code.google.com/p/scalalab/wiki/SciLabInScalaLab >> >> https://code.google.com/p/jlabgroovy/wiki/SciLabInGroovyLab >> >> >> Best Regards >> Stergios >> _______________________________________________ >> 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 Tue Jun 3 12:07:23 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Tue, 03 Jun 2014 12:07:23 +0200 Subject: [Scilab-users] ScalaLab and GroovyLab support for caling SciLab In-Reply-To: <538D9986.1080909@teikav.edu.gr> References: <538B2A9F.1030108@teikav.edu.gr> <1401788083.2431.13.camel@paros> <538D9986.1080909@teikav.edu.gr> Message-ID: <1401790043.2431.23.camel@paros> JIMS is part of Scilab 5.5.0 and documented in the Scilab help. There is also a getting started page : http://help.scilab.org/docs/5.5.0/en_US/jims-getting-started.html -- Cl?ment Le mardi 03 juin 2014 ? 12:46 +0300, Stergios Papadimitriou a ?crit : > Hi Cl?ment, > > I need to study the Scilab-Java interface in order to improve also the > interface of Scilab - ScalaLab/GroovyLab. > (... this week I'm very pressed, perhaps I can do more next week!) > > Do you know (or somebody else) where to find good staff on SciLab-Java > interfacing? > > I do not know JIMS, a possible link for it? > > > .. perhaps, it is good to have both the compactness of SciLab and the > general purpose scripting of Scala/Groovy, > this is why I like to improve on that interfacing work ... > > Best Regards > Stergios > > > > ???? 3/6/2014 12:34 ??, ?/? Cl?ment David ??????: > > Hi Stergios, > > > > Great job ! I did not known scalalab nor jlabgroovy before. > > > > Do you know that Scilab can call any Java library using JIMS ? Do you > > think that it may be easily possible to use groovy are scala instead ? > > > > -- > > Cl?ment > > > > Le dimanche 01 juin 2014 ? 16:29 +0300, Stergios Papadimitriou a ?crit : > >> Hi all, > >> > >> I like to inform the group (mostly concerns the Java funs of it) > >> that I interfaced my projects ScalaLab and GroovyLab to > >> cooperate with SciLab scripts: > >> > >> https://code.google.com/p/scalalab/wiki/SciLabInScalaLab > >> > >> https://code.google.com/p/jlabgroovy/wiki/SciLabInGroovyLab > >> > >> > >> Best Regards > >> Stergios > >> _______________________________________________ > >> 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 bmbouter at gmail.com Tue Jun 3 12:21:08 2014 From: bmbouter at gmail.com (Brian Bouterse) Date: Tue, 3 Jun 2014 06:21:08 -0400 Subject: [Scilab-users] GROCER VAR-HMM Example Output Different from Chapter 23 In-Reply-To: References: Message-ID: Hi Eric, Thanks for the quick and helpful reply! I adjusted the first parameter from 'cte' to 'all' and that produces the correct form of output, and the values look close to the expected values. This is great! I've read the help documents several times, but I didn't realize the impact of adjusting the value of the first parameter. The doc example on the ms_var help page describes the parameter, but the only example on that page uses 'cte'. There is no example with 'all' on the ms_var() help page. Also including some sample output may also be good there. I would like to contribute an example like this to the docs, manual, or help files. I think that would have caused me to solve my own problems instead of writing long e-mails. I also tried using a larger data file 10,000 observations in a column format, and then impexc2bd transformed the data file correctly. That is also solved thanks to your suggestion. I know the code is available for download, but would you ever consider putting a copy of the GROCER code in a repository on github.com? If it were in a place like that, I, and others, would have an easier time of contributing to the codebase instead of e-mailing patches. What do you think about an idea like this? Just a friendly suggestion based on what has worked well for me in the past. I can help put it up there, and configure it if that is helpful. Thanks again for making this great software, and helping me use it! Best, Brian On Thu, May 29, 2014 at 1:41 PM, Eric Dubois wrote: > Hello Brian > > You did not recover the intial parameters because you did not estimate the > same model as the one you generated. > > This is because you used in ms_var the option 'cte' which means that only > the constant is allowed to switch: and this is why the AR coeffcients are > only given for 'all regimes'. I f you rather use the option 'all', then all > parameters are allowed to switch and the result is now consistent with the > parameters of the generated model, that is the true values are within the > confidence interval of the estimated ones. This is explained in the help > command (see help ms_var): may I ask you why you did not infer it from the > help files? > > As for the big file, I suspect it is truncated into Scilab. Can you try > with data in column instead of in rows? > > ?ric. > > > 2014-05-29 13:24 GMT+02:00 Brian Bouterse : > > Hi Eric, >> >> Thanks again for the response. What you told me to do worked. If I closed >> Scilab, and rerun the example from my earlier e-mail, I do not see any >> warnings, which is good. The output I received is not what I expected in >> several ways, so I'll ask some more questions here. There are a lot of >> different questions here, but they all relate to the example below, and the >> use of ms_var and GROCER. Thanks in advance for any insight into the issues >> I am experiencing. >> >> The example: >> >> I have a file arhmm_example_280.dat [0] that was generated and loaded >> from arhmm_example_280.csv [1] using the following commands: >> >> impexc2bd('arhmm_example_280.csv', ';', 'arhmm_example_280.dat') >> load('arhmm_example_280.dat') >> >> I then run, what I expect to be a single variable, 2-regime, >> >> nb_states=2 >> switch_var=2 >> var_opt=3 >> >> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >> >> I then receive this output: http://ur1.ca/hej98 >> >> 1) In the "coefficients" section of the output, the output doesn't seem >> to contain enough numbers to describe the expected number of outputs. The >> model has 2 regimes, and it is autoregressive with order 3, so I expect >> there to be 3 coefficients for EACH regime (regime 1 and regime 2). >> Instead, I see a single coefficient for Regime 1 (0.9313736) and a single >> coefficient for Regime 2 (8.2685973). Where are the three autoregressive >> parameters for regime 1 and the other three autoregressive parameters for >> regime 2? >> >> 2) There is a section in the output called "All Regimes" which is in the >> "Coefficients" section contains 3 coefficients. Do these correspond to the >> autoregressive aspect of the model [0.2258037, 0.4757391, -0.3161730). I >> expected to see a section like this for Regime 1, and one for Regime 2, but >> I did not. I'm not sure how these 3 coefficients for "All Regimes" fit in. >> The Markov Switching Model switches the autoregressive coefficients with >> each Regime, so I'm not sure how to interpret coefficients that are for >> "both regimes". >> >> 3) I believe each regime will have its error term estimated, and I do see >> that in the output. For example under Regime 2, the variance of the error >> term is 2.3164941. This is listed in the section "Variance-covariance >> matrix of residuals" which is a subsection under "Regime 2". Am I >> interpreting the output correctly? Where is the mean value for the regime? >> >> 4) My overall goal is to validate that I can hand a list of observations >> to GROCER+scilab, along with the order of the autoregressive model, and the >> number of regimes, and have it estimate the Markov hidden state transition >> probabilities, the estimated regime that is in effect at time t, and >> coefficients for each regime, and the error term that can be used for >> prediction of the t+1 observation. I don't know of a "well known" data set >> that allows for this type of validation, so I wrote a simple AR-HMM >> generator named koa marhn [2], which outputs AR-HMM data that follows a set >> of parameters. This is the program that generated the 280 point data set in >> the arhmm_example_280.csv [1] file. You can see the parameters used to >> generate this file here: >> >> https://github.com/bmbouter/koa_marhn/blob/master/arhmm_example.py >> >> A[i][j] is the probability transition from regime i to regime j. >> C[x][y] is the autoregressive coefficient for regime x used for >> generation of the t+1 element from the t, t-1, t-2, ..., t-y summation. >> Each y is a coefficient for regime x. >> R[q] contains mean and std_deviation values for regime q >> pi[m] contains the starting state probabilities that the regime is in >> state m at time t=0 >> >> The output from http://ur1.ca/hej98 does not correspond at all with the >> values of these inputs, which means something unexpected is going on! Do >> you have insight into why scilab doesn't produce these expected numbers? >> >> 5) One possibility is that maybe it is "working", but there isn't enough >> data for it to estimate the parameters correctly. I had koa marhn [2] >> output 10,000 observations, but I ran into trouble transforming it from a >> csv file to a dat file using impexc2bd. My input file is >> arhmm_example_10000.csv [3]. I tried to convert to .dat using: >> >> impexc2bd('arhmm_example_10000.csv', ';', 'arhmm_example_10000.dat') >> >> I receive this output, and the .dat file in NOT created >> >> !--error 10000 >> dates are entered neither in chronlogical order nor in reverseerse >> chronological order >> at line 87 of function read_dates called by : >> at line 188 of function impexc2bd called by : >> arhmm_example_10000.dat' >> >> How do you transform a large .csv file to .dat? Is there some problem >> with the data file, the way I am making it? >> >> [0]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_280.dat >> [1]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_280.csv >> [2]: https://github.com/bmbouter/koa_marhn >> [3]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_10000.csv >> [4]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_1000.csv >> [5]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_1000.dat >> >> This e-mail is way to long, but I wanted to fully recap the issues I've >> considered as I use ms_var(). Thanks for any help you can provide. You can >> also find me in the #scilab channel on the OFTC IRC servers. >> >> Thanks, >> Brian >> >> >> >> On Fri, May 23, 2014 at 3:48 PM, Eric Dubois >> wrote: >> >>> Dear Brian >>> >>> 1) For me it works fine; I suspect that you have run ms_var with your >>> data after having run with the manual example. I expect that if you reopen >>> Scialb and run your problem, starting from load('arhmm_example.dat'), >>> it will work. A good advice anyway is to set the bounds before each >>> estimation or to run: >>> --> bounds() >>> if you want to use the greatest available time span with your data >>> 2) I agree; indeed the example in the manual as well as yours are >>> univariate, but ms_var also works with multivariate series (I have run some >>> tests which worked well) >>> 3) If you want to contribute, do not hesitate to send me code (at >>> grocer.toolbox at gmail.com or grocer.toolbox at free.fr); add your >>> copyright; if you can create a help file it would still be better (I have >>> some tools do help doing that if you want them); and if you can add to the >>> manual, it would be marvellous! >>> If you want to imporve the docs you are also welcome; I can sned you the >>> OpenOffice files you need if you find it suitable >>> >>> ?ric. >>> >>> >>> >>> 2014-05-23 14:13 GMT+02:00 Brian Bouterse : >>> >>> Thanks for the reply Eric. It is great to get hints and suggestions from >>>> the author directly! >>>> >>>> I used the commands that you outlined, and they were able to reproduce >>>> the expected output verbatim, which is great. Thanks for clearing that up >>>> for me. I've gotten further towards my goal. >>>> >>>> I've now got three questions: >>>> >>>> 1) I get an unexpected result when I run ms_var() on my own data. I >>>> run these commands: >>>> >>>> load('arhmm_example.dat') >>>> nb_states=2 >>>> switch_var=2 // variances are switching >>>> var_opt=3 // unrestricted var-cov matrix >>>> >>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>> >>>> I receive this output: >>>> >>>> WARNING: in overlay, series number 2 has been ignored because of a bad >>>> frequency >>>> !--error 10000 >>>> series ends before the end date of the bounds >>>> at line 39 of function ts2vec0 called by : >>>> at line 101 of function explone called by : >>>> at line 253 of function ms_var called by : >>>> >>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>> >>>> The arhmm_example.dat file is available here[0], and it was made by >>>> running.the following command on the original csv[1] file arhmm_example.csv. >>>> >>>> impexc2bd('arhmm_example.csv', ';', 'arhmm_example.dat') >>>> >>>> I believe I either don't have the dates configured correctly, or it >>>> requires a specific number of data points to match the frequency value, >>>> which also may be wrong. Do you have some insight into this error message? >>>> I've been reading the docs on the ts structure, and I will continue to try >>>> to solve this roadblock. >>>> >>>> >>>> 2) My goal in doing all this is to analyze Autoregressive Hidden Markov >>>> Models. As I understand it, the VAR-HMM that ms_var provides is a >>>> multivariate case of an Autoregressive Hidden Markov Model. The terms >>>> Markov Switching Model, and Hidden Markov Model refer to the same thing. >>>> Using a single variable with ms_var() as I show above in the example, will >>>> simulate an AR-HMM(3). I would like to check if these statements agree with >>>> your understanding. >>>> >>>> >>>> 3) How could I contribute to the grocer code. At the very least I could >>>> improve the docs some. >>>> >>>> >>>> [0]: >>>> https://s3.amazonaws.com/dfsklfdsklfds/fdsjkfsdjkfds/arhmm_example.dat >>>> [1]: >>>> https://s3.amazonaws.com/dfsklfdsklfds/fdsjkfsdjkfds/arhmm_example.csv >>>> >>>> Thanks, >>>> Brian >>>> >>>> >>>> On Tue, May 20, 2014 at 3:55 PM, Eric Dubois >>>> wrote: >>>> >>>>> Hello Brian >>>>> >>>>> Sorry for this late answer, but I have been quite busy these days. >>>>> >>>>> I did not notice the problem with this version of the MS programs. >>>>> Indeed I have changed the optimization device of all GROCER programs and I >>>>> have not adapted the defaults for the MS programs. >>>>> >>>>> If you run: >>>>> --> global GROCERDIR; >>>>> --> load(GROCERDIR+'\data\us_revu.dat') >>>>> --> bounds('1967m4','2004m2') >>>>> >>>>> --> nb_states=2 >>>>> --> switch_var=2 // variances are switching >>>>> --> var_opt=3 // unrestricted var-cov matrix >>>>> >>>>> --> >>>>> r=ms_var('cte',3,'100*(log(us_revu)-lagts(2,log(us_revu)))',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>> (see chapter 6 of the manual for explanations) >>>>> >>>>> Then the results of the ms_var demo is restaured. >>>>> >>>>> I will change the default in Grocer next version . >>>>> >>>>> Regards. >>>>> >>>>> ?ric. >>>>> >>>>> >>>>> 2014-05-19 12:56 GMT+02:00 Brian Bouterse : >>>>> >>>>>> Hi Scilab community! >>>>>> >>>>>> I'm new to Scilab, and the AR-HMM and VAR-HMM solving capabilities of >>>>>> GROCER are what interest me. >>>>>> >>>>>> I have a question relating to Chapter 23 from the GROCER manual[0]. >>>>>> This is the univariate MS-AR(3) solved using the function ms_reg_d() on >>>>>> the us_revu.dat data included with GROCER. I have made no adjustment from >>>>>> the example statements in Chapter 23. >>>>>> >>>>>> The example output is shown on pages 4 and 5 of the Chapter 23 >>>>>> module. Compare that against the output I receive. >>>>>> >>>>>> http://fpaste.org/102978/14004958/ >>>>>> >>>>>> Here are my questions: >>>>>> >>>>>> 1. The numerical output is completely different. I expected it to >>>>>> be the same since the data is provided by GROCER, and I've done the example >>>>>> exactly as shown in Chapter 23. Is there some explanation to why the >>>>>> solved solution I receive is different than the example output in the >>>>>> chapter? >>>>>> >>>>>> 2. I see output like %i*8.4469016 which seems like an error because >>>>>> %i looks like a variable that yet needs to be replaced, and then multiplied >>>>>> to get to its final value. Is this some kind of bug or error? >>>>>> >>>>>> Thanks for any help the community can provide. We'll be using this >>>>>> for a seminar on HMM, AR-HMM, and VAR-HMM at North Carolina State >>>>>> University. I'm also a developer, so I really appreciate all the effort >>>>>> that has been put into scilab and GROCER. >>>>>> >>>>>> Thanks, >>>>>> Brian >>>>>> >>>>>> >>>>>> [0]: http://dubois.ensae.net/Grocer_manual_v1.6.zip >>>>>> >>>>>> >>>>>> -- >>>>>> Brian Bouterse >>>>>> >>>>>> _______________________________________________ >>>>>> 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 >>>>> >>>>> >>>> >>>> >>>> -- >>>> Brian Bouterse >>>> >>>> _______________________________________________ >>>> 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 >>> >>> >> >> >> -- >> Brian Bouterse >> >> _______________________________________________ >> 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 > > -- Brian Bouterse -------------- next part -------------- An HTML attachment was scrubbed... URL: From sterg at teikav.edu.gr Tue Jun 3 12:37:01 2014 From: sterg at teikav.edu.gr (Stergios Papadimitriou) Date: Tue, 03 Jun 2014 13:37:01 +0300 Subject: [Scilab-users] ScalaLab and GroovyLab support for caling SciLab In-Reply-To: <1401790043.2431.23.camel@paros> References: <538B2A9F.1030108@teikav.edu.gr> <1401788083.2431.13.camel@paros> <538D9986.1080909@teikav.edu.gr> <1401790043.2431.23.camel@paros> Message-ID: <538DA54D.4070107@teikav.edu.gr> Good stuff!! with my academic and research interests been displaced from Computational Intelligence towards Compilers and Operating Systems, I consider as very interested to study the SciLab sources, and perhaps to be involved in parts of the SciLab core, for the Linux platform that I like and I know more!! .. I need to admit however, that Windows SciLab installation is much easier! Stergios ???? 3/6/2014 1:07 ??, ?/? Cl?ment David ??????: > JIMS is part of Scilab 5.5.0 and documented in the Scilab help. There is > also a getting started page : > http://help.scilab.org/docs/5.5.0/en_US/jims-getting-started.html > > -- > Cl?ment > > Le mardi 03 juin 2014 ? 12:46 +0300, Stergios Papadimitriou a ?crit : >> Hi Cl?ment, >> >> I need to study the Scilab-Java interface in order to improve also the >> interface of Scilab - ScalaLab/GroovyLab. >> (... this week I'm very pressed, perhaps I can do more next week!) >> >> Do you know (or somebody else) where to find good staff on SciLab-Java >> interfacing? >> >> I do not know JIMS, a possible link for it? >> >> >> .. perhaps, it is good to have both the compactness of SciLab and the >> general purpose scripting of Scala/Groovy, >> this is why I like to improve on that interfacing work ... >> >> Best Regards >> Stergios >> >> >> >> ???? 3/6/2014 12:34 ??, ?/? Cl?ment David ??????: >>> Hi Stergios, >>> >>> Great job ! I did not known scalalab nor jlabgroovy before. >>> >>> Do you know that Scilab can call any Java library using JIMS ? Do you >>> think that it may be easily possible to use groovy are scala instead ? >>> >>> -- >>> Cl?ment >>> >>> Le dimanche 01 juin 2014 ? 16:29 +0300, Stergios Papadimitriou a ?crit : >>>> Hi all, >>>> >>>> I like to inform the group (mostly concerns the Java funs of it) >>>> that I interfaced my projects ScalaLab and GroovyLab to >>>> cooperate with SciLab scripts: >>>> >>>> https://code.google.com/p/scalalab/wiki/SciLabInScalaLab >>>> >>>> https://code.google.com/p/jlabgroovy/wiki/SciLabInGroovyLab >>>> >>>> >>>> Best Regards >>>> Stergios >>>> _______________________________________________ >>>> 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 > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From sdr at durietz.se Tue Jun 3 15:29:54 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Tue, 03 Jun 2014 15:29:54 +0200 Subject: [Scilab-users] Figure In-Reply-To: <538CBCF7.4010601@durietz.se> References: <53838A7A.2080405@durietz.se> <53839BB2.3020005@free.fr> <538639E7.7080102@durietz.se> <538C62F4.8020109@free.fr> <538CAE3A.5090904@durietz.se> <538CBCF7.4010601@durietz.se> Message-ID: <538DCDD2.10504@durietz.se> On 2014-06-02 20:05, Stefan Du Rietz wrote: -------------------- > On 2014-06-02 19:02, Stefan Du Rietz wrote: > -------------------- >> On 2014-06-02 13:41, Samuel Gougeon wrote: >> -------------------- >>> Le 28/05/2014 21:32, Stefan Du Rietz a ?crit : >>>> On 2014-05-26 21:53, Samuel Gougeon wrote: >>>>> Le 26/05/2014 20:39, Stefan Du Rietz a ?crit : >>>>>> Hello, >>>>>> I am trying to get my figures from 5.4.1 in 5.5.0. >>>>>> >>>>>> Why is everything red? (no help found and no answer to an equal >>>>>> question on the Internet) >>>>> Like here?: >>>>> http://bugzilla.scilab.org/attachment.cgi?id=3561 >>>>> http://bugzilla.scilab.org/show_bug.cgi?id=12788 >>>>> >>>> Thank you, Samuel. I have made a comment to your bug report. >>>> >>>> Another question: >>>> Is it possible to get rid of the extra title bar in a figure? It >>>> just repeats the text in the main title bar and makes the axes space >>>> lower, to no avail whatsoever. >>> It is used as handle for the docking system. When a window is set as >>> .dockable="off", this bar could be masked (since then it becomes >>> useless). >>> Scilab 5.5.0 allows to create a window without this bar, through the >>> new function createWindow(). >>> >> Excellent! Thank you very much! >> >> Is there any help on docking? >> > And how can I plot anything in this figure window? gcf() returns a > handle to this figure, but newaxes() and plot() creates a new window ... > > However, I can put GUI:s in it. > Thanks Antoine for your reply on 2014-06-03 08:53 to Re: [Scilab-users] createWindow and uicontrol That also explains why I couldn't plot anything. Apparently I must have "default_axes","on" Stefan From grocer.toolbox at gmail.com Tue Jun 3 22:40:24 2014 From: grocer.toolbox at gmail.com (Eric Dubois) Date: Tue, 3 Jun 2014 22:40:24 +0200 Subject: [Scilab-users] GROCER VAR-HMM Example Output Different from Chapter 23 In-Reply-To: References: Message-ID: Hello Brian. Thank you for your feedback. I will try to make the distinction between 'cte' and 'all' more precise in the manual. I will also think about your suggestion to put the code on a hub. Regards. ?ric. 2014-06-03 12:21 GMT+02:00 Brian Bouterse : > Hi Eric, > > Thanks for the quick and helpful reply! I adjusted the first parameter > from 'cte' to 'all' and that produces the correct form of output, and the > values look close to the expected values. This is great! > > I've read the help documents several times, but I didn't realize the > impact of adjusting the value of the first parameter. The doc example on > the ms_var help page describes the parameter, but the only example on that > page uses 'cte'. There is no example with 'all' on the ms_var() help page. > Also including some sample output may also be good there. I would like to > contribute an example like this to the docs, manual, or help files. I think > that would have caused me to solve my own problems instead of writing long > e-mails. > > I also tried using a larger data file 10,000 observations in a column > format, and then impexc2bd transformed the data file correctly. That is > also solved thanks to your suggestion. > > I know the code is available for download, but would you ever consider > putting a copy of the GROCER code in a repository on github.com? If it > were in a place like that, I, and others, would have an easier time of > contributing to the codebase instead of e-mailing patches. What do you > think about an idea like this? Just a friendly suggestion based on what has > worked well for me in the past. I can help put it up there, and configure > it if that is helpful. > > Thanks again for making this great software, and helping me use it! > > Best, > Brian > > > > > > On Thu, May 29, 2014 at 1:41 PM, Eric Dubois > wrote: > >> Hello Brian >> >> You did not recover the intial parameters because you did not estimate >> the same model as the one you generated. >> >> This is because you used in ms_var the option 'cte' which means that only >> the constant is allowed to switch: and this is why the AR coeffcients are >> only given for 'all regimes'. I f you rather use the option 'all', then all >> parameters are allowed to switch and the result is now consistent with the >> parameters of the generated model, that is the true values are within the >> confidence interval of the estimated ones. This is explained in the help >> command (see help ms_var): may I ask you why you did not infer it from the >> help files? >> >> As for the big file, I suspect it is truncated into Scilab. Can you try >> with data in column instead of in rows? >> >> ?ric. >> >> >> 2014-05-29 13:24 GMT+02:00 Brian Bouterse : >> >> Hi Eric, >>> >>> Thanks again for the response. What you told me to do worked. If I >>> closed Scilab, and rerun the example from my earlier e-mail, I do not see >>> any warnings, which is good. The output I received is not what I expected >>> in several ways, so I'll ask some more questions here. There are a lot of >>> different questions here, but they all relate to the example below, and the >>> use of ms_var and GROCER. Thanks in advance for any insight into the issues >>> I am experiencing. >>> >>> The example: >>> >>> I have a file arhmm_example_280.dat [0] that was generated and loaded >>> from arhmm_example_280.csv [1] using the following commands: >>> >>> impexc2bd('arhmm_example_280.csv', ';', 'arhmm_example_280.dat') >>> load('arhmm_example_280.dat') >>> >>> I then run, what I expect to be a single variable, 2-regime, >>> >>> nb_states=2 >>> switch_var=2 >>> var_opt=3 >>> >>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>> >>> I then receive this output: http://ur1.ca/hej98 >>> >>> 1) In the "coefficients" section of the output, the output doesn't seem >>> to contain enough numbers to describe the expected number of outputs. The >>> model has 2 regimes, and it is autoregressive with order 3, so I expect >>> there to be 3 coefficients for EACH regime (regime 1 and regime 2). >>> Instead, I see a single coefficient for Regime 1 (0.9313736) and a single >>> coefficient for Regime 2 (8.2685973). Where are the three autoregressive >>> parameters for regime 1 and the other three autoregressive parameters for >>> regime 2? >>> >>> 2) There is a section in the output called "All Regimes" which is in the >>> "Coefficients" section contains 3 coefficients. Do these correspond to the >>> autoregressive aspect of the model [0.2258037, 0.4757391, -0.3161730). I >>> expected to see a section like this for Regime 1, and one for Regime 2, but >>> I did not. I'm not sure how these 3 coefficients for "All Regimes" fit in. >>> The Markov Switching Model switches the autoregressive coefficients with >>> each Regime, so I'm not sure how to interpret coefficients that are for >>> "both regimes". >>> >>> 3) I believe each regime will have its error term estimated, and I do >>> see that in the output. For example under Regime 2, the variance of the >>> error term is 2.3164941. This is listed in the section "Variance-covariance >>> matrix of residuals" which is a subsection under "Regime 2". Am I >>> interpreting the output correctly? Where is the mean value for the regime? >>> >>> 4) My overall goal is to validate that I can hand a list of >>> observations to GROCER+scilab, along with the order of the autoregressive >>> model, and the number of regimes, and have it estimate the Markov hidden >>> state transition probabilities, the estimated regime that is in effect at >>> time t, and coefficients for each regime, and the error term that can be >>> used for prediction of the t+1 observation. I don't know of a "well known" >>> data set that allows for this type of validation, so I wrote a simple >>> AR-HMM generator named koa marhn [2], which outputs AR-HMM data that >>> follows a set of parameters. This is the program that generated the 280 >>> point data set in the arhmm_example_280.csv [1] file. You can see the >>> parameters used to generate this file here: >>> >>> https://github.com/bmbouter/koa_marhn/blob/master/arhmm_example.py >>> >>> A[i][j] is the probability transition from regime i to regime j. >>> C[x][y] is the autoregressive coefficient for regime x used for >>> generation of the t+1 element from the t, t-1, t-2, ..., t-y summation. >>> Each y is a coefficient for regime x. >>> R[q] contains mean and std_deviation values for regime q >>> pi[m] contains the starting state probabilities that the regime is in >>> state m at time t=0 >>> >>> The output from http://ur1.ca/hej98 does not correspond at all with the >>> values of these inputs, which means something unexpected is going on! Do >>> you have insight into why scilab doesn't produce these expected numbers? >>> >>> 5) One possibility is that maybe it is "working", but there isn't >>> enough data for it to estimate the parameters correctly. I had koa marhn >>> [2] output 10,000 observations, but I ran into trouble transforming it from >>> a csv file to a dat file using impexc2bd. My input file is >>> arhmm_example_10000.csv [3]. I tried to convert to .dat using: >>> >>> impexc2bd('arhmm_example_10000.csv', ';', 'arhmm_example_10000.dat') >>> >>> I receive this output, and the .dat file in NOT created >>> >>> !--error 10000 >>> dates are entered neither in chronlogical order nor in reverseerse >>> chronological order >>> at line 87 of function read_dates called by : >>> at line 188 of function impexc2bd called by : >>> arhmm_example_10000.dat' >>> >>> How do you transform a large .csv file to .dat? Is there some problem >>> with the data file, the way I am making it? >>> >>> [0]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_280.dat >>> [1]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_280.csv >>> [2]: https://github.com/bmbouter/koa_marhn >>> [3]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_10000.csv >>> [4]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_1000.csv >>> [5]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_1000.dat >>> >>> This e-mail is way to long, but I wanted to fully recap the issues I've >>> considered as I use ms_var(). Thanks for any help you can provide. You can >>> also find me in the #scilab channel on the OFTC IRC servers. >>> >>> Thanks, >>> Brian >>> >>> >>> >>> On Fri, May 23, 2014 at 3:48 PM, Eric Dubois >>> wrote: >>> >>>> Dear Brian >>>> >>>> 1) For me it works fine; I suspect that you have run ms_var with your >>>> data after having run with the manual example. I expect that if you reopen >>>> Scialb and run your problem, starting from load('arhmm_example.dat'), >>>> it will work. A good advice anyway is to set the bounds before each >>>> estimation or to run: >>>> --> bounds() >>>> if you want to use the greatest available time span with your data >>>> 2) I agree; indeed the example in the manual as well as yours are >>>> univariate, but ms_var also works with multivariate series (I have run some >>>> tests which worked well) >>>> 3) If you want to contribute, do not hesitate to send me code (at >>>> grocer.toolbox at gmail.com or grocer.toolbox at free.fr); add your >>>> copyright; if you can create a help file it would still be better (I have >>>> some tools do help doing that if you want them); and if you can add to the >>>> manual, it would be marvellous! >>>> If you want to imporve the docs you are also welcome; I can sned you >>>> the OpenOffice files you need if you find it suitable >>>> >>>> ?ric. >>>> >>>> >>>> >>>> 2014-05-23 14:13 GMT+02:00 Brian Bouterse : >>>> >>>> Thanks for the reply Eric. It is great to get hints and suggestions >>>>> from the author directly! >>>>> >>>>> I used the commands that you outlined, and they were able to reproduce >>>>> the expected output verbatim, which is great. Thanks for clearing that up >>>>> for me. I've gotten further towards my goal. >>>>> >>>>> I've now got three questions: >>>>> >>>>> 1) I get an unexpected result when I run ms_var() on my own data. I >>>>> run these commands: >>>>> >>>>> load('arhmm_example.dat') >>>>> nb_states=2 >>>>> switch_var=2 // variances are switching >>>>> var_opt=3 // unrestricted var-cov matrix >>>>> >>>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>> >>>>> I receive this output: >>>>> >>>>> WARNING: in overlay, series number 2 has been ignored because of a bad >>>>> frequency >>>>> !--error 10000 >>>>> series ends before the end date of the bounds >>>>> at line 39 of function ts2vec0 called by : >>>>> at line 101 of function explone called by : >>>>> at line 253 of function ms_var called by : >>>>> >>>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>> >>>>> The arhmm_example.dat file is available here[0], and it was made by >>>>> running.the following command on the original csv[1] file arhmm_example.csv. >>>>> >>>>> impexc2bd('arhmm_example.csv', ';', 'arhmm_example.dat') >>>>> >>>>> I believe I either don't have the dates configured correctly, or it >>>>> requires a specific number of data points to match the frequency value, >>>>> which also may be wrong. Do you have some insight into this error message? >>>>> I've been reading the docs on the ts structure, and I will continue to try >>>>> to solve this roadblock. >>>>> >>>>> >>>>> 2) My goal in doing all this is to analyze Autoregressive Hidden >>>>> Markov Models. As I understand it, the VAR-HMM that ms_var provides is a >>>>> multivariate case of an Autoregressive Hidden Markov Model. The terms >>>>> Markov Switching Model, and Hidden Markov Model refer to the same thing. >>>>> Using a single variable with ms_var() as I show above in the example, will >>>>> simulate an AR-HMM(3). I would like to check if these statements agree with >>>>> your understanding. >>>>> >>>>> >>>>> 3) How could I contribute to the grocer code. At the very least I >>>>> could improve the docs some. >>>>> >>>>> >>>>> [0]: >>>>> https://s3.amazonaws.com/dfsklfdsklfds/fdsjkfsdjkfds/arhmm_example.dat >>>>> [1]: >>>>> https://s3.amazonaws.com/dfsklfdsklfds/fdsjkfsdjkfds/arhmm_example.csv >>>>> >>>>> Thanks, >>>>> Brian >>>>> >>>>> >>>>> On Tue, May 20, 2014 at 3:55 PM, Eric Dubois >>>> > wrote: >>>>> >>>>>> Hello Brian >>>>>> >>>>>> Sorry for this late answer, but I have been quite busy these days. >>>>>> >>>>>> I did not notice the problem with this version of the MS programs. >>>>>> Indeed I have changed the optimization device of all GROCER programs and I >>>>>> have not adapted the defaults for the MS programs. >>>>>> >>>>>> If you run: >>>>>> --> global GROCERDIR; >>>>>> --> load(GROCERDIR+'\data\us_revu.dat') >>>>>> --> bounds('1967m4','2004m2') >>>>>> >>>>>> --> nb_states=2 >>>>>> --> switch_var=2 // variances are switching >>>>>> --> var_opt=3 // unrestricted var-cov matrix >>>>>> >>>>>> --> >>>>>> r=ms_var('cte',3,'100*(log(us_revu)-lagts(2,log(us_revu)))',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>>> (see chapter 6 of the manual for explanations) >>>>>> >>>>>> Then the results of the ms_var demo is restaured. >>>>>> >>>>>> I will change the default in Grocer next version . >>>>>> >>>>>> Regards. >>>>>> >>>>>> ?ric. >>>>>> >>>>>> >>>>>> 2014-05-19 12:56 GMT+02:00 Brian Bouterse : >>>>>> >>>>>>> Hi Scilab community! >>>>>>> >>>>>>> I'm new to Scilab, and the AR-HMM and VAR-HMM solving capabilities >>>>>>> of GROCER are what interest me. >>>>>>> >>>>>>> I have a question relating to Chapter 23 from the GROCER manual[0]. >>>>>>> This is the univariate MS-AR(3) solved using the function ms_reg_d() on >>>>>>> the us_revu.dat data included with GROCER. I have made no adjustment from >>>>>>> the example statements in Chapter 23. >>>>>>> >>>>>>> The example output is shown on pages 4 and 5 of the Chapter 23 >>>>>>> module. Compare that against the output I receive. >>>>>>> >>>>>>> http://fpaste.org/102978/14004958/ >>>>>>> >>>>>>> Here are my questions: >>>>>>> >>>>>>> 1. The numerical output is completely different. I expected it to >>>>>>> be the same since the data is provided by GROCER, and I've done the example >>>>>>> exactly as shown in Chapter 23. Is there some explanation to why the >>>>>>> solved solution I receive is different than the example output in the >>>>>>> chapter? >>>>>>> >>>>>>> 2. I see output like %i*8.4469016 which seems like an error because >>>>>>> %i looks like a variable that yet needs to be replaced, and then multiplied >>>>>>> to get to its final value. Is this some kind of bug or error? >>>>>>> >>>>>>> Thanks for any help the community can provide. We'll be using this >>>>>>> for a seminar on HMM, AR-HMM, and VAR-HMM at North Carolina State >>>>>>> University. I'm also a developer, so I really appreciate all the effort >>>>>>> that has been put into scilab and GROCER. >>>>>>> >>>>>>> Thanks, >>>>>>> Brian >>>>>>> >>>>>>> >>>>>>> [0]: http://dubois.ensae.net/Grocer_manual_v1.6.zip >>>>>>> >>>>>>> >>>>>>> -- >>>>>>> Brian Bouterse >>>>>>> >>>>>>> _______________________________________________ >>>>>>> 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 >>>>>> >>>>>> >>>>> >>>>> >>>>> -- >>>>> Brian Bouterse >>>>> >>>>> _______________________________________________ >>>>> 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 >>>> >>>> >>> >>> >>> -- >>> Brian Bouterse >>> >>> _______________________________________________ >>> 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 >> >> > > > -- > Brian Bouterse > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From tanthiamhuat at yahoo.com Tue Jun 3 10:29:51 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Tue, 3 Jun 2014 01:29:51 -0700 (PDT) Subject: [Scilab-users] Scilab 5.5.0 integrate with Symbolic Toolbox Message-ID: <1401784191878-4030657.post@n3.nabble.com> Can I know how Scilab 5.5.0 in Windows 7 can be integrated with Symbolic Toolbox? for example, http://es.wikibooks.org/wiki/Manual_de_Scilab/Xcos/C%C3%A1lculo_simb%C3%B3lico_con_SciMax -- View this message in context: http://mailinglists.scilab.org/Scilab-5-5-0-integrate-with-Symbolic-Toolbox-tp4030657.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From mladen.meduric at gmail.com Tue Jun 3 13:17:07 2014 From: mladen.meduric at gmail.com (mladen) Date: Tue, 3 Jun 2014 04:17:07 -0700 (PDT) Subject: [Scilab-users] plot x,y,z point with color In-Reply-To: <1401585474365-4030640.post@n3.nabble.com> References: <1401585474365-4030640.post@n3.nabble.com> Message-ID: <1401794227560-4030664.post@n3.nabble.com> Anyone? It's a 4d plot query. I've tried plot3d(x,y,list(z,p),) to no avail. Kind regards. -- View this message in context: http://mailinglists.scilab.org/plot-x-y-z-point-with-color-tp4030640p4030664.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sergio.pocas at energytop.pt Tue Jun 3 20:02:13 2014 From: sergio.pocas at energytop.pt (=?iso-8859-1?Q?S=E9rgio_Po=E7as?=) Date: Tue, 3 Jun 2014 19:02:13 +0100 Subject: [Scilab-users] Eletrical diagram with XCOS Message-ID: Anyone can help me with this project? I have one error that I can?t solve. S?rgio Po?as Departamento de Automa??o Tlm: 937 825 257 sergio.pocas at energytop.pt cid:image002.jpg at 01CDD9E4.165D0C20 ?h? mais de 10 anos ao servi?o da Energia? cid:image002.png at 01CF2E62.6FFC6B80 Rua Sto Ant?nio, 710/742 Armaz?m J Tel: 229 021 086 | Fax: 229 023 249 4475 ? 611- Sta Maria de Avioso ? Maia www.energytop.pt | energytop at energytop.pt cid:image003.jpg at 01CDD9E3.119874A0 Se puder, n?o imprima! A Natureza agradece. ________________________________________________________________________ Toda a informa??o contida nesta mensagem ? confidencial e dirigida apenas aos destinat?rios. Se n?o ? destinat?rio, n?o fa?a uso da informa??o nela contida. Por favor apague-a e informe o remetente. O uso n?o autorizado, divulga??o ou altera??o desta mensagem ? estritamente proibido. A EnergyTop n?o aceita responsabilidade por danos causados pela rece??o incorreta desta mensagem. -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image001.jpg Type: image/jpeg Size: 3188 bytes Desc: image001.jpg URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image002.png Type: image/png Size: 96865 bytes Desc: image002.png URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image003.jpg Type: image/jpeg Size: 941 bytes Desc: image003.jpg URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Esquema RL RC Sergio.zcos Type: application/octet-stream Size: 14442 bytes Desc: Esquema RL RC Sergio.zcos URL: From Christophe.Dang at sidel.com Wed Jun 4 09:22:54 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 4 Jun 2014 09:22:54 +0200 Subject: [Scilab-users] plot x,y,z point with color In-Reply-To: <1401794227560-4030664.post@n3.nabble.com> References: <1401585474365-4030640.post@n3.nabble.com> <1401794227560-4030664.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32B4170A0@301EX00100.sidel.com> Hello, > De : mladen > Envoy? : mardi 3 juin 2014 13:17 > > It's a 4d plot query. Do you mean you have a (x, y, z, t) vector, and you want a colour scale for t? If so, maybe param3d1 with the "list(z, colours)" option http://help.scilab.org/docs/5.5.0/en_US/param3d1.html Probably a bit painful... I think there was an example given on the mailing list, I'll try to find it on the archive (or if someone remembers). -- 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 Wed Jun 4 09:25:48 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 4 Jun 2014 09:25:48 +0200 Subject: [Scilab-users] plot x,y,z point with color In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F32B4170A0@301EX00100.sidel.com> References: <1401585474365-4030640.post@n3.nabble.com> <1401794227560-4030664.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32B4170A0@301EX00100.sidel.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32B4170B1@301EX00100.sidel.com> > I think there was an example given on the mailing list Got it: http://mailinglists.scilab.org/Scilab-users-RGB-color-plot-trying-to-color-each-mark-with-its-respective-color-td4027586.html HTH 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 clement.david at scilab-enterprises.com Wed Jun 4 09:31:16 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Wed, 04 Jun 2014 09:31:16 +0200 Subject: [Scilab-users] ScalaLab and GroovyLab support for caling SciLab In-Reply-To: <538DA54D.4070107@teikav.edu.gr> References: <538B2A9F.1030108@teikav.edu.gr> <1401788083.2431.13.camel@paros> <538D9986.1080909@teikav.edu.gr> <1401790043.2431.23.camel@paros> <538DA54D.4070107@teikav.edu.gr> Message-ID: <1401867076.2241.4.camel@paros> Hello, Feel free to join the low-traffic dev mailing-list and ask question about scilab internals there. If you have issue compiling Scilab (regardless of your operating system), just ask. -- Cl?ment Le mardi 03 juin 2014 ? 13:37 +0300, Stergios Papadimitriou a ?crit : > Good stuff!! > > with my academic and research interests been displaced from > Computational Intelligence > towards Compilers and Operating Systems, > I consider as very interested to study the SciLab sources, > and perhaps to be involved in parts of the SciLab core, > for the Linux platform that I like and I know more!! > > .. I need to admit however, that Windows SciLab installation is much > easier! > > > Stergios > > ???? 3/6/2014 1:07 ??, ?/? Cl?ment David ??????: > > JIMS is part of Scilab 5.5.0 and documented in the Scilab help. There is > > also a getting started page : > > http://help.scilab.org/docs/5.5.0/en_US/jims-getting-started.html > > > > -- > > Cl?ment > > > > Le mardi 03 juin 2014 ? 12:46 +0300, Stergios Papadimitriou a ?crit : > >> Hi Cl?ment, > >> > >> I need to study the Scilab-Java interface in order to improve also the > >> interface of Scilab - ScalaLab/GroovyLab. > >> (... this week I'm very pressed, perhaps I can do more next week!) > >> > >> Do you know (or somebody else) where to find good staff on SciLab-Java > >> interfacing? > >> > >> I do not know JIMS, a possible link for it? > >> > >> > >> .. perhaps, it is good to have both the compactness of SciLab and the > >> general purpose scripting of Scala/Groovy, > >> this is why I like to improve on that interfacing work ... > >> > >> Best Regards > >> Stergios > >> > >> > >> > >> ???? 3/6/2014 12:34 ??, ?/? Cl?ment David ??????: > >>> Hi Stergios, > >>> > >>> Great job ! I did not known scalalab nor jlabgroovy before. > >>> > >>> Do you know that Scilab can call any Java library using JIMS ? Do you > >>> think that it may be easily possible to use groovy are scala instead ? > >>> > >>> -- > >>> Cl?ment > >>> > >>> Le dimanche 01 juin 2014 ? 16:29 +0300, Stergios Papadimitriou a ?crit : > >>>> Hi all, > >>>> > >>>> I like to inform the group (mostly concerns the Java funs of it) > >>>> that I interfaced my projects ScalaLab and GroovyLab to > >>>> cooperate with SciLab scripts: > >>>> > >>>> https://code.google.com/p/scalalab/wiki/SciLabInScalaLab > >>>> > >>>> https://code.google.com/p/jlabgroovy/wiki/SciLabInGroovyLab > >>>> > >>>> > >>>> Best Regards > >>>> Stergios > >>>> _______________________________________________ > >>>> 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 > > > > _______________________________________________ > > 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 varun.krk at gmail.com Wed Jun 4 09:26:32 2014 From: varun.krk at gmail.com (Varun Singh) Date: Wed, 4 Jun 2014 12:56:32 +0530 Subject: [Scilab-users] Eucledian Distance between two colors in an image. Message-ID: Hello This is another mail i am sending with regard to the problem which i have been facing from a long time now.I did'nt got the reply of the previous message but if u can please reply this time then it will be very generous. My problem is that I want to find the distance between the green and red color in an image but don't know how to proceed with that. I browsed the web and got to know about the L*a*b* CIE_LAB option to find the eucledian distance and i converted the image in L*a*b* but i dnt know how to proceed next. How to find the eucledian distance after converting the image in that form. Please reply asap. Thankyou Regards Varun Singh -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Wed Jun 4 21:09:39 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Wed, 04 Jun 2014 21:09:39 +0200 Subject: [Scilab-users] Is FPPO available ? Message-ID: <538F6EF3.4000108@gmail.com> Hi there I have generated some 60 seconds of (pink) noise, sampled at 48 kHz, which I can fft to get a linear-frequency representation (predominantly interested in 20 - 20k Hz). If I plot this data on a log-frequency plot, then of course high frequency spectrum plots look very noisy. I need a fixed-point-per-octave (FPPO) transformation of my data so I can study the spectrum better at higher frequencies with a more readable trace of my frequency response data. I would be interested in for example a 1/3 octave moving-average smoothing. In effect, the FPPO technique utilizes a measurement time window that varies as a function of frequency, utilizing a long time window at low frequencies (for narrow frequency resolution) and a successively shorter time window at high frequencies (but averaged through the entire 60 seconds). Is FPPO available in Scilab? (Signal Processing Toolbox maybe?) ... I cannot find it. /Claus P.S. You could see this PDF file for some pictures (go to page 6, Figure 4): http://www.rationalacoustics.com/files/FFT_Fundamentals.pdf --- This email is free from viruses and malware because avast! Antivirus protection is active. http://www.avast.com From sdr at durietz.se Wed Jun 4 21:17:55 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Wed, 04 Jun 2014 21:17:55 +0200 Subject: [Scilab-users] Is FPPO available ? In-Reply-To: <538F6EF3.4000108@gmail.com> References: <538F6EF3.4000108@gmail.com> Message-ID: <538F70E3.2000308@durietz.se> On 2014-06-04 21:09, Claus Futtrup wrote: -------------------- > Hi there > > I have generated some 60 seconds of (pink) noise, sampled at 48 kHz, > which I can fft to get a linear-frequency representation > (predominantly interested in 20 - 20k Hz). > > If I plot this data on a log-frequency plot, then of course high > frequency spectrum plots look very noisy. > > I need a fixed-point-per-octave (FPPO) transformation of my data so I > can study the spectrum better at higher frequencies with a more > readable trace of my frequency response data. I would be interested in > for example a 1/3 octave moving-average smoothing. > > In effect, the FPPO technique utilizes a measurement time window that > varies as a function of frequency, utilizing a long time window at low > frequencies (for narrow frequency resolution) and a successively > shorter time window at high frequencies (but averaged through the > entire 60 seconds). > > Is FPPO available in Scilab? (Signal Processing Toolbox maybe?) ... I > cannot find it. > > /Claus > > P.S. You could see this PDF file for some pictures (go to page 6, > Figure 4): > http://www.rationalacoustics.com/files/FFT_Fundamentals.pdf > What about wavelets? Stefan From tanthiamhuat at yahoo.com Thu Jun 5 04:47:24 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Wed, 4 Jun 2014 19:47:24 -0700 (PDT) Subject: [Scilab-users] Using Scilab Genetic Algorithm to solve Simple Card Problem Message-ID: <1401936444905-4030679.post@n3.nabble.com> Could this Simple Card Problem be included in Scilab Genetic Algorithm (GA), such that users can have a quickstart of using GA in Scilab? E.g, how start the first random initial population, how this population goes through GA operators like Crossover (Recombination) and Mutation, etc. And lastly how the solution can be presented. The Simple Card Problem is very easy to understand, and using GA to solve it is a very good example. It can be referred from http://www.codeproject.com/Articles/16286/AI-Simple-Genetic-Algorithm-GA-to-solve-a-card-pro Prompt actions and efforts would be greatly appreciated. thanks. -- View this message in context: http://mailinglists.scilab.org/Using-Scilab-Genetic-Algorithm-to-solve-Simple-Card-Problem-tp4030679.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From edmund.okoroigwe at unn.edu.ng Thu Jun 5 11:29:03 2014 From: edmund.okoroigwe at unn.edu.ng (Edmund Okoroigwe) Date: Thu, 5 Jun 2014 02:29:03 -0700 Subject: [Scilab-users] Plotting time (hours) in scilab Message-ID: I want to plot time (8:00,9:00,10:00,11:00, 12:00, 13:00,14:00) as x-axis. How can I achieve this. The data is saved in excel file. Thanks Edmund -------------- next part -------------- An HTML attachment was scrubbed... URL: From sdr at durietz.se Thu Jun 5 12:41:54 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Thu, 05 Jun 2014 12:41:54 +0200 Subject: [Scilab-users] Plotting time (hours) in scilab In-Reply-To: References: Message-ID: <53904972.6000609@durietz.se> On 2014-06-05 11:29, Edmund Okoroigwe wrote: -------------------- > I want to plot time (8:00,9:00,10:00,11:00, 12:00, 13:00,14:00) as > x-axis. How can I achieve this. The data is saved in excel file. > Thanks > > Edmund > Maybe like this: t = 8:14; a = gca(); // Set the x-axis limits to the first and last element of t a.tight_limits = "on"; a.data_bounds(:, 1) = t([1 $])'; // Set the ticks and format the labels a.x_ticks = tlist(["ticks","locations","labels"], ... t, msprintf("%02i:00\n", t')); Here the x-axis is scaled in hours. If you have each time (date and time) as a vector [year month day hour minute second] you can convert them with datenum() to get numbers with one day = 1.0 (and one hour = 1/24). See the help on datenum() and msprintf(). Regards Stefan From sdr at durietz.se Thu Jun 5 13:00:13 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Thu, 05 Jun 2014 13:00:13 +0200 Subject: [Scilab-users] Plotting time (hours) in scilab In-Reply-To: <53904972.6000609@durietz.se> References: <53904972.6000609@durietz.se> Message-ID: <53904DBD.9010508@durietz.se> On 2014-06-05 12:41, Stefan Du Rietz wrote: -------------------- > On 2014-06-05 11:29, Edmund Okoroigwe wrote: > -------------------- >> I want to plot time (8:00,9:00,10:00,11:00, 12:00, 13:00,14:00) as >> x-axis. How can I achieve this. The data is saved in excel file. >> Thanks >> >> Edmund >> > Maybe like this: > > t = 8:14; > a = gca(); > // Set the x-axis limits to the first and last element of t > a.tight_limits = "on"; > a.data_bounds(:, 1) = t([1 $])'; > // Set the ticks and format the labels > a.x_ticks = tlist(["ticks","locations","labels"], ... > t, msprintf("%02i:00\n", t')); > > Here the x-axis is scaled in hours. If you have each time (date and > time) as a vector [year month day hour minute second] you can convert > them with datenum() to get numbers with one day = 1.0 (and one hour = > 1/24). See the help on datenum() and msprintf(). Sorry, I mishandled the transpose operators! t = (8:14)'; a = gca(); // Set the x-axis limits to the first and last element of t a.tight_limits = "on"; a.data_bounds(:, 1) = t([1 $])'; // Set the ticks and format the labels a.x_ticks = tlist(["ticks","locations","labels"], ... t, msprintf("%02i:00\n", t)); > Regards > Stefan From samuel.enibe at unn.edu.ng Thu Jun 5 13:46:42 2014 From: samuel.enibe at unn.edu.ng (Samuel Enibe) Date: Thu, 5 Jun 2014 12:46:42 +0100 Subject: [Scilab-users] Plotting time (hours) in scilab In-Reply-To: References: Message-ID: The solution by Stefan is adequate to my understanding. Samuel Enibe On Thu, Jun 5, 2014 at 10:29 AM, Edmund Okoroigwe < edmund.okoroigwe at unn.edu.ng> wrote: > I want to plot time (8:00,9:00,10:00,11:00, 12:00, 13:00,14:00) as > x-axis. How can I achieve this. The data is saved in excel file. > Thanks > > Edmund > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From sdr at durietz.se Thu Jun 5 14:04:56 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Thu, 05 Jun 2014 14:04:56 +0200 Subject: [Scilab-users] Plotting time (hours) in scilab In-Reply-To: <53904DBD.9010508@durietz.se> References: <53904972.6000609@durietz.se> <53904DBD.9010508@durietz.se> Message-ID: <53905CE8.6040401@durietz.se> On 2014-06-05 13:00, Stefan Du Rietz wrote: > On 2014-06-05 12:41, Stefan Du Rietz wrote: > -------------------- >> On 2014-06-05 11:29, Edmund Okoroigwe wrote: >> -------------------- >>> I want to plot time (8:00,9:00,10:00,11:00, 12:00, 13:00,14:00) as >>> x-axis. How can I achieve this. The data is saved in excel file. >>> Thanks >>> >>> Edmund >>> >> Maybe like this: >> >> t = 8:14; >> a = gca(); >> // Set the x-axis limits to the first and last element of t >> a.tight_limits = "on"; >> a.data_bounds(:, 1) = t([1 $])'; >> // Set the ticks and format the labels >> a.x_ticks = tlist(["ticks","locations","labels"], ... >> t, msprintf("%02i:00\n", t')); >> >> Here the x-axis is scaled in hours. If you have each time (date and >> time) as a vector [year month day hour minute second] you can convert >> them with datenum() to get numbers with one day = 1.0 (and one hour = >> 1/24). See the help on datenum() and msprintf(). > > Sorry, I mishandled the transpose operators! > > t = (8:14)'; > a = gca(); > // Set the x-axis limits to the first and last element of t > a.tight_limits = "on"; > a.data_bounds(:, 1) = t([1 $])'; > // Set the ticks and format the labels > a.x_ticks = tlist(["ticks","locations","labels"], ... > t, msprintf("%02i:00\n", t)); > Once again ... After fixing t, I forgot to change to this: a.data_bounds(:, 1) = t([1 $]); >> Regards >> Stefan From cfuttrup at gmail.com Thu Jun 5 19:14:41 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Thu, 05 Jun 2014 19:14:41 +0200 Subject: [Scilab-users] Is FPPO available ? In-Reply-To: <538F70E3.2000308@durietz.se> References: <538F6EF3.4000108@gmail.com> <538F70E3.2000308@durietz.se> Message-ID: <5390A581.7000306@gmail.com> Hi Stefan et al. >What about wavelets? Maybe wavelet toolbox? I am not so familiar with it - which wavelet function can do the per octave smoothing I ask for? /Claus On 04-Jun-14 21:17, Stefan Du Rietz wrote: > On 2014-06-04 21:09, Claus Futtrup wrote: > -------------------- >> Hi there >> >> I have generated some 60 seconds of (pink) noise, sampled at 48 kHz, >> which I can fft to get a linear-frequency representation >> (predominantly interested in 20 - 20k Hz). >> >> If I plot this data on a log-frequency plot, then of course high >> frequency spectrum plots look very noisy. >> >> I need a fixed-point-per-octave (FPPO) transformation of my data so I >> can study the spectrum better at higher frequencies with a more >> readable trace of my frequency response data. I would be interested in >> for example a 1/3 octave moving-average smoothing. >> >> In effect, the FPPO technique utilizes a measurement time window that >> varies as a function of frequency, utilizing a long time window at low >> frequencies (for narrow frequency resolution) and a successively >> shorter time window at high frequencies (but averaged through the >> entire 60 seconds). >> >> Is FPPO available in Scilab? (Signal Processing Toolbox maybe?) ... I >> cannot find it. >> >> /Claus >> >> P.S. You could see this PDF file for some pictures (go to page 6, >> Figure 4): >> http://www.rationalacoustics.com/files/FFT_Fundamentals.pdf >> > What about wavelets? > > Stefan > > > _______________________________________________ > 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 sdr at durietz.se Thu Jun 5 19:58:43 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Thu, 05 Jun 2014 19:58:43 +0200 Subject: [Scilab-users] 5.5.0 axes data_bounds Message-ID: <5390AFD3.7060006@durietz.se> Hello, it seems like there is a new bug in Scilab 5.5.0. If I set a.data_bounds and then plot something, even with xpoly(), they are changed, despite that a.auto_scale=="off" and a.tight_limits=="on". It worked in 5.4.1. Shall I report a bug? Regards Stefan From amonmayr at laas.fr Thu Jun 5 23:54:10 2014 From: amonmayr at laas.fr (Antoine Monmayrant) Date: Thu, 05 Jun 2014 23:54:10 +0200 Subject: [Scilab-users] =?utf-8?q?5=2E5=2E0_axes_data=5Fbounds?= In-Reply-To: <5390AFD3.7060006@durietz.se> Message-ID: <1a63-5390e700-5-b088840@114759739> Le Jeudi 5 Juin 2014 19:58 CEST, Stefan Du Rietz a ?crit: > Hello, > it seems like there is a new bug in Scilab 5.5.0. > > If I set a.data_bounds and then plot something, even with xpoly(), they > are changed, despite that a.auto_scale=="off" and a.tight_limits=="on". > It worked in 5.4.1. > > Shall I report a bug? Hi Stefan, If you are sure that the behaviour has changed between 5.5.0 and the previous version, I think you should report this as a bug! Cheers, Antoine > > Regards > Stefan > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From mladen.meduric at gmail.com Thu Jun 5 12:24:15 2014 From: mladen.meduric at gmail.com (mladen) Date: Thu, 5 Jun 2014 03:24:15 -0700 (PDT) Subject: [Scilab-users] plot x,y,z point with color In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F32B4170B1@301EX00100.sidel.com> References: <1401585474365-4030640.post@n3.nabble.com> <1401794227560-4030664.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32B4170A0@301EX00100.sidel.com> <3B5FFC67498DFF49AE7271A584867D16F32B4170B1@301EX00100.sidel.com> Message-ID: <1401963855121-4030681.post@n3.nabble.com> Thank you Christophe! I think this will work. I just need to overcome following issue (applying your example within my routine): /Invalid index. at line 39 of function generic_i_h called by : at line 2 of function %c_i_h called by : e.children.mark_mode = "on"; at line 46 of exec file called by : exec('D:\work\Scilab\CFD data\q3d point plot_3.sce', -1)/ where line 46 is that command: e.children.mark_mode = "on" Kind regards. -- View this message in context: http://mailinglists.scilab.org/plot-x-y-z-point-with-color-tp4030640p4030681.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From henri.pages at banque-france.fr Fri Jun 6 14:53:51 2014 From: henri.pages at banque-france.fr (henri.pages at banque-france.fr) Date: Fri, 6 Jun 2014 14:53:51 +0200 Subject: [Scilab-users] please unsubscribe Message-ID: ************************************************************** Ce courrier ?lectronique, y compris les pi?ces jointes, est ? l'attention exclusive des destinataires d?sign?s et rev?t un caract?re confidentiel. Si vous recevez ce courrier ?lectronique par erreur, merci d'en informer sans d?lai l'exp?diteur et de supprimer son contenu et ses pi?ces jointes. Le contenu de ce courrier ?lectronique ne pourrait engager la responsabilit? de la Banque de France que s'il a ?t? ?mis par une personne d?ment habilit?e agissant dans le strict cadre des fonctions auxquelles elle est employ?e et ? des fins non ?trang?res ? ses attributions. L'exp?diteur de ce courrier ?lectronique ne peut pas garantir la s?curit? de l'acheminement par voie ?lectronique et ne saurait d?s lors encourir une quelconque responsabilit? en cas d'alt?ration de son contenu. ************************************************************** This e-mail, attachments included, is intended solely for the addressees and should be considered as confidential. Should you receive this message by error, please notify the sender immediately and destroy this e-mail and its attachments. Banque de France cannot be considered as liable for the content of this message unless the sender has been duly authorized and has acted strictly in the course of his/her tasks as an employee. The sender of this e-mail cannot ensure the security of its electronic transmission and consequently will not be liable should its content be altered. ************************************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: From sdr at durietz.se Fri Jun 6 18:43:20 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Fri, 06 Jun 2014 18:43:20 +0200 Subject: [Scilab-users] 5.5.0 axes data_bounds In-Reply-To: <1a63-5390e700-5-b088840@114759739> References: <1a63-5390e700-5-b088840@114759739> Message-ID: <5391EFA8.2040406@durietz.se> On 2014-06-05 23:54, Antoine Monmayrant wrote: > > Le Jeudi 5 Juin 2014 19:58 CEST, Stefan Du Rietz a ?crit: > >> Hello, >> it seems like there is a new bug in Scilab 5.5.0. >> >> If I set a.data_bounds and then plot something, even with xpoly(), they >> are changed, despite that a.auto_scale=="off" and a.tight_limits=="on". >> It worked in 5.4.1. >> >> Shall I report a bug? > > Hi Stefan, > > If you are sure that the behaviour has changed between 5.5.0 and the previous version, I think you should report this as a bug! > > Cheers, > > Antoine > >> OK, Bug # 13462 >> Regards >> Stefan From chuox.bece at gmail.com Sat Jun 7 02:06:33 2014 From: chuox.bece at gmail.com (Chuox) Date: Fri, 6 Jun 2014 17:06:33 -0700 (PDT) Subject: [Scilab-users] Get random number from a matrix Message-ID: <1402099593613-4030691.post@n3.nabble.com> Good evening. I need to get a ramdom value form a set Matrix, for example: Matrix=[25 ,32 ,40, 50] I need to get a random number form that matrix. Thanks -- View this message in context: http://mailinglists.scilab.org/Get-random-number-from-a-matrix-tp4030691.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From chuox.bece at gmail.com Sat Jun 7 02:30:34 2014 From: chuox.bece at gmail.com (Chuox) Date: Fri, 6 Jun 2014 17:30:34 -0700 (PDT) Subject: [Scilab-users] Get random number from a matrix In-Reply-To: <1402099593613-4030691.post@n3.nabble.com> References: <1402099593613-4030691.post@n3.nabble.com> Message-ID: <1402101034175-4030692.post@n3.nabble.com> i found this solution, but i'm guessing there must be a more easy and efficient solution. The scrip is this: Matrix=[25 ,32 ,40, 50] ; tam=size(Matrix); number=grand(1,1,"uin",1,tam(2)); random=Matrix(number) Hope its helps someone. -- View this message in context: http://mailinglists.scilab.org/Get-random-number-from-a-matrix-tp4030691p4030692.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From mail2kristophorus at gmail.com Sat Jun 7 07:41:47 2014 From: mail2kristophorus at gmail.com (hadiono) Date: Fri, 6 Jun 2014 22:41:47 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic Message-ID: <1402119707393-4030693.post@n3.nabble.com> Hello friends.. I'm new to Scilab and doesn't have any experience with Matlab also. Still read and learn the Scilab tutorial from Scilab website and others. I'm doing a research based on interval arithmetic (http://en.wikipedia.org/wiki/Interval_arithmetic) and want to do the interval arithmetic calculation using scilab. Searching with Google and found there is a toolbox/library or similar things for interval arithmetic but some how the development is stop. Can somebody here help me for doing interval arithmetic calculation with scilab? Any link to websites or else are welcome. Regards, Hadiono. -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From mail2kristophorus at gmail.com Sat Jun 7 07:45:31 2014 From: mail2kristophorus at gmail.com (Kristophorus Hadiono) Date: Sat, 7 Jun 2014 12:45:31 +0700 Subject: [Scilab-users] posting Message-ID: -- Regards, Kristophorus -------------- next part -------------- An HTML attachment was scrubbed... URL: From hourglass at free.fr Sat Jun 7 11:17:38 2014 From: hourglass at free.fr (hourglass) Date: Sat, 7 Jun 2014 02:17:38 -0700 (PDT) Subject: [Scilab-users] WScilex.exe does not start anymore In-Reply-To: <1387188954311-4028047.post@n3.nabble.com> References: <1387112657318-4028045.post@n3.nabble.com> <52ADF9CE.8060308@centre-cired.fr> <1387188954311-4028047.post@n3.nabble.com> Message-ID: <1402132658309-4030695.post@n3.nabble.com> Hello, I have a similar problem here. I have been using scilab from a long time and it is not starting anymore when: - I double click sci/sce or - launch from the shortcut (scilab 5.5.0, win 7 64). - double click the Wscilex or scilex executable However, it starts fine from the command line and even from a .bat where the command is "C:\Program Files\scilab-5.5.0\bin\WScilex.exe" I tried deleting all files from the Appdata folder but it did not solve the issue... Any idea? Thanks Philippe -- View this message in context: http://mailinglists.scilab.org/WScilex-exe-does-not-start-anymore-tp4028045p4030695.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Sat Jun 7 13:31:31 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sat, 07 Jun 2014 13:31:31 +0200 Subject: [Scilab-users] Get random number from a matrix In-Reply-To: <1402101034175-4030692.post@n3.nabble.com> References: <1402099593613-4030691.post@n3.nabble.com> <1402101034175-4030692.post@n3.nabble.com> Message-ID: <5392F813.2010104@free.fr> Le 07/06/2014 02:30, Chuox a ?crit : > i found this solution, but i'm guessing there must be a more easy and > efficient solution. not really, may be just a little shorter: Matrix = [25 ,32 ,40, 50] ; random = Matrix(grand(1,1,"uin",1,length(Matrix))) If Matrix is rather of text, length(Matrix) must be replaced with size(Matrix,"*") Regards From akhorshidi at live.com Sun Jun 8 13:20:17 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 8 Jun 2014 04:20:17 -0700 (PDT) Subject: [Scilab-users] Get random number from a matrix In-Reply-To: <5392F813.2010104@free.fr> References: <1402099593613-4030691.post@n3.nabble.com> <1402101034175-4030692.post@n3.nabble.com> <5392F813.2010104@free.fr> Message-ID: <1402226417419-4030697.post@n3.nabble.com> Samuel GOUGEON wrote > Le 07/06/2014 02:30, Chuox a ?crit : >> i found this solution, but i'm guessing there must be a more easy and >> efficient solution. > not really, may be just a little shorter: > > Matrix = [25 ,32 ,40, 50] ; > random = Matrix(grand(1,1,"uin",1,length(Matrix))) > > If Matrix is rather of text, length(Matrix) must be replaced with > size(Matrix,"*") > > Regards Hi all; Samuel, I think it will be handy if we have a function which can extract randomly one, two or more entities of a matrix as well as can also be used to return a random row/column of a matrix. Attached you can see a preliminary version of that function. Regards, Mehran random_extraction.sce _ -- View this message in context: http://mailinglists.scilab.org/Get-random-number-from-a-matrix-tp4030691p4030697.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From akhorshidi at live.com Sun Jun 8 13:29:28 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 8 Jun 2014 04:29:28 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic In-Reply-To: <1402119707393-4030693.post@n3.nabble.com> References: <1402119707393-4030693.post@n3.nabble.com> Message-ID: <1402226968732-4030698.post@n3.nabble.com> Hi; There is a Scilab toolbox for interval analysis but, the toolbox homepage is not available at this time and return 403 access denied error. Unfortunately! :( http://www-sop.inria.fr/coprin/logiciels/Int4Sci/ Could anyone get an alternative link to download this toolbox? BTW, I refer you to the book "Introduction to interval analysis" written by Ramon E. Moore et al. This book has been written based on the use of Matlab's INtLAB toolbox. If one has good knowledge regarding interval analysis as well as Scilab programming, it wouldn't be hard to write her/his own Scilab functions to implement interval arithmetic. Regards, Mehran _ -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030698.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From mail2kristophorus at gmail.com Sun Jun 8 13:54:16 2014 From: mail2kristophorus at gmail.com (hadiono) Date: Sun, 8 Jun 2014 04:54:16 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic In-Reply-To: <1402226968732-4030698.post@n3.nabble.com> References: <1402119707393-4030693.post@n3.nabble.com> <1402226968732-4030698.post@n3.nabble.com> Message-ID: That what I found from google and stuck, even I get the toolbox still won't execute in Scilab. The cause is the version of Scilab and changes in Scilab version. I'm joined to this mailinglist is because of that and hope someone here has a better solution for me. I just want to implement the simple interval processing like adding, multiplication, dividing, subtracting. My adviser told me to use Mathematica or Matlab, but I'm using Linux as my main operating system and even they had Linux version I'm not sure they will implement the same tools like the paid version in windows. Anyway, your advise is nice and I'm looking that book now. On Sun, Jun 8, 2014 at 6:29 PM, A Khorshidi [via Scilab / Xcos - Mailing Lists Archives] wrote: > > Hi; > > There is a Scilab toolbox for interval analysis but, the toolbox homepage > is not available at this time and return 403 access denied error. > Unfortunately! :( > http://www-sop.inria.fr/coprin/logiciels/Int4Sci/ > > Could anyone get an alternative link to download this toolbox? > > BTW, I refer you to the book "Introduction to interval analysis" written > by Ramon E. Moore et al. This book has been written based on the use of > Matlab's INtLAB toolbox. > If one has good knowledge regarding interval analysis as well as Scilab > programming, it wouldn't be hard to write her/his own Scilab functions to > implement interval arithmetic. > > Regards, > Mehran > _ > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030698.html > To unsubscribe from Interval arithmetic, click here > > . > NAML > > -- Regards, Kristophorus ----- Best regards, Hadiono -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030699.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From akhorshidi at live.com Sun Jun 8 15:26:00 2014 From: akhorshidi at live.com (A Khorshidi) Date: Sun, 8 Jun 2014 06:26:00 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic In-Reply-To: References: <1402119707393-4030693.post@n3.nabble.com> <1402226968732-4030698.post@n3.nabble.com> Message-ID: <1402233960929-4030700.post@n3.nabble.com> So, alternatively, you can use the following libraries for interval arithmetic in C++: PROFIL/BIAS: http://www.ti3.tuhh.de/keil/profil/index_e.html FILIB++: http://www2.math.uni-wuppertal.de/~xsc/software/filib.html _ -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030700.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From cfuttrup at gmail.com Sun Jun 8 16:46:13 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sun, 08 Jun 2014 16:46:13 +0200 Subject: [Scilab-users] Interval arithmetic In-Reply-To: <1402119707393-4030693.post@n3.nabble.com> References: <1402119707393-4030693.post@n3.nabble.com> Message-ID: <53947735.8070709@gmail.com> Hello Hadiono, et al. Gareth Hargreaves 50 pages on Interval Analysis in Matlab (December 2002) is available here: http://www.ti3.tuhh.de/rump/intlab/narep416.pdf Besides I received the latest (never published) Int4Sci package from the original author (2 years ago). Int4ScibinWin32.zip is about 1.8 MB. I can send in private mail, who/where shall I send it to? I never had the package working for me. I think it links to a C-library (and I don't have a C-compiler on my machine). The package includes htm documentation, examples, test script, etc. ... many of the files are from 2008. /Claus On 07-Jun-14 07:41, hadiono wrote: > Hello friends.. > > I'm new to Scilab and doesn't have any experience with Matlab also. Still > read and learn the Scilab tutorial from Scilab website and others. > I'm doing a research based on interval arithmetic > (http://en.wikipedia.org/wiki/Interval_arithmetic) and want to do the > interval arithmetic calculation using scilab. > Searching with Google and found there is a toolbox/library or similar things > for interval arithmetic but some how the development is stop. > Can somebody here help me for doing interval arithmetic calculation with > scilab? > Any link to websites or else are welcome. > > Regards, > Hadiono. > > > > -- > View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693.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 mail2kristophorus at gmail.com Sun Jun 8 17:06:45 2014 From: mail2kristophorus at gmail.com (hadiono) Date: Sun, 8 Jun 2014 08:06:45 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic In-Reply-To: <1402233960929-4030700.post@n3.nabble.com> References: <1402119707393-4030693.post@n3.nabble.com> <1402226968732-4030698.post@n3.nabble.com> <1402233960929-4030700.post@n3.nabble.com> Message-ID: Okay... thank you, let me see first if it can works or not.. On Sun, Jun 8, 2014 at 8:26 PM, A Khorshidi [via Scilab / Xcos - Mailing Lists Archives] wrote: > > So, alternatively, you can use the following libraries for interval > arithmetic in C++: > PROFIL/BIAS: http://www.ti3.tuhh.de/keil/profil/index_e.html > FILIB++: http://www2.math.uni-wuppertal.de/~xsc/software/filib.html > > _ > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030700.html > To unsubscribe from Interval arithmetic, click here > > . > NAML > > -- Regards, Kristophorus ----- Best regards, Hadiono -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030702.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From mail2kristophorus at gmail.com Sun Jun 8 17:09:57 2014 From: mail2kristophorus at gmail.com (hadiono) Date: Sun, 8 Jun 2014 08:09:57 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic In-Reply-To: <53947735.8070709@gmail.com> References: <1402119707393-4030693.post@n3.nabble.com> <53947735.8070709@gmail.com> Message-ID: Thank you my friend, I download the files from your web site. Let me check first and hope it will works... On Sun, Jun 8, 2014 at 9:46 PM, Claus Futtrup [via Scilab / Xcos - Mailing Lists Archives] wrote: > Hello Hadiono, et al. > > Gareth Hargreaves 50 pages on Interval Analysis in Matlab (December > 2002) is available here: > > http://www.ti3.tuhh.de/rump/intlab/narep416.pdf > > Besides I received the latest (never published) Int4Sci package from the > original author (2 years ago). > > Int4ScibinWin32.zip is about 1.8 MB. I can send in private mail, > who/where shall I send it to? > > I never had the package working for me. I think it links to a C-library > (and I don't have a C-compiler on my machine). > > The package includes htm documentation, examples, test script, etc. ... > many of the files are from 2008. > > /Claus > > On 07-Jun-14 07:41, hadiono wrote: > > > Hello friends.. > > > > I'm new to Scilab and doesn't have any experience with Matlab also. > Still > > read and learn the Scilab tutorial from Scilab website and others. > > I'm doing a research based on interval arithmetic > > (http://en.wikipedia.org/wiki/Interval_arithmetic) and want to do the > > interval arithmetic calculation using scilab. > > Searching with Google and found there is a toolbox/library or similar > things > > for interval arithmetic but some how the development is stop. > > Can somebody here help me for doing interval arithmetic calculation with > > scilab? > > Any link to websites or else are welcome. > > > > Regards, > > Hadiono. > > > > > > > > -- > > View this message in context: > http://mailinglists.scilab.org/Interval-arithmetic-tp4030693.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 > > > > > --- > This email is free from viruses and malware because avast! Antivirus > protection is active. > http://www.avast.com > > _______________________________________________ > 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/Interval-arithmetic-tp4030693p4030701.html > To unsubscribe from Interval arithmetic, click here > > . > NAML > > -- Regards, Kristophorus ----- Best regards, Hadiono -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030703.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Sun Jun 8 22:53:06 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sun, 08 Jun 2014 22:53:06 +0200 Subject: [Scilab-users] Interval arithmetic In-Reply-To: References: <1402119707393-4030693.post@n3.nabble.com> <53947735.8070709@gmail.com> Message-ID: <5394CD32.1050909@gmail.com> Hi there The person whom I received the Int4Sci package from is David Daney. I found it here: http://hal.inria.fr/hal-00987461 You might still be able to reach him at the INRIA address at the bottom of above link. /Claus On 08-Jun-14 17:09, hadiono wrote: > Thank you my friend, I download the files from your web site. Let me > check first and hope it will works... > > > On Sun, Jun 8, 2014 at 9:46 PM, Claus Futtrup [via Scilab / Xcos - > Mailing Lists Archives] <[hidden email] > > wrote: > > Hello Hadiono, et al. > > Gareth Hargreaves 50 pages on Interval Analysis in Matlab (December > 2002) is available here: > > http://www.ti3.tuhh.de/rump/intlab/narep416.pdf > > Besides I received the latest (never published) Int4Sci package > from the > original author (2 years ago). > > Int4ScibinWin32.zip is about 1.8 MB. I can send in private mail, > who/where shall I send it to? > > I never had the package working for me. I think it links to a > C-library > (and I don't have a C-compiler on my machine). > > The package includes htm documentation, examples, test script, > etc. ... > many of the files are from 2008. > > /Claus > > On 07-Jun-14 07:41, hadiono wrote: > > > Hello friends.. > > > > I'm new to Scilab and doesn't have any experience with Matlab > also. Still > > read and learn the Scilab tutorial from Scilab website and others. > > I'm doing a research based on interval arithmetic > > (http://en.wikipedia.org/wiki/Interval_arithmetic) and want to > do the > > interval arithmetic calculation using scilab. > > Searching with Google and found there is a toolbox/library or > similar things > > for interval arithmetic but some how the development is stop. > > Can somebody here help me for doing interval arithmetic > calculation with > > scilab? > > Any link to websites or else are welcome. > > > > Regards, > > Hadiono. > > > > > > > > -- > > View this message in context: > http://mailinglists.scilab.org/Interval-arithmetic-tp4030693.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 > > > > > --- > This email is free from viruses and malware because avast! > Antivirus protection is active. > http://www.avast.com > > _______________________________________________ > 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/Interval-arithmetic-tp4030693p4030701.html > > To unsubscribe from Interval arithmetic, click here. > NAML > > > > > > > -- > Regards, > Kristophorus > Best regards, > Hadiono > > ------------------------------------------------------------------------ > View this message in context: Re: Interval arithmetic > > 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 -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Sun Jun 8 23:04:11 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sun, 08 Jun 2014 23:04:11 +0200 Subject: [Scilab-users] Interval arithmetic In-Reply-To: References: <1402119707393-4030693.post@n3.nabble.com> <53947735.8070709@gmail.com> Message-ID: <5394CFCB.5050201@gmail.com> Hi there David Daney wrote this presentation in 2010: http://scan2010.ens-lyon.fr/SCAN2010_Daney.pdf /Claus On 08-Jun-14 17:09, hadiono wrote: > Thank you my friend, I download the files from your web site. Let me > check first and hope it will works... > > > On Sun, Jun 8, 2014 at 9:46 PM, Claus Futtrup [via Scilab / Xcos - > Mailing Lists Archives] <[hidden email] > > wrote: > > Hello Hadiono, et al. > > Gareth Hargreaves 50 pages on Interval Analysis in Matlab (December > 2002) is available here: > > http://www.ti3.tuhh.de/rump/intlab/narep416.pdf > > Besides I received the latest (never published) Int4Sci package > from the > original author (2 years ago). > > Int4ScibinWin32.zip is about 1.8 MB. I can send in private mail, > who/where shall I send it to? > > I never had the package working for me. I think it links to a > C-library > (and I don't have a C-compiler on my machine). > > The package includes htm documentation, examples, test script, > etc. ... > many of the files are from 2008. > > /Claus > > On 07-Jun-14 07:41, hadiono wrote: > > > Hello friends.. > > > > I'm new to Scilab and doesn't have any experience with Matlab > also. Still > > read and learn the Scilab tutorial from Scilab website and others. > > I'm doing a research based on interval arithmetic > > (http://en.wikipedia.org/wiki/Interval_arithmetic) and want to > do the > > interval arithmetic calculation using scilab. > > Searching with Google and found there is a toolbox/library or > similar things > > for interval arithmetic but some how the development is stop. > > Can somebody here help me for doing interval arithmetic > calculation with > > scilab? > > Any link to websites or else are welcome. > > > > Regards, > > Hadiono. > > > > > > > > -- > > View this message in context: > http://mailinglists.scilab.org/Interval-arithmetic-tp4030693.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 > > > > > --- > This email is free from viruses and malware because avast! > Antivirus protection is active. > http://www.avast.com > > _______________________________________________ > 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/Interval-arithmetic-tp4030693p4030701.html > > To unsubscribe from Interval arithmetic, click here. > NAML > > > > > > > -- > Regards, > Kristophorus > Best regards, > Hadiono > > ------------------------------------------------------------------------ > View this message in context: Re: Interval arithmetic > > 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 -------------- next part -------------- An HTML attachment was scrubbed... URL: From pierre-aime.agnel at scilab-enterprises.com Mon Jun 9 08:40:09 2014 From: pierre-aime.agnel at scilab-enterprises.com (=?UTF-8?B?UGllcnJlLUFpbcOpIEFnbmVs?=) Date: Mon, 09 Jun 2014 08:40:09 +0200 Subject: [Scilab-users] Get random number from a matrix In-Reply-To: <1402226417419-4030697.post@n3.nabble.com> References: <1402099593613-4030691.post@n3.nabble.com> <1402101034175-4030692.post@n3.nabble.com> <5392F813.2010104@free.fr> <1402226417419-4030697.post@n3.nabble.com> Message-ID: <539556C9.7090207@scilab-enterprises.com> Hi, The functions already exists for different kind of sampling (based on grand) https://help.scilab.org/docs/5.5.0/en_US/sample.html https://help.scilab.org/docs/5.5.0/en_US/samplef.html https://help.scilab.org/docs/5.5.0/en_US/samwr.html They also sample from rows/columns Regards On 06/08/2014 01:20 PM, A Khorshidi wrote: > Samuel GOUGEON wrote >> Le 07/06/2014 02:30, Chuox a ?crit : >>> i found this solution, but i'm guessing there must be a more easy and >>> efficient solution. >> not really, may be just a little shorter: >> >> Matrix = [25 ,32 ,40, 50] ; >> random = Matrix(grand(1,1,"uin",1,length(Matrix))) >> >> If Matrix is rather of text, length(Matrix) must be replaced with >> size(Matrix,"*") >> >> Regards > Hi all; > > Samuel, I think it will be handy if we have a function which can extract > randomly one, two or more entities of a matrix as well as can also be used > to return a random row/column of a matrix. > Attached you can see a preliminary version of that function. > > Regards, > Mehran > > random_extraction.sce > > > _ > > > > > -- > View this message in context: http://mailinglists.scilab.org/Get-random-number-from-a-matrix-tp4030691p4030697.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 -- Pierre-Aim? Agnel Development Team ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.68 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From mail2kristophorus at gmail.com Mon Jun 9 10:37:28 2014 From: mail2kristophorus at gmail.com (hadiono) Date: Mon, 9 Jun 2014 01:37:28 -0700 (PDT) Subject: [Scilab-users] Interval arithmetic In-Reply-To: <5394CFCB.5050201@gmail.com> References: <1402119707393-4030693.post@n3.nabble.com> <53947735.8070709@gmail.com> <5394CFCB.5050201@gmail.com> Message-ID: yes... I got all what you suggest me. I'm trying to contact David but not sure he can solve my problem, Thank you.. On Mon, Jun 9, 2014 at 4:04 AM, Claus Futtrup [via Scilab / Xcos - Mailing Lists Archives] wrote: > Hi there > > David Daney wrote this presentation in 2010: > http://scan2010.ens-lyon.fr/SCAN2010_Daney.pdf > > > /Claus > > On 08-Jun-14 17:09, hadiono wrote: > > Thank you my friend, I download the files from your web site. Let me check > first and hope it will works... > > > > On Sun, Jun 8, 2014 at 9:46 PM, Claus Futtrup [via Scilab / Xcos - Mailing > Lists Archives] <[hidden email] > > wrote: > >> Hello Hadiono, et al. >> >> Gareth Hargreaves 50 pages on Interval Analysis in Matlab (December >> 2002) is available here: >> >> http://www.ti3.tuhh.de/rump/intlab/narep416.pdf >> >> Besides I received the latest (never published) Int4Sci package from the >> original author (2 years ago). >> >> Int4ScibinWin32.zip is about 1.8 MB. I can send in private mail, >> who/where shall I send it to? >> >> I never had the package working for me. I think it links to a C-library >> (and I don't have a C-compiler on my machine). >> >> The package includes htm documentation, examples, test script, etc. ... >> many of the files are from 2008. >> >> /Claus >> >> On 07-Jun-14 07:41, hadiono wrote: >> >> > Hello friends.. >> > >> > I'm new to Scilab and doesn't have any experience with Matlab also. >> Still >> > read and learn the Scilab tutorial from Scilab website and others. >> > I'm doing a research based on interval arithmetic >> > (http://en.wikipedia.org/wiki/Interval_arithmetic) and want to do the >> > interval arithmetic calculation using scilab. >> > Searching with Google and found there is a toolbox/library or similar >> things >> > for interval arithmetic but some how the development is stop. >> > Can somebody here help me for doing interval arithmetic calculation >> with >> > scilab? >> > Any link to websites or else are welcome. >> > >> > Regards, >> > Hadiono. >> > >> > >> > >> > -- >> > View this message in context: >> http://mailinglists.scilab.org/Interval-arithmetic-tp4030693.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 >> > >> >> >> --- >> This email is free from viruses and malware because avast! Antivirus >> protection is active. >> http://www.avast.com >> >> _______________________________________________ >> 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/Interval-arithmetic-tp4030693p4030701.html >> To unsubscribe from Interval arithmetic, click here. >> NAML >> >> > > > > -- > Regards, > Kristophorus > Best regards, > Hadiono > > ------------------------------ > View this message in context: Re: Interval arithmetic > > > 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 > > > > > ------------------------------ > > > This email is free from viruses and malware because avast! Antivirus > protection is active. > > > _______________________________________________ > 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/Interval-arithmetic-tp4030693p4030705.html > To unsubscribe from Interval arithmetic, click here > > . > NAML > > -- Regards, Kristophorus ----- Best regards, Hadiono -- View this message in context: http://mailinglists.scilab.org/Interval-arithmetic-tp4030693p4030707.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From tanthiamhuat at yahoo.com Tue Jun 10 00:40:43 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Mon, 9 Jun 2014 15:40:43 -0700 (PDT) Subject: [Scilab-users] Using Scilab Genetic Algorithm to solve Simple Card Problem In-Reply-To: <1401936444905-4030679.post@n3.nabble.com> References: <1401936444905-4030679.post@n3.nabble.com> Message-ID: <1402353643850-4030708.post@n3.nabble.com> Could anyone tell me how to formulate the objective functions (to be minimized) such that 1) The sum of the first pile is as close as possible to 36. 2) And the product of all in the second pile is as close as possible to 360. See http://www.codeproject.com/Articles/16286/AI-Simple-Genetic-Algorithm-GA-to-solve-a-card-pro thank you. -- View this message in context: http://mailinglists.scilab.org/Using-Scilab-Genetic-Algorithm-to-solve-Simple-Card-Problem-tp4030679p4030708.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From lfrodrigues at gmail.com Tue Jun 10 16:57:35 2014 From: lfrodrigues at gmail.com (Luis Rodrigues) Date: Tue, 10 Jun 2014 15:57:35 +0100 Subject: [Scilab-users] interp implementation Message-ID: Hi all, I'm using a cubic spline to calibrate a sensor, the code in scilab is as follows. n=10; x=linspace(min(dists(:,1)), max(dists(:,1)), n); [y, d] = lsq_splin(dists(:,1), dists(:,2), x); ys = interp(dists(:,1), x, y, d,'natural'); I would like to store x, y, d in an eeprom and calculate the calibrated value (basically running interp) Could someone please tell me where I can find the implementation of interp so that I can implement it in Arduino? Best regards, Luis -------------- next part -------------- An HTML attachment was scrubbed... URL: From antoine.monmayrant at laas.fr Tue Jun 10 17:06:17 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Tue, 10 Jun 2014 17:06:17 +0200 Subject: [Scilab-users] interp implementation In-Reply-To: References: Message-ID: <53971EE9.2080309@laas.fr> On 06/10/2014 04:57 PM, Luis Rodrigues wrote: > Hi all, > > I'm using a cubic spline to calibrate a sensor, the code in scilab is > as follows. > > n=10; > x=linspace(min(dists(:,1)), max(dists(:,1)), n); > [y, d] = lsq_splin(dists(:,1), dists(:,2), x); > ys = interp(dists(:,1), x, y, d,'natural'); > > > I would like to store x, y, d in an eeprom and calculate the > calibrated value (basically running interp) > Could someone please tell me where I can find the implementation of > interp so that I can implement it in Arduino? I think the source code is here: http://doxygen.scilab.org/master/d7/d3b/sci__interp_8c.html Hope it helps, Antoine > > Best regards, > Luis > > > > _______________________________________________ > 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 lfrodrigues at gmail.com Tue Jun 10 22:16:17 2014 From: lfrodrigues at gmail.com (Luis Rodrigues) Date: Tue, 10 Jun 2014 21:16:17 +0100 Subject: [Scilab-users] interp implementation In-Reply-To: <53971EE9.2080309@laas.fr> References: <53971EE9.2080309@laas.fr> Message-ID: Thanks Antoine, I've also found this: http://doxygen.scilab.org/master/d1/d31/someinterp_8c.html Do you know what is the relation? Or which one is actually used? Best, Luis On 10 June 2014 16:06, Antoine Monmayrant wrote: > On 06/10/2014 04:57 PM, Luis Rodrigues wrote: > > Hi all, > > I'm using a cubic spline to calibrate a sensor, the code in scilab is as > follows. > > n=10; > x=linspace(min(dists(:,1)), max(dists(:,1)), n); > [y, d] = lsq_splin(dists(:,1), dists(:,2), x); > ys = interp(dists(:,1), x, y, d,'natural'); > > > I would like to store x, y, d in an eeprom and calculate the calibrated > value (basically running interp) > Could someone please tell me where I can find the implementation of > interp so that I can implement it in Arduino? > > > I think the source code is here: > > http://doxygen.scilab.org/master/d7/d3b/sci__interp_8c.html > > Hope it helps, > > Antoine > > > Best regards, > Luis > > > > _______________________________________________ > users mailing listusers at lists.scilab.orghttp://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 > > +++++++++++++++++++++++++++++++++++++++++++++++++++++++ > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From tanthiamhuat at yahoo.com Wed Jun 11 05:34:08 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Tue, 10 Jun 2014 20:34:08 -0700 (PDT) Subject: [Scilab-users] Generate 5 unique random integer number between 1 and 10 Message-ID: <1402457648518-4030712.post@n3.nabble.com> how do I generate 5 UNIQUE, DISTINCT random integer between 1 and 10? I try this, but it is not unique. Random_5=round(rand(1,5)*9)+1 -- View this message in context: http://mailinglists.scilab.org/Generate-5-unique-random-integer-number-between-1-and-10-tp4030712.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From tanthiamhuat at yahoo.com Wed Jun 11 07:35:55 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Tue, 10 Jun 2014 22:35:55 -0700 (PDT) Subject: [Scilab-users] Using Scilab Genetic Algorithm to solve Simple Card Problem In-Reply-To: <1401936444905-4030679.post@n3.nabble.com> References: <1401936444905-4030679.post@n3.nabble.com> Message-ID: <1402464955580-4030713.post@n3.nabble.com> managed to write a objective function for this GA problem. How do I get the final solution for Card_Stack1 and Card_Stack2? which variables are they stored in? clear;clc;clf //function y=f(x) // y = sum(x.^2) //endfunction function out = objectivefunction(AllCards) AllCards = 1:1:10; Random_10 = grand(1, "prm", (1:10)')' // random unique number 1 to 10 Cards_Stack1 = AllCards(Random_10(1:5)) Cards_Stack2 = AllCards(Random_10(6:10)) err1 = sum(Cards_Stack1) - 36; err2 = prod(Cards_Stack2) - 360; out = err1^2 + err2^2; endfunction PopSize = 100; Proba_cross = 0.7; Proba_mut = 0.1; NbGen = 10; NbCouples = 110; Log = %T; ga_params = init_param(); // Parameters to control the initial population. //ga_params = add_param(ga_params,"dimension",1); [pop_opt, fobj_pop_opt] = .. optim_ga(objectivefunction, PopSize, NbGen, Proba_mut, Proba_cross, Log, ga_params); // Display basic statistics // min, mean and max function values of the population. disp([min(fobj_pop_opt) mean(fobj_pop_opt) max(fobj_pop_opt)]) // Get the best x (i.e. the one which achieves the minimum function value) [fmin ,k] = min(fobj_pop_opt); xmin = pop_opt(k); // Get the worst x [fmax ,k] = max(fobj_pop_opt); xmax = pop_opt(k); -- View this message in context: http://mailinglists.scilab.org/Using-Scilab-Genetic-Algorithm-to-solve-Simple-Card-Problem-tp4030679p4030713.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.bignier at scilab-enterprises.com Wed Jun 11 08:41:37 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Wed, 11 Jun 2014 08:41:37 +0200 Subject: [Scilab-users] Generate 5 unique random integer number between 1 and 10 In-Reply-To: <1402457648518-4030712.post@n3.nabble.com> References: <1402457648518-4030712.post@n3.nabble.com> Message-ID: <5397FA21.4000000@scilab-enterprises.com> Hello, How about taking the five first elements of a permutation of vector 1:10 ? Random_10 = grand(1, "prm", 1:10); Random_5 = Random_10(1:5) http://help.scilab.org/docs/5.5.0/en_US/grand.html Regards, Paul On 06/11/2014 05:34 AM, tanthiamhuat wrote: > how do I generate 5 UNIQUE, DISTINCT random integer between 1 and 10? > > I try this, but it is not unique. > > Random_5=round(rand(1,5)*9)+1 > > > > -- > View this message in context: http://mailinglists.scilab.org/Generate-5-unique-random-integer-number-between-1-and-10-tp4030712.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 tanthiamhuat at yahoo.com Wed Jun 11 08:44:05 2014 From: tanthiamhuat at yahoo.com (TAN TH) Date: Tue, 10 Jun 2014 23:44:05 -0700 (PDT) Subject: [Scilab-users] Generate 5 unique random integer number between 1 and 10 In-Reply-To: <5397FA21.4000000@scilab-enterprises.com> References: <1402457648518-4030712.post@n3.nabble.com> <5397FA21.4000000@scilab-enterprises.com> Message-ID: <1402469045.58110.YahooMailNeo@web125504.mail.ne1.yahoo.com> yes, I got it already, thanks a lot. On Wednesday, 11 June 2014, 14:42, Paul Bignier wrote: Hello, How about taking the five first elements of a permutation of vector 1:10 ? Random_10 = grand(1, "prm", 1:10); Random_5? = Random_10(1:5) http://help.scilab.org/docs/5.5.0/en_US/grand.html Regards, Paul On 06/11/2014 05:34 AM, tanthiamhuat wrote: > how do I generate 5 UNIQUE, DISTINCT random integer between 1 and 10? > > I try this, but it is not unique. > > Random_5=round(rand(1,5)*9)+1 > > > > -- > View this message in context: http://mailinglists.scilab.org/Generate-5-unique-random-integer-number-between-1-and-10-tp4030712.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 _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From raghu.char at gmail.com Wed Jun 11 08:17:26 2014 From: raghu.char at gmail.com (Raghu Char) Date: Wed, 11 Jun 2014 11:47:26 +0530 Subject: [Scilab-users] bode plot Message-ID: List, I'd like to increase the font size of the bode plot( jpeg attached), i have a function to do this but as there are 2 figures it works only on 1, could not find any help on this, how do i do this? code snippet: function y=set_gr_txt(fnt_sz) // sets textize in plots to larger (4) // [out, inp] = argn(0); // //if inp < 1 then, x = 1, end //set default font size to 1 a = gca(); if ~exists("fnt_sz","local") then // default size fnt_sz = 4; end a.font_size = 4; // axis font size x_label=a.x_label; x_label.font_size=fnt_sz; y_label=a.y_label; y_label.font_size=fnt_sz; t = a.title; t.font_size=fnt_sz; //disp(is_handle_valid(a),typeof(h),'typeof handle'); endfunction thanx in advance. -- Regards raghu -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: bode_plt.jpg Type: image/jpeg Size: 102638 bytes Desc: not available URL: From paul.bignier at scilab-enterprises.com Wed Jun 11 10:02:13 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Wed, 11 Jun 2014 10:02:13 +0200 Subject: [Scilab-users] bode plot In-Reply-To: References: Message-ID: <53980D05.2000909@scilab-enterprises.com> Hi Raghu, You need to retrieve the axes of the second plot thanks to the children of gcf(). Here is how I did it, based on your snippet: function set_gr_txt(fnt_sz) // sets textize in plots to larger (4) // [out, inp] = argn(0); // //if inp < 1 then, x = 1, end //set default font size to 1 f = gcf(); a1 = f.children(1); if ~exists("fnt_sz","local") then // default size fnt_sz = 4; end a1.font_size = 4; // axis font size x_label=a1.x_label; x_label.font_size=fnt_sz; y_label=a1.y_label; y_label.font_size=fnt_sz; t = a1.title; t.font_size=fnt_sz; //disp(is_handle_valid(a1),typeof(h),'typeof handle'); a2 = f.children(2); if ~exists("fnt_sz","local") then // default size fnt_sz = 4; end a2.font_size = 4; // axis font size x_label=a2.x_label; x_label.font_size=fnt_sz; y_label=a2.y_label; y_label.font_size=fnt_sz; t = a2.title; t.font_size=fnt_sz; //disp(is_handle_valid(a2),typeof(h),'typeof handle'); endfunction // Execution: bode() set_gr_txt(4) Regards, Paul On 06/11/2014 08:17 AM, Raghu Char wrote: > List, > > I'd like to increase the font size of the bode plot( jpeg attached), i > have a function to do this > but as there are 2 figures it works only on 1, could not find any > help on this, how do i do this? > > > code snippet: > function y=set_gr_txt(fnt_sz) > // sets textize in plots to larger (4) > // [out, inp] = argn(0); > // //if inp< 1 then, x = 1, end //set default font size to 1 > > a = gca(); > if ~exists("fnt_sz","local") then // default size > fnt_sz = 4; > end > > a.font_size = 4; // axis font size > x_label=a.x_label; > x_label.font_size=fnt_sz; > y_label=a.y_label; > y_label.font_size=fnt_sz; > t = a.title; > t.font_size=fnt_sz; > //disp(is_handle_valid(a),typeof(h),'typeof handle'); > > > endfunction > > thanx in advance. > > > -- > Regards raghu > > > _______________________________________________ > 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 julian.estevez at ehu.es Wed Jun 11 16:29:39 2014 From: julian.estevez at ehu.es (Jeibros) Date: Wed, 11 Jun 2014 07:29:39 -0700 (PDT) Subject: [Scilab-users] Draw a figure in specific points Message-ID: <1402496979188-4030718.post@n3.nabble.com> Hi all, I would like to plot a cross made with two lines in specific points. Those crosses must have equal length. I am using right now param3d command to try to draw the lines if the crosses. anybody knows about the most efficient way to programme this? Thanks -- View this message in context: http://mailinglists.scilab.org/Draw-a-figure-in-specific-points-tp4030718.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Wed Jun 11 16:53:18 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 11 Jun 2014 16:53:18 +0200 Subject: [Scilab-users] Draw a figure in specific points In-Reply-To: <1402496979188-4030718.post@n3.nabble.com> References: <1402496979188-4030718.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32B5FB62E@301EX00100.sidel.com> Hello, > De Jeibros > Envoy? : mercredi 11 juin 2014 16:30 > > I would like to plot a cross made with two lines in specific points. > [...] > anybody knows about the most efficient way to programme this? Not sure what you want exactly, but could plot(xvalue, yvalue, "+") do the trick? Or maybe xpoly([xvalue-hw ; xvalue+hw], [yvalue ; yvalue]) xpoly([xvalue ; xvalue], [yvalue-hw ; yvalue+hw]) beware that the half width hw might be different according to the scales. HTH -- 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 tanthiamhuat at yahoo.com Thu Jun 12 03:26:21 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Wed, 11 Jun 2014 18:26:21 -0700 (PDT) Subject: [Scilab-users] Genetic Algorithm (GA) in Scilab: Retrieve the final solution Message-ID: <1402536381784-4030721.post@n3.nabble.com> Hi Scilab experts in GA, I would like to find out how can I retrieve the final GA solutions in Scilab variables. Firstly, I am trying to solve a Simple Card Problem, as define here: http://www.codeproject.com/Articles/16286/AI-Simple-Genetic-Algorithm-GA-to-solve-a-card-pro My Scilab code is as below: But I could not retrieve its final solution for Cards_Stack1 and Cards_Stack2. How should I change my code, such that those variables (Cards_Stack1, Cards_Stack2) are available? clear;clc;clf function out = objectivefunction(AllCards) AllCards = 1:1:10; Cards_Stack1 = samwr(5,1,AllCards)' Cards_Stack2 = setdiff(AllCards,Cards_Stack1) err1 = sum(Cards_Stack1) - 36; err2 = prod(Cards_Stack2) - 360; out = err1^2 + err2^2; endfunction PopSize = 100; Proba_cross = 0.7; Proba_mut = 0.1; NbGen = 10; NbCouples = 110; Log = %T; ga_params = init_param(); // Parameters to control the initial population. //ga_params = add_param(ga_params,"dimension",1); [pop_opt, fobj_pop_opt] = .. optim_ga(objectivefunction, PopSize, NbGen, Proba_mut, Proba_cross, Log, ga_params); // Display basic statistics // min, mean and max function values of the population. disp([min(fobj_pop_opt) mean(fobj_pop_opt) max(fobj_pop_opt)]) // Get the best x (i.e. the one which achieves the minimum function value) [fmin ,k] = min(fobj_pop_opt); xmin = pop_opt(k); // Get the worst x [fmax ,k] = max(fobj_pop_opt); xmax = pop_opt(k); -- View this message in context: http://mailinglists.scilab.org/Genetic-Algorithm-GA-in-Scilab-Retrieve-the-final-solution-tp4030721.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From julian.estevez at ehu.es Wed Jun 11 18:19:47 2014 From: julian.estevez at ehu.es (Jeibros) Date: Wed, 11 Jun 2014 09:19:47 -0700 (PDT) Subject: [Scilab-users] Draw a figure in specific points In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F32B5FB62E@301EX00100.sidel.com> References: <1402496979188-4030718.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32B5FB62E@301EX00100.sidel.com> Message-ID: <1402503587289-4030720.post@n3.nabble.com> not exactly. Let's assume I have a list of x, y, z coordinates. In each of these coordinates, I would like to plot a specific figure. Imagine I want to plot a circle centered in each one of those coordinates with a 2 units radius and oriented somehow. -- View this message in context: http://mailinglists.scilab.org/Draw-a-figure-in-specific-points-tp4030718p4030720.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Thu Jun 12 09:31:56 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Thu, 12 Jun 2014 09:31:56 +0200 Subject: [Scilab-users] Draw a figure in specific points In-Reply-To: <1402503587289-4030720.post@n3.nabble.com> References: <1402496979188-4030718.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32B5FB62E@301EX00100.sidel.com> <1402503587289-4030720.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F32B5FBA21@301EX00100.sidel.com> Hello, > De Jeibros > Envoy? : mercredi 11 juin 2014 18:20 > > Let's assume I have a list of x, y, z coordinates. Sorry. I should have guessed it was 3D when I read "param3d", but I answered too quickly. > In each of these coordinates, I would like to plot a specific figure. So you have a figure, and you want to translate it to a specific (x, y, z) position ? Let's go back to the crosses. You can have : param3d(xvalue, yvalue, zvalue) h = gce() h.mark_style=1 The cross will be exactly the same whatever the point of view. Now if you want to draw anything else than a mark, the shape will change with the point of view. So, you can define your figure in a function, and provide the position as a parameter, something like function []=myfigure(x, y, z) param3d([x-1 x+1], [y y], [z z]) param3d([x x], [y-1 y+1], [z z]) param3d([x x], [y y], [z-1 z+1]) endfunction myfigure(0, 0, 0) myfigure(2, 2, 2) You can also try the move() function, but after a few tests, It semms the handle only refers to one single polyline. -- 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 bhark at kkumail.com Thu Jun 12 10:11:09 2014 From: bhark at kkumail.com (bhark) Date: Thu, 12 Jun 2014 01:11:09 -0700 (PDT) Subject: [Scilab-users] transform image contour to coordinate system Message-ID: <1402560669064-4030724.post@n3.nabble.com> I have the contour image of jellyfish, so i'm trying to quantify bell surface area. How can i transform contour of bell into points on coordinate system? -- View this message in context: http://mailinglists.scilab.org/transform-image-contour-to-coordinate-system-tp4030724.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sakiss123second at gmail.com Thu Jun 12 18:36:14 2014 From: sakiss123second at gmail.com (Stephen Kiss) Date: Thu, 12 Jun 2014 12:36:14 -0400 Subject: [Scilab-users] please email me my password - I tried to reregister with my other email sak123second@gmail.com but that password is not working either. Thanks. In-Reply-To: <53991255.7050800@att.net> References: <53991255.7050800@att.net> Message-ID: I made a typo on my email address that I want to register for: sakiss123second at gmail.com Is the correct email address. Thanks. On Wed, Jun 11, 2014 at 10:37 PM, Stephen Kiss wrote: > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From bhark at kkumail.com Fri Jun 13 03:54:21 2014 From: bhark at kkumail.com (bhark) Date: Thu, 12 Jun 2014 18:54:21 -0700 (PDT) Subject: [Scilab-users] transform image contour to coordinate system In-Reply-To: <1402560669064-4030724.post@n3.nabble.com> References: <1402560669064-4030724.post@n3.nabble.com> Message-ID: <1402624461996-4030736.post@n3.nabble.com> -- View this message in context: http://mailinglists.scilab.org/transform-image-contour-to-coordinate-system-tp4030724p4030736.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From aweeks at hidglobal.com Fri Jun 13 12:07:47 2014 From: aweeks at hidglobal.com (aweeks at hidglobal.com) Date: Fri, 13 Jun 2014 11:07:47 +0100 Subject: [Scilab-users] SciLab: possible bug in graphics. Message-ID: Hi Everyone, One of my programs has a 'terminal' where incoming text is shown. I've done this by writing text to a graphics figure. The text is stored as a global vector ('displaytext') consisting of lines (strings) which are written to the figure by means of the handle to the text entity (texthdl.text = displaytext). The code for this is shown below. The display is updated around four times per second. Notice also an eventhandler in the figure, which allows the looping which generates the text to be stopped. All too often SciLab crashes and the error messages which appear in the console are shown below. I've tried to work around this by closing and re-opening the figure after #disptot lines have been written - this is the purpose of variables 'disptot' and 'restart'. However, #disptot can be quite low, 40 say, and SciLab will still crash. It looks to me like memory is not being freed after use and so runs out. I haven't been able to work out whether it's the number of lines or the number of characters which determines how soon it crashes. Any suggestions ? This looks like a bug to me. Are there any fixes ? Alternatively, if I knew what made it crash I might be able to make the software adaptive so it only re-opens the figure as infrequently as necessary - re-opening the figure is slightly annoying, having to continually close and re-launch SciLab is a real pain. I'm using SciLab 5.4.0 in Windows XP SP3 32 bit. Many thanks in anticipation, Adrian. #ldisplay = 50 #disptot = 100 function StartDisplay(restart) global figrhdl; global texthdl; global displaytext; global disptot; global stophdl; disptot = 0; figrhdl = scf(1); figrhdl.figure_position = [#fig1_xTLC, #fig1_yTLC]; figrhdl.figure_size = [#fig1_width, #fig1_height]; figrhdl.figure_name = "Message window."; xstring(-0.1,1.1,"LEFT CLICK HERE TO STOP LOOPING"); stophdl = gce(); stophdl.font_style = 0; stophdl.font_size = 2; stophdl.box = 'on'; seteventhandler('StopLoop'); figrhdl.event_handler = 'StopLoop'; figrhdl.event_handler_enable = 'on'; if ~restart then // if restart = %t the displaytext is not cleared displaytext(1: #ldisplay) = " "; end xstring(-0.1, -0.1, displaytext); // this writes the text initially texthdl = gce(); texthdl.font_style = 0; texthdl.font_size = 2; endfunction function DisplayText(textlines, clrfirst) global figrhdl; global texthdl; global displaytext; global disptot; if clrfirst then displaytext(1: #ldisplay) = " "; end // if clrfirst = %t then displaytext is cleared nlines = size(textlines, 'r'); // these lines, from here offset = #ldisplay - nlines; if nlines < #ldisplay then i = 1: offset; displaytext(i) = displaytext(i + nlines); i = 1: nlines; displaytext(offset + i) = textlines(i) + ' '; else i = 1: #ldisplay; displaytext(i) = textlines(i - offset) + ' '; end // to here update displaytext by replacing lines with the new 'textlines' if disptot >= #disptot then // if disptot lines have been written the figure is closed and re-opened close(figrhdl); StartDisplay(%t); end texthdl.text = displaytext; // this writes the text to the figure disptot = disptot + 1; endfunction Exception in thread "Thread-1782" javax.media.opengl.GLException: java.lang.reflect.InvocationTargetException at javax.media.opengl.awt.GLJPanel.display(GLJPanel.java:265) at org.scilab.forge.scirenderer.implementation.jogl.JoGLCanvas $CanvasAnimator.run(Unknown Source) at java.lang.Thread.run(Unknown Source) Caused by: java.lang.reflect.InvocationTargetException at java.awt.EventQueue.invokeAndWait(Unknown Source) at javax.media.opengl.awt.GLJPanel.display(GLJPanel.java:263) ... 2 more Caused by: java.lang.OutOfMemoryError: Java heap space at java.awt.image.DataBufferInt.(Unknown Source) at java.awt.image.Raster.createPackedRaster(Unknown Source) at java.awt.image.DirectColorModel.createCompatibleWritableRaster (Unknown Source) at java.awt.image.BufferedImage.(Unknown Source) at org.scilab.forge.scirenderer.texture.TextureBufferedImage. (Unknown Source) at org.scilab.forge.scirenderer.texture.BufferedImageTextureDrawingTools. (Unknown Source) at org.scilab.forge.scirenderer.texture.DrawnTextureDataProvider.reDraw (Unknown Source) at org.scilab.forge.scirenderer.texture.DrawnTextureDataProvider.getImage (Unknown Source) at org.scilab.forge.scirenderer.texture.DrawnTextureDataProvider.getSubImage (Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.texture.JoGLTextureManager $JoGLTexture.checkData(Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.texture.JoGLTextureManager $JoGLTexture.preDraw(Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.texture.JoGLTextureManager.draw (Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.JoGLDrawingTools.draw (Unknown Source) at org.scilab.modules.renderer.JoGLView.text.TextManager.draw(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.visit(Unknown Source) at org.scilab.modules.graphic_objects.textObject.Text.accept(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.askAcceptVisitor(Unknown Source) at org.scilab.modules.renderer.JoGLView.axes.AxesDrawer.draw(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.visit(Unknown Source) at org.scilab.modules.graphic_objects.axes.Axes.accept(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.askAcceptVisitor(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.visit(Unknown Source) at org.scilab.modules.graphic_objects.figure.Figure.accept(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.draw(Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.JoGLCanvas.display(Unknown Source) at jogamp.opengl.GLDrawableHelper.displayImpl (GLDrawableHelper.java:189) at jogamp.opengl.GLDrawableHelper.display(GLDrawableHelper.java:177) at javax.media.opengl.awt.GLJPanel$Updater.display(GLJPanel.java:655) at jogamp.opengl.GLDrawableHelper.displayImpl (GLDrawableHelper.java:189) at jogamp.opengl.GLDrawableHelper.display(GLDrawableHelper.java:177) at jogamp.opengl.GLPbufferImpl$DisplayAction.run (GLPbufferImpl.java:300) at jogamp.opengl.GLDrawableHelper.invokeGLImpl (GLDrawableHelper.java:425) Exception in thread "AWT-EventQueue-0" java.lang.OutOfMemoryError: Java heap space at java.awt.image.DataBufferInt.(Unknown Source) at java.awt.image.Raster.createPackedRaster(Unknown Source) at java.awt.image.DirectColorModel.createCompatibleWritableRaster (Unknown Source) at java.awt.image.BufferedImage.(Unknown Source) at org.scilab.forge.scirenderer.texture.TextureBufferedImage. (Unknown Source) at org.scilab.forge.scirenderer.texture.BufferedImageTextureDrawingTools. (Unknown Source) at org.scilab.forge.scirenderer.texture.DrawnTextureDataProvider.reDraw (Unknown Source) at org.scilab.forge.scirenderer.texture.DrawnTextureDataProvider.getImage (Unknown Source) at org.scilab.forge.scirenderer.texture.DrawnTextureDataProvider.getSubImage (Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.texture.JoGLTextureManager $JoGLTexture.checkData(Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.texture.JoGLTextureManager $JoGLTexture.preDraw(Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.texture.JoGLTextureManager.draw (Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.JoGLDrawingTools.draw (Unknown Source) at org.scilab.modules.renderer.JoGLView.text.TextManager.draw(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.visit(Unknown Source) at org.scilab.modules.graphic_objects.textObject.Text.accept(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.askAcceptVisitor(Unknown Source) at org.scilab.modules.renderer.JoGLView.axes.AxesDrawer.draw(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.visit(Unknown Source) at org.scilab.modules.graphic_objects.axes.Axes.accept(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.askAcceptVisitor(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.visit(Unknown Source) at org.scilab.modules.graphic_objects.figure.Figure.accept(Unknown Source) at org.scilab.modules.renderer.JoGLView.DrawerVisitor.draw(Unknown Source) at org.scilab.forge.scirenderer.implementation.jogl.JoGLCanvas.display(Unknown Source) at jogamp.opengl.GLDrawableHelper.displayImpl (GLDrawableHelper.java:189) at jogamp.opengl.GLDrawableHelper.display(GLDrawableHelper.java:177) at javax.media.opengl.awt.GLJPanel$Updater.display(GLJPanel.java:655) at jogamp.opengl.GLDrawableHelper.displayImpl (GLDrawableHelper.java:189) at jogamp.opengl.GLDrawableHelper.display(GLDrawableHelper.java:177) at jogamp.opengl.GLPbufferImpl$DisplayAction.run (GLPbufferImpl.java:300) at jogamp.opengl.GLDrawableHelper.invokeGLImpl (GLDrawableHelper.java:425) 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: 12696196.gif Type: image/gif Size: 4559 bytes Desc: not available URL: From julian.estevez at ehu.es Fri Jun 13 11:53:00 2014 From: julian.estevez at ehu.es (Jeibros) Date: Fri, 13 Jun 2014 02:53:00 -0700 (PDT) Subject: [Scilab-users] Draw a figure in specific points In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F32B5FBA21@301EX00100.sidel.com> References: <1402496979188-4030718.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32B5FB62E@301EX00100.sidel.com> <1402503587289-4030720.post@n3.nabble.com> <3B5FFC67498DFF49AE7271A584867D16F32B5FBA21@301EX00100.sidel.com> Message-ID: <1402653180976-4030737.post@n3.nabble.com> Thanks, that got me things a little bit more clearer. -- View this message in context: http://mailinglists.scilab.org/Draw-a-figure-in-specific-points-tp4030718p4030737.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From saber.s.gamei at chmail.ir Fri Jun 13 15:26:36 2014 From: saber.s.gamei at chmail.ir (saber_s) Date: Fri, 13 Jun 2014 06:26:36 -0700 (PDT) Subject: [Scilab-users] optim (fminbnd) and alpha_star Message-ID: <1402665996959-4030739.post@n3.nabble.com> salam hi! 1-f=x^2 - 2*x; f & x optim in scilab? (fminbnd?) plz -- View this message in context: http://mailinglists.scilab.org/optim-fminbnd-and-alpha-star-tp4030739.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sdr at durietz.se Sat Jun 14 13:35:55 2014 From: sdr at durietz.se (Stefan Du Rietz) Date: Sat, 14 Jun 2014 13:35:55 +0200 Subject: [Scilab-users] Bug in csim in Scilab 5.5.0? Message-ID: <539C339B.4070500@durietz.se> Hello, I have a function using csim that worked OK in 5.4.1 but not in 5.5.0. However, if I use the csim.sci from 5.4.1 it also works in 5.5.0. I stop the simulation, read the states, insert an impulse by adding to the first state, and continue with the the states as initial states (x0) in a new run. So it seems like the new csim does something different with the initial states? Regards Stefan From amoeller1371 at gmail.com Fri Jun 13 18:47:24 2014 From: amoeller1371 at gmail.com (A. Moeller) Date: Fri, 13 Jun 2014 18:47:24 +0200 Subject: [Scilab-users] Problems with colors when exporting figures in scilab Message-ID: Hi all, when I try to export figures in scilab 5.5.0 under windows 7, I habe the problem that the colors in the graphics windows and in the exported .pdf-file are not the same. Is there some option to tell scilab to use the same colors in the exported .pdf-files as in the graphics windows? (In the graphics window any colors palette seems to work fine.) Thanks in advance A. Moeller From Christophe.Dang at sidel.com Mon Jun 16 10:12:10 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 16 Jun 2014 10:12:10 +0200 Subject: [Scilab-users] Problems with colors when exporting figures in scilab In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F407F20CCD@301EX00100.sidel.com> Hello, > De : A. Moeller > Envoy? : vendredi 13 juin 2014 18:47 > > the colors in the graphics windows and in the exported .pdf-file are not the same. Maybe it's the same as the bug 11610 comment #5? http://bugzilla.scilab.org/show_bug.cgi?id=11610#c5 Samule, you told on comment #8 that you posted this as a separate bug. Could we have the reference for it? 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 calixte.denizet at scilab-enterprises.com Mon Jun 16 11:10:43 2014 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Mon, 16 Jun 2014 11:10:43 +0200 Subject: [Scilab-users] Problems with colors when exporting figures in scilab In-Reply-To: References: Message-ID: <539EB493.4070103@scilab-enterprises.com> Hi, Could you provide a test case ? a screenshot ? the generated pdf ? It is difficult to help you without any informations. Best regards Calixte On 13/06/2014 18:47, A. Moeller wrote: > Hi all, > > when I try to export figures in scilab 5.5.0 under windows 7, I habe > the problem that the colors in the graphics windows and in the > exported .pdf-file are not the same. > > Is there some option to tell scilab to use the same colors in the > exported .pdf-files as in the graphics windows? (In the graphics > window any colors palette seems to work fine.) > > Thanks in advance > A. Moeller > _______________________________________________ > 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 From ashishjaen at gmail.com Thu Jun 19 09:44:50 2014 From: ashishjaen at gmail.com (Ashish Jaen) Date: Thu, 19 Jun 2014 00:44:50 -0700 Subject: [Scilab-users] Large scale linear optimization Message-ID: Hi, I am new to this mailing list, so please pardon me if this question has been asked before. Has anyone used Scilab for large scale linear optimization ? The problem size I am looking at is about 2000 variables and 5000 constraints. Will be great if you can also share some performance numbers. I have used Matlab in the past but no longer have the license as I am not in university any more :( Regards. -------------- next part -------------- An HTML attachment was scrubbed... URL: From bmbouter at gmail.com Fri Jun 20 12:52:53 2014 From: bmbouter at gmail.com (Brian Bouterse) Date: Fri, 20 Jun 2014 06:52:53 -0400 Subject: [Scilab-users] GROCER VAR-HMM Example Output Different from Chapter 23 In-Reply-To: References: Message-ID: Hi Eric, I want to reach out and see what your thoughts are on putting a copy of grocer on github? There is already a grocer on github that is a different software project so github.com/grocer/grocer is already taken. Maybe it could be housed at github.com//grocer for now. I could also post it on github ( github.com/bmbouter ), but having a canonical version from the author would be better. I would be willing to submit some documentation patches through github. Thanks again for the great software and the help; it's working for us! Best, Brian On Tue, Jun 3, 2014 at 4:40 PM, Eric Dubois wrote: > Hello Brian. > > Thank you for your feedback. I will try to make the distinction between > 'cte' and 'all' more precise in the manual. > > I will also think about your suggestion to put the code on a hub. > > Regards. > > ?ric. > > > 2014-06-03 12:21 GMT+02:00 Brian Bouterse : > > Hi Eric, >> >> Thanks for the quick and helpful reply! I adjusted the first parameter >> from 'cte' to 'all' and that produces the correct form of output, and the >> values look close to the expected values. This is great! >> >> I've read the help documents several times, but I didn't realize the >> impact of adjusting the value of the first parameter. The doc example on >> the ms_var help page describes the parameter, but the only example on that >> page uses 'cte'. There is no example with 'all' on the ms_var() help page. >> Also including some sample output may also be good there. I would like to >> contribute an example like this to the docs, manual, or help files. I think >> that would have caused me to solve my own problems instead of writing long >> e-mails. >> >> I also tried using a larger data file 10,000 observations in a column >> format, and then impexc2bd transformed the data file correctly. That is >> also solved thanks to your suggestion. >> >> I know the code is available for download, but would you ever consider >> putting a copy of the GROCER code in a repository on github.com? If it >> were in a place like that, I, and others, would have an easier time of >> contributing to the codebase instead of e-mailing patches. What do you >> think about an idea like this? Just a friendly suggestion based on what has >> worked well for me in the past. I can help put it up there, and configure >> it if that is helpful. >> >> Thanks again for making this great software, and helping me use it! >> >> Best, >> Brian >> >> >> >> >> >> On Thu, May 29, 2014 at 1:41 PM, Eric Dubois >> wrote: >> >>> Hello Brian >>> >>> You did not recover the intial parameters because you did not estimate >>> the same model as the one you generated. >>> >>> This is because you used in ms_var the option 'cte' which means that >>> only the constant is allowed to switch: and this is why the AR coeffcients >>> are only given for 'all regimes'. I f you rather use the option 'all', then >>> all parameters are allowed to switch and the result is now consistent with >>> the parameters of the generated model, that is the true values are within >>> the confidence interval of the estimated ones. This is explained in the >>> help command (see help ms_var): may I ask you why you did not infer it from >>> the help files? >>> >>> As for the big file, I suspect it is truncated into Scilab. Can you try >>> with data in column instead of in rows? >>> >>> ?ric. >>> >>> >>> 2014-05-29 13:24 GMT+02:00 Brian Bouterse : >>> >>> Hi Eric, >>>> >>>> Thanks again for the response. What you told me to do worked. If I >>>> closed Scilab, and rerun the example from my earlier e-mail, I do not see >>>> any warnings, which is good. The output I received is not what I expected >>>> in several ways, so I'll ask some more questions here. There are a lot of >>>> different questions here, but they all relate to the example below, and the >>>> use of ms_var and GROCER. Thanks in advance for any insight into the issues >>>> I am experiencing. >>>> >>>> The example: >>>> >>>> I have a file arhmm_example_280.dat [0] that was generated and loaded >>>> from arhmm_example_280.csv [1] using the following commands: >>>> >>>> impexc2bd('arhmm_example_280.csv', ';', 'arhmm_example_280.dat') >>>> load('arhmm_example_280.dat') >>>> >>>> I then run, what I expect to be a single variable, 2-regime, >>>> >>>> nb_states=2 >>>> switch_var=2 >>>> var_opt=3 >>>> >>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>> >>>> I then receive this output: http://ur1.ca/hej98 >>>> >>>> 1) In the "coefficients" section of the output, the output doesn't seem >>>> to contain enough numbers to describe the expected number of outputs. The >>>> model has 2 regimes, and it is autoregressive with order 3, so I expect >>>> there to be 3 coefficients for EACH regime (regime 1 and regime 2). >>>> Instead, I see a single coefficient for Regime 1 (0.9313736) and a single >>>> coefficient for Regime 2 (8.2685973). Where are the three autoregressive >>>> parameters for regime 1 and the other three autoregressive parameters for >>>> regime 2? >>>> >>>> 2) There is a section in the output called "All Regimes" which is in >>>> the "Coefficients" section contains 3 coefficients. Do these correspond to >>>> the autoregressive aspect of the model [0.2258037, 0.4757391, -0.3161730). >>>> I expected to see a section like this for Regime 1, and one for Regime 2, >>>> but I did not. I'm not sure how these 3 coefficients for "All Regimes" fit >>>> in. The Markov Switching Model switches the autoregressive coefficients >>>> with each Regime, so I'm not sure how to interpret coefficients that are >>>> for "both regimes". >>>> >>>> 3) I believe each regime will have its error term estimated, and I do >>>> see that in the output. For example under Regime 2, the variance of the >>>> error term is 2.3164941. This is listed in the section "Variance-covariance >>>> matrix of residuals" which is a subsection under "Regime 2". Am I >>>> interpreting the output correctly? Where is the mean value for the regime? >>>> >>>> 4) My overall goal is to validate that I can hand a list of >>>> observations to GROCER+scilab, along with the order of the autoregressive >>>> model, and the number of regimes, and have it estimate the Markov hidden >>>> state transition probabilities, the estimated regime that is in effect at >>>> time t, and coefficients for each regime, and the error term that can be >>>> used for prediction of the t+1 observation. I don't know of a "well known" >>>> data set that allows for this type of validation, so I wrote a simple >>>> AR-HMM generator named koa marhn [2], which outputs AR-HMM data that >>>> follows a set of parameters. This is the program that generated the 280 >>>> point data set in the arhmm_example_280.csv [1] file. You can see the >>>> parameters used to generate this file here: >>>> >>>> https://github.com/bmbouter/koa_marhn/blob/master/arhmm_example.py >>>> >>>> A[i][j] is the probability transition from regime i to regime j. >>>> C[x][y] is the autoregressive coefficient for regime x used for >>>> generation of the t+1 element from the t, t-1, t-2, ..., t-y summation. >>>> Each y is a coefficient for regime x. >>>> R[q] contains mean and std_deviation values for regime q >>>> pi[m] contains the starting state probabilities that the regime is in >>>> state m at time t=0 >>>> >>>> The output from http://ur1.ca/hej98 does not correspond at all with >>>> the values of these inputs, which means something unexpected is going on! >>>> Do you have insight into why scilab doesn't produce these expected numbers? >>>> >>>> 5) One possibility is that maybe it is "working", but there isn't >>>> enough data for it to estimate the parameters correctly. I had koa marhn >>>> [2] output 10,000 observations, but I ran into trouble transforming it from >>>> a csv file to a dat file using impexc2bd. My input file is >>>> arhmm_example_10000.csv [3]. I tried to convert to .dat using: >>>> >>>> impexc2bd('arhmm_example_10000.csv', ';', 'arhmm_example_10000.dat') >>>> >>>> I receive this output, and the .dat file in NOT created >>>> >>>> !--error 10000 >>>> dates are entered neither in chronlogical order nor in reverseerse >>>> chronological order >>>> at line 87 of function read_dates called by : >>>> at line 188 of function impexc2bd called by : >>>> arhmm_example_10000.dat' >>>> >>>> How do you transform a large .csv file to .dat? Is there some problem >>>> with the data file, the way I am making it? >>>> >>>> [0]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_280.dat >>>> [1]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_280.csv >>>> [2]: https://github.com/bmbouter/koa_marhn >>>> [3]: >>>> https://s3.amazonaws.com/scilab_data_files/arhmm_example_10000.csv >>>> [4]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_1000.csv >>>> [5]: https://s3.amazonaws.com/scilab_data_files/arhmm_example_1000.dat >>>> >>>> This e-mail is way to long, but I wanted to fully recap the issues I've >>>> considered as I use ms_var(). Thanks for any help you can provide. You can >>>> also find me in the #scilab channel on the OFTC IRC servers. >>>> >>>> Thanks, >>>> Brian >>>> >>>> >>>> >>>> On Fri, May 23, 2014 at 3:48 PM, Eric Dubois >>>> wrote: >>>> >>>>> Dear Brian >>>>> >>>>> 1) For me it works fine; I suspect that you have run ms_var with your >>>>> data after having run with the manual example. I expect that if you reopen >>>>> Scialb and run your problem, starting from load('arhmm_example.dat'), >>>>> it will work. A good advice anyway is to set the bounds before each >>>>> estimation or to run: >>>>> --> bounds() >>>>> if you want to use the greatest available time span with your data >>>>> 2) I agree; indeed the example in the manual as well as yours are >>>>> univariate, but ms_var also works with multivariate series (I have run some >>>>> tests which worked well) >>>>> 3) If you want to contribute, do not hesitate to send me code (at >>>>> grocer.toolbox at gmail.com or grocer.toolbox at free.fr); add your >>>>> copyright; if you can create a help file it would still be better (I have >>>>> some tools do help doing that if you want them); and if you can add to the >>>>> manual, it would be marvellous! >>>>> If you want to imporve the docs you are also welcome; I can sned you >>>>> the OpenOffice files you need if you find it suitable >>>>> >>>>> ?ric. >>>>> >>>>> >>>>> >>>>> 2014-05-23 14:13 GMT+02:00 Brian Bouterse : >>>>> >>>>> Thanks for the reply Eric. It is great to get hints and suggestions >>>>>> from the author directly! >>>>>> >>>>>> I used the commands that you outlined, and they were able to >>>>>> reproduce the expected output verbatim, which is great. Thanks for clearing >>>>>> that up for me. I've gotten further towards my goal. >>>>>> >>>>>> I've now got three questions: >>>>>> >>>>>> 1) I get an unexpected result when I run ms_var() on my own data. I >>>>>> run these commands: >>>>>> >>>>>> load('arhmm_example.dat') >>>>>> nb_states=2 >>>>>> switch_var=2 // variances are switching >>>>>> var_opt=3 // unrestricted var-cov matrix >>>>>> >>>>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>>> >>>>>> I receive this output: >>>>>> >>>>>> WARNING: in overlay, series number 2 has been ignored because of a >>>>>> bad frequency >>>>>> !--error 10000 >>>>>> series ends before the end date of the bounds >>>>>> at line 39 of function ts2vec0 called by : >>>>>> at line 101 of function explone called by : >>>>>> at line 253 of function ms_var called by : >>>>>> >>>>>> r=ms_var('cte',3,'ardata',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>>> >>>>>> The arhmm_example.dat file is available here[0], and it was made by >>>>>> running.the following command on the original csv[1] file arhmm_example.csv. >>>>>> >>>>>> impexc2bd('arhmm_example.csv', ';', 'arhmm_example.dat') >>>>>> >>>>>> I believe I either don't have the dates configured correctly, or it >>>>>> requires a specific number of data points to match the frequency value, >>>>>> which also may be wrong. Do you have some insight into this error message? >>>>>> I've been reading the docs on the ts structure, and I will continue to try >>>>>> to solve this roadblock. >>>>>> >>>>>> >>>>>> 2) My goal in doing all this is to analyze Autoregressive Hidden >>>>>> Markov Models. As I understand it, the VAR-HMM that ms_var provides is a >>>>>> multivariate case of an Autoregressive Hidden Markov Model. The terms >>>>>> Markov Switching Model, and Hidden Markov Model refer to the same thing. >>>>>> Using a single variable with ms_var() as I show above in the example, will >>>>>> simulate an AR-HMM(3). I would like to check if these statements agree with >>>>>> your understanding. >>>>>> >>>>>> >>>>>> 3) How could I contribute to the grocer code. At the very least I >>>>>> could improve the docs some. >>>>>> >>>>>> >>>>>> [0]: >>>>>> https://s3.amazonaws.com/dfsklfdsklfds/fdsjkfsdjkfds/arhmm_example.dat >>>>>> [1]: >>>>>> https://s3.amazonaws.com/dfsklfdsklfds/fdsjkfsdjkfds/arhmm_example.csv >>>>>> >>>>>> Thanks, >>>>>> Brian >>>>>> >>>>>> >>>>>> On Tue, May 20, 2014 at 3:55 PM, Eric Dubois < >>>>>> grocer.toolbox at gmail.com> wrote: >>>>>> >>>>>>> Hello Brian >>>>>>> >>>>>>> Sorry for this late answer, but I have been quite busy these days. >>>>>>> >>>>>>> I did not notice the problem with this version of the MS programs. >>>>>>> Indeed I have changed the optimization device of all GROCER programs and I >>>>>>> have not adapted the defaults for the MS programs. >>>>>>> >>>>>>> If you run: >>>>>>> --> global GROCERDIR; >>>>>>> --> load(GROCERDIR+'\data\us_revu.dat') >>>>>>> --> bounds('1967m4','2004m2') >>>>>>> >>>>>>> --> nb_states=2 >>>>>>> --> switch_var=2 // variances are switching >>>>>>> --> var_opt=3 // unrestricted var-cov matrix >>>>>>> >>>>>>> --> >>>>>>> r=ms_var('cte',3,'100*(log(us_revu)-lagts(2,log(us_revu)))',nb_states,switch_var,var_opt,'prt=initial;final','opt_convg=0') >>>>>>> (see chapter 6 of the manual for explanations) >>>>>>> >>>>>>> Then the results of the ms_var demo is restaured. >>>>>>> >>>>>>> I will change the default in Grocer next version . >>>>>>> >>>>>>> Regards. >>>>>>> >>>>>>> ?ric. >>>>>>> >>>>>>> >>>>>>> 2014-05-19 12:56 GMT+02:00 Brian Bouterse : >>>>>>> >>>>>>>> Hi Scilab community! >>>>>>>> >>>>>>>> I'm new to Scilab, and the AR-HMM and VAR-HMM solving capabilities >>>>>>>> of GROCER are what interest me. >>>>>>>> >>>>>>>> I have a question relating to Chapter 23 from the GROCER manual[0]. >>>>>>>> This is the univariate MS-AR(3) solved using the function ms_reg_d() on >>>>>>>> the us_revu.dat data included with GROCER. I have made no adjustment from >>>>>>>> the example statements in Chapter 23. >>>>>>>> >>>>>>>> The example output is shown on pages 4 and 5 of the Chapter 23 >>>>>>>> module. Compare that against the output I receive. >>>>>>>> >>>>>>>> http://fpaste.org/102978/14004958/ >>>>>>>> >>>>>>>> Here are my questions: >>>>>>>> >>>>>>>> 1. The numerical output is completely different. I expected it to >>>>>>>> be the same since the data is provided by GROCER, and I've done the example >>>>>>>> exactly as shown in Chapter 23. Is there some explanation to why the >>>>>>>> solved solution I receive is different than the example output in the >>>>>>>> chapter? >>>>>>>> >>>>>>>> 2. I see output like %i*8.4469016 which seems like an error >>>>>>>> because %i looks like a variable that yet needs to be replaced, and then >>>>>>>> multiplied to get to its final value. Is this some kind of bug or error? >>>>>>>> >>>>>>>> Thanks for any help the community can provide. We'll be using this >>>>>>>> for a seminar on HMM, AR-HMM, and VAR-HMM at North Carolina State >>>>>>>> University. I'm also a developer, so I really appreciate all the effort >>>>>>>> that has been put into scilab and GROCER. >>>>>>>> >>>>>>>> Thanks, >>>>>>>> Brian >>>>>>>> >>>>>>>> >>>>>>>> [0]: http://dubois.ensae.net/Grocer_manual_v1.6.zip >>>>>>>> >>>>>>>> >>>>>>>> -- >>>>>>>> Brian Bouterse >>>>>>>> >>>>>>>> _______________________________________________ >>>>>>>> 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 >>>>>>> >>>>>>> >>>>>> >>>>>> >>>>>> -- >>>>>> Brian Bouterse >>>>>> >>>>>> _______________________________________________ >>>>>> 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 >>>>> >>>>> >>>> >>>> >>>> -- >>>> Brian Bouterse >>>> >>>> _______________________________________________ >>>> 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 >>> >>> >> >> >> -- >> Brian Bouterse >> >> _______________________________________________ >> 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 > > -- Brian Bouterse -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Sat Jun 21 17:22:59 2014 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sat, 21 Jun 2014 17:22:59 +0200 Subject: [Scilab-users] Large scale linear optimization In-Reply-To: References: Message-ID: <53A5A353.4080902@gmail.com> Hi Ashish Jaen I of course cannot know about your problem explicitly, but I'd say that 2000 values (variables) is no problem. I've worked with vectors up to 2 million in size. Also the 5000 constraints, well, it sounds like it can be done. There is a problem with Scilab in its current form having limits to managing memory. You should execute the command stacksize('max'); to use as much available memory as possible. I believe that Scilab is working on removing this limit in a future edition of Scilab (v. 6). Whether you will run into a memory problem I can't say. I guess, the only way to find out is to try. Best regards, Claus On 19-Jun-14 09:44, Ashish Jaen wrote: > Hi, > I am new to this mailing list, so please pardon me if this question > has been asked before. > > Has anyone used Scilab for large scale linear optimization ? The > problem size I am looking at is about 2000 variables and 5000 constraints. > > Will be great if you can also share some performance numbers. > > I have used Matlab in the past but no longer have the license as I am > not in university any more :( > > Regards. > > > _______________________________________________ > 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 -------------- next part -------------- An HTML attachment was scrubbed... URL: From tanthiamhuat at yahoo.com Sun Jun 22 06:20:33 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Sat, 21 Jun 2014 21:20:33 -0700 (PDT) Subject: [Scilab-users] Partial Fraction: how to display the variable in x? Message-ID: <1403410833956-4030748.post@n3.nabble.com> -->x=poly(0,'x') x = x -->pfss((3*x+5)/(x^2+3*x+2)) ans = ans(1) 2 ----- 1 + s ans(2) 1 ----- 2 + s ans(1) and ans(2) are displayed in the variable 's'. How do I make them to display as variable 'x'? -- View this message in context: http://mailinglists.scilab.org/Partial-Fraction-how-to-display-the-variable-in-x-tp4030748.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From tanthiamhuat at yahoo.com Sun Jun 22 06:59:10 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Sat, 21 Jun 2014 21:59:10 -0700 (PDT) Subject: [Scilab-users] pfss question... Message-ID: <1403413150618-4030749.post@n3.nabble.com> this returns this answer... -->x=poly(0,'x') x = x -->pfss((x^2)/(x^3-3*x+2)) ans = ans(1) 0.4444444 --------- 2 + s ans(2) - 0.2222222 + 0.5555556s ---------------------- 2 1 - 2s + s But I expect this answer: 5/(9*(x - 1)) + 1/(3*(x - 1)^2) + 4/(9*(x + 2)) how can Scilab return above answer? -- View this message in context: http://mailinglists.scilab.org/pfss-question-tp4030749.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Sun Jun 22 13:34:13 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 22 Jun 2014 13:34:13 +0200 Subject: [Scilab-users] pfss question... In-Reply-To: <1403413150618-4030749.post@n3.nabble.com> References: <1403413150618-4030749.post@n3.nabble.com> Message-ID: <53A6BF35.2030406@free.fr> Hello, Le 22/06/2014 06:59, tanthiamhuat a ?crit : > > -->pfss((x^2)/(x^3-3*x+2)) > ans = > ans(1) > > 0.4444444 > --------- > 2 + s > > ans(2) > > - 0.2222222 + 0.5555556s > ---------------------- > 2 > 1 - 2s + s > > But I expect this answer: > 5/(9*(x - 1)) + 1/(3*(x - 1)^2) + 4/(9*(x + 2)) > > how can Scilab return above answer? Apart the x=>s bug, the uncomplete decomposition is due to the uncomplete factorization/simplification. Such a limitation has been reported here: http://bugzilla.scilab.org/show_bug.cgi?id=4191 This bug has been tagged FIXED instead of WONT (be fixed), while IMO the algorithm could post-process the initial set of raw roots, agregating subsets of roots near enough from each others to be considered as the same (average over the subset). In this way, i think that roots of high multiplicity would be better handled. Regards Samuel PS : a similar effect is reported there for complex coefficients: http://bugzilla.scilab.org/show_bug.cgi?id=8099 From sgougeon at free.fr Sun Jun 22 13:41:51 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 22 Jun 2014 13:41:51 +0200 Subject: [Scilab-users] pfss question... In-Reply-To: <53A6BF35.2030406@free.fr> References: <1403413150618-4030749.post@n3.nabble.com> <53A6BF35.2030406@free.fr> Message-ID: <53A6C0FF.4030504@free.fr> Le 22/06/2014 13:34, Samuel Gougeon a ?crit : > Apart the x=>s bug, the uncomplete decomposition is due to the > uncomplete factorization/simplification. Illustration for the case you have submitted: -->f = polfact(x^3-3*x+2) f = 1 2 + x - 1 + x - 1.0000000 + x -->f(3)==f(4) ans = F From tanthiamhuat at yahoo.com Mon Jun 23 01:25:19 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Sun, 22 Jun 2014 16:25:19 -0700 (PDT) Subject: [Scilab-users] pfss question... In-Reply-To: <53A6C0FF.4030504@free.fr> References: <1403413150618-4030749.post@n3.nabble.com> <53A6BF35.2030406@free.fr> <53A6C0FF.4030504@free.fr> Message-ID: <1403479519573-4030752.post@n3.nabble.com> so would those bugs be fixed in the next version of Scilab? -- View this message in context: http://mailinglists.scilab.org/pfss-question-tp4030749p4030752.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Mon Jun 23 07:14:39 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 23 Jun 2014 07:14:39 +0200 Subject: [Scilab-users] pfss question... In-Reply-To: <1403479519573-4030752.post@n3.nabble.com> References: <1403413150618-4030749.post@n3.nabble.com> <53A6BF35.2030406@free.fr> <53A6C0FF.4030504@free.fr> <1403479519573-4030752.post@n3.nabble.com> Message-ID: <53A7B7BF.7000108@free.fr> Le 23/06/2014 01:25, tanthiamhuat a ?crit : > so would those bugs be fixed in the next version of Scilab? For the "x"=>"s" one, i reported it and posted its fix yesterday: http://bugzilla.scilab.org/13482 http://bugzilla.scilab.org/13481 So i hope so. For the roots agregation issue: i wrote a cleanroots() function in 2010 doing the basic job, to be inserted in polfact() etc. I will review it against Scilab 5.5.0 and will post it ASAP. Regards Samuel From jrafaelbguerra at hotmail.com Mon Jun 23 16:59:43 2014 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Mon, 23 Jun 2014 15:59:43 +0100 Subject: [Scilab-users] Scilab colorbar titles? Message-ID: Hello, Does anybody know how to write titles above colorbars (say: "series1","series2") as in example below: clf;X=1:10;Y=1:20;surf(X,Y,-X.*.Y'); colorbar(-20,-1);X=21:30;surf(X,Y,X.*.Y'); colorbar(1,20); Thanks and regards, Rafael G. -------------- next part -------------- An HTML attachment was scrubbed... URL: From jrafaelbguerra at hotmail.com Mon Jun 23 18:10:45 2014 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Mon, 23 Jun 2014 17:10:45 +0100 Subject: [Scilab-users] Scilab colorbar titles? In-Reply-To: References: Message-ID: And a second question: - How to make the colorbar numbers more legible? In the example below the minus sign overlays (in my PC) the colorbar tick marks and cannot be read Thanks and regards Rafael G. From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Rafael Guerra Sent: Monday, June 23, 2014 4:00 PM To: 'International users mailing list for Scilab.' Subject: [Scilab-users] Scilab colorbar titles? Hello, Does anybody know how to write titles above colorbars (say: "series1","series2") as in example below: clf;X=1:10;Y=1:20;surf(X,Y,-X.*.Y'); colorbar(-20,-1);X=21:30;surf(X,Y,X.*.Y'); colorbar(1,20); Thanks and regards, Rafael G. -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Mon Jun 23 22:12:23 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 23 Jun 2014 22:12:23 +0200 Subject: [Scilab-users] Scilab colorbar titles? In-Reply-To: References: Message-ID: <53A88A27.6050204@free.fr> Le 23/06/2014 18:10, Rafael Guerra a ?crit : > > And a second question: > > -How to make the colorbar numbers more legible? In the example below > the minus sign overlays (in my PC) the colorbar tick marks and cannot > be read > This bug of 5.5.0 has been recently fixed. -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Mon Jun 23 22:18:00 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 23 Jun 2014 22:18:00 +0200 Subject: [Scilab-users] Scilab colorbar titles? In-Reply-To: References: Message-ID: <53A88B78.7030108@free.fr> Hello Rafael, Le 23/06/2014 16:59, Rafael Guerra a ?crit : > > Hello, > > Does anybody know how to write titles above colorbars (say: > "series1","series2") as in example below: > clf;X=1:10;Y=1:20;surf(X,Y,-X.*.Y'); colorbar(-20,-1); e = gce(); e.parent.title.text="Series 1"; X=21:30; surf(X,Y,X.*.Y'); colorbar(1,20); e = gce(); e.parent.title.text="Series 2"; HTH Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From rei.listas at yahoo.com Mon Jun 23 22:56:33 2014 From: rei.listas at yahoo.com (Reinaldo) Date: Mon, 23 Jun 2014 13:56:33 -0700 Subject: [Scilab-users] I have problem to plot a graph of the trigonometric function. Message-ID: <1403556993.75985.YahooMailNeo@web121602.mail.ne1.yahoo.com> Hi Scilab users, I have a problem to plot a graph of the trigonometric function, as follows: // Close all opened figures and clear workspace //xdel(winsid()); clear;clc;function f=myfunction(x) f = x*sin(x) endfunction // Plot scf(1);clf(1);xdata= linspace ( -10 ,10 ,100 );ydata= myfunction ( xdata);plot ( xdata,ydata);xlabel("$-5\le x\le 5$","fontsize",4,"color","red");ylabel("y = x . sin(x)","fontsize",4,"color","red");title("Function","red","fontsize",4);legend("Function evaluation"); Could someone tell me what's wrong? Thank you in advance. All best, Reinaldo. -------------- next part -------------- An HTML attachment was scrubbed... URL: From jrafaelbguerra at hotmail.com Mon Jun 23 23:16:09 2014 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Mon, 23 Jun 2014 22:16:09 +0100 Subject: [Scilab-users] Scilab colorbar titles? In-Reply-To: <53A88B78.7030108@free.fr> References: <53A88B78.7030108@free.fr> Message-ID: Hi Samuel, That is brilliant! Thank you. Rafael From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Gougeon Sent: Monday, June 23, 2014 9:18 PM To: International users mailing list for Scilab. Subject: Re: [Scilab-users] Scilab colorbar titles? Hello Rafael, Le 23/06/2014 16:59, Rafael Guerra a ?crit : Hello, Does anybody know how to write titles above colorbars (say: ?series1?,?series2?) as in example below: clf;X=1:10;Y=1:20;surf(X,Y,-X.*.Y'); colorbar(-20,-1); e = gce(); e.parent.title.text = "Series 1"; X = 21:30; surf(X,Y,X.*.Y'); colorbar(1,20); e = gce(); e.parent.title.text = "Series 2"; HTH Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From tanthiamhuat at yahoo.com Tue Jun 24 03:09:02 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Mon, 23 Jun 2014 18:09:02 -0700 (PDT) Subject: [Scilab-users] What is D in Scilab? what is its use? Message-ID: <1403572142661-4030760.post@n3.nabble.com> -->1D ans = 1. -->2D ans = 2. -->3D ans = 3. -->3D-2 ans = 0.03 -->4D-2 ans = 0.04 -->5D-2 ans = 0.05 -->5D-3 ans = 0.005 -->5D-4 ans = 0.0005 -->5D-6 ans = 0.000005 -- View this message in context: http://mailinglists.scilab.org/What-is-D-in-Scilab-what-is-its-use-tp4030760.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From ashishjaen at gmail.com Tue Jun 24 08:31:20 2014 From: ashishjaen at gmail.com (Ashish Jaen) Date: Mon, 23 Jun 2014 23:31:20 -0700 Subject: [Scilab-users] Large scale linear optimization In-Reply-To: <53A5A353.4080902@gmail.com> References: <53A5A353.4080902@gmail.com> Message-ID: Hi Claus, Thanks for your answer. If possible, can you share your performance numbers ? Details about Runtime / Hardware / RAM / Number of constraints will help. Thanks again, Ashish On Sat, Jun 21, 2014 at 8:22 AM, Claus Futtrup wrote: > Hi Ashish Jaen > > I of course cannot know about your problem explicitly, but I'd say that > 2000 values (variables) is no problem. I've worked with vectors up to 2 > million in size. Also the 5000 constraints, well, it sounds like it can be > done. > > There is a problem with Scilab in its current form having limits to > managing memory. You should execute the command stacksize('max'); to use as > much available memory as possible. I believe that Scilab is working on > removing this limit in a future edition of Scilab (v. 6). > > Whether you will run into a memory problem I can't say. I guess, the only > way to find out is to try. > > Best regards, > Claus > > > On 19-Jun-14 09:44, Ashish Jaen wrote: > > Hi, > I am new to this mailing list, so please pardon me if this question has > been asked before. > > Has anyone used Scilab for large scale linear optimization ? The problem > size I am looking at is about 2000 variables and 5000 constraints. > > Will be great if you can also share some performance numbers. > > I have used Matlab in the past but no longer have the license as I am not > in university any more :( > > Regards. > > > _______________________________________________ > users mailing listusers at lists.scilab.orghttp://lists.scilab.org/mailman/listinfo/users > > > > > ------------------------------ > > > This email is free from viruses and malware because avast! Antivirus > protection is active. > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.bignier at scilab-enterprises.com Tue Jun 24 08:52:45 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Tue, 24 Jun 2014 08:52:45 +0200 Subject: [Scilab-users] What is D in Scilab? what is its use? In-Reply-To: <1403572142661-4030760.post@n3.nabble.com> References: <1403572142661-4030760.post@n3.nabble.com> Message-ID: <53A9203D.9060904@scilab-enterprises.com> Hello, You may replace "XdY" and "XeY" with "X*10^Y". So 3d2 = 3*10^2 = 300. Regards, Paul On 06/24/2014 03:09 AM, tanthiamhuat wrote: > -->1D > ans = > > 1. > > -->2D > ans = > > 2. > > -->3D > ans = > > 3. > > -->3D-2 > ans = > > 0.03 > > -->4D-2 > ans = > > 0.04 > > -->5D-2 > ans = > > 0.05 > > -->5D-3 > ans = > > 0.005 > > -->5D-4 > ans = > > 0.0005 > > -->5D-6 > ans = > > 0.000005 > > > > -- > View this message in context: http://mailinglists.scilab.org/What-is-D-in-Scilab-what-is-its-use-tp4030760.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 tanthiamhuat at yahoo.com Tue Jun 24 09:25:05 2014 From: tanthiamhuat at yahoo.com (tanthiamhuat) Date: Tue, 24 Jun 2014 00:25:05 -0700 (PDT) Subject: [Scilab-users] What is D in Scilab? what is its use? In-Reply-To: <53A9203D.9060904@scilab-enterprises.com> References: <1403572142661-4030760.post@n3.nabble.com> <53A9203D.9060904@scilab-enterprises.com> Message-ID: <1403594705951-4030764.post@n3.nabble.com> thanks. now I see... but I still prefer -->3*10^-3 ans = 0.003 to be more readable than -->3d-3 ans = 0.003 -- View this message in context: http://mailinglists.scilab.org/What-is-D-in-Scilab-what-is-its-use-tp4030760p4030764.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From pierre-aime.agnel at scilab-enterprises.com Tue Jun 24 09:36:13 2014 From: pierre-aime.agnel at scilab-enterprises.com (=?ISO-8859-1?Q?Pierre-Aim=E9_Agnel?=) Date: Tue, 24 Jun 2014 09:36:13 +0200 Subject: [Scilab-users] I have problem to plot a graph of the trigonometric function. In-Reply-To: <1403556993.75985.YahooMailNeo@web121602.mail.ne1.yahoo.com> References: <1403556993.75985.YahooMailNeo@web121602.mail.ne1.yahoo.com> Message-ID: <53A92A6D.1040402@scilab-enterprises.com> Hi, You should use the element-wise multiplication .*instead of the matrix multiplication * in your function f = x * sin(x) You also had a missing property "color" in the function title (and I modified your xdata to be between -5 and 5) // Close all opened figures and clear workspace //xdel(winsid()); clear; clc; function f=myfunction(x) f = x.*sin(x) endfunction // Plot scf(1); clf(1); xdata = linspace ( -5 , 5 , 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); xlabel("$-5\le x\le 5$","fontsize",4,"color","red"); ylabel("y = x . sin(x)","fontsize",4,"color","red"); title("Function","color", "red","fontsize",4); legend("Function evaluation"); Best On 06/23/2014 10:56 PM, Reinaldo wrote: > Hi Scilab users, > > I have a problem to plot a graph of the trigonometric function, as > follows: > > // Close all opened figures and clear workspace > //xdel(winsid()); > clear; > clc; > > function f=myfunction(x) > f = x*sin(x) > endfunction > > // Plot > scf(1); > clf(1); > > xdata= linspace ( -10 ,10 ,100 ); > ydata= myfunction ( xdata); > plot ( xdata , ydata); > xlabel("$-5\le x\le 5$","fontsize",4,"color","red"); > ylabel("y = x . sin(x)","fontsize",4,"color","red"); > title("Function","red","fontsize",4); > legend("Function evaluation"); > Could someone tell me what's wrong? > Thank you in advance. > All best, > Reinaldo. > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Pierre-Aim? Agnel Development Team ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.68 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Tue Jun 24 09:46:03 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Tue, 24 Jun 2014 09:46:03 +0200 Subject: [Scilab-users] What is D in Scilab? what is its use? In-Reply-To: <1403594705951-4030764.post@n3.nabble.com> References: <1403572142661-4030760.post@n3.nabble.com> <53A9203D.9060904@scilab-enterprises.com> <1403594705951-4030764.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F4087751A8@301EX00100.sidel.com> Hello, > De : tanthiamhuat > Envoy? : mardi 24 juin 2014 09:25 > > but I still prefer > -->3*10^-3 > [...] to be more readable than > -->3d-3 I think this is quite different. When you write "3e-3", you define a number; "0.003" is just another way to write it. When you write "3*10^-3", you define an operation which result is 3e-3. So this takes more time and more resource. You can see it if you type "10^10". The difference is negligible for one number, but imagine you handle a matrix with millions of reals or complexes... 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 abhi24394 at gmail.com Mon Jun 23 17:57:25 2014 From: abhi24394 at gmail.com (abhishek sharma) Date: Mon, 23 Jun 2014 21:27:25 +0530 Subject: [Scilab-users] Queries about vector and complex numbers Message-ID: I want to know how to perform divergence and curl of a vector in scilab,also how to change a cimplex no. From one form to another. I am using scilab version 5.4.1 Thanks -------------- next part -------------- An HTML attachment was scrubbed... URL: From digitalchkn at yahoo.com Tue Jun 24 03:18:29 2014 From: digitalchkn at yahoo.com (digitalchkn) Date: Mon, 23 Jun 2014 18:18:29 -0700 (PDT) Subject: [Scilab-users] Pausing script execution for custom GUI Message-ID: <1403572709715-4030761.post@n3.nabble.com> Hi folks,I am a bit new to GUI programming and have a question about structuring code execution when involving a custom GUI figure. What I have is roughly the following*Top level file* T1: execute some initial code (inits global variables, reads in data from some files, some initial code..) T2: execute *GUI script* using the "exec" T3: remainder of main code *GUI script* G1: function handlers for the GUI buttons below G2: creates a simple GUI with "uicontrols" that includes custom buttons. Buttons call the handlers of step G1 when pressed to set a global variableWhen I execute the top level file, the whole thing executes sequentially:T1->T2->G2->T3and user has no opportunity to make a selection. What I want to do is have it execute steps:T1 -> T2 -> G2 but *hold* execution of step T3 until user pressed the buttons in the GUI and G1 is finished executingI tried stragetically inserting the "pause" after G2 but inserting resume in each of the functions in G1, but what that does is just prevent the button handlers of G1 to execute altogether. Anyone has a good suggestion on how to deal with this? -- View this message in context: http://mailinglists.scilab.org/Pausing-script-execution-for-custom-GUI-tp4030761.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From rei.listas at yahoo.com Tue Jun 24 17:07:53 2014 From: rei.listas at yahoo.com (Reinaldo) Date: Tue, 24 Jun 2014 08:07:53 -0700 Subject: [Scilab-users] I have problem to plot a graph of the trigonometric function. In-Reply-To: <53A92A6D.1040402@scilab-enterprises.com> References: <1403556993.75985.YahooMailNeo@web121602.mail.ne1.yahoo.com> <53A92A6D.1040402@scilab-enterprises.com> Message-ID: <1403622473.37198.YahooMailNeo@web121601.mail.ne1.yahoo.com> Hi Pierre-Aim?, Thank you for your help! All best, Reinaldo.? ________________________________ De: Pierre-Aim? Agnel Para: Reinaldo ; International users mailing list for Scilab. Enviadas: Ter?a-feira, 24 de Junho de 2014 4:36 Assunto: Re: [Scilab-users] I have problem to plot a graph of the trigonometric function. Hi, You should use the element-wise multiplication .*instead of the matrix multiplication * in your function f = x * sin(x) You also had a missing property "color" in the function title (and I modified your xdata to be between -5 and 5) // Close all opened figures and clear workspace //xdel(winsid()); clear; clc; function f=myfunction(x) f = x.*sin(x) endfunction // Plot scf(1); clf(1); xdata = linspace ( -5 , 5 , 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); xlabel("$-5\le x\le 5$","fontsize",4,"color","red"); ylabel("y = x . sin(x)","fontsize",4,"color","red"); title("Function","color", "red","fontsize",4); legend("Function evaluation"); Best On 06/23/2014 10:56 PM, Reinaldo wrote: Hi Scilab users, > > >I have a problem to plot a graph of the trigonometric function, as follows: > > >// Close all opened figures and clear workspace //xdel(winsid()); clear; clc; function f=myfunction(x) f = x*sin(x) endfunction // Plot scf(1); clf(1); xdata = linspace ( -10 , 10 , 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); xlabel("$-5\le x\le 5$","fontsize",4,"color","red"); ylabel("y = x . sin(x)","fontsize",4,"color","red"); title("Function","red","fontsize",4); legend("Function evaluation"); >Could someone tell me what's wrong? >Thank you in advance. >All best, >Reinaldo. > > >_______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -- Pierre-Aim? Agnel Development Team ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.68 http://www.scilab-enterprises.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Wed Jun 25 09:09:30 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Wed, 25 Jun 2014 09:09:30 +0200 Subject: [Scilab-users] Queries about vector and complex numbers In-Reply-To: References: Message-ID: <3B5FFC67498DFF49AE7271A584867D16F408800B58@301EX00100.sidel.com> Hello, > De : abhishek sharma > Envoy? : lundi 23 juin 2014 17:57 > > I want to know how to perform divergence and curl of a vector in scilab, > also how to change a cimplex no. From one form to another. I'm afraid you have to do it yourself. For the divergence, you might use something like sum(numderivative(f, x)) For the curl, you might as well use the approximate jacobian matrix calculated by numderivative(). Sorry, I don't know what you mean by "cimplex no.". 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 paul.bignier at scilab-enterprises.com Wed Jun 25 10:02:56 2014 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Wed, 25 Jun 2014 10:02:56 +0200 Subject: [Scilab-users] Queries about vector and complex numbers In-Reply-To: <3B5FFC67498DFF49AE7271A584867D16F408800B58@301EX00100.sidel.com> References: <3B5FFC67498DFF49AE7271A584867D16F408800B58@301EX00100.sidel.com> Message-ID: <53AA8230.4010407@scilab-enterprises.com> also how to change a cimplex no. From one form to another. Hello, "how to change a cimplex no. From one form to another" Here is a link to the available functions to manipulate complex numbers. Note that the transpose operator is special for complex numbers (returns the conjugate transpose of the input). Regards, Paul On 06/25/2014 09:09 AM, Dang, Christophe wrote: > Hello, > >> De : abhishek sharma >> Envoy? : lundi 23 juin 2014 17:57 >> >> I want to know how to perform divergence and curl of a vector in scilab, >> also how to change a cimplex no. From one form to another. > I'm afraid you have to do it yourself. > > For the divergence, you might use something like > > sum(numderivative(f, x)) > > For the curl, you might as well use the approximate jacobian matrix calculated by numderivative(). > > Sorry, I don't know what you mean by "cimplex no.". > > Hope this helps > -- 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 sgougeon at free.fr Wed Jun 25 10:46:12 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Wed, 25 Jun 2014 10:46:12 +0200 Subject: [Scilab-users] Queries about vector and complex numbers In-Reply-To: References: Message-ID: <53AA8C54.2090002@free.fr> Hello, Le 23/06/2014 17:57, abhishek sharma a ?crit : > > ...also how to change a cimplex no. From one form to another. I am > using scilab version 5.4.1 > modulus (from cartesian) : abs(c) phase from cartesian: atan(imag(c), real(c)) real part from polar (m, phase): re = m.*cos(phase) imaginary part : im = m.*sin(phase) Beware that [m, phase] = polar(c) works as expected here only if c is a scalar complex. Regards Samuel From grivet at cnrs-orleans.fr Wed Jun 25 14:05:43 2014 From: grivet at cnrs-orleans.fr (grivet) Date: Wed, 25 Jun 2014 14:05:43 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <1402464955580-4030713.post@n3.nabble.com> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> Message-ID: <53AABB17.1080407@cnrs-orleans.fr> Hello, I have noticed the following small bug in Scilab 5.5.0 When exporting a graph to an .eps file, a pop-up window proposes to choose between portrait or landscape orientation. This option is for me without ffect, the orientation of the resulting figure is always landcape. Cheers, JP Grivet From calixte.denizet at scilab-enterprises.com Wed Jun 25 18:45:46 2014 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Wed, 25 Jun 2014 18:45:46 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53AABB17.1080407@cnrs-orleans.fr> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> Message-ID: <53AAFCBA.4030701@scilab-enterprises.com> Hi JP, I just tryed with a simple plot2d and that worked for me. So could you provide me a test case and the generated eps file ? Best regards, Calixte On 25/06/2014 14:05, grivet wrote: > Hello, > I have noticed the following small bug in Scilab 5.5.0 > When exporting a graph to an .eps file, a pop-up window proposes to > choose between portrait or landscape orientation. > This option is for me without ffect, the orientation of the resulting > figure is always landcape. > Cheers, > JP Grivet > _______________________________________________ > 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 From fahn at zhaw.ch Wed Jun 25 22:08:34 2014 From: fahn at zhaw.ch (Fahrni Pascal (fahn)) Date: Wed, 25 Jun 2014 20:08:34 +0000 Subject: [Scilab-users] user terminate program (exit while loop) Message-ID: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> Dear scilab users, I'm using a scilab script to monitor a serial port, read from it and do some calculations based on that data. The program is wrapped by a while loop and contains an xpause statement. What is the best way to allow the user to terminate the program by user input? Best wishes Pascal -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Wed Jun 25 23:01:05 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Wed, 25 Jun 2014 23:01:05 +0200 Subject: [Scilab-users] user terminate program (exit while loop) In-Reply-To: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> References: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> Message-ID: <53AB3891.4030703@free.fr> Hello, Le 25/06/2014 22:08, Fahrni Pascal (fahn) a ?crit : > > Dear scilab users, > > I'm using a scilab script to monitor a serial port, read from it and > do some calculations based on that data. The program is wrapped by a > while loop and contains an xpause statement. What is the best way to > allow the user to terminate the program by user input? > This is an untested proposal: create a virtual figure and use its event handler sensitive to on-the-flight keyboard events: * create a figure * set its .visible property to "off" * set properly its .event_handler and its event_handler_enable properties * in the event handler function: catch the code of the event; if it's a keyboard event (or a specific one that you require from the user), issue an abort instruction; HTH Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From michael.ungelehrt at th-deg.de Wed Jun 25 11:00:47 2014 From: michael.ungelehrt at th-deg.de (Michael Ungelehrt) Date: Wed, 25 Jun 2014 11:00:47 +0200 Subject: [Scilab-users] Expression block - Initialisation problem: singularity in a block Message-ID: <006901cf9053$f5cf69d0$e16e3d70$@th-deg.de> Dear users, I copied a simulation which I built up at a Windows 7 PC 64bit with Scilab 5.5.0 to a new Windows 7 PC 64bit with Scilab 5.5.0 installed just a view days ago. On the old PC the simulation run without any problems, on the new PC I get an error message (Initialisation problem: singularity in a block) at the Expression block. I reduced the simulation to a simple model (see in the attachment) to reduce other possible error factors. The simulation you find in the attachment run on one PC, on the other not. I have absolutely no idea where a mistake is. Simulation settings are the same, too. Can you please help me? Best regards, Michael Ungelehrt Logo_TC_Cham Michael Ungelehrt, B.Eng. Laboratory engineer DIT ? Deggendorf Institute of Technology Technologie Campus Cham Badstra?e 21 D-93413 Cham Tel: +49 9971 99673-29 Fax: +49 9971 99673-69 E-Mail: michael.ungelehrt at th-deg.de Web: www.tc-cham.th-deg.de -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: image002.jpg Type: image/jpeg Size: 4175 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: Expression block.zcos Type: application/octet-stream Size: 4162 bytes Desc: not available URL: From grivet at cnrs-orleans.fr Thu Jun 26 11:52:03 2014 From: grivet at cnrs-orleans.fr (grivet) Date: Thu, 26 Jun 2014 11:52:03 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53AAFCBA.4030701@scilab-enterprises.com> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> Message-ID: <53ABED43.4020803@cnrs-orleans.fr> x = 0:200; y = sin(2*%pi*x/50); plot(x,y) I forgot to mention that I work under Win7. Regards, JP Grivet From antoine.monmayrant at laas.fr Thu Jun 26 12:28:35 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Thu, 26 Jun 2014 12:28:35 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53ABED43.4020803@cnrs-orleans.fr> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> <53ABED43.4020803@cnrs-orleans.fr> Message-ID: <53ABF5D3.7070007@laas.fr> On 06/26/2014 11:52 AM, grivet wrote: > x = 0:200; > y = sin(2*%pi*x/50); > plot(x,y) it works for me: I got two different orientations. (under linux ubuntu 12.04 64bits) From calixte.denizet at scilab-enterprises.com Thu Jun 26 12:30:34 2014 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Thu, 26 Jun 2014 12:30:34 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53ABED43.4020803@cnrs-orleans.fr> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> <53ABED43.4020803@cnrs-orleans.fr> Message-ID: <53ABF64A.2050504@scilab-enterprises.com> Could you attach the generated EPS file please ? Calixte On 26/06/2014 11:52, grivet wrote: > x = 0:200; > y = sin(2*%pi*x/50); > plot(x,y) > > I forgot to mention that I work under Win7. > > Regards, > JP Grivet > > _______________________________________________ > 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 From antoine.monmayrant at laas.fr Thu Jun 26 13:56:13 2014 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Thu, 26 Jun 2014 13:56:13 +0200 Subject: [Scilab-users] Writing & reading sod file from another software (Labview) Message-ID: <53AC0A5D.9030701@laas.fr> Hi everyone, Is there a way to write *.sod files (and eventually read them also) from outside Scilab, more specifically from Labview? I need to analyse with Scilab fairly big datasets (ideally a uint16 1280*1024*4096 hypermatrix) acquired with Labview. For the moment I save/open *csv files, but it would be faster and more convenient using Scilab native file format... Antoine -- +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 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 david.cheze at cea.fr Thu Jun 26 14:23:50 2014 From: david.cheze at cea.fr (=?UTF-8?Q?David_Ch=C3=A8ze?=) Date: Thu, 26 Jun 2014 05:23:50 -0700 (PDT) Subject: [Scilab-users] Writing & reading sod file from another software (Labview) In-Reply-To: <53AC0A5D.9030701@laas.fr> References: <53AC0A5D.9030701@laas.fr> Message-ID: <1403785429986-4030785.post@n3.nabble.com> Hi Antoine, I've never used it until now but there's a gateway Scilab/Labview that should do it easily: https://www.scilab.org/fr/scilab/interoperability/connection_applications/labview David -- View this message in context: http://mailinglists.scilab.org/Scilab-users-Writing-reading-sod-file-from-another-software-Labview-tp4030783p4030785.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From grivet at cnrs-orleans.fr Thu Jun 26 15:43:16 2014 From: grivet at cnrs-orleans.fr (grivet) Date: Thu, 26 Jun 2014 15:43:16 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53ABF64A.2050504@scilab-enterprises.com> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> <53ABED43.4020803@cnrs-orleans.fr> <53ABF64A.2050504@scilab-enterprises.com> Message-ID: <53AC2374.9020808@cnrs-orleans.fr> Le 26/06/2014 12:30, Calixte Denizet a ?crit : > Could you attach the generated EPS file please ? > > Calixte > > On 26/06/2014 11:52, grivet wrote: >> x = 0:200; >> y = sin(2*%pi*x/50); >> plot(x,y) >> >> I forgot to mention that I work under Win7. >> >> Regards, >> JP Grivet >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > Here they are. Opening both file with Numpad, I notice that the ounding box is idential in th two cases, and wider than the page. Cordialement JP Grivet -------------- next part -------------- A non-text attachment was scrubbed... Name: paysage.eps Type: application/postscript Size: 14206 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: portrait.eps Type: application/postscript Size: 14205 bytes Desc: not available URL: From clement.david at scilab-enterprises.com Thu Jun 26 17:10:05 2014 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Thu, 26 Jun 2014 17:10:05 +0200 Subject: [Scilab-users] Expression block - Initialisation problem: singularity in a block In-Reply-To: <006901cf9053$f5cf69d0$e16e3d70$@th-deg.de> References: <006901cf9053$f5cf69d0$e16e3d70$@th-deg.de> Message-ID: <1403795405.2024.52.camel@paros> Hello Michael, It seems to be a scilab bug, please report to bugzilla.scilab.org. Thanks in advance, -- Cl?ment Le mercredi 25 juin 2014 ? 11:00 +0200, Michael Ungelehrt a ?crit : > Dear users, > > > > I copied a simulation which I built up at a Windows 7 PC 64bit with > Scilab 5.5.0 to a new Windows 7 PC 64bit with Scilab 5.5.0 installed > just a view days ago. On the old PC the simulation run without any > problems, on the new PC I get an error message (Initialisation > problem: singularity in a block) at the Expression block. I reduced > the simulation to a simple model (see in the attachment) to reduce > other possible error factors. > > > > The simulation you find in the attachment run on one PC, on the other > not. I have absolutely no idea where a mistake is. Simulation settings > are the same, too. > > > > Can you please help me? > > > > Best regards, > > Michael Ungelehrt > > > > > > Logo_TC_Cham > > > > Michael Ungelehrt, B.Eng. > > Laboratory engineer > > > > DIT ? Deggendorf Institute of Technology > Technologie Campus Cham > > Badstra?e 21 > > D-93413 Cham > > > > Tel: +49 9971 99673-29 > > Fax: +49 9971 99673-69 > > E-Mail: michael.ungelehrt at th-deg.de > > Web: www.tc-cham.th-deg.de > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From calixte.denizet at scilab-enterprises.com Thu Jun 26 17:40:02 2014 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Thu, 26 Jun 2014 17:40:02 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53AC2374.9020808@cnrs-orleans.fr> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> <53ABED43.4020803@cnrs-orleans.fr> <53ABF64A.2050504@scilab-enterprises.com> <53AC2374.9020808@cnrs-orleans.fr> Message-ID: <53AC3ED2.6050704@scilab-enterprises.com> For me the EPS files are correct: %%PageOrientation is correct in each case. There is probably a problem with your postscript viewer. But maybe the exporter could rotate the figure itself (as it is done with PDF exporter)... Please report a bug at bugzilla.scilab.org Best regards Calixte On 26/06/2014 15:43, grivet wrote: > Le 26/06/2014 12:30, Calixte Denizet a ?crit : >> Could you attach the generated EPS file please ? >> >> Calixte >> >> On 26/06/2014 11:52, grivet wrote: >>> x = 0:200; >>> y = sin(2*%pi*x/50); >>> plot(x,y) >>> >>> I forgot to mention that I work under Win7. >>> >>> Regards, >>> JP Grivet >>> >>> _______________________________________________ >>> users mailing list >>> users at lists.scilab.org >>> http://lists.scilab.org/mailman/listinfo/users >> >> > Here they are. Opening both file with Numpad, I notice that the > ounding box is idential in th two cases, and wider than the page. > Cordialement > JP Grivet > > > > _______________________________________________ > 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 grivet at cnrs-orleans.fr Thu Jun 26 18:24:54 2014 From: grivet at cnrs-orleans.fr (grivet) Date: Thu, 26 Jun 2014 18:24:54 +0200 Subject: [Scilab-users] orientation of encapsulated postscript In-Reply-To: <53AC3ED2.6050704@scilab-enterprises.com> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> <53ABED43.4020803@cnrs-orleans.fr> <53ABF64A.2050504@scilab-enterprises.com> <53AC2374.9020808@cnrs-orleans.fr> <53AC3ED2.6050704@scilab-enterprises.com> Message-ID: <53AC4956.6040904@cnrs-orleans.fr> Le 26/06/2014 17:40, Calixte Denizet a ?crit : > For me the EPS files are correct: %%PageOrientation is correct in each > case. > There is probably a problem with your postscript viewer. > But maybe the exporter could rotate the figure itself (as it is done > with PDF exporter)... Please report a bug at bugzilla.scilab.org I use gsview32.exe, version 5.0 of 2012-01-17, which seems to perform quite well in other cases. I will report the bug. Cordialement JP Grivet -------------- next part -------------- An HTML attachment was scrubbed... URL: From grivet at cnrs-orleans.fr Thu Jun 26 18:33:34 2014 From: grivet at cnrs-orleans.fr (grivet) Date: Thu, 26 Jun 2014 18:33:34 +0200 Subject: [Scilab-users] bug in fftw In-Reply-To: <53AC3ED2.6050704@scilab-enterprises.com> References: <1401936444905-4030679.post@n3.nabble.com> <1402464955580-4030713.post@n3.nabble.com> <53AABB17.1080407@cnrs-orleans.fr> <53AAFCBA.4030701@scilab-enterprises.com> <53ABED43.4020803@cnrs-orleans.fr> <53ABF64A.2050504@scilab-enterprises.com> <53AC2374.9020808@cnrs-orleans.fr> <53AC3ED2.6050704@scilab-enterprises.com> Message-ID: <53AC4B5E.8070900@cnrs-orleans.fr> Hi, The bug that affected the dct transform in Scilab 5.4.0 is still present in Scilab 5.5.0. To evidence the bug, please run the following lines: x = 0:127; xx = 0:255; y = cos(2*%pi*x/32).*exp(-x/32); yy = cos(2*%pi*xx/64).*exp(-xx/64); w = dct(y); z = fft(y); //ww = dct(yy,"dct1"); where I define a damped sinusoid and attempt to compute an fft *AND* a dct from it. The output is: Warning !!! Scilab has found a critical error (EXCEPTION_ACCESS_VIOLATION) with "fftw" function. Save your data and restart Scilab. On the other hand, if I compute 'ww' instead of 'w', i.e. dct and fft on vectors of different lengths, the problem disappears. It would seem that one of the numerous fftw parameters is not reset between two calls. JP Grivet From sgougeon at free.fr Fri Jun 27 10:07:48 2014 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 27 Jun 2014 10:07:48 +0200 Subject: [Scilab-users] user terminate program (exit while loop) Message-ID: <53AD2654.8080805@free.fr> Hello, Le 25/06/2014 22:08, Fahrni Pascal (fahn) a ?crit : > > Dear scilab users, > > I'm using a scilab script to monitor a serial port, read from it and > do some calculations based on that data. The program is wrapped by a > while loop and contains an xpause statement. What is the best way to > allow the user to terminate the program by user input? > This is an untested proposal: create a virtual figure and use its event handler sensitive to on-the-flight keyboard events: * create a figure * set its .visible property to "off" * set properly its .event_handler and its event_handler_enable properties * in the event handler function: catch the code of the event; if it's a keyboard event (or a specific one that you require from the user), issue an abort instruction; HTH Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From aweeks at hidglobal.com Fri Jun 27 15:45:22 2014 From: aweeks at hidglobal.com (aweeks at hidglobal.com) Date: Fri, 27 Jun 2014 14:45:22 +0100 Subject: [Scilab-users] user terminate program (exit while loop) In-Reply-To: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> References: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> Message-ID: Hi Pascal, I've wanted to do this too. You might like to look at link http://compgroups.net/comp.soft-sys.math.scilab/non-blocking-keyboard-input/2083213 which is very similar to Samuel Gougeon's suggestion. Samuel says his proposal is untested but I use this method and it works for me. I've pasted the basis of my code below and apologise if it doesn't work in your setup because I position the windows and text boxes to suit the layout on my monitor. If you have a question feel free to ask. global figrhdl; global stophdl; global stoploop; function StopLoop(win, x, y, ibut) global stoploop; if ibut == -1000 then stoploop = %t; end if ibut == 3 then stoploop = %t; mprintf("Loop stopped.\n"); end endfunction function StartDisplay() global figrhdl; global stophdl; figrhdl = scf(1); xstring(-0.1,1.1,"LEFT CLICK HERE TO STOP LOOPING"); stophdl = gce(); stophdl.font_style = 0; stophdl.font_size = 2; stophdl.box = 'on'; seteventhandler('StopLoop'); figrhdl.event_handler = 'StopLoop'; figrhdl.event_handler_enable = 'on'; endfunction stoploop = %f; // The endless loop which is to be stopped. while ~stoploop do .... whatever end 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 From: "Fahrni Pascal (fahn)" To: "users at lists.scilab.org" Date: 25/06/2014 21:15 Subject: [Scilab-users] user terminate program (exit while loop) Sent by: "users" Dear scilab users, I?m using a scilab script to monitor a serial port, read from it and do some calculations based on that data. The program is wrapped by a while loop and contains an xpause statement. What is the best way to allow the user to terminate the program by user input? Best wishes Pascal _______________________________________________ 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: 12445500.gif Type: image/gif Size: 4559 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: graycol.gif Type: image/gif Size: 105 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: ecblank.gif Type: image/gif Size: 45 bytes Desc: not available URL: From aweeks at hidglobal.com Fri Jun 27 19:23:47 2014 From: aweeks at hidglobal.com (aweeks at hidglobal.com) Date: Fri, 27 Jun 2014 18:23:47 +0100 Subject: [Scilab-users] user terminate program (exit while loop) In-Reply-To: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> References: <83D0B391F6E5B6479C08794764094D743BC5D8E2@srv-mail-111.zhaw.ch> Message-ID: Hi Pascal, Further to my previous e-mail, I tested a few ideas. If you add a print statement to the StopLoop function you can see what values come back as the mouse is moved. function StopLoop(win, x, y, ibut) global stoploop; mprintf(" win: %d x: %d y: %d ibut: %d\n", win, x, y, ibut) if ibut == -1000 then stoploop = %t; end if ibut == 3 then stoploop = %t; mprintf("Loop stopped.\n"); end endfunction Values are only returned when the mouse cursor is over the graphics window to which the event handler is attached. This is true even if the 'focus' is not on that window. However, no value seems to be returned when a keyboard key is pressed unless the 'focus' is on the graphics window and the mouse cursor is in the window. This is a little disappointing - it would have been nice to be able to 'press any key to stop loop'. To anyone else who may be reading this: Is there a way to switch the 'windows focus' onto a particular window from within the running SciLab program? Is there any chance that a function equivalent to the Borland Turbo C function 'kbhit' could be added to future versions of SciLab? Best wishes, Adrian. 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 From: "Fahrni Pascal (fahn)" To: "users at lists.scilab.org" Date: 25/06/2014 21:15 Subject: [Scilab-users] user terminate program (exit while loop) Sent by: "users" Dear scilab users, I?m using a scilab script to monitor a serial port, read from it and do some calculations based on that data. The program is wrapped by a while loop and contains an xpause statement. What is the best way to allow the user to terminate the program by user input? Best wishes Pascal _______________________________________________ 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: 32475108.gif Type: image/gif Size: 4559 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: graycol.gif Type: image/gif Size: 105 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: ecblank.gif Type: image/gif Size: 45 bytes Desc: not available URL: From breitu at arcor.de Sat Jun 28 11:01:35 2014 From: breitu at arcor.de (breitu) Date: Sat, 28 Jun 2014 02:01:35 -0700 (PDT) Subject: [Scilab-users] CLSS in XCOS: wrong use of input? Message-ID: <1403946095341-4030797.post@n3.nabble.com> Hi all, If I use CLSS to simulate a simple linear system, for example: x'=x and I have the matrices: A=[1]; B=[]; C=[1]; D=[], then I get the exponential function. If I now add a little bit normal noise from the random generator (mean 0 and variance 1) changing B to B=[0.01], I do not get what I expect a slightly blurred exponential function, but practically white noise. ONLY if I set also D=[0], I get the blurred exponential. There is no error message that with B =[0.01] one now has to give a value for D too. It looks as the input from the random generator is used as if D=[1]. Question: Is my observation correct and is this a bug or a feature? Ciao Karl -- View this message in context: http://mailinglists.scilab.org/CLSS-in-XCOS-wrong-use-of-input-tp4030797.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From mi.jaritz at gmail.com Sun Jun 29 00:26:24 2014 From: mi.jaritz at gmail.com (iiyama) Date: Sat, 28 Jun 2014 15:26:24 -0700 (PDT) Subject: [Scilab-users] Calculate electrical field with Scifreefem Message-ID: <1403994384113-4030798.post@n3.nabble.com> Dear all I want to use Scifreefem for calculating the electrical field between two wires. I managed to solve the laplace equation for the voltage but I was not able to implement the right code for the electrical field (E=-grad(V)). I would be grateful for your information.Thx Code: // Standard formulation problem, Equations are described in FreeFeem language. clear close clc // User can enable / disable graphics by setting Graphics variable to %t or %f if ~exists('Graphics') then Graphics = %t; end if Graphics then xset('wdim',700,700); xset('colormap',jetcolormap(254)); colorbar(-100,100) end // Definition of the border deff('[x,y] = f1(t)','x = 1+0.5*cos(t); y = 1+0.5*sin(t)'); //Linker Kreis deff('[x,y] = f2(t)','x = 3+0.5*cos(t); y = 1+0.5*sin(t)'); //Rechter Kreis bord2 = tlist(['border';'a';'b';'c';'d';'e';'f'],... list('x = t; y = 0',0,5,40,1),... list('x = 5; y = t',0,2,20, 1),... list('x = 5 - t; y = 2',0,5,40,1),... list('x = 0; y = t',2,0,20,1),... list(valf(f1, 0:2*%pi/100:2*%pi), 1),... list(valf(f2, 0:2*%pi/100:2*%pi), -1)); // buildMesh(bord2,'th'); [noeul,trianl] = getffResult(); // Get mesh description in Scilab erreur = 0.1; coef = 0.1^(1./5.); // Error level will be divided by 10 every 5 iterations. [noeul,trianl] = getffResult(); // Get mesh description in Scilab u = defvar('u'); // Set the unknown f = defvar('0'); // Set the right hand side of the equation ff_problem(strcat(['solve(u) {',... pde_sol(u,div(grad(u)),f),... 'on(a,b,c,d) u = 0;... on(e) u = 100;... on(f) u = -100;};'])); // 'on(a,b,c,d,e,f) u = 0; };'])); // Get the result as a Scilab variable [noeul,trianl,tFunc] = getffResult('u'); clf(); //xsetech([0,0,1,0.5]); xsetech([0,0,1,1]); resultvisu(noeul,trianl,tFunc); colorbar(-100,100) ff_end(); // Destroy Fem interpretor -- View this message in context: http://mailinglists.scilab.org/Calculate-electrical-field-with-Scifreefem-tp4030798.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From amalia39ar at yahoo.com.ar Sun Jun 29 05:12:56 2014 From: amalia39ar at yahoo.com.ar (Amalia Susana) Date: Sat, 28 Jun 2014 20:12:56 -0700 (PDT) Subject: [Scilab-users] Polynomic regression Message-ID: <1404011576794-4030799.post@n3.nabble.com> Hello: I?m new in Scilab and need some help in calculate a regression analysis between the coefficient of asimmetry and the shape parameter. The equation is very complicated, and I need to search a solution of the equation in polynomic form. The equation is: function y=f(x); y= (x/abs(x))*(-gamma(1+3*x)+3*gamma(1+x)*gamma(1+2*x)-2*gamma(1+x)^3)/ /(gamma(1+2*x)-gamma(1+x)^2)^(3/2) and I want to get a polynomic with x as the dependent variable of y (x is shape parameter and y is the coefficient of assimettry). if x < 0 and (1.14 < y<10) one polynomic expression if x > 0 and (-2< y < 1.14) another polynomic expression if x < 0 and (-10 < y < 0) another how can I shall do something like this? Thanks!! -- View this message in context: http://mailinglists.scilab.org/Polynomic-regression-tp4030799.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Mon Jun 30 10:50:53 2014 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 30 Jun 2014 10:50:53 +0200 Subject: [Scilab-users] Polynomic regression In-Reply-To: <1404011576794-4030799.post@n3.nabble.com> References: <1404011576794-4030799.post@n3.nabble.com> Message-ID: <3B5FFC67498DFF49AE7271A584867D16F4088DE8B9@301EX00100.sidel.com> Hello, I'm afaraid I don't understand very well your problem. > De : Amalia Susana > Envoy? : dimanche 29 juin 2014 05:13 > > need some help in calculate a regression analysis between the coefficient of asimmetry > and the shape parameter. > [...] The equation is: > function y=f(x); y=[...] You give us a function, not an equation. > and I want to get a polynomic with x as the dependent variable of y (x is shape > parameter and y is the coefficient of assimettry). You have a collection of values, and you expect them to match a parametric model. The regression is a way to adjust the parameters of the model. >From what I understand, you have a collection of (x, y) values, and want to find a polynomial relationship x = p(y) is it right? A kind of inverse form of the f function? > if x < 0 and (1.14 < y<10) one polynomic expression if x > 0 and (-2< y < 1.14) another > polynomic expression if x < 0 and (-10 < y < 0) another -- 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. ______________________________________________________________________