From Christophe.Dang at sidel.com Mon Jun 1 09:51:53 2015 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 1 Jun 2015 07:51:53 +0000 Subject: [Scilab-users] Previous Scilab version. In-Reply-To: <32C7C66D72F14F2C9979C801C21638FF@giganti> References: <32C7C66D72F14F2C9979C801C21638FF@giganti> Message-ID: Hello, > De Mario Giganti > Envoy? : mercredi 27 mai 2015 17:20 > > I've bought the book "CURVES AND SURFACES How to draw different shapes through SCILAB" > edited by Pitagora in Bologna, Italy. > > Now I learn that I must download an old version of SCILAB, namely 5.2.1. What do you mean by "I must"? Maybe you could just try with the latest version and see if things go right. > However I got stuck with the question: What kind of archive do you need ? > Binary, Sources, Prerequirements? > [...] > The OS here is 32-bit Windows XP. You probably have got a "standard" configuration, so "binary" should be OK. 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 n.strelkov at gmail.com Mon Jun 1 10:07:31 2015 From: n.strelkov at gmail.com (Nikolay Strelkov) Date: Mon, 1 Jun 2015 11:07:31 +0300 Subject: [Scilab-users] Previous Scilab version. In-Reply-To: <32C7C66D72F14F2C9979C801C21638FF@giganti> References: <32C7C66D72F14F2C9979C801C21638FF@giganti> Message-ID: Dear, Mario! You can download Scilab 5.2.1 for example from here ( http://uni-smr.ac.ru/archive/win/science/math/scilab/5x/win/ ). But I think you do not need exact version, because of stable Scilab syntax. Moreover modern versions of Scilab are more user-friendly. With best regards, maintainer and developer of Mathieu functions toolbox for Scilab, IEEE member, Ph.D., Nikolay Strelkov. 2015-05-27 18:20 GMT+03:00 Mario Giganti : > To Scilab team > > I've bought the book "CURVES AND > SURFACES How to draw different shapes through SCILAB" edited by Pitagora in > Bologna, Italy. > > Now I learn that I must download an old version of SCILAB, namely 5.2.1. > > I was going to mail the "Previous Scilab Version" form. > > However I got stuck with the question: What kind of archive do you need ? > Binary, Sources, Prerequirements? > > I ask kindly help from you as far as *what is the right answer* to the > above request, when I want to download the SCILAB program as is normally > done from the Scilab.org internet site. > > The OS here is 32-bit Windows XP. > > I intend to use the Scilab program merely for personal educational > purpose. No commercial usage is foreseen. > > Thanks for the help. > > GIGANTI MARIO > > Country: Italy > > giganti.mario-wolit at tiscali.it > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From communication at scilab-enterprises.com Mon Jun 1 15:23:26 2015 From: communication at scilab-enterprises.com (Scilab Communications) Date: Mon, 01 Jun 2015 15:23:26 +0200 Subject: [Scilab-users] ScilabTEC 2015 : Exchanges with passion Message-ID: <556C5CCE.7010804@scilab-enterprises.com> Dear Scilab Users, The 7th edition of the International Scilab Users Conference, hosted at Mozilla Space in Paris on May 21 & 22, has been a real success. ScilabTEC attendees have been able to enjoy two days of presentations about innovative applications using Scilab and of fruitful networking moments. We, the Scilab team, want to express our sincere thanks to the community, the speakers and all the attendees for coming together to network, exchange ideas, learn and also have fun with us! Presentations will be available shortly. Meanwhile, you can consult the photos of the event at https://plus.google.com/+scilab/posts/DaDTztgcZtb See you at ScilabTEC 2016! -- Communication Department, Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles (France) http://www.scilab-enterprises.com - http://www.scilab.org -------------- next part -------------- An HTML attachment was scrubbed... URL: From francois.le.berre at free.fr Mon Jun 1 21:26:24 2015 From: francois.le.berre at free.fr (francoislb) Date: Mon, 1 Jun 2015 12:26:24 -0700 (MST) Subject: [Scilab-users] Error: Scilab cannot create Scilab Java Main-class In-Reply-To: <1429714213332-4032122.post@n3.nabble.com> References: <1429714213332-4032122.post@n3.nabble.com> Message-ID: <1433186784804-4032373.post@n3.nabble.com> Hello, i had exactly the same problem as the one described by Floyd. impossible to run scilab V5.5.2 (X64), under windows 7. The error message is the same (in French version : "Scilab ne peut pas acc?der ? la classe principale Java scilab..."). the V5.5.2 was supposed to fix this bug. Fran?ois -- View this message in context: http://mailinglists.scilab.org/Error-Scilab-cannot-create-Scilab-Java-Main-class-tp4032122p4032373.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From brentonhorne77 at gmail.com Thu Jun 4 01:38:48 2015 From: brentonhorne77 at gmail.com (Brenton Horne) Date: Thu, 04 Jun 2015 09:38:48 +1000 Subject: [Scilab-users] Building on 32 bit Ubuntu 15.04 In-Reply-To: <556F8F0C.6070207@ymail.com> References: <556F8F0C.6070207@ymail.com> Message-ID: <556F9008.8000006@gmail.com> Hi, I have tried to install Scilab 5.5.2 on my Ubuntu machine from source code. Now I know there's a binary but I want an installation that creates a /usr/bin/scilab file so that the command 'scilab' launches the Scilab GUI. I have tried creating a symbolic link (via sudo ln -s /path/to/scilab-5.5.2/bin /usr/bin/scilab) to allow this but for whatever reason running 'scilab' from the command line won't seem to launch the Scilab GUI after I create such a link. My build attempt may be summarized by this bash script (before this script was when I downloaded the source tarball and the prerequisite tarball to Downloads using Chrome):| | |||cd ~/Downloads for i in *.tar.gz do sudo tar xvzf $i -C /opt; done sudo chmod 777 -R /opt/scilab-5.5.2 cd /opt/scilab-5.5.2 sudo ./configure sudo make all| ? no errors were encountered up until the last line, the last six lines of output were: |Found dependency from scilab-cli-bin to libstdc++ found Makefile:2092: recipe for target 'check-libstdcpp-dep' failed make[1]: *** [check-libstdcpp-dep] Error 2 make[1]: Leaving directory '/opt/scilab-5.5.2' Makefile:1541: recipe for target 'all-recursive' failed make: *** [all-recursive] Error 1| ? If it is important I have been trying to follow the "Compiling Scilab 5.x under GNU-Linux Unix" Wiki article and to get any prerequisites not in the tarball I ran: |sudo apt-get build-dep scilab| ? The version of Scilab found in the vivid repos is 5.5.1 so I thought that the dependencies for 5.5.2 would be similar enough to allow Scilab 5.5.2 to be built using 5.5.1's dependencies. Thanks for your time, Brenton -------------- next part -------------- An HTML attachment was scrubbed... URL: From tim at wescottdesign.com Thu Jun 4 02:07:53 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Wed, 03 Jun 2015 17:07:53 -0700 Subject: [Scilab-users] Building on 32 bit Ubuntu 15.04 In-Reply-To: <556F9008.8000006@gmail.com> References: <556F8F0C.6070207@ymail.com> <556F9008.8000006@gmail.com> Message-ID: <1433376473.3900.0.camel@servo> For my locally-installed Scilab I have a file at /usr/local/bin/scilab which is marked executable and which reads: #!/bin/bash /usr/local/share/scilab-5.5.1/bin/scilab $* It isn't pretty, but it works. On Thu, 2015-06-04 at 09:38 +1000, Brenton Horne wrote: > Hi, > > I have tried to install Scilab 5.5.2 on my Ubuntu machine from source > code. Now I know there's a binary but I want an installation that > creates a /usr/bin/scilab file so that the command 'scilab' launches > the Scilab GUI. I have tried creating a symbolic link (via sudo ln > -s /path/to/scilab-5.5.2/bin /usr/bin/scilab) to allow this but for > whatever reason running 'scilab' from the command line won't seem to > launch the Scilab GUI after I create such a link. > > My build attempt may be summarized by this bash script (before this > script was when I downloaded the source tarball and the prerequisite > tarball to Downloads using Chrome): > cd ~/Downloads for i in *.tar.gz do sudo tar xvzf $i -C /opt; > done sudo chmod 777 -R /opt/scilab-5.5.2 cd /opt/scilab-5.5.2 > sudo ./configure sudo make all > ? > no errors were encountered up until the last line, the last six lines > of output were: > Found dependency from scilab-cli-bin to libstdc++ found > Makefile:2092: recipe for target 'check-libstdcpp-dep' failed > make[1]: *** [check-libstdcpp-dep] Error 2 make[1]: Leaving > directory '/opt/scilab-5.5.2' Makefile:1541: recipe for target > 'all-recursive' failed make: *** [all-recursive] Error 1 > ? > If it is important I have been trying to follow the "Compiling Scilab > 5.x under GNU-Linux Unix" Wiki article and to get any prerequisites > not in the tarball I ran: > sudo apt-get build-dep scilab > > ? > The version of Scilab found in the vivid repos is 5.5.1 so I thought > that the dependencies for 5.5.2 would be similar enough to allow > Scilab 5.5.2 to be built using 5.5.1's dependencies. > > Thanks for your time, > Brenton > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From tim at wescottdesign.com Thu Jun 4 07:30:27 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Wed, 03 Jun 2015 22:30:27 -0700 Subject: [Scilab-users] New computer, graphics(?) crashes In-Reply-To: References: <1432682395.30072.2.camel@servo> Message-ID: <1433395827.11119.1.camel@Servo> I just installed Xubuntu, 14.04, 64 bit. Even without the fglrx drivers, it's faster and it doesn't seem to crash (at least, it's 570 seconds into a simulation, with tons of updates to the monitor graph, without problem). Thank you all. On Wed, 2015-05-27 at 08:15 +0300, Nikolay Strelkov wrote: > Dear Tim! > > For me it seems that you need to use proprietary AMD drivers - fglrx. > Open drivers are buggy and show low-performance. > > If you want to use open driver on Ubuntu 12.04 LTS (precise) you > should upgrade its HWE to 14.04 LTS (trusty). > > With best regards, > Nikolay. > > > 27 ??? 2015 ?. 3:22 ???????????? "Tim Wescott" > ???????: > Running a Scilab algorithm that puts up a lot of graphs, > Scilab crashes > with the report below. It's only happened with Scilab, but it > appears > to be the video driver that's complaining, so I don't know if > this is a > problem with Linux or with Scilab. The processor is an AMD > A10 with > on-board video processing which is considerably different from > my > 12-year-old Dell that this computer is replacing. > > Does anyone have clues for the clueless? Particularly > suggestions for > isolating, fixing, or working around the problem? > > radeon: Failed to allocate a buffer: > radeon: size : 3328 bytes > radeon: alignment : 256 bytes > radeon: domains : 2 > EE ../../../../../../src/gallium/drivers/radeon/r600_texture.c:1012 > r600_texture_transfer_map - failed to create temporary texture > to hold > untiled copy > radeon: Failed to allocate a buffer: > radeon: size : 1280 bytes > radeon: alignment : 256 bytes > radeon: domains : 4 > radeon: Failed to allocate a buffer: > radeon: size : 1 bytes > radeon: alignment : 1 bytes > radeon: domains : 2 > radeon: Failed to allocate a buffer: > radeon: size : 1 bytes > radeon: alignment : 1 bytes > radeon: domains : 2 > Segmentation fault (core dumped) > > > -- > > Tim Wescott > www.wescottdesign.com > Control & Communications systems, circuit & software design. > Phone: 503.631.7815 > Cell: 503.349.8432 > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From fujimoto2005 at gmail.com Fri Jun 5 12:36:58 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Fri, 5 Jun 2015 03:36:58 -0700 (MST) Subject: [Scilab-users] for-loop does't work Message-ID: <1433500618236-4032382.post@n3.nabble.com> I run the following script. It works at 4th line But It failes from 5th line to 7th line. It seems to be right Where is this script wrong? ******************* A=zeros(151,3); x=[0.0 5.0 10.0 15.0 19.0] i=5 A(1,:)=x(i) for i=1:5 A(1,:)=x(i)? end -- View this message in context: http://mailinglists.scilab.org/for-loop-does-t-work-tp4032382.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.carrico at free.fr Fri Jun 5 13:40:20 2015 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Fri, 5 Jun 2015 13:40:20 +0200 (CEST) Subject: [Scilab-users] =?utf-8?q?Re=C2=A0=3A__for-loop_does=27t_work?= In-Reply-To: <1433500618236-4032382.post@n3.nabble.com> Message-ID: <97025184.510425714.1433504420188.JavaMail.root@zimbra5-e1.priv.proxad.net> hi Maybe you should first explain what you would like to do ... A is a matrix 151x3 x is a matrix 1x5 I do not understand the prupose of the loop???? Depending on what you are expecting for, a good solution will be found Paul ----- Mail d'origine ----- De: fujimoto2005 ?: users at lists.scilab.org Envoy?: Fri, 05 Jun 2015 12:36:58 +0200 (CEST) Objet: [Scilab-users] for-loop does't work I run the following script. It works at 4th line But It failes from 5th line to 7th line. It seems to be right Where is this script wrong? ******************* A=zeros(151,3); x=[0.0 5.0 10.0 15.0 19.0] i=5 A(1,:)=x(i) for i=1:5 A(1,:)=x(i)? end -- View this message in context: http://mailinglists.scilab.org/for-loop-does-t-work-tp4032382.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users From fujimoto2005 at gmail.com Fri Jun 5 14:01:31 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Fri, 5 Jun 2015 05:01:31 -0700 (MST) Subject: [Scilab-users] for-loop does't work In-Reply-To: <1433500618236-4032382.post@n3.nabble.com> References: <1433500618236-4032382.post@n3.nabble.com> Message-ID: <1433505691505-4032384.post@n3.nabble.com> It is a part of such as the following script. A=zeros(151,3); x=[0.0 5.0 10.0 15.0 19.0] i=5 A(1,:)=x(i) B=[] for i=1:5 A(1,:)=x(i)? A(2:$,:) =f(i) B=[B,A] end But I found the cause of the error. '?'?after x(i) is mistyped as two bite character '?' instead of ';'. -- View this message in context: http://mailinglists.scilab.org/for-loop-does-t-work-tp4032382p4032384.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From pierre-aime.agnel at scilab-enterprises.com Fri Jun 5 14:02:39 2015 From: pierre-aime.agnel at scilab-enterprises.com (=?UTF-8?B?UGllcnJlLUFpbcOpIEFnbmVs?=) Date: Fri, 05 Jun 2015 14:02:39 +0200 Subject: [Scilab-users] =?utf-8?q?Re=C2=A0=3A__for-loop_does=27t_work?= In-Reply-To: <97025184.510425714.1433504420188.JavaMail.root@zimbra5-e1.priv.proxad.net> References: <97025184.510425714.1433504420188.JavaMail.root@zimbra5-e1.priv.proxad.net> Message-ID: <55718FDF.9090000@scilab-enterprises.com> Hi, I believe the problem is on the encoding of the text A=zeros(151,3); x=[0.0 5.0 10.0 15.0 19.0] i=5 A(1,:)=x(i) for i=1:5 A(1,:)=x(i); end This worked for me without any problems compare ascii("A(1,:)=x(i)?") // Text copied from your original mail and ascii("A(1,:)=x(i);") // Text copied from my console Best, Le 05/06/2015 13:40, paul.carrico at free.fr a ?crit : > hi > > Maybe you should first explain what you would like to do ... > > A is a matrix 151x3 > x is a matrix 1x5 > > I do not understand the prupose of the loop???? > > Depending on what you are expecting for, a good solution will be found > > Paul > > ----- Mail d'origine ----- > De: fujimoto2005 > ?: users at lists.scilab.org > Envoy?: Fri, 05 Jun 2015 12:36:58 +0200 (CEST) > Objet: [Scilab-users] for-loop does't work > > I run the following script. > It works at 4th line > But It failes from 5th line to 7th line. > It seems to be right > Where is this script wrong? > > ******************* > A=zeros(151,3); > x=[0.0 5.0 10.0 15.0 19.0] > i=5 > A(1,:)=x(i) > for i=1:5 > A(1,:)=x(i)? > end > > > > -- > View this message in context: http://mailinglists.scilab.org/for-loop-does-t-work-tp4032382.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Pierre-Aim? Agnel R&D Projects Manager Phone: +33.1.80.77.04.67 Mobile: +33.6.82.49.35.23 ----------------------------------------------------------- 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 amonmayr at laas.fr Fri Jun 5 23:17:51 2015 From: amonmayr at laas.fr (Antoine Monmayrant) Date: Fri, 05 Jun 2015 23:17:51 +0200 Subject: [Scilab-users] =?utf-8?q?Re=C2=A0=3A__for-loop_does=27t_work?= In-Reply-To: <55718FDF.9090000@scilab-enterprises.com> Message-ID: <4784-55721200-7-2ea5a0c0@133674753> Hi, I think Pierre Aim? is right about a possible issue with encoding: your ";" looks like a semi-colon but seems to be a completely different beast on my computer: ... ascii("A(1,:)=x(i)???") // Text copied from your original mail and ascii("A(1,:)=x(i);") // Text copied from my console ... "???" is definitely not a ";" Cheers, Antoine Le Vendredi 5 Juin 2015 14:02 CEST, Pierre-Aim? Agnel a ?crit: > Hi, > > I believe the problem is on the encoding of the text > > A=zeros(151,3); > x=[0.0 5.0 10.0 15.0 19.0] > i=5 > A(1,:)=x(i) > for i=1:5 > A(1,:)=x(i); > end > > This worked for me without any problems > compare > > ascii("A(1,:)=x(i)?") // Text copied from your original mail > and > ascii("A(1,:)=x(i);") // Text copied from my console > > Best, > > Le 05/06/2015 13:40, paul.carrico at free.fr a ?crit : > > hi > > > > Maybe you should first explain what you would like to do ... > > > > A is a matrix 151x3 > > x is a matrix 1x5 > > > > I do not understand the prupose of the loop???? > > > > Depending on what you are expecting for, a good solution will be found > > > > Paul > > > > ----- Mail d'origine ----- > > De: fujimoto2005 > > ?: users at lists.scilab.org > > Envoy?: Fri, 05 Jun 2015 12:36:58 +0200 (CEST) > > Objet: [Scilab-users] for-loop does't work > > > > I run the following script. > > It works at 4th line > > But It failes from 5th line to 7th line. > > It seems to be right > > Where is this script wrong? > > > > ******************* > > A=zeros(151,3); > > x=[0.0 5.0 10.0 15.0 19.0] > > i=5 > > A(1,:)=x(i) > > for i=1:5 > > A(1,:)=x(i)? > > end > > > > > > > > -- > > View this message in context: http://mailinglists.scilab.org/for-loop-does-t-work-tp4032382.html > > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > http://lists.scilab.org/mailman/listinfo/users > > -- > Pierre-Aim? Agnel > R&D Projects Manager > Phone: +33.1.80.77.04.67 > Mobile: +33.6.82.49.35.23 > ----------------------------------------------------------- > Scilab Enterprises > 143bis rue Yves Le Coz - 78000 Versailles, France > Phone: +33.1.80.77.04.68 > http://www.scilab-enterprises.com > From fujimoto2005 at gmail.com Sat Jun 6 00:49:55 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Fri, 5 Jun 2015 15:49:55 -0700 (MST) Subject: [Scilab-users] =?utf-8?q?Re=C2=A0=3A__for-loop_does=27t_work?= In-Reply-To: <4784-55721200-7-2ea5a0c0@133674753> References: <1433500618236-4032382.post@n3.nabble.com> <97025184.510425714.1433504420188.JavaMail.root@zimbra5-e1.priv.proxad.net> <55718FDF.9090000@scilab-enterprises.com> <4784-55721200-7-2ea5a0c0@133674753> Message-ID: <1433544595196-4032393.post@n3.nabble.com> Hi,Pierre-Aim? Agne and Antoine Monmayrant Yes, you are right. That's exactly due to my mistype ';' The scinote is difficult for me to see whethe IME mode is on or off. I did't have such kind of the experience when I used MATLAB or VISUAL BASIC. Is there any other editor which you recommend? Best regards -- View this message in context: http://mailinglists.scilab.org/for-loop-does-t-work-tp4032382p4032393.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From fujimoto2005 at gmail.com Sun Jun 7 10:14:02 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Sun, 7 Jun 2015 01:14:02 -0700 (MST) Subject: [Scilab-users] rename many variables Message-ID: <1433664842623-4032394.post@n3.nabble.com> When there are so many variable to be rebamed, finding an old one and rename it to a new name one by one is easy to mistake. If there is any methods to read a list of old names with corresponding new names and rename automatically, please teach me. Best regards -- View this message in context: http://mailinglists.scilab.org/rename-many-variables-tp4032394.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Sun Jun 7 12:44:05 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 07 Jun 2015 12:44:05 +0200 Subject: [Scilab-users] rename many variables In-Reply-To: <1433664842623-4032394.post@n3.nabble.com> References: <1433664842623-4032394.post@n3.nabble.com> Message-ID: <55742075.5000807@free.fr> Le 07/06/2015 10:14, fujimoto2005 a ?crit : > When there are so many variable to be rebamed, finding an old one and rename > it to a new name one by one is easy to mistake. > If there is any methods to read a list of old names with corresponding new > names and rename automatically, please teach me. According to which criterium do you decide/know that a name is an old or a new one? From fujimoto2005 at gmail.com Sun Jun 7 13:35:42 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Sun, 7 Jun 2015 04:35:42 -0700 (MST) Subject: [Scilab-users] rename many variables In-Reply-To: <55742075.5000807@free.fr> References: <1433664842623-4032394.post@n3.nabble.com> <55742075.5000807@free.fr> Message-ID: <1433676942906-4032396.post@n3.nabble.com> It's on refactoring the existing script. Some of the variable name are not good enough in the sense that I do not understand quickly when I see them. I am sure that I wil feel it difficult to understand my script quickly after 6 month. So I want to change those variable-names into more understandable names. Old(existing) names are those I want to change. New names are those I feel understandable quickly. old names new names a1 b1 a2 b2 ????????????? an bn I want to rename 'ai' into 'bi'. Best regards -- View this message in context: http://mailinglists.scilab.org/rename-many-variables-tp4032394p4032396.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From adelson.oliveira at gmail.com Sun Jun 7 16:26:31 2015 From: adelson.oliveira at gmail.com (Adelson) Date: Sun, 7 Jun 2015 07:26:31 -0700 (MST) Subject: [Scilab-users] interp_color_vector Message-ID: <1433687191704-4032397.post@n3.nabble.com> Could someone please post here a simple and clear example of plotting a curve using a specific interp_color_vector? -- View this message in context: http://mailinglists.scilab.org/interp-color-vector-tp4032397.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From jrafaelbguerra at hotmail.com Sun Jun 7 18:07:53 2015 From: jrafaelbguerra at hotmail.com (Rafael Guera) Date: Sun, 7 Jun 2015 17:07:53 +0100 Subject: [Scilab-users] Scinotes special keywords Message-ID: Dear Scilabers, Do you know if there any other keywords pre-defined in Scinotes which change the colours of the comments? (similar to the behaviour of keyword " TODO: " which turns comments into red colour). Thanks and regards, Rafael From calixte.denizet at scilab-enterprises.com Sun Jun 7 18:55:50 2015 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Sun, 07 Jun 2015 18:55:50 +0200 Subject: [Scilab-users] Scinotes special keywords In-Reply-To: References: Message-ID: <55747796.6010308@scilab-enterprises.com> Hi, copy/paste these comments: // http://www.scilab.org // TODO: a todo comment // mailto:foo.bar at scilab-enterprises.com // $$\sum_{n=1}^{+\infty}\frac1{n^2} = \frac{\pi^2}6$$ you can put your mouse over the LaTeX expression to visualize it. Regards Calixte On 07/06/2015 18:07, Rafael Guera wrote: > Dear Scilabers, > > Do you know if there any other keywords pre-defined in Scinotes which change the > colours of the comments? > (similar to the behaviour of keyword " TODO: " which turns comments into red > colour). > > Thanks and regards, > > Rafael > > > _______________________________________________ > 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 sgougeon at free.fr Mon Jun 8 04:43:31 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 08 Jun 2015 04:43:31 +0200 Subject: [Scilab-users] interp_color_vector In-Reply-To: <1433687191704-4032397.post@n3.nabble.com> References: <1433687191704-4032397.post@n3.nabble.com> Message-ID: <55750153.702@free.fr> Le 07/06/2015 16:26, Adelson a ?crit : > Could someone please post here a simple and clear example of plotting a curve > using a specific interp_color_vector? Here it is: clf drawlater f = gcf(); Nc = 100; // Number of colors used f.color_map = jetcolormap(Nc); x = [0 0.6 1 0.2]; y = [0 0.4 1 0.7]; xpoly(x,y, "lines",1) e = gce(); e.fill_mode = "on"; // the background must be filled e.interp_color_vector = [Nc Nc 1 1]; // defines distibution of graded colors e.interp_color_mode = "on"; // actives graded colors drawnow The result of a more elaborated example (could be added as a demo) : -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dcehgfcj.png Type: image/png Size: 50501 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dajfgici.png Type: image/png Size: 9512 bytes Desc: not available URL: From herve.bellot at irstea.fr Thu Jun 4 16:47:30 2015 From: herve.bellot at irstea.fr (Herve B.) Date: Thu, 4 Jun 2015 07:47:30 -0700 (MST) Subject: [Scilab-users] mget 64 32 bit Message-ID: <1433429250648-4032379.post@n3.nabble.com> I used Scilab 5.3.3 to open digital data format lb using the "mget". But I can not open them with Scilab 5.5.2. There was apparently made changes of this function has to be compatible with 64-bit data but I do not understand how to fix it and open my data seems to be 32 bit. thank you Herv? -- View this message in context: http://mailinglists.scilab.org/mget-64-32-bit-tp4032379.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From Christophe.Dang at sidel.com Mon Jun 8 12:55:46 2015 From: Christophe.Dang at sidel.com (Dang, Christophe) Date: Mon, 8 Jun 2015 10:55:46 +0000 Subject: [Scilab-users] =?utf-8?q?Re=C2=A0=3A__for-loop_does=27t_work?= In-Reply-To: <1433544595196-4032393.post@n3.nabble.com> References: <1433500618236-4032382.post@n3.nabble.com> <97025184.510425714.1433504420188.JavaMail.root@zimbra5-e1.priv.proxad.net> <55718FDF.9090000@scilab-enterprises.com> <4784-55721200-7-2ea5a0c0@133674753> <1433544595196-4032393.post@n3.nabble.com> Message-ID: Hello, > De : fujimoto2005 > Envoy? : samedi 6 juin 2015 00:50 > > The scinote is difficult for me to see whethe IME mode is on or off. Maybe you could report a bug in Bugzilla to ask such an improvement http://bugzilla.scilab.org/ 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 adelson.oliveira at gmail.com Mon Jun 8 18:56:32 2015 From: adelson.oliveira at gmail.com (Adelson Oliveira) Date: Mon, 8 Jun 2015 13:56:32 -0300 Subject: [Scilab-users] interp_color_vector In-Reply-To: <55750153.702@free.fr> References: <1433687191704-4032397.post@n3.nabble.com> <55750153.702@free.fr> Message-ID: Thanks for the answer, I intend to plot a curve that changes its color. Following your code I've tried, clfdrawlaterf = gcf();Nc = 100; // Number of colors usedf.color_map = jetcolormap(Nc); x=[0:0.1:4*%pi]; y=sin(x); plot2d(x,y,axesflag=0) e = gce();e.fill_mode = "on"; // the background must be fillede.interp_color_vector = [Nc Nc 1 1]; // defines distibution of graded colorse.interp_color_mode = "on"; // actives graded colorsdrawnow And got errors like, "... the property 'fill_mode' or 'data' does not exist for this handler ..." (a tranlation from portuguese) Is it possible to have a curve changing color as it proceeds? Thanks again 2015-06-07 23:43 GMT-03:00 Samuel Gougeon : > Le 07/06/2015 16:26, Adelson a ?crit : > > Could someone please post here a simple and clear example of plotting a curve > using a specific interp_color_vector? > > Here it is: > > clfdrawlaterf = gcf();Nc = 100; // Number of colors usedf.color_map = jetcolormap(Nc);x = [0 0.6 1 0.2];y = [0 0.4 1 0.7];xpoly(x,y, "lines",1)e = gce();e.fill_mode = "on"; // the background must be fillede.interp_color_vector = [Nc Nc 1 1]; // defines distibution of graded colorse.interp_color_mode = "on"; // actives graded colorsdrawnow > > > The result of a more elaborated example (could be added as a demo) : > > > > _______________________________________________ > 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: dcehgfcj.png Type: image/png Size: 50501 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dajfgici.png Type: image/png Size: 9512 bytes Desc: not available URL: From Serge.Steer at inria.fr Mon Jun 8 18:56:59 2015 From: Serge.Steer at inria.fr (Serge Steer) Date: Mon, 08 Jun 2015 18:56:59 +0200 Subject: [Scilab-users] Questions about handling cursor from Scilab Message-ID: <5575C95B.2090509@inria.fr> Hello, Does anyone knows how to * modify the pointeur coordinates into a Scilab graphic window with Scilab instruction? * modify the pointeur icon with Scilab instruction? Thanks Serge Steer -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Mon Jun 8 21:45:11 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 08 Jun 2015 21:45:11 +0200 Subject: [Scilab-users] uman : help in the console, with language switch, automatic redirections, related online pages. Message-ID: <5575F0C7.3080207@free.fr> Hello, *uman* is a *U*nified display command of the Scilab User *MAN*ual. It is a general purpose function that might be useful for every Scilab user and developer. *uman* 1.2 is just released on ATOMS (web and Scilab interface).Install it and enjoy! *Pointers* : * ATOMS page : http://atoms.scilab.org/toolboxes/uman * Help (en): http://sgougeon.free.fr/scilab/help/uman/scilab_en_US_help/uman.html * Help (fr) : http://sgougeon.free.fr/scilab/help/uman/scilab_fr_FR_help/uman.html *Main features* : * A single function to see help pages from Scilab, from installed ATOMS modules (no need to load them), from other external modules, and heading comments in functions. Around 50 external references available on the Scilab FileExchange can also be addressed. * A single function to target and display contents in text mode in the console, or as usual in the help browser, or as well online (help, bugzilla, atoms, forge, or even external websites). *uman* aims to embrace as many Scilab ressources as possible. * Just give a language code en | de | fr | ja | pt | ru | zh in option, and you get the right version of the page in the console or online. No need to change the session's language. Watching the reference en_US version of the page is now straightforward, without leaving your locales. * Just display main infos in brief: path in help, calling sequences, parameters, See also. And be free to display more on request: description, examples, history, parent ToC... Or finally display the whole page with the "a" option. It is so easy to copy/past pieces of examples from the console... to the console! * Are you used to use another scientific language? Almost 200 automatic redirections will target the right or nearest Scilab equivalence for you. Other handy shortcuts are also defined for all users. * Stucked with a bug? The "wb" option lets you list documented ones related to the faulty entry, for Scilab and many ATOMS packages. Online users comments have never been so easily reachable. * Last but not least: get a smart text layout in console. Nested lists of items, tables, pieces of codes, notes and warnings, text styles.. are supported. The code style of calling sequences and examples is also improved. *Supported OS* : Windows, Linux and Mac OS, with Scilab 5.5 and Scilab 6. *uman* has been extensively tested for Windows, a few tested for Linux, and is almost untested for Mac OS. Issues may be reported in comments on the ATOMS page. Please have a look to the list of known limitations before posting. Try *uman*, and adopt it ;) Best regards Samuel Gougeon -------------- next part -------------- An HTML attachment was scrubbed... URL: From antoine.monmayrant at laas.fr Tue Jun 9 09:27:29 2015 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Tue, 09 Jun 2015 09:27:29 +0200 Subject: [Scilab-users] uman : help in the console, with language switch, automatic redirections, related online pages. In-Reply-To: <5575F0C7.3080207@free.fr> References: <5575F0C7.3080207@free.fr> Message-ID: <55769561.8050605@laas.fr> I just started to play with it and it looks good! It provides a feature that scilab is lacking: a decent command line documentation. Thanks, Antoine Le 06/08/2015 09:45 PM, Samuel Gougeon a ?crit : > Hello, > > *uman* is a *U*nified display command of the Scilab User *MAN*ual. It > is a general purpose > function that might be useful for every Scilab user and developer. > *uman* 1.2 is just released on ATOMS (web and Scilab > interface).Install it and enjoy! > > *Pointers* : > * ATOMS page : http://atoms.scilab.org/toolboxes/uman > * Help (en): > http://sgougeon.free.fr/scilab/help/uman/scilab_en_US_help/uman.html > * Help (fr) : > http://sgougeon.free.fr/scilab/help/uman/scilab_fr_FR_help/uman.html > > *Main features* : > > * A single function to see help pages from Scilab, from installed > ATOMS modules (no need > to load them), from other external modules, and heading comments in > functions. Around > 50 external references available on the Scilab FileExchange can > also be addressed. > > * A single function to target and display contents in text mode in > the console, or as usual > in the help browser, or as well online (help, bugzilla, atoms, > forge, or even external websites). > *uman* aims to embrace as many Scilab ressources as possible. > > * Just give a language code en | de | fr | ja | pt | ru | zh in > option, and you get the > right version of the page in the console or online. No need to > change the session's > language. Watching the reference en_US version of the page is now > straightforward, > without leaving your locales. > > * Just display main infos in brief: path in help, calling sequences, > parameters, See > also. And be free to display more on request: description, > examples, history, parent > ToC... Or finally display the whole page with the "a" option. It is > so easy to copy/past > pieces of examples from the console... to the console! > > * Are you used to use another scientific language? Almost 200 > automatic redirections will > target the right or nearest Scilab equivalence for you. Other handy > shortcuts are also > defined for all users. > > * Stucked with a bug? The "wb" option lets you list documented ones > related to the > faulty entry, for Scilab and many ATOMS packages. Online users > comments have > never been so easily reachable. > > * Last but not least: get a smart text layout in console. Nested > lists of items, tables, > pieces of codes, notes and warnings, text styles.. are supported. > The code style of > calling sequences and examples is also improved. > > *Supported OS* : Windows, Linux and Mac OS, with Scilab 5.5 and Scilab 6. > > *uman* has been extensively tested for Windows, a few tested for > Linux, and is > almost untested for Mac OS. Issues may be reported in comments on > the ATOMS > page. Please have a look to the list of known limitations before > posting. > > Try *uman*, and adopt it ;) > > Best regards > Samuel Gougeon > > > > _______________________________________________ > 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 Tue Jun 9 12:41:59 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 09 Jun 2015 12:41:59 +0200 Subject: [Scilab-users] interp_color_vector In-Reply-To: References: <1433687191704-4032397.post@n3.nabble.com> <55750153.702@free.fr> Message-ID: <5576C2F7.4050604@free.fr> Hi, Le 08/06/2015 18:56, Adelson Oliveira a ?crit : > ../.. > Is it possible to have a curve changing color as it proceeds? Are you talking about the line color or the filling color (of the underlying polygon)? Is the same varying color applied to the whole curve, or are they applied segment per segment? What do you mean by "as it proceed"? is the curve moving? is it lengthening with previous segments fixed? is it reshaping? is it in 3D with nodes varying in z? SG -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Tue Jun 9 13:13:43 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 09 Jun 2015 13:13:43 +0200 Subject: [Scilab-users] Getting names of attributes of an XML tag In-Reply-To: <556B7FA8.6010501@free.fr> References: <556B7FA8.6010501@free.fr> Message-ID: <5576CA67.1090002@free.fr> Le 31/05/2015 23:39, Samuel Gougeon a ?crit : > Hello, > > I am trying to get the list of _names_ of attributes of an XML tag. > Getting their values is OK, but i fail catching related names! > Is this possible? > Yes! I finally got them: > s = ""+.. > ""+.. > "Hello world" > doc = xmlReadStr (s); > Then: > *a= xmlXPath**(**doc.root.children(1)**,** **"@*"**)* and finally: > > -->*a.name* > ans = > > !att rib ! > as well as > > -->*a.content* > ans = > > !foo bar ! . I would have expected that .name and .content be usable on e, but previous trials showed that it is/was not the case . > -->e = doc.root.children(1).attributes > e = > > XML Attributes > att --> foo > rib --> bar > > -->e.name > ans = > [] > > -->e.content > ans = > [] It would be quite intuitive and more straightforward. Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Tue Jun 9 13:32:08 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 09 Jun 2015 13:32:08 +0200 Subject: [Scilab-users] Getting names of attributes of an XML tag In-Reply-To: <5576CA67.1090002@free.fr> References: <556B7FA8.6010501@free.fr> <5576CA67.1090002@free.fr> Message-ID: <5576CEB8.5040200@free.fr> Le 09/06/2015 13:13, Samuel Gougeon a ?crit : > .../... > I would have expected that .name and .content be usable on e, but > previous trials showed that it is/was not the case > . >> -->e = doc.root.children(1).attributes >> e = >> >> XML Attributes >> att --> foo >> rib --> bar >> >> -->e.name >> ans = >> [] >> >> -->e.content >> ans = >> [] > It would be quite intuitive and more straightforward. . OK, i have just caught why this is not implemented as this: e being an XML attribute object, e.mystring points to the value of the attribute named myString. So, e.name and e.content point to the value of attributes named "name" and "content", as in bla bla I don't know well XML conventions for attributes names. I am pretty sure that attributes names should exclude some special chars like "&" or some others. So, addressing the names could be done through e.name& and values through e.value& ... without having to use xmlXPath(). A proposal. Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From calixte.denizet at scilab-enterprises.com Tue Jun 9 14:00:18 2015 From: calixte.denizet at scilab-enterprises.com (Calixte Denizet) Date: Tue, 09 Jun 2015 14:00:18 +0200 Subject: [Scilab-users] Getting names of attributes of an XML tag In-Reply-To: <5576CEB8.5040200@free.fr> References: <556B7FA8.6010501@free.fr> <5576CA67.1090002@free.fr> <5576CEB8.5040200@free.fr> Message-ID: <5576D552.5090809@scilab-enterprises.com> On 09/06/2015 13:32, Samuel Gougeon wrote: > Le 09/06/2015 13:13, Samuel Gougeon a ?crit : >> .../... >> I would have expected that .name and .content be usable on e, but >> previous trials showed that it is/was not the case >> . >>> -->e = doc.root.children(1).attributes >>> e = >>> >>> XML Attributes >>> att --> foo >>> rib --> bar >>> >>> -->e.name >>> ans = >>> [] >>> >>> -->e.content >>> ans = >>> [] >> It would be quite intuitive and more straightforward. > . > OK, i have just caught why this is not implemented as this: > e being an XML attribute object, e.mystring points to the value of > the attribute named myString. So, e.name and e.content point to the > value of attributes named "name" and "content", as in name="test" content="value">bla bla > You're right ! > I don't know well XML conventions for attributes names. I am pretty > sure that attributes names should exclude some special chars like "&" > or some others. So, addressing the names could be done through e.name& > and values through e.value& ... without having to use xmlXPath(). A > proposal. Bof... Maybe xmlName is what you expect: http://help.scilab.org/docs/5.5.2/en_US/xmlName.html Calixte > > Samuel > > > > _______________________________________________ > 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 sgougeon at free.fr Tue Jun 9 14:19:57 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 09 Jun 2015 14:19:57 +0200 Subject: [Scilab-users] Getting names of attributes of an XML tag In-Reply-To: <5576D552.5090809@scilab-enterprises.com> References: <556B7FA8.6010501@free.fr> <5576CA67.1090002@free.fr> <5576CEB8.5040200@free.fr> <5576D552.5090809@scilab-enterprises.com> Message-ID: <5576D9ED.1010409@free.fr> Le 09/06/2015 14:00, Calixte Denizet a ?crit : > .../... > Maybe xmlName is what you expect: > http://help.scilab.org/docs/5.5.2/en_US/xmlName.html Exactly! The right function's name, for the proper goal.. and i missed it! :( Thanks -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Tue Jun 9 14:44:52 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 09 Jun 2015 14:44:52 +0200 Subject: [Scilab-users] xsegs() <= Re: interp_color_vector In-Reply-To: References: <1433687191704-4032397.post@n3.nabble.com> <55750153.702@free.fr> Message-ID: <5576DFC4.2090000@free.fr> Le 08/06/2015 18:56, Adelson Oliveira a ?crit : > Thanks for the answer, > > I intend to plot a curve that changes its color. Following your code > I've tried, . Here is an example with a 1D curve laying in the 2D plane, using xsegs(): f = gcf(); f.color_map = jetcolormap(100); // with Scilab 6, gcf().colormap = jetcolormap(100); // will work. Cool :) clf x = linspace(0,10,101); y = sin(x); xs = [ x(1:$-1) ; x(2:$)]; ys = [ y(1:$-1) ; y(2:$)]; xsegs(xs,ys,1:100) e = gce(); e.thickness = 5; You may customize it. Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: ajeahfbc.png Type: image/png Size: 2606 bytes Desc: not available URL: From clement.david at scilab-enterprises.com Tue Jun 9 15:29:18 2015 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Tue, 09 Jun 2015 15:29:18 +0200 Subject: [Scilab-users] FROMWSB Scilab In-Reply-To: <1432699815234-4032337.post@n3.nabble.com> References: <1432699815234-4032337.post@n3.nabble.com> Message-ID: <1433856558.4859.27.camel@scilab-enterprises.com> Hello, The parameters are fixed at diagram compilation time so you cannot modify them at simulation time however advanced usage of scicos_simulate can simulate this case (passing and modifying the Info list). -- Cl?ment Le mardi 26 mai 2015 ? 21:10 -0700, ANDREY_PSTU a ?crit : > Hello, guys. I have a problem with FROMWSB block. I simulate a simple > model > which use a workspace parameter (V.time and V.values) and draw a > simple > graph.When i use scicosim function i want to change a parameters of > FROMWSB > block (new V.time and V.values) and continue simulation. But i don't > know > how to change a parameters of this block. Is it possible? Please, I > need > your help. Thank you. > > > > -- > View this message in context: > http://mailinglists.scilab.org/FROMWSB-Scilab-tp4032337.html > Sent from the Scilab users - Mailing Lists Archives mailing list > archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From clement.david at scilab-enterprises.com Tue Jun 9 15:54:13 2015 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Tue, 09 Jun 2015 15:54:13 +0200 Subject: [Scilab-users] Questions about handling cursor from Scilab In-Reply-To: <5575C95B.2090509@inria.fr> References: <5575C95B.2090509@inria.fr> Message-ID: <1433858053.4859.30.camel@scilab-enterprises.com> Hello Serge, Out of any Java specific code, it might probably be impossible. -- Cl?ment Le lundi 08 juin 2015 ? 18:56 +0200, Serge Steer a ?crit : > Hello, > > Does anyone knows how to > modify the pointeur coordinates into a Scilab graphic window with > Scilab instruction? > modify the pointeur icon with Scilab instruction? > Thanks > Serge Steer > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From paul.carrico at esterline.com Tue Jun 9 16:36:17 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Tue, 9 Jun 2015 14:36:17 +0000 Subject: [Scilab-users] advice on a pseudo-"bandwidth" calculation Message-ID: <3A6B7233274DB449A2A0053A47684F952B7A7121@BGS-EX01.auxitrol.ad> Dear all I'm currently thinking in a way to compare experimental have a bell shape, composed of points and not connected to any parametric curve (the bell shape assumption could answer to the need in a first step). A way I can imagine is to calculate a kind of " bandwidth" (-3 db). I add a look in the Scilab help but nothing obvious appears what is the best way to proceed to: - calculate/determine the intersection points with the "pseudo"-curves, - the number of points may change and may have different abscissa ? - the points could not necessarily by expressed by a parametric curve (to fit on before calculating the intersection), Any suggestion are welcome NB: please find hereafter a basic example using a gauss curve (just to illustrate the purpose) Thanks Paul ###################################################################### mode(0) // gauss curve coef = 10; x = [-6 : 0.1 : 6]'; [nl,nc]=size(x); y1 = coef*(1/sqrt(2*%pi))*exp(-0.5*x.^2); // -3db line top_ = max(y1); three_db = top_/sqrt(2); y2=three_db*ones(nl,1); plot2d(x,y1); // gausse curve plot2d(x,y2); // -3db line EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Tue Jun 9 20:04:49 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Tue, 09 Jun 2015 20:04:49 +0200 Subject: [Scilab-users] advice on a pseudo-"bandwidth" calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7A7121@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7A7121@BGS-EX01.auxitrol.ad> Message-ID: <55772AC1.6000107@free.fr> Hi Paul, I am not sure to well understand the purpose and constrains. Why not fitting experimental points with a model? -- yes, the gaussian law could be a good start, or a lorentzian one, or other, depending on the type of observed phenomena. Then, you would get a set of modeling parameters (max amplitude, mean position, standard deviation), and you could compare fitting parameters from different sets of points, among which the dispersion of points. This would be a standard way to proceed. Samuel From tim at wescottdesign.com Tue Jun 9 20:30:08 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Tue, 09 Jun 2015 11:30:08 -0700 Subject: [Scilab-users] advice on a pseudo-"bandwidth" calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7A7121@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7A7121@BGS-EX01.auxitrol.ad> Message-ID: <1433874608.2388.62.camel@Servo> On Tue, 2015-06-09 at 14:36 +0000, Carrico, Paul wrote: > Dear all > > > > I?m currently thinking in a way to compare experimental have a bell > shape, composed of points and not connected to any parametric curve > (the bell shape assumption could answer to the need in a first step). > > > > A way I can imagine is to calculate a kind of ? bandwidth? (-3 db). > > > > I add a look in the Scilab help but nothing obvious appears > > > > what is the best way to proceed to: > > - calculate/determine the intersection points with the > ?pseudo?-curves, > > - the number of points may change and may have different abscissa ? > > - the points could not necessarily by expressed by a parametric curve > (to fit on before calculating the intersection), > > > > Any suggestion are welcome > > > > > > NB: please find hereafter a basic example using a gauss curve (just to > illustrate the purpose) > To the extent that you can do so without giving away any proprietary information, could you please tell us what you're really trying to do? Citing a Gaussian distribution (a "bell curve") implies that you're trying to find a probability density, but citing a "bandwidth" implies that you're trying to characterize some sort of a signal processing system. Those two concepts don't exactly fit well together when you're talking of one set of data. If you're taking a bunch of measurements and you're assuming a Gaussian distribution of errors, and you want to estimate the mean and variance of the distribution, then you want to use mean(x) and var(x). This is kind of the quick & dirty way to do things, but with a large enough sample it's not going to be very far off the mark (if you have less than a couple of dozen samples you may want to investigate the Student t distribution, but I don't like recommending it if you don't have a solid background in statistics). -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From adelson.oliveira at gmail.com Wed Jun 10 04:34:59 2015 From: adelson.oliveira at gmail.com (Adelson) Date: Tue, 9 Jun 2015 19:34:59 -0700 (MST) Subject: [Scilab-users] xsegs() <= Re: interp_color_vector In-Reply-To: <5576DFC4.2090000@free.fr> References: <1433687191704-4032397.post@n3.nabble.com> <55750153.702@free.fr> <5576DFC4.2090000@free.fr> Message-ID: <1433903699373-4032421.post@n3.nabble.com> Thanks Samuel, I had already tried changing colors of individual segments of a line using a trick with plot2d. This example of yours using xsegs is much better than my trick because it generates just "one element" of graphics as opposed to a bunch of polylines made with plot2d. I have just one more question: Is it possible to assign a rounded border to all line segments in the curve? Thanks again -- View this message in context: http://mailinglists.scilab.org/interp-color-vector-tp4032397p4032421.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From manasdas17 at gmail.com Wed Jun 10 07:24:17 2015 From: manasdas17 at gmail.com (manasdas17) Date: Tue, 9 Jun 2015 22:24:17 -0700 (MST) Subject: [Scilab-users] Source folder of GUI of Xcos palette browser Xcos GUI Menu bar Message-ID: <1433913857686-4032422.post@n3.nabble.com> HI, I need to make some changes to GUI of Xcos palette browser menubar and also in the menubar of Xcos editor.But i couldn't find the scilab files regarding this in the share/scilab/macros folder. GUI for Xcos is also n't find in the gui folder in macros folder.So,plz let me know where all the files are there in which the gui of Xcos palette browser and xcos editor is defined. Regards, Manas -- View this message in context: http://mailinglists.scilab.org/Source-folder-of-GUI-of-Xcos-palette-browser-Xcos-GUI-Menu-bar-tp4032422.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From clement.david at scilab-enterprises.com Wed Jun 10 09:00:54 2015 From: clement.david at scilab-enterprises.com (=?ISO-8859-1?Q?Cl=E9ment?= David) Date: Wed, 10 Jun 2015 09:00:54 +0200 Subject: [Scilab-users] Source folder of GUI of Xcos palette browser Xcos GUI Menu bar In-Reply-To: <1433913857686-4032422.post@n3.nabble.com> References: <1433913857686-4032422.post@n3.nabble.com> Message-ID: <1433919654.4643.2.camel@scilab-enterprises.com> Hello, Xcos is a java implemented module so most of the source code is in the modules/xcos/src/java directory. Especially for the palette there is a java package named 'palette' which handle the GUI. For Xcos, Scilab callbacks can be added to the 'Tools' menu but to change the others you have to take a look at the Java source code. Regards, -- Cl?ment Le mardi 09 juin 2015 ? 22:24 -0700, manasdas17 a ?crit : > HI, > I need to make some changes to GUI of Xcos palette browser > menubar and > also in the menubar of Xcos editor.But i couldn't find the scilab > files > regarding this in the share/scilab/macros folder. GUI for Xcos is > also n't > find in the gui folder in macros folder.So,plz let me know where all > the > files are there in which the gui of Xcos palette browser and xcos > editor is > defined. > > Regards, > Manas > > > > -- > View this message in context: http://mailinglists.scilab.org/Source > -folder-of-GUI-of-Xcos-palette-browser-Xcos-GUI-Menu-bar > -tp4032422.html > Sent from the Scilab users - Mailing Lists Archives mailing list > archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From paul.carrico at free.fr Wed Jun 10 23:11:03 2015 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Wed, 10 Jun 2015 23:11:03 +0200 (CEST) Subject: [Scilab-users] equivalent to eigs in scilab Message-ID: <1821059326.532022691.1433970663276.JavaMail.root@zimbra5-e1.priv.proxad.net> Dear All I get a Matlab file and I would like to know what is the equivalent to Matlab: n = 6; [u,v]=eigs(K,M,n,'SM'); // where K and U and 2 sparse matrix //'SM' means Smallest for my understanting -> looking for the 6th lower eigen values, isn't it ? Scilab: Spec or bdiag seems to be equivalent, but I failed in using them How should I to proceed ? Thanks Paul From tim at wescottdesign.com Wed Jun 10 23:45:28 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Wed, 10 Jun 2015 14:45:28 -0700 Subject: [Scilab-users] equivalent to eigs in scilab In-Reply-To: <1821059326.532022691.1433970663276.JavaMail.root@zimbra5-e1.priv.proxad.net> References: <1821059326.532022691.1433970663276.JavaMail.root@zimbra5-e1.priv.proxad.net> Message-ID: <1433972728.2388.85.camel@Servo> On Wed, 2015-06-10 at 23:11 +0200, paul.carrico at free.fr wrote: > Dear All > > I get a Matlab file and I would like to know what is the equivalent to > > Matlab: > n = 6; > [u,v]=eigs(K,M,n,'SM'); // where K and U and 2 sparse matrix > //'SM' means Smallest for my understanting -> looking for the 6th lower eigen values, isn't it ? > > Scilab: > Spec or bdiag seems to be equivalent, but I failed in using them > > How should I to proceed ? > > Thanks > > Paul spec (short for "spectrum"): -->A = rand(3,3,"n") A = 0.3501626 - 1.4061926 0.5546874 1.0478272 - 1.0384734 - 0.2143931 - 1.3218008 - 1.7350313 - 2.0735088 -->spec(A) ans = - 0.2010690 + 1.3816266i - 0.2010690 - 1.3816266i - 2.3596815 -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From pierre-aime.agnel at scilab-enterprises.com Thu Jun 11 01:13:42 2015 From: pierre-aime.agnel at scilab-enterprises.com (=?windows-1252?Q?Pierre-Aim=E9_Agnel?=) Date: Thu, 11 Jun 2015 01:13:42 +0200 Subject: [Scilab-users] equivalent to eigs in scilab In-Reply-To: <1821059326.532022691.1433970663276.JavaMail.root@zimbra5-e1.priv.proxad.net> References: <1821059326.532022691.1433970663276.JavaMail.root@zimbra5-e1.priv.proxad.net> Message-ID: <5578C4A6.2060609@scilab-enterprises.com> Hello, You can use the eigs routine from ARPACK http://help.scilab.org/docs/5.5.2/en_US/eigs.html An example A = sprand(1000, 1000, 0.02, "normal"); // creates a random sparse matrix of size 1000 x 1000 and density 0.02 [d, v] = eigs(A, [], 10, "SM"); // computes the 10 smallest eigenvalues of A Hope it helps, Le 10/06/2015 23:11, paul.carrico at free.fr a ?crit : > Dear All > > I get a Matlab file and I would like to know what is the equivalent to > > Matlab: > n = 6; > [u,v]=eigs(K,M,n,'SM'); // where K and U and 2 sparse matrix > //'SM' means Smallest for my understanting -> looking for the 6th lower eigen values, isn't it ? > > Scilab: > Spec or bdiag seems to be equivalent, but I failed in using them > > How should I to proceed ? > > Thanks > > Paul > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -- Pierre-Aim? Agnel R&D Projects Manager Phone: +33.1.80.77.04.67 Mobile: +33.6.82.49.35.23 ----------------------------------------------------------- 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 kiabhilasha08 at gmail.com Sat Jun 13 20:22:19 2015 From: kiabhilasha08 at gmail.com (kiabhilasha08) Date: Sat, 13 Jun 2015 11:22:19 -0700 (MST) Subject: [Scilab-users] Regarding compilation of scilab 5.5.2 in windows Message-ID: Sir, I am trying to compile the source code of scilab version 5.5.2 in windows 7 (64-bit) according to the steps given in http://wiki.scilab.org/Visual%20Studio%20C%2B%2B%202008%20Pro%20%26%20Intel%20Fortran%2011 but when I tried to build the scilab_f -- View this message in context: http://mailinglists.scilab.org/Regarding-compilation-of-scilab-5-5-2-in-windows-tp4032444.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From kiabhilasha08 at gmail.com Sat Jun 13 20:29:41 2015 From: kiabhilasha08 at gmail.com (kiabhilasha08) Date: Sat, 13 Jun 2015 11:29:41 -0700 (MST) Subject: [Scilab-users] Compilation of scilab 5.5.2 in windows Message-ID: Sir, I am trying to compile scilab 5.5.2 source using the instructions given in this link http://wiki.scilab.org/Visual%20Studio%20C%2B%2B%202008%20Express but when i tried to build scilab_f2c.sln , it shows 'this solution is made to build with the newer version of the application' . when I tried to build it with Visual Studio 2012 ,the error message thrown was 'Cannot find V120 toolset' , I also tried it on VS 2013 but erros were thrown. Please help me in solving the issue as i am in an urgent need of building the source code. Regards Abhilasha Sancheti #rd year Engineering student -- View this message in context: http://mailinglists.scilab.org/Compilation-of-scilab-5-5-2-in-windows-tp4032445.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From c_rabanal0 at yahoo.com.ar Sat Jun 13 01:53:30 2015 From: c_rabanal0 at yahoo.com.ar (Cristian) Date: Fri, 12 Jun 2015 16:53:30 -0700 (MST) Subject: [Scilab-users] Grocer - Markov Switching model estimation Message-ID: <1434153210108-4032440.post@n3.nabble.com> Hello everyone, I try to estimate a MS with Grocer . I have *annual data*. I have read the chapter 23 (manual), where I found an perfect example for quartely data. Even though, I have alredy a doubt. The chapter 23 gives us an example to etimate a MS-AR(3). The specification is the following: r=ms_reg('100*(lpbi-lagts(2,lpbi))',['100*(lagts(1,lpbi)-lagts(3,lpbi))';... -->'100*(lagts(2,lpbi)-lagts(4,lpbi))';'100*(lagts(3,lpbi)-lagts(5,lpbi))'],'cte',nb_states,switch_var,var_opt,'transf=stud',... -->'prt=initial;final') So, the question is: How should I modify that expression in order to estimate a MS-AR(1) with anual data? Thanks in advance Cristian -- View this message in context: http://mailinglists.scilab.org/Grocer-Markov-Switching-model-estimation-tp4032440.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From grocer.toolbox at gmail.com Mon Jun 15 11:49:28 2015 From: grocer.toolbox at gmail.com (=?utf-8?Q?=C3=89ric_Dubois?=) Date: Mon, 15 Jun 2015 11:49:28 +0200 Subject: [Scilab-users] Grocer - Markov Switching model estimation In-Reply-To: <1434153210108-4032440.post@n3.nabble.com> References: <1434153210108-4032440.post@n3.nabble.com> Message-ID: <10AD7A2D-20E1-41C2-8FC7-3CD6E308112C@gmail.com> Hello Christian. The fact that your data are annuel does not have any influence except for the definition of the bounds if your data are Grocer Time series and you do not want to use the maximum time span. If your data are called mydata and whethrer they are Grocer Time series or a vector, then you may run: --> r=ms_var('all',1,mydata,2,2,1) (Run : Help ms_var for d?tails) ?ric Envoy? de mon iPhone > Le 13 juin 2015 ? 01:53, Cristian a ?crit : > > Hello everyone, I try to estimate a MS with Grocer . I have *annual data*. I > have read the chapter 23 (manual), where I found an perfect example for > quartely data. Even though, I have alredy a doubt. > > The chapter 23 gives us an example to etimate a MS-AR(3). The specification > is the following: > > r=ms_reg('100*(lpbi-lagts(2,lpbi))',['100*(lagts(1,lpbi)-lagts(3,lpbi))';... > -->'100*(lagts(2,lpbi)-lagts(4,lpbi))';'100*(lagts(3,lpbi)-lagts(5,lpbi))'],'cte',nb_states,switch_var,var_opt,'transf=stud',... > -->'prt=initial;final') > > So, the question is: How should I modify that expression in order to > estimate a MS-AR(1) with anual data? > > Thanks in advance > > Cristian > > > > -- > View this message in context: http://mailinglists.scilab.org/Grocer-Markov-Switching-model-estimation-tp4032440.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From rei.listas at yahoo.com Tue Jun 16 02:12:12 2015 From: rei.listas at yahoo.com (Reinaldo) Date: Tue, 16 Jun 2015 00:12:12 +0000 (UTC) Subject: [Scilab-users] How can I change the x_ticks and y_ticks? Message-ID: <807388956.3533964.1434413532097.JavaMail.yahoo@mail.yahoo.com> Hi Scilab members, How could I change the x_ticks and y_ticks? For example: In this following script, the x_ticks displays each 2 units on X axis. I'd like each 0.5; and the y_ticks displays each 0.2 units on Y axis. I'd like each 0.1. How could I change those parameters? function f=myfunction(x) f = sin(x)./x endfunction // Plot scf(1); clf(1); xdata = linspace ( -12, 12, 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); a=gca(); // obtendo o manipulador dos eixos correntes a.grid=[1 1]; xlabel("X","fontsize",4,"color","red"); ylabel("Y","fontsize",4,"color","red");Thank you in advance. All best,Reinaldo. -------------- next part -------------- An HTML attachment was scrubbed... URL: From rei.listas at yahoo.com Tue Jun 16 02:23:43 2015 From: rei.listas at yahoo.com (Reinaldo) Date: Tue, 16 Jun 2015 00:23:43 +0000 (UTC) Subject: [Scilab-users] How can I change the x_ticks and y_ticks? In-Reply-To: <807388956.3533964.1434413532097.JavaMail.yahoo@mail.yahoo.com> References: <807388956.3533964.1434413532097.JavaMail.yahoo@mail.yahoo.com> Message-ID: <1513463060.3578230.1434414223738.JavaMail.yahoo@mail.yahoo.com> Hi Scilab members, How could I change the x_ticks and y_ticks? For example: In this following script, the x_ticks displays each 2 units on X axis. I'd like each 0.5; and the y_ticks displays each 0.2 units on Y axis. I'd like each 0.1. How could I change those parameters? function f=myfunction(x) f = sin(x)./x endfunction // Plot scf(1); clf(1); xdata = linspace ( -12, 12, 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); a=gca(); // obtendo o manipulador dos eixos correntes a.grid=[1 1]; xlabel("X","fontsize",4,"color","red"); ylabel("Y","fontsize",4,"color","red");Thank you in advance. All best,Reinaldo. -------------- next part -------------- An HTML attachment was scrubbed... URL: From Jean-Yves.Baudais at insa-rennes.fr Tue Jun 16 09:52:34 2015 From: Jean-Yves.Baudais at insa-rennes.fr (jbaudais@insa-rennes.fr) Date: Tue, 16 Jun 2015 09:52:34 +0200 Subject: [Scilab-users] How can I change the x_ticks and y_ticks? In-Reply-To: <807388956.3533964.1434413532097.JavaMail.yahoo@mail.yahoo.com> References: <807388956.3533964.1434413532097.JavaMail.yahoo@mail.yahoo.com> Message-ID: <557FD5C2.3080008@insa-rennes.fr> Hello, Le 16/06/2015 02:12, Reinaldo a ?crit : > How could I change the x_ticks and y_ticks? > For example: In this following script, the x_ticks displays each 2 units on X axis. I'd like each 0.5; and the y_ticks displays each 0.2 units on Y axis. I'd like each 0.1. How could I change those parameters? > function f=myfunction(x) > f = sin(x)./x > endfunction > > // Plot > scf(1); > clf(1); > xdata = linspace ( -12, 12, 100 ); > ydata = myfunction ( xdata ); > plot ( xdata , ydata ); > a=gca(); // obtendo o manipulador dos eixos correntes > a.grid=[1 1]; > xlabel("X","fontsize",4,"color","red"); > ylabel("Y","fontsize",4,"color","red");Thank you in advance. > All best,Reinaldo. For example, you could use function f=myfunction(x) f = sin(x)./x endfunction // Plot scf(1); clf(1); xdata = linspace ( -12, 12, 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); a=gca(); // obtendo o manipulador dos eixos correntes a.tight_limits = "on"; a.grid=[1 1]; xlabel("X","fontsize",4,"color","red"); ylabel("Y","fontsize",4,"color","red"); tmp=a.x_ticks; tmp.locations=[-12:12]'; tmp.labels=string(tmp.locations); a.x_ticks=tmp; tmp=a.y_ticks; tmp.locations=[-0.4:.1:1]'; tmp.labels=string(tmp.locations); a.y_ticks=tmp; See https://groups.google.com/forum/#!topic/comp.soft-sys.math.scilab/W2NwNbnys_I/discussion but I'm not sure that's the good link... -- Jean-Yves From paul.carrico at free.fr Wed Jun 17 09:05:51 2015 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Wed, 17 Jun 2015 09:05:51 +0200 (CEST) Subject: [Scilab-users] opposite to 'diag' keyword Message-ID: <2091806424.556909895.1434524751720.JavaMail.root@zimbra5-e1.priv.proxad.net> Dear All, The question is quite basic: I'm trying to opposite to 'diag' keyword in order to build a diagonal matrix from a vector, but I fail so far ... what is the trick? Thanks Paul From Jean-Yves.Baudais at insa-rennes.fr Wed Jun 17 09:14:39 2015 From: Jean-Yves.Baudais at insa-rennes.fr (jbaudais@insa-rennes.fr) Date: Wed, 17 Jun 2015 09:14:39 +0200 Subject: [Scilab-users] opposite to 'diag' keyword In-Reply-To: <2091806424.556909895.1434524751720.JavaMail.root@zimbra5-e1.priv.proxad.net> References: <2091806424.556909895.1434524751720.JavaMail.root@zimbra5-e1.priv.proxad.net> Message-ID: <55811E5F.6080406@insa-rennes.fr> Hello, Le 17/06/2015 09:05, paul.carrico at free.fr a ?crit : > The question is quite basic: I'm trying to opposite to 'diag' keyword in order to build a diagonal matrix from a vector, but I fail so far ... what is the trick? The diag help says "diagonal including or extracting" > diag([1 2]) gives the matrix [1 0;0 2] Is it the "opposite" you want? --Jean-Yves From paul.carrico at free.fr Wed Jun 17 09:22:12 2015 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Wed, 17 Jun 2015 09:22:12 +0200 (CEST) Subject: [Scilab-users] =?utf-8?b?UmXCoDogUmU6ICBvcHBvc2l0ZSB0byAnZGlhZycg?= =?utf-8?q?keyword?= In-Reply-To: <55811E5F.6080406@insa-rennes.fr> Message-ID: <1806214808.556968627.1434525732642.JavaMail.root@zimbra5-e1.priv.proxad.net> No A basic example: -->A = [1 0 0 0 ; 0 6 0 0 ; 0 0 -5 0 ; 0 0 0 31] A = 1. 0. 0. 0. 0. 6. 0. 0. 0. 0. - 5. 0. 0. 0. 0. 31. -->B = diag(A) B = 1. 6. - 5. 31. how a rebuild A from B .... NB: - the purpose is to find a basic but fast way to build a symmetrical matrix using the following syntax for example : K = K +Kd + K' - where K is the triangle upper (or lower) matrix and Kd is the diagonal - I want to decrease the amount of memory in avoiding to use a (nxn) matrix Paul ----- Mail d'origine ----- De: jbaudais at insa-rennes.fr ?: users at lists.scilab.org Envoy?: Wed, 17 Jun 2015 09:14:39 +0200 (CEST) Objet: Re: [Scilab-users] opposite to 'diag' keyword Hello, Le 17/06/2015 09:05, paul.carrico at free.fr a ?crit : > The question is quite basic: I'm trying to opposite to 'diag' keyword in order to build a diagonal matrix from a vector, but I fail so far ... what is the trick? The diag help says "diagonal including or extracting" > diag([1 2]) gives the matrix [1 0;0 2] Is it the "opposite" you want? --Jean-Yves _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users From pierre-aime.agnel at scilab-enterprises.com Wed Jun 17 09:28:38 2015 From: pierre-aime.agnel at scilab-enterprises.com (=?UTF-8?B?UGllcnJlLUFpbcOpIEFnbmVs?=) Date: Wed, 17 Jun 2015 09:28:38 +0200 Subject: [Scilab-users] =?utf-8?b?UmXCoDogUmU6ICBvcHBvc2l0ZSB0byAnZGlhZycg?= =?utf-8?q?keyword?= In-Reply-To: <1806214808.556968627.1434525732642.JavaMail.root@zimbra5-e1.priv.proxad.net> References: <1806214808.556968627.1434525732642.JavaMail.root@zimbra5-e1.priv.proxad.net> Message-ID: <558121A6.1090402@scilab-enterprises.com> Hello, diag works both ways: Le 17/06/2015 09:22, paul.carrico at free.fr a ?crit : > No > > A basic example: > > -->A = [1 0 0 0 ; 0 6 0 0 ; 0 0 -5 0 ; 0 0 0 31] > A = > > 1. 0. 0. 0. > 0. 6. 0. 0. > 0. 0. - 5. 0. > 0. 0. 0. 31. > > -->B = diag(A) > B = > > 1. > 6. > - 5. > 31. > -->A = [1 0 0 0 ; 0 6 0 0 ; 0 0 -5 0 ; 0 0 0 31] A = 1. 0. 0. 0. 0. 6. 0. 0. 0. 0. - 5. 0. 0. 0. 0. 31. -->B = diag(A) B = 1. 6. - 5. 31. -->C = diag(B) C = 1. 0. 0. 0. 0. 6. 0. 0. 0. 0. - 5. 0. 0. 0. 0. 31. > > how a rebuild A from B .... > > NB: > - the purpose is to find a basic but fast way to build a symmetrical matrix using the following syntax for example : K = K +Kd + K' > - where K is the triangle upper (or lower) matrix and Kd is the diagonal > - I want to decrease the amount of memory in avoiding to use a (nxn) matrix > > > Paul > > ----- Mail d'origine ----- > De: jbaudais at insa-rennes.fr > ?: users at lists.scilab.org > Envoy?: Wed, 17 Jun 2015 09:14:39 +0200 (CEST) > Objet: Re: [Scilab-users] opposite to 'diag' keyword > > Hello, > > Le 17/06/2015 09:05, paul.carrico at free.fr a ?crit : >> The question is quite basic: I'm trying to opposite to 'diag' keyword in order to build a diagonal matrix from a vector, but I fail so far ... what is the trick? > The diag help says "diagonal including or extracting" > > diag([1 2]) > gives the matrix [1 0;0 2] > Is it the "opposite" you want? > > --Jean-Yves > _______________________________________________ > 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 -- Pierre-Aim? Agnel R&D Projects Manager Phone: +33.1.80.77.04.67 Mobile: +33.6.82.49.35.23 ----------------------------------------------------------- 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 antoine.monmayrant at laas.fr Wed Jun 17 09:30:31 2015 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Wed, 17 Jun 2015 09:30:31 +0200 Subject: [Scilab-users] =?utf-8?b?UmXCoDogUmU6ICBvcHBvc2l0ZSB0byAnZGlhZycg?= =?utf-8?q?keyword?= In-Reply-To: <1806214808.556968627.1434525732642.JavaMail.root@zimbra5-e1.priv.proxad.net> References: <1806214808.556968627.1434525732642.JavaMail.root@zimbra5-e1.priv.proxad.net> Message-ID: <55812217.4000603@laas.fr> Le 06/17/2015 09:22 AM, paul.carrico at free.fr a ?crit : > No > > A basic example: > > -->A = [1 0 0 0 ; 0 6 0 0 ; 0 0 -5 0 ; 0 0 0 31] > A = > > 1. 0. 0. 0. > 0. 6. 0. 0. > 0. 0. - 5. 0. > 0. 0. 0. 31. > > -->B = diag(A) > B = > > 1. > 6. > - 5. > 31. > > > how a rebuild A from B .... easy: diag(B) Antoine > > NB: > - the purpose is to find a basic but fast way to build a symmetrical matrix using the following syntax for example : K = K +Kd + K' > - where K is the triangle upper (or lower) matrix and Kd is the diagonal > - I want to decrease the amount of memory in avoiding to use a (nxn) matrix > > > Paul > > ----- Mail d'origine ----- > De: jbaudais at insa-rennes.fr > ?: users at lists.scilab.org > Envoy?: Wed, 17 Jun 2015 09:14:39 +0200 (CEST) > Objet: Re: [Scilab-users] opposite to 'diag' keyword > > Hello, > > Le 17/06/2015 09:05, paul.carrico at free.fr a ?crit : >> The question is quite basic: I'm trying to opposite to 'diag' keyword in order to build a diagonal matrix from a vector, but I fail so far ... what is the trick? > The diag help says "diagonal including or extracting" > > diag([1 2]) > gives the matrix [1 0;0 2] > Is it the "opposite" you want? > > --Jean-Yves > _______________________________________________ > 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 paul.carrico at free.fr Wed Jun 17 09:32:10 2015 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Wed, 17 Jun 2015 09:32:10 +0200 (CEST) Subject: [Scilab-users] =?utf-8?b?UmXCoDogUmU6ICBSZcKgOiBSZTogIG9wcG9zaXRl?= =?utf-8?q?_to_=27diag=27_keyword?= In-Reply-To: <558121A6.1090402@scilab-enterprises.com> Message-ID: <2052513025.557000364.1434526330696.JavaMail.root@zimbra5-e1.priv.proxad.net> oups I probably did something wrong when I tested it ... Thanks Paul ----- Mail d'origine ----- De: Pierre-Aim? Agnel ?: users at lists.scilab.org Envoy?: Wed, 17 Jun 2015 09:28:38 +0200 (CEST) Objet: Re: [Scilab-users] Re?: Re: opposite to 'diag' keyword Hello, diag works both ways: Le 17/06/2015 09:22, paul.carrico at free.fr a ?crit : > No > > A basic example: > > -->A = [1 0 0 0 ; 0 6 0 0 ; 0 0 -5 0 ; 0 0 0 31] > A = > > 1. 0. 0. 0. > 0. 6. 0. 0. > 0. 0. - 5. 0. > 0. 0. 0. 31. > > -->B = diag(A) > B = > > 1. > 6. > - 5. > 31. > -->A = [1 0 0 0 ; 0 6 0 0 ; 0 0 -5 0 ; 0 0 0 31] A = 1. 0. 0. 0. 0. 6. 0. 0. 0. 0. - 5. 0. 0. 0. 0. 31. -->B = diag(A) B = 1. 6. - 5. 31. -->C = diag(B) C = 1. 0. 0. 0. 0. 6. 0. 0. 0. 0. - 5. 0. 0. 0. 0. 31. > > how a rebuild A from B .... > > NB: > - the purpose is to find a basic but fast way to build a symmetrical matrix using the following syntax for example : K = K +Kd + K' > - where K is the triangle upper (or lower) matrix and Kd is the diagonal > - I want to decrease the amount of memory in avoiding to use a (nxn) matrix > > > Paul > > ----- Mail d'origine ----- > De: jbaudais at insa-rennes.fr > ?: users at lists.scilab.org > Envoy?: Wed, 17 Jun 2015 09:14:39 +0200 (CEST) > Objet: Re: [Scilab-users] opposite to 'diag' keyword > > Hello, > > Le 17/06/2015 09:05, paul.carrico at free.fr a ?crit : >> The question is quite basic: I'm trying to opposite to 'diag' keyword in order to build a diagonal matrix from a vector, but I fail so far ... what is the trick? > The diag help says "diagonal including or extracting" > > diag([1 2]) > gives the matrix [1 0;0 2] > Is it the "opposite" you want? > > --Jean-Yves > _______________________________________________ > 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 -- Pierre-Aim? Agnel R&D Projects Manager Phone: +33.1.80.77.04.67 Mobile: +33.6.82.49.35.23 ----------------------------------------------------------- Scilab Enterprises 143bis rue Yves Le Coz - 78000 Versailles, France Phone: +33.1.80.77.04.68 http://www.scilab-enterprises.com From dominic.ernst at bystronic.com Wed Jun 17 12:25:50 2015 From: dominic.ernst at bystronic.com (DominicErnst) Date: Wed, 17 Jun 2015 03:25:50 -0700 (MST) Subject: [Scilab-users] Interactive zoom on plot with multiple axes Message-ID: <1434536750422-4032460.post@n3.nabble.com> Dear all, I have successfully created a plot with two y-axes (left / right), including some customization on axis labeling as well as legends. Now, if I try to interactively zoom the plot, only the last active axis is zoom-able. I tried to activate the other axis by issuing the gca(handle) command, but it did not change anything. Source: t = linspace(0,10,1000); P = 0.5*t; alpha = 2.5*sin(0.2*t); // Processed raw data figure f_proc = scf(2); f_proc.figure_name = "Processed raw Values"; f_proc.color_map = [get(sdf(),"color_map");[0,0.5,0]]; // add dark-green at end of default color table clf(f_proc); plot(t,alpha,'b'); //first plot alpha blue plot(t,P); //dummy plot -> only for working legend ax_Ang = gca(); ax_Ang.axes_visible = ['off','on','off']; // hide 2nd x-axis ax_Ang.y_location = 'left'; // alpha axis left ax_Ang.children.children(1).foreground = 33; // assign green to dummy plot -> correct display in legend ax_Ang.children.children(1).visible = 'off'; // hide dummy plot ax_Ang.filled = 'off'; // no background for this plot ax_Ang.tight_limits = 'on'; ax_Ang.auto_scale = 'off'; ax_Ang.data_bounds = [min(t),min(alpha);max(t),max(alpha)]; // i dont quite get why this line is necessary, since tight_limits is on ax_Ang.font_size = 2; ticks = ax_Ang.y_ticks; // change y-axis ticks to readable display ticks.locations = [0;%pi/4;%pi/2;3*%pi/4;%pi;5*%pi/4]; ticks.labels = ['$0$';'$\frac{\pi}{4}$';'$\frac{\pi}{2}$';'$\frac{3\pi}{4}$';'$\pi$';'$\frac{5\pi}{4}$']; ax_Ang.y_ticks = ticks; ax_Y = newaxes(); // create secondary axis plot(t,P); // plot second plot ax_Y.axes_visible = ['on','on','off']; // show x and y-axis ax_Y.y_location = 'right'; // y-axis right ax_Y.children.children.foreground = ax_Ang.children.children(1).foreground; // apply previously defined colors ax_Y.filled = ax_Ang.filled; // apply previously defined background ax_Y.tight_limits = ax_Ang.tight_limits; ax_Y.font_size = ax_Ang.font_size; title('Angle and Position vs Time','font_size',3); xlabel(ax_Y,'Time [ms]','font_size',2); ylabel(ax_Y,'Position [mm]','font_size',2); ylabel(ax_Ang,'Angle [-]','font_size',2); legend(ax_Ang,['Angle','Position']); Has anyone some experience with that. Would be a great help to explore my data a little easier. Any help is greatly appreciated Regards Dominic ----- -- Dominic Ernst CNC-PLC Software Developer @ Bystronic Laser AG -- View this message in context: http://mailinglists.scilab.org/Interactive-zoom-on-plot-with-multiple-axes-tp4032460.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.carrico at esterline.com Wed Jun 17 15:50:02 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Wed, 17 Jun 2015 13:50:02 +0000 Subject: [Scilab-users] eigs calculation Message-ID: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> Dear All, I'm performing a (mechanical) calculation using the eigs and I've been noticing that the results are strongly sensitive on the unit system I'm using; I can understand that high numbers can lead to some numerical "issues" ... Is there a way to increase the accuracy ? Paul PS: the 2 types of results NB: 1 (MPa) = 1E6 (Pa) 1 (mm) = 1E-3 (m) 1 (Kg/m^3) = 1E12 (T/mm^3) [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); a) calculation 1 in Pa, m, Kg/m^3 Natural frequency calculation: - Resonance 1 : 497.956 Hz - Resonance 2 : 3120.64 Hz - Resonance 3 : 5277.8 Hz - Resonance 4 : 6948.69 Hz - Resonance 5 : 8737.88 Hz - Resonance 6 : 15832.1 Hz - Resonance 7 : 17122.8 Hz - Resonance 8 : 20847.8 Hz - Resonance 9 : 26382.5 Hz - Resonance 10 : 28305.1 Hz - Resonance 11 : 34752 Hz - Resonance 12 : 36926.4 Hz b) Calculation in MPa, mm, T/mm^3 .... Natural frequency calculation: - Resonance 1 : 497.955 Hz - Resonance 2 : 497.956 Hz - Resonance 3 : 3120.59 Hz - Resonance 4 : 3120.64 Hz - Resonance 5 : 6948.69 Hz - Resonance 6 : 7463.93 Hz - Resonance 7 : 8737.56 Hz - Resonance 8 : 8737.88 Hz - Resonance 9 : 17121.6 Hz - Resonance 10 : 17122.8 Hz - Resonance 11 : 20847.8 Hz - Resonance 12 : 22390 Hz EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Wed Jun 17 16:30:17 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Wed, 17 Jun 2015 07:30:17 -0700 (MST) Subject: [Scilab-users] console character trouble Message-ID: <1434551417223-4032462.post@n3.nabble.com> Suddenly my console screen does't work.Please look at As well Scinote doesn't work I reinstalled Scilab but same trouble happens again. -- View this message in context: http://mailinglists.scilab.org/console-character-trouble-tp4032462.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From tim at wescottdesign.com Wed Jun 17 17:55:50 2015 From: tim at wescottdesign.com (tim at wescottdesign.com) Date: Wed, 17 Jun 2015 08:55:50 -0700 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> Message-ID: On 2015-06-17 06:50, Carrico, Paul wrote: > Dear All, > > I'm performing a (mechanical) calculation using the eigs and I've been > noticing that the results are strongly sensitive on the unit system > I'm using; I can understand that high numbers can lead to some > numerical "issues" ? > > Is there a way to increase the accuracy ? > > Paul > > PS: the 2 types of results > > _NB_: > > 1 (MPa) = 1E6 (Pa) > > 1 (mm) = 1E-3 (m) > > 1 (Kg/m^3) = 1E12 (T/mm^3) > > [u,v] = > eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); > > > a) calculation 1 in Pa, m, Kg/m^3 > > Natural frequency calculation: > > - Resonance 1 : 497.956 Hz > > - Resonance 2 : 3120.64 Hz > > - Resonance 3 : 5277.8 Hz > > - Resonance 4 : 6948.69 Hz > > - Resonance 5 : 8737.88 Hz > > - Resonance 6 : 15832.1 Hz > > - Resonance 7 : 17122.8 Hz > > - Resonance 8 : 20847.8 Hz > > - Resonance 9 : 26382.5 Hz > > - Resonance 10 : 28305.1 Hz > > - Resonance 11 : 34752 Hz > > - Resonance 12 : 36926.4 Hz > > b) Calculation in MPa, mm, T/mm^3 ?. > > Natural frequency calculation: > > - Resonance 1 : 497.955 Hz > > - Resonance 2 : 497.956 Hz > > - Resonance 3 : 3120.59 Hz > > - Resonance 4 : 3120.64 Hz > > - Resonance 5 : 6948.69 Hz > > - Resonance 6 : 7463.93 Hz > > - Resonance 7 : 8737.56 Hz > > - Resonance 8 : 8737.88 Hz > > - Resonance 9 : 17121.6 Hz > > - Resonance 10 : 17122.8 Hz > > - Resonance 11 : 20847.8 Hz > > - Resonance 12 : 22390 Hz Hi Paul: I can't tell you about the innards of Scilab specifically, but eigenvalue calculation in general can be very sensitive to numerical issues. If you're entering the data by hand or otherwise truncating the source data your entire difference in results may just be from rounding error in your source data. If you're starting from one set of source data and multiplying by conversion constants, then you can try changing the tolerance (if it's not in the function then there's a global one, called, I think, %TOL). There are ways to make the matrices more numerically stable. I am absolutely positively not an expert on this, but I think that the more you can make your matrix into something with a band of non-zero numbers around the main diagonal and zeros elsewhere, the more stable the problem will be. If mechanical systems are like control systems, then the numerical stability of the matrix that describes the system dynamics is just a reflection of the real sensitivity of the real system to manufacturing variations -- it may be that, in a group of several physical units all assembled to the same specification, you'll find that much variation in the real world! From roger.cormier at ncf.ca Wed Jun 17 19:50:13 2015 From: roger.cormier at ncf.ca (roger.cormier at ncf.ca) Date: Wed, 17 Jun 2015 13:50:13 -0400 Subject: [Scilab-users] eigs calculation In-Reply-To: References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> Message-ID: Paul: Did you consider checking the condition number of the matrix? What I sometimes do, for problems that are ill conditioned by their very nature, is normalise the entries in an attempt to be as far away from singularity (i.e. to keep the rows and columns as orthogonal as possible) and then do a check on the condition number to see if there is an improvement. The lower the condition number the better. Once your calculations are done you can convert back to the un-normalised entries, does improve a bit when inverting matrices of ill-conditioned problems such as target-motion analysis. I'm willing to bet that the condition number of the second matrix is a bit lower than of the first one. Also, I'm noticing that there are four double poles, possibly six in total (resonnances 5 and 6, and 11 and 12), in the second calculation. You knowing your system, that could provide you with some hints as to which calculation to trust. Let me know how you make out. Regards, Roger. ___________________________ Dr. Roger Cormier, P.Eng. Le mer. 17 juin 2015 ? 11:55, tim at wescottdesign.com a ?crit : > On 2015-06-17 06:50, Carrico, Paul wrote: >> Dear All, >> I'm performing a (mechanical) calculation using the eigs and I've been >> noticing that the results are strongly sensitive on the unit system >> I'm using; I can understand that high numbers can lead to some >> numerical "issues" ? >> Is there a way to increase the accuracy ? >> Paul >> PS: the 2 types of results >> _NB_: >> 1 (MPa) = 1E6 (Pa) >> 1 (mm) = 1E-3 (m) >> 1 (Kg/m^3) = 1E12 (T/mm^3) >> [u,v] = >> eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); >> a) calculation 1 in Pa, m, Kg/m^3 >> Natural frequency calculation: >> - Resonance 1 : 497.956 Hz >> - Resonance 2 : 3120.64 Hz >> - Resonance 3 : 5277.8 Hz >> - Resonance 4 : 6948.69 Hz >> - Resonance 5 : 8737.88 Hz >> - Resonance 6 : 15832.1 Hz >> - Resonance 7 : 17122.8 Hz >> - Resonance 8 : 20847.8 Hz >> - Resonance 9 : 26382.5 Hz >> - Resonance 10 : 28305.1 Hz >> - Resonance 11 : 34752 Hz >> - Resonance 12 : 36926.4 Hz >> b) Calculation in MPa, mm, T/mm^3 ?. >> Natural frequency calculation: >> - Resonance 1 : 497.955 Hz >> - Resonance 2 : 497.956 Hz >> - Resonance 3 : 3120.59 Hz >> - Resonance 4 : 3120.64 Hz >> - Resonance 5 : 6948.69 Hz >> - Resonance 6 : 7463.93 Hz >> - Resonance 7 : 8737.56 Hz >> - Resonance 8 : 8737.88 Hz >> - Resonance 9 : 17121.6 Hz >> - Resonance 10 : 17122.8 Hz >> - Resonance 11 : 20847.8 Hz >> - Resonance 12 : 22390 Hz > > Hi Paul: > > I can't tell you about the innards of Scilab specifically, but eigenvalue calculation in general can be very sensitive to numerical issues. If you're entering the data by hand or otherwise truncating the source data your entire difference in results may just be from rounding error in your source data. > > If you're starting from one set of source data and multiplying by conversion constants, then you can try changing the tolerance (if it's not in the function then there's a global one, called, I think, %TOL). > > There are ways to make the matrices more numerically stable. I am absolutely positively not an expert on this, but I think that the more you can make your matrix into something with a band of non-zero numbers around the main diagonal and zeros elsewhere, the more stable the problem will be. > > If mechanical systems are like control systems, then the numerical stability of the matrix that describes the system dynamics is just a reflection of the real sensitivity of the real system to manufacturing variations -- it may be that, in a group of several physical units all assembled to the same specification, you'll find that much variation in the real world! > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From sgougeon at free.fr Wed Jun 17 20:32:15 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Wed, 17 Jun 2015 20:32:15 +0200 Subject: [Scilab-users] Interactive zoom on plot with multiple axes In-Reply-To: <1434536750422-4032460.post@n3.nabble.com> References: <1434536750422-4032460.post@n3.nabble.com> Message-ID: <5581BD2F.8050504@free.fr> Hi Dominic, This is a noticeable graphical regression arrived with Scilab 5.4. It is reported here http://bugzilla.scilab.org/show_bug.cgi?id=12110 and an even more disturbing and closely related one, because almost always met: http://bugzilla.scilab.org/show_bug.cgi?id=11959 You may awake these threads on Bugzilla. There is no known workaround. Regards Samuel Le 17/06/2015 12:25, DominicErnst a ?crit : > Dear all, > > I have successfully created a plot with two y-axes (left / right), including > some customization on axis labeling as well as legends. Now, if I try to > interactively zoom the plot, only the last active axis is zoom-able. I tried > to activate the other axis by issuing the gca(handle) command, but it did > not change anything. > > Source: > .../... > > Has anyone some experience with that. Would be a great help to explore my > data a little easier. > Any help is greatly appreciated > > Regards > > Dominic > From paul.carrico at esterline.com Wed Jun 17 22:18:01 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Wed, 17 Jun 2015 20:18:01 +0000 Subject: [Scilab-users] eigs calculation In-Reply-To: References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> Message-ID: <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> Dear All Thanks for the answers. To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); Nota: ddl means dof I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... In any way, thanks for the debate Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de roger.cormier at ncf.ca Envoy??: 17 June 2015 19:50 ??: International users mailing list for Scilab. Objet?: Re: [Scilab-users] eigs calculation Paul: Did you consider checking the condition number of the matrix? What I sometimes do, for problems that are ill conditioned by their very nature, is normalise the entries in an attempt to be as far away from singularity (i.e. to keep the rows and columns as orthogonal as possible) and then do a check on the condition number to see if there is an improvement. The lower the condition number the better. Once your calculations are done you can convert back to the un-normalised entries, does improve a bit when inverting matrices of ill-conditioned problems such as target-motion analysis. I'm willing to bet that the condition number of the second matrix is a bit lower than of the first one. Also, I'm noticing that there are four double poles, possibly six in total (resonnances 5 and 6, and 11 and 12), in the second calculation. You knowing your system, that could provide you with some hints as to which calculation to trust. Let me know how you make out. Regards, Roger. ___________________________ Dr. Roger Cormier, P.Eng. Le mer. 17 juin 2015 ? 11:55, tim at wescottdesign.com a ?crit : > On 2015-06-17 06:50, Carrico, Paul wrote: >> Dear All, >> I'm performing a (mechanical) calculation using the eigs and I've been >> noticing that the results are strongly sensitive on the unit system >> I'm using; I can understand that high numbers can lead to some >> numerical "issues" . >> Is there a way to increase the accuracy ? >> Paul >> PS: the 2 types of results >> _NB_: >> 1 (MPa) = 1E6 (Pa) >> 1 (mm) = 1E-3 (m) >> 1 (Kg/m^3) = 1E12 (T/mm^3) >> [u,v] = >> eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); >> a) calculation 1 in Pa, m, Kg/m^3 >> Natural frequency calculation: >> - Resonance 1 : 497.956 Hz >> - Resonance 2 : 3120.64 Hz >> - Resonance 3 : 5277.8 Hz >> - Resonance 4 : 6948.69 Hz >> - Resonance 5 : 8737.88 Hz >> - Resonance 6 : 15832.1 Hz >> - Resonance 7 : 17122.8 Hz >> - Resonance 8 : 20847.8 Hz >> - Resonance 9 : 26382.5 Hz >> - Resonance 10 : 28305.1 Hz >> - Resonance 11 : 34752 Hz >> - Resonance 12 : 36926.4 Hz >> b) Calculation in MPa, mm, T/mm^3 .. >> Natural frequency calculation: >> - Resonance 1 : 497.955 Hz >> - Resonance 2 : 497.956 Hz >> - Resonance 3 : 3120.59 Hz >> - Resonance 4 : 3120.64 Hz >> - Resonance 5 : 6948.69 Hz >> - Resonance 6 : 7463.93 Hz >> - Resonance 7 : 8737.56 Hz >> - Resonance 8 : 8737.88 Hz >> - Resonance 9 : 17121.6 Hz >> - Resonance 10 : 17122.8 Hz >> - Resonance 11 : 20847.8 Hz >> - Resonance 12 : 22390 Hz > > Hi Paul: > > I can't tell you about the innards of Scilab specifically, but eigenvalue calculation in general can be very sensitive to numerical issues. If you're entering the data by hand or otherwise truncating the source data your entire difference in results may just be from rounding error in your source data. > > If you're starting from one set of source data and multiplying by conversion constants, then you can try changing the tolerance (if it's not in the function then there's a global one, called, I think, %TOL). > > There are ways to make the matrices more numerically stable. I am absolutely positively not an expert on this, but I think that the more you can make your matrix into something with a band of non-zero numbers around the main diagonal and zeros elsewhere, the more stable the problem will be. > > If mechanical systems are like control systems, then the numerical stability of the matrix that describes the system dynamics is just a reflection of the real sensitivity of the real system to manufacturing variations -- it may be that, in a group of several physical units all assembled to the same specification, you'll find that much variation in the real world! > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIF-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=X5rGSmMzj_mrDPhlxODiYEfNco7jW0285o_vP15-sDI&s=AUkkoXcGhr1-zgfahFwMmtOJuuCDaaCeKtS2-qHM8NQ&e= _______________________________________________ users mailing list users at lists.scilab.org https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIF-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=X5rGSmMzj_mrDPhlxODiYEfNco7jW0285o_vP15-sDI&s=AUkkoXcGhr1-zgfahFwMmtOJuuCDaaCeKtS2-qHM8NQ&e= EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data From fujimoto2005 at gmail.com Thu Jun 18 02:30:28 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Wed, 17 Jun 2015 17:30:28 -0700 (MST) Subject: [Scilab-users] console character trouble In-Reply-To: <1434551417223-4032462.post@n3.nabble.com> References: <1434551417223-4032462.post@n3.nabble.com> Message-ID: <1434587428116-4032468.post@n3.nabble.com> I used 'System Restore' of windows to restore my system files and settings to an point in time when scilab worked with no trouble. The troubles have disappeared. There are some automatically updated programs for securities before the trouble happened and they might be the cause of trouble. Are there anyone who have same experiences? -- View this message in context: http://mailinglists.scilab.org/console-character-trouble-tp4032462p4032468.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From pgsmntry at gmail.com Thu Jun 18 09:15:59 2015 From: pgsmntry at gmail.com (Parthageet) Date: Thu, 18 Jun 2015 00:15:59 -0700 (MST) Subject: [Scilab-users] Mathematical Modelling using Xcos Message-ID: <1434611759685-4032469.post@n3.nabble.com> Dear All, I am doing a project and I am supposed to do the mathematical modelling for a Phase Locked Loop. Can anyone please help me in this regard? Actually, I am unable to figure out how to device a mathematical model for PLL using XCOS in Scilab. -- View this message in context: http://mailinglists.scilab.org/Mathematical-Modelling-using-Xcos-tp4032469.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.bignier at scilab-enterprises.com Thu Jun 18 09:28:51 2015 From: paul.bignier at scilab-enterprises.com (Paul Bignier) Date: Thu, 18 Jun 2015 09:28:51 +0200 Subject: [Scilab-users] Mathematical Modelling using Xcos In-Reply-To: <1434611759685-4032469.post@n3.nabble.com> References: <1434611759685-4032469.post@n3.nabble.com> Message-ID: <55827333.5070608@scilab-enterprises.com> Hello, Did you check all the Xcos demos? There are several filtering/controlling examples that may approach your need. Regards, Paul On 06/18/2015 09:15 AM, Parthageet wrote: > Dear All, > > I am doing a project and I am supposed to do the mathematical modelling for > a Phase Locked Loop. Can anyone please help me in this regard? Actually, I > am unable to figure out how to device a mathematical model for PLL using > XCOS in Scilab. > > > > -- > View this message in context: http://mailinglists.scilab.org/Mathematical-Modelling-using-Xcos-tp4032469.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 Serge.Steer at inria.fr Thu Jun 18 10:08:02 2015 From: Serge.Steer at inria.fr (Serge Steer) Date: Thu, 18 Jun 2015 10:08:02 +0200 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> Message-ID: <55827C62.3010605@inria.fr> Le 17/06/2015 22:18, Carrico, Paul a ?crit : > Dear All > > Thanks for the answers. > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); > > Nota: ddl means dof > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... There is no way to improve accurary with an option. In general the numerical algorithms try to return the best solution. But it should be possible to improve accuracy by a well suited normalisation (unit change for example!) The condition number of u gives a measure of the numerical difficulty to solve the problem Note that as stated the eigenvalue computation may be a ill conditionned problem in particular when there are clusters of eigenvalues Please find below a little script which illustrate a typicall case A=zeros(10,10)+diag(ones(1,9),1);A(10,1)=%eps; s=spec(A); clf;plot(real(s),imag(s),'+'); Serge Steer > In any way, thanks for the debate > > Paul From antoine.monmayrant at laas.fr Thu Jun 18 10:38:02 2015 From: antoine.monmayrant at laas.fr (Antoine Monmayrant) Date: Thu, 18 Jun 2015 10:38:02 +0200 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> Message-ID: <5582836A.2040604@laas.fr> Le 06/17/2015 10:18 PM, Carrico, Paul a ?crit : > Dear All > > Thanks for the answers. > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); > > Nota: ddl means dof > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... > > In any way, thanks for the debate Hi Paul, I won't go into the ill-conditionned matrix inversion, I don't know much about it. What I say below might be out of topic but I might also be useful. Here is a rule: whenever you use floating point arithmetic: stay close to one! That is, apply whatever normalization makes sense for your problem so that the numbers your are crunching are around 1 (your matrix elements in your case). Never, never use "units" in your actual calculation, just apply the proper renormalization to the result of your number crunching. The rational behind this rule is that the distance to the next float (or double) you can represent scales with the float value (see attached plot). Matlab has a function eps to measure this distance: http://fr.mathworks.com/help/matlab/ref/eps.html (no equivalent in scilab) You should have a look here for more fun with floats: http://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html Cheers, Antoine > > Paul > > -----Message d'origine----- > De : users [mailto:users-bounces at lists.scilab.org] De la part de roger.cormier at ncf.ca > Envoy? : 17 June 2015 19:50 > ? : International users mailing list for Scilab. > Objet : Re: [Scilab-users] eigs calculation > > Paul: > > Did you consider checking the condition number of the matrix? > > What I sometimes do, for problems that are ill conditioned by their very nature, is normalise the entries in an attempt to be as far away from singularity (i.e. to keep the rows and columns as orthogonal as possible) and then do a check on the condition number to see if there is an improvement. The lower the condition number the better. Once your calculations are done you can convert back to the un-normalised entries, does improve a bit when inverting matrices of ill-conditioned problems such as target-motion analysis. > > I'm willing to bet that the condition number of the second matrix is a bit lower than of the first one. Also, I'm noticing that there are four double poles, possibly six in total (resonnances 5 and 6, and 11 and 12), in the second calculation. You knowing your system, that could provide you with some hints as to which calculation to trust. > > Let me know how you make out. > > Regards, > > > Roger. > > > ___________________________ > Dr. Roger Cormier, P.Eng. > > > Le mer. 17 juin 2015 ? 11:55, tim at wescottdesign.com a ?crit : > >> On 2015-06-17 06:50, Carrico, Paul wrote: >>> Dear All, >>> I'm performing a (mechanical) calculation using the eigs and I've been >>> noticing that the results are strongly sensitive on the unit system >>> I'm using; I can understand that high numbers can lead to some >>> numerical "issues" . >>> Is there a way to increase the accuracy ? >>> Paul >>> PS: the 2 types of results >>> _NB_: >>> 1 (MPa) = 1E6 (Pa) >>> 1 (mm) = 1E-3 (m) >>> 1 (Kg/m^3) = 1E12 (T/mm^3) >>> [u,v] = >>> eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); >>> a) calculation 1 in Pa, m, Kg/m^3 >>> Natural frequency calculation: >>> - Resonance 1 : 497.956 Hz >>> - Resonance 2 : 3120.64 Hz >>> - Resonance 3 : 5277.8 Hz >>> - Resonance 4 : 6948.69 Hz >>> - Resonance 5 : 8737.88 Hz >>> - Resonance 6 : 15832.1 Hz >>> - Resonance 7 : 17122.8 Hz >>> - Resonance 8 : 20847.8 Hz >>> - Resonance 9 : 26382.5 Hz >>> - Resonance 10 : 28305.1 Hz >>> - Resonance 11 : 34752 Hz >>> - Resonance 12 : 36926.4 Hz >>> b) Calculation in MPa, mm, T/mm^3 .. >>> Natural frequency calculation: >>> - Resonance 1 : 497.955 Hz >>> - Resonance 2 : 497.956 Hz >>> - Resonance 3 : 3120.59 Hz >>> - Resonance 4 : 3120.64 Hz >>> - Resonance 5 : 6948.69 Hz >>> - Resonance 6 : 7463.93 Hz >>> - Resonance 7 : 8737.56 Hz >>> - Resonance 8 : 8737.88 Hz >>> - Resonance 9 : 17121.6 Hz >>> - Resonance 10 : 17122.8 Hz >>> - Resonance 11 : 20847.8 Hz >>> - Resonance 12 : 22390 Hz >> Hi Paul: >> >> I can't tell you about the innards of Scilab specifically, but eigenvalue calculation in general can be very sensitive to numerical issues. If you're entering the data by hand or otherwise truncating the source data your entire difference in results may just be from rounding error in your source data. >> >> If you're starting from one set of source data and multiplying by conversion constants, then you can try changing the tolerance (if it's not in the function then there's a global one, called, I think, %TOL). >> >> There are ways to make the matrices more numerically stable. I am absolutely positively not an expert on this, but I think that the more you can make your matrix into something with a band of non-zero numbers around the main diagonal and zeros elsewhere, the more stable the problem will be. >> >> If mechanical systems are like control systems, then the numerical stability of the matrix that describes the system dynamics is just a reflection of the real sensitivity of the real system to manufacturing variations -- it may be that, in a group of several physical units all assembled to the same specification, you'll find that much variation in the real world! >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIF-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=X5rGSmMzj_mrDPhlxODiYEfNco7jW0285o_vP15-sDI&s=AUkkoXcGhr1-zgfahFwMmtOJuuCDaaCeKtS2-qHM8NQ&e= > _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIF-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=X5rGSmMzj_mrDPhlxODiYEfNco7jW0285o_vP15-sDI&s=AUkkoXcGhr1-zgfahFwMmtOJuuCDaaCeKtS2-qHM8NQ&e= > > EXPORT CONTROL : > Cet email ne contient pas de donn?es techniques > This email does not contain technical data > _______________________________________________ > 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 -------------- A non-text attachment was scrubbed... Name: untitled.png Type: image/png Size: 16440 bytes Desc: not available URL: From dominic.ernst at bystronic.com Thu Jun 18 10:46:04 2015 From: dominic.ernst at bystronic.com (DominicErnst) Date: Thu, 18 Jun 2015 01:46:04 -0700 (MST) Subject: [Scilab-users] Interactive zoom on plot with multiple axes In-Reply-To: <5581BD2F.8050504@free.fr> References: <1434536750422-4032460.post@n3.nabble.com> <5581BD2F.8050504@free.fr> Message-ID: <1434617164670-4032473.post@n3.nabble.com> Hi Samuel, Thanks for pointing out these two bug reports. I have been looking into them but was not successful finding a report discussing my issue. I'll give it a shot reactivating them. Best, Dominic ----- -- Dominic Ernst CNC-PLC Software Developer @ Bystronic Laser AG -- View this message in context: http://mailinglists.scilab.org/Interactive-zoom-on-plot-with-multiple-axes-tp4032460p4032473.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From rei.listas at yahoo.com Thu Jun 18 16:58:05 2015 From: rei.listas at yahoo.com (Reinaldo) Date: Thu, 18 Jun 2015 07:58:05 -0700 (MST) Subject: [Scilab-users] How can I change the x_ticks and y_ticks? In-Reply-To: <557FD5C2.3080008@insa-rennes.fr> References: <807388956.3533964.1434413532097.JavaMail.yahoo@mail.yahoo.com> <557FD5C2.3080008@insa-rennes.fr> Message-ID: <1148562337.1176700.1434639458988.JavaMail.yahoo@mail.yahoo.com> Hi Jean, Thanks for your answer. This solves my problem! I reckon Scilab is an excellent and powerful mathematical software. However, it requires a medium/deep knowledge of programming to solve a simple problem. For example: To change some basic parameters such as: x- and y-ticks, the Scilab users need to program some tens of (complex) code instructions: tmp=a.x_ticks;? tmp.locations=[-12:12]';? tmp.labels=string(tmp.locations);? a.x_ticks=tmp;? tmp=a.y_ticks;? tmp.locations=[-0.4:.1:1]';? tmp.labels=string(tmp.locations);? a.y_ticks=tmp;? Why hasn't Scilab got a single code to change x- and y- parameters (e.g., a.x_ticks=[-12:1:12] and a.y_ticks=[-0.4:.1:1])? This could be easier for Scilab users and makes the script be compact! Furthermore, if I don't use the code:?a.tight_limits = "on" the y-axis range is from -0.4 to 1; however, if I use it, the y-axis range changes to -0.2 to 1. Why?? In addition, how could I change the y-axis range (i.e., y_min and y_max)? Thank you in advance. All best,Reinaldo. De: jbaudais [via Scilab / Xcos - Mailing Lists Archives] Para: Reinaldo Enviadas: Ter?a-feira, 16 de Junho de 2015 5:04 Assunto: Re: How can I change the x_ticks and y_ticks? Hello, Le 16/06/2015 02:12, Reinaldo a ?crit : > How could I change the x_ticks and y_ticks? > For example: In this following script, the x_ticks displays each 2 units on X axis. I'd like each 0.5; and the y_ticks displays each 0.2 units on Y axis. I'd like each 0.1. How could I change those parameters? > function f=myfunction(x) > f = sin(x)./x > endfunction > > // Plot > scf(1); > clf(1); > xdata = linspace ( -12, 12, 100 ); > ydata = myfunction ( xdata ); > plot ( xdata , ydata ); > a=gca(); // obtendo o manipulador dos eixos correntes > a.grid=[1 1]; > xlabel("X","fontsize",4,"color","red"); > ylabel("Y","fontsize",4,"color","red");Thank you in advance. > All best,Reinaldo. For example, you could use function f=myfunction(x) f = sin(x)./x endfunction // Plot scf(1); clf(1); xdata = linspace ( -12, 12, 100 ); ydata = myfunction ( xdata ); plot ( xdata , ydata ); a=gca(); // obtendo o manipulador dos eixos correntes a.tight_limits = "on"; a.grid=[1 1]; xlabel("X","fontsize",4,"color","red"); ylabel("Y","fontsize",4,"color","red"); tmp=a.x_ticks; tmp.locations=[-12:12]'; tmp.labels=string(tmp.locations); a.x_ticks=tmp; tmp=a.y_ticks; tmp.locations=[-0.4:.1:1]'; tmp.labels=string(tmp.locations); a.y_ticks=tmp; See https://groups.google.com/forum/#!topic/comp.soft-sys.math.scilab/W2NwNbnys_I/discussion but I'm not sure that's the good link... -- Jean-Yves _______________________________________________ users mailing list [hidden email] http://lists.scilab.org/mailman/listinfo/users If you reply to this email, your message will be added to the discussion below: http://mailinglists.scilab.org/Scilab-users-How-can-I-change-the-x-ticks-and-y-ticks-tp4032449p4032451.html To start a new topic under Scilab users - Mailing Lists Archives, email ml-node+s994242n2602246h53 at n3.nabble.com To unsubscribe from Scilab users - Mailing Lists Archives, click here. NAML -- View this message in context: http://mailinglists.scilab.org/Scilab-users-How-can-I-change-the-x-ticks-and-y-ticks-tp4032449p4032474.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From tim at wescottdesign.com Thu Jun 18 18:20:58 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Thu, 18 Jun 2015 09:20:58 -0700 Subject: [Scilab-users] Mathematical Modelling using Xcos In-Reply-To: <1434611759685-4032469.post@n3.nabble.com> References: <1434611759685-4032469.post@n3.nabble.com> Message-ID: <1434644458.2689.94.camel@Servo> On Thu, 2015-06-18 at 00:15 -0700, Parthageet wrote: > Dear All, > > I am doing a project and I am supposed to do the mathematical modelling for > a Phase Locked Loop. Can anyone please help me in this regard? Actually, I > am unable to figure out how to device a mathematical model for PLL using > XCOS in Scilab. Are you having trouble with the mathematical modeling part, or the XCos part? The biggest issue that I find in modeling phase-locked loops is first in deciding what sort of model I want (i.e., do you want a full-on continuous-time circuit model, a fully linearized model of the loop in lock, or something in between), and then, if I'm modeling something more abstract (e.g. a fully linearized model), how I want to represent everything (this is particularly challenging if the fully linearized model is in sampled time, and uses something non-traditional for the phase detect, "NCO", or both). -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From sgougeon at free.fr Thu Jun 18 21:00:45 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Thu, 18 Jun 2015 21:00:45 +0200 Subject: [Scilab-users] eigs calculation In-Reply-To: <5582836A.2040604@laas.fr> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> <5582836A.2040604@laas.fr> Message-ID: <5583155D.9070804@free.fr> Hello, Le 18/06/2015 10:38, Antoine Monmayrant a ?crit : > .../... > Matlab has a function eps to measure this distance: > http://fr.mathworks.com/help/matlab/ref/eps.html (no equivalent in > scilab) Beside Scilab's %eps, does Matlab have an equivalent to nearfloat()?: http://help.scilab.org/docs/5.5.2/en_US/nearfloat.html See also number_properties: http://help.scilab.org/docs/5.5.2/en_US/number_properties.html Samuel From sgougeon at free.fr Thu Jun 18 21:07:39 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Thu, 18 Jun 2015 21:07:39 +0200 Subject: [Scilab-users] Interactive zoom on plot with multiple axes In-Reply-To: <1434617164670-4032473.post@n3.nabble.com> References: <1434536750422-4032460.post@n3.nabble.com> <5581BD2F.8050504@free.fr> <1434617164670-4032473.post@n3.nabble.com> Message-ID: <558316FB.1080109@free.fr> Le 18/06/2015 10:46, DominicErnst a ?crit : > Thanks for pointing out these two bug reports. I have been looking into them > but was not successful finding a report discussing my issue. ?? http://bugzilla.scilab.org/show_bug.cgi?id=12110 /=> Bug description: Since the release 5.4 it's impossible to zoom on multiple plots. This is necessary, at least as an option, for plots with multiple axes (plotyy).// / Weren't you describing this case? If this is something else, you might open a new bug report. Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From dominic.ernst at bystronic.com Fri Jun 19 09:17:25 2015 From: dominic.ernst at bystronic.com (DominicErnst) Date: Fri, 19 Jun 2015 00:17:25 -0700 (MST) Subject: [Scilab-users] Interactive zoom on plot with multiple axes In-Reply-To: <558316FB.1080109@free.fr> References: <1434536750422-4032460.post@n3.nabble.com> <5581BD2F.8050504@free.fr> <1434617164670-4032473.post@n3.nabble.com> <558316FB.1080109@free.fr> Message-ID: <143FE62907A02945864AE5D3F4FC4E81BD71E5@CHNIS210.bystronic.metal.firm> Hei Samuel Sorry, misunderstanding! I was locking into bug-reports before asking the question on the mailinglist, but was not successful then. But what you pointed out is exactly the error I am experiencing. Thanks a lot again, and sorry for the confusion. Best, Dominic Ernst CNC-PLC Software Developer _______________________________________________ Bystronic Laser AG CH-3362 Niederoenz, Switzerland Phone +41 62 956 35 68, Fax +41 62 956 36 79 dominic.ernst at bystronic.com; http://www.bystronic.com Bystronic: Best choice. Laser | Bending | Waterjet Important Note This e-mail may contain confidential information that is for the sole attention of the intended recipient. Please notify the sender if you are not the intended recipient and delete the message. All unauthorized use of the information contained in this message is expressly prohibited. Wichtiger Hinweis Der Inhalt dieses E-Mails kann vertrauliche Angaben enthalten, die nur f?r den/die namentlich genannten Empf?nger bestimmt sind. Falls Sie nicht der Adressat dieses E-Mails sind, nehmen Sie Verbindung mit dem Absender auf und l?schen Sie diese Mitteilung. Jede unbefugte Verwendung der in dieser Mitteilung enthaltenen Informationen ist untersagt. Von: Samuel GOUGEON [via Scilab / Xcos - Mailing Lists Archives] [mailto:ml-node+s994242n4032477h8 at n3.nabble.com] Gesendet: Donnerstag, 18. Juni 2015 21:08 An: Ernst Dominic Betreff: Re: Interactive zoom on plot with multiple axes Le 18/06/2015 10:46, DominicErnst a ?crit : Thanks for pointing out these two bug reports. I have been looking into them but was not successful finding a report discussing my issue. ?? http://bugzilla.scilab.org/show_bug.cgi?id=12110 => Bug description: Since the release 5.4 it's impossible to zoom on multiple plots. This is necessary, at least as an option, for plots with multiple axes (plotyy). Weren't you describing this case? If this is something else, you might open a new bug report. Samuel _______________________________________________ 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/Interactive-zoom-on-plot-with-multiple-axes-tp4032460p4032477.html To unsubscribe from Interactive zoom on plot with multiple axes, click here. NAML ----- -- Dominic Ernst CNC-PLC Software Developer @ Bystronic Laser AG -- View this message in context: http://mailinglists.scilab.org/Interactive-zoom-on-plot-with-multiple-axes-tp4032460p4032478.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.carrico at esterline.com Fri Jun 19 10:49:57 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Fri, 19 Jun 2015 08:49:57 +0000 Subject: [Scilab-users] eigs calculation In-Reply-To: <55827C62.3010605@inria.fr> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> <55827C62.3010605@inria.fr> Message-ID: <3A6B7233274DB449A2A0053A47684F952B7B9360@BGS-EX01.auxitrol.ad> Dear I've been reading all the advices that (all of them) speak about the need to normalize the data in order to avoid the ill-conditioned issues (and to be the closest as possible to the value 1.0) ; the example shared by Serge Steer is quite interesting and highlights the ill-conditioning issue. Nevertheless I'm a bit confuse on how to do it in order to rebuild the results afterward: [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); pulse_ =sqrt(diag(u)) //pulse - circular frequency freq_ = pulse_/(2*%pi); //natural frequency Naively I thought that the following would have fix the problem: - using max(K) and max(M) to normalize the sparse matrix, - removing the smallest values (i.e. bellow a threshold value) ... but the matrix becomes singular (why ?). So is there any additional advice ? Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Serge Steer Envoy??: jeudi 18 juin 2015 10:08 ??: users at lists.scilab.org Objet?: Re: [Scilab-users] eigs calculation Le 17/06/2015 22:18, Carrico, Paul a ?crit : > Dear All > > Thanks for the answers. > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > [u,v] = > eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$) > ,n,"SM"); > > Nota: ddl means dof > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... There is no way to improve accurary with an option. In general the numerical algorithms try to return the best solution. But it should be possible to improve accuracy by a well suited normalisation (unit change for example!) The condition number of u gives a measure of the numerical difficulty to solve the problem Note that as stated the eigenvalue computation may be a ill conditionned problem in particular when there are clusters of eigenvalues Please find below a little script which illustrate a typicall case A=zeros(10,10)+diag(ones(1,9),1);A(10,1)=%eps; s=spec(A); clf;plot(real(s),imag(s),'+'); Serge Steer > In any way, thanks for the debate > > Paul _______________________________________________ users mailing list users at lists.scilab.org https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIF-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=S2zme-g7wK4GeuiYqaod2UALQUgCMA5H4YAmnwGMSjE&s=yrKeTLe9iyrw3ACr0KvvHnUusYyMaXo-RoPpv8yhipg&e= EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data From tim at wescottdesign.com Fri Jun 19 18:38:29 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Fri, 19 Jun 2015 09:38:29 -0700 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7B9360@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> <55827C62.3010605@inria.fr> <3A6B7233274DB449A2A0053A47684F952B7B9360@BGS-EX01.auxitrol.ad> Message-ID: <1434731909.3271.11.camel@Servo> On Fri, 2015-06-19 at 08:49 +0000, Carrico, Paul wrote: (top posting fixed) > -----Message d'origine----- > De : users [mailto:users-bounces at lists.scilab.org] De la part de Serge Steer Envoy? : jeudi 18 juin 2015 10:08 ? : users at lists.scilab.org Objet : Re: [Scilab-users] eigs calculation > > Le 17/06/2015 22:18, Carrico, Paul a ?crit : > > Dear All > > > > Thanks for the answers. > > > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > > [u,v] = > > eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$) > > ,n,"SM"); > > > > Nota: ddl means dof > > > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... > There is no way to improve accurary with an option. In general the numerical algorithms try to return the best solution. > But it should be possible to improve accuracy by a well suited normalisation (unit change for example!) > > The condition number of u gives a measure of the numerical difficulty to solve the problem > > > Note that as stated the eigenvalue computation may be a ill conditionned problem in particular when there are clusters of eigenvalues > > Please find below a little script which illustrate a typicall case A=zeros(10,10)+diag(ones(1,9),1);A(10,1)=%eps; > s=spec(A); > clf;plot(real(s),imag(s),'+'); > > Serge Steer > > In any way, thanks for the debate > > > > Paul Dear > > I've been reading all the advices that (all of them) speak about the need to normalize the data in order to avoid the ill-conditioned issues (and to be the closest as possible to the value 1.0) ; the example shared by Serge Steer is quite interesting and highlights the ill-conditioning issue. > > Nevertheless I'm a bit confuse on how to do it in order to rebuild the results afterward: > [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1): $,(ddl_elem+1):$),n,"SM"); > > pulse_ =sqrt(diag(u)) //pulse - circular frequency > freq_ = pulse_/(2*%pi); //natural frequency > > Naively I thought that the following would have fix the problem: > - using max(K) and max(M) to normalize the sparse matrix, > - removing the smallest values (i.e. bellow a threshold value) > > ... but the matrix becomes singular (why ?). > > So is there any additional advice ? > > Paul > "Normalization" refers to jiggering the numbers around in a way that does not change the problem in a strictly mathematical sense, but which makes it more tractable. d = eigs(A, B) computes the solutions to A * v = lambda * B * v So any nonsingular square matrix N won't change the problem if it's multiplied in: N * A * v = N * lambda * B * v Because lambda is a scalar (well, I hope I'm getting that right) you can change the problem to A_ = N * A, B_ = N * B, and d = eigs(A_, B_) will, theoretically, get the same answers as with the original matrices, but possibly with better numerical conditioning. This is what you're doing when you change units. If you don't mind the meaning of your eigenvectors changing, you can do a similarity transform. Start with N * A * N^(-1) * N * v = N * lambda * B * N^(-1) * N * v Now set A_ = N * A * N^(-1) B_ = N * B * N^(-1) v_ = N * v [d, v_] = eigs(A_, B_); v = N^(-1) * v_; will, again, theoretically give you the same numbers as before, but it may be better conditioned numerically. Actually _finding_ N, or giving you advise on how to do so, is beyond my powers -- but maybe this will set you on a better road. -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From dominic.ernst at bystronic.com Mon Jun 22 11:31:06 2015 From: dominic.ernst at bystronic.com (DominicErnst) Date: Mon, 22 Jun 2015 02:31:06 -0700 (MST) Subject: [Scilab-users] Redefined function does not remain redefined in macro Message-ID: <1434965466660-4032484.post@n3.nabble.com> Dear all, I stumbled across a strange behavior I don't understand correctly. Please also refer to the simplified code extract below. Apparently beta(x,y) refers to a built in SciLab function (beta function (Euler integral of the first kind)). If I redefine the function by a simple vector variable, I obviously get the warning that the function was redefined (funcprot(1)). Which is what I would expect. Now, if I try to write-access a particular value in this now initialized vector I get the Warning "Warning: obsolete use of '=' instead of '=='.". To my understanding, this happens as scilab still interprets the function as the default function beta and not the redefined vector. Is there any way to circumvent this? Or am I forced to choose another variable name? Example code snippet: clc; clear; funcprot(1); b = 2; p= linspace(0,20,10)'; beta= zeros(1:10)'; for i = 1:length(p)-1 do a(i)= b*cos(beta(i)/2); beta(i+1)= p(i+1)*a(i); end Any help is greatly appreciated Best, Dominic ----- -- Dominic Ernst CNC-PLC Software Developer @ Bystronic Laser AG -- View this message in context: http://mailinglists.scilab.org/Redefined-function-does-not-remain-redefined-in-macro-tp4032484.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.carrico at esterline.com Mon Jun 22 12:08:15 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Mon, 22 Jun 2015 10:08:15 +0000 Subject: [Scilab-users] eigs calculation In-Reply-To: <1434731909.3271.11.camel@Servo> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> <55827C62.3010605@inria.fr> <3A6B7233274DB449A2A0053A47684F952B7B9360@BGS-EX01.auxitrol.ad> <1434731909.3271.11.camel@Servo> Message-ID: <3A6B7233274DB449A2A0053A47684F952B7B96A0@BGS-EX01.auxitrol.ad> Dear Tim, Thanks for the enlightenment, but unfortunately I confess I don't know how to proceed now; I feel that I need to fix such issue; another example if needed (see screenshot in attachment), the product herebellow should lead to the Identity matrix Verif1 = u' * M * u (vhere u is the eigenvector matrix) Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Tim Wescott Envoy??: vendredi 19 juin 2015 18:38 ??: International users mailing list for Scilab. Objet?: Re: [Scilab-users] eigs calculation "Normalization" refers to jiggering the numbers around in a way that does not change the problem in a strictly mathematical sense, but which makes it more tractable. d = eigs(A, B) computes the solutions to A * v = lambda * B * v So any nonsingular square matrix N won't change the problem if it's multiplied in: N * A * v = N * lambda * B * v Because lambda is a scalar (well, I hope I'm getting that right) you can change the problem to A_ = N * A, B_ = N * B, and d = eigs(A_, B_) will, theoretically, get the same answers as with the original matrices, but possibly with better numerical conditioning. This is what you're doing when you change units. If you don't mind the meaning of your eigenvectors changing, you can do a similarity transform. Start with N * A * N^(-1) * N * v = N * lambda * B * N^(-1) * N * v Now set A_ = N * A * N^(-1) B_ = N * B * N^(-1) v_ = N * v [d, v_] = eigs(A_, B_); v = N^(-1) * v_; will, again, theoretically give you the same numbers as before, but it may be better conditioned numerically. Actually _finding_ N, or giving you advise on how to do so, is beyond my powers -- but maybe this will set you on a better road. -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 _______________________________________________ users mailing list users at lists.scilab.org https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=qhzR7ofb8i_twSz_CIbw-knVke_-DuytzJ4yEXImbIo&s=pz1Y6NzKAzCPEEaszY2aGRyCKq2L1x-sRi5OiI2QgKM&e= EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data -------------- next part -------------- A non-text attachment was scrubbed... Name: screenshot.PNG Type: image/png Size: 26943 bytes Desc: screenshot.PNG URL: From fujimoto2005 at gmail.com Mon Jun 22 16:22:41 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Mon, 22 Jun 2015 07:22:41 -0700 (MST) Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix Message-ID: <1434982961621-4032488.post@n3.nabble.com> Let x=zeros(m,n,l) if m==1, x is generated as 2 dimension matrix, not 3 dim hyper-matrix. In that case the specification of elemnt like x(i,j,k) fails. But the value of m is unknown in advance. Is there any code which I can deal with such a case in general way? I don't like such as if m==1 then ??else??, because the actual number of the dimensions of hyper matrix i'm using is more than 4 and the nest of if-cluase is troublesome. Best ragards -- View this message in context: http://mailinglists.scilab.org/how-to-respond-the-shrinkage-of-hyper-matrix-tp4032488.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From pierre-aime.agnel at scilab-enterprises.com Mon Jun 22 16:43:56 2015 From: pierre-aime.agnel at scilab-enterprises.com (=?UTF-8?B?UGllcnJlLUFpbcOpIEFnbmVs?=) Date: Mon, 22 Jun 2015 16:43:56 +0200 Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix In-Reply-To: <1434982961621-4032488.post@n3.nabble.com> References: <1434982961621-4032488.post@n3.nabble.com> Message-ID: <55881F2C.8050408@scilab-enterprises.com> Hi, I don't think there is a problem (on my 5.5.2 version) : function hm_display(m,n,l) a = rand(m, n, l); for i = 1:m for j = 1:n for k = 1:l mprintf("value of index [%d, %d, %d] = %f \n", i, j, k, a(i, j, k)) end end end endfunction -->hm_display(1, 2, 1) value of index [1, 1, 1] = 0.614538 value of index [1, 2, 1] = 0.925896 -->hm_display(2, 2, 2) value of index [1, 1, 1] = 0.099382 value of index [1, 1, 2] = 0.921327 value of index [1, 2, 1] = 0.943183 value of index [1, 2, 2] = 0.900707 value of index [2, 1, 1] = 0.428058 value of index [2, 1, 2] = 0.944902 value of index [2, 2, 1] = 0.032740 value of index [2, 2, 2] = 0.809432 Le 22/06/2015 16:22, fujimoto2005 a ?crit : > Let x=zeros(m,n,l) > if m==1, x is generated as 2 dimension matrix, not 3 dim hyper-matrix. > In that case the specification of elemnt like x(i,j,k) fails. > But the value of m is unknown in advance. > Is there any code which I can deal with such a case in general way? > > I don't like such as if m==1 then ??else??, because the actual number of the > dimensions of hyper matrix i'm using is more than 4 and the nest of > if-cluase is troublesome. > > Best ragards > > > > -- > View this message in context: http://mailinglists.scilab.org/how-to-respond-the-shrinkage-of-hyper-matrix-tp4032488.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 R&D Projects Manager Phone: +33.1.80.77.04.67 Mobile: +33.6.82.49.35.23 ----------------------------------------------------------- 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 tim at wescottdesign.com Mon Jun 22 16:47:41 2015 From: tim at wescottdesign.com (tim at wescottdesign.com) Date: Mon, 22 Jun 2015 07:47:41 -0700 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7B96A0@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> <55827C62.3010605@inria.fr> <3A6B7233274DB449A2A0053A47684F952B7B9360@BGS-EX01.auxitrol.ad> <1434731909.3271.11.camel@Servo> <3A6B7233274DB449A2A0053A47684F952B7B96A0@BGS-EX01.auxitrol.ad> Message-ID: That's not at all atypical, and it's what you expect from a numerical process to arrive at an eigenvector matrix. In the context of ones on the main diagonal, 1e-16 is "mostly zero". On 2015-06-22 03:08, Carrico, Paul wrote: > Dear Tim, > > Thanks for the enlightenment, but unfortunately I confess I don't know > how to proceed now; > > I feel that I need to fix such issue; another example if needed (see > screenshot in attachment), the product herebellow should lead to the > Identity matrix > > Verif1 = u' * M * u > > (vhere u is the eigenvector matrix) > > > > Paul > > > -----Message d'origine----- > De?: users [mailto:users-bounces at lists.scilab.org] De la part de Tim > Wescott Envoy??: vendredi 19 juin 2015 18:38 ??: International users > mailing list for Scilab. > Objet?: Re: [Scilab-users] eigs calculation > > "Normalization" refers to jiggering the numbers around in a way that > does not change the problem in a strictly mathematical sense, but > which makes it more tractable. > > d = eigs(A, B) computes the solutions to > > A * v = lambda * B * v > > So any nonsingular square matrix N won't change the problem if it's > multiplied in: > > N * A * v = N * lambda * B * v > > Because lambda is a scalar (well, I hope I'm getting that right) you > can change the problem to > > A_ = N * A, B_ = N * B, and > > d = eigs(A_, B_) > > will, theoretically, get the same answers as with the original > matrices, but possibly with better numerical conditioning. This is > what you're doing when you change units. > > If you don't mind the meaning of your eigenvectors changing, you can > do a similarity transform. Start with > > N * A * N^(-1) * N * v = N * lambda * B * N^(-1) * N * v > > Now set > > A_ = N * A * N^(-1) > B_ = N * B * N^(-1) > v_ = N * v > > [d, v_] = eigs(A_, B_); v = N^(-1) * v_; > > will, again, theoretically give you the same numbers as before, but it > may be better conditioned numerically. > > Actually _finding_ N, or giving you advise on how to do so, is beyond > my powers -- but maybe this will set you on a better road. > > -- > > Tim Wescott > www.wescottdesign.com > Control & Communications systems, circuit & software design. > Phone: 503.631.7815 > Cell: 503.349.8432 > > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=qhzR7ofb8i_twSz_CIbw-knVke_-DuytzJ4yEXImbIo&s=pz1Y6NzKAzCPEEaszY2aGRyCKq2L1x-sRi5OiI2QgKM&e= > > EXPORT CONTROL : > Cet email ne contient pas de donn?es techniques > This email does not contain technical data From sgougeon at free.fr Mon Jun 22 22:05:02 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 22 Jun 2015 22:05:02 +0200 Subject: [Scilab-users] Redefined function does not remain redefined in macro In-Reply-To: <1434965466660-4032484.post@n3.nabble.com> References: <1434965466660-4032484.post@n3.nabble.com> Message-ID: <55886A6E.8020300@free.fr> Hello Dominic, Le 22/06/2015 11:31, DominicErnst a ?crit : > .../... > Now, if I try to write-access a particular value in this now initialized > vector I get the Warning "Warning: obsolete use of '=' instead of '=='.". It looks like a bug. After the beta=zeros(1:10)'; overwriting assignment, funptr("beta") still returns the beta()'s id. > .../... > > Is there any way to circumvent this? Or am I forced to choose another > variable name? You may use "Beta" instead of "beta". Samuel From dominic.ernst at bystronic.com Tue Jun 23 07:57:55 2015 From: dominic.ernst at bystronic.com (DominicErnst) Date: Mon, 22 Jun 2015 22:57:55 -0700 (MST) Subject: [Scilab-users] Redefined function does not remain redefined in macro In-Reply-To: <55886A6E.8020300@free.fr> References: <1434965466660-4032484.post@n3.nabble.com> <55886A6E.8020300@free.fr> Message-ID: <1435039075601-4032493.post@n3.nabble.com> Hi Samuel Thanks for your reply. That's what I thought, but as a beginner I'm always careful saying 'This is a bug', as I am not always sure I got everything right. Using Beta instead of beta is a good idea. Thanks a lot. Best, Dominic ----- -- Dominic Ernst CNC-PLC Software Developer @ Bystronic Laser AG -- View this message in context: http://mailinglists.scilab.org/Redefined-function-does-not-remain-redefined-in-macro-tp4032484p4032493.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From fujimoto2005 at gmail.com Wed Jun 24 05:45:49 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Tue, 23 Jun 2015 20:45:49 -0700 (MST) Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix In-Reply-To: <55881F2C.8050408@scilab-enterprises.com> References: <1434982961621-4032488.post@n3.nabble.com> <55881F2C.8050408@scilab-enterprises.com> Message-ID: <1435117549376-4032494.post@n3.nabble.com> Thanks for reply. Let me take another example. The following script fails because x is made as a 1?2?3?hyper matrix and I can't use 4th and 5th indexes. In my actual code, the number of 4th and 5th dimensions are variable and I cant't know in advance. x=zeros(1,2,3,1,1); disp(size(x)) a=x(1,1,1,1,1) Best regards -- View this message in context: http://mailinglists.scilab.org/how-to-respond-the-shrinkage-of-hyper-matrix-tp4032488p4032494.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Wed Jun 24 10:27:36 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Wed, 24 Jun 2015 10:27:36 +0200 Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix In-Reply-To: <1435117549376-4032494.post@n3.nabble.com> References: <1434982961621-4032488.post@n3.nabble.com> <55881F2C.8050408@scilab-enterprises.com> <1435117549376-4032494.post@n3.nabble.com> Message-ID: <558A69F8.9060608@free.fr> Le 24/06/2015 05:45, fujimoto2005 a ?crit : > Thanks for reply. > > Let me take another example. > The following script fails because x is made as a 1?2?3?hyper matrix and I > can't use 4th and 5th indexes. > In my actual code, the number of 4th and 5th dimensions are variable and I > cant't know in advance. > > x=zeros(1,2,3,1,1); > disp(size(x)) > a=x(1,1,1,1,1) It is a known behavior, reported for instance here for Scilab 5.5 : http://bugzilla.scilab.org/13237 It is fixed in Scilab 6 (so-called YaSp). Regards From ritikastone at gmail.com Mon Jun 22 07:34:26 2015 From: ritikastone at gmail.com (Ritika Jaiswal) Date: Sun, 21 Jun 2015 22:34:26 -0700 (MST) Subject: [Scilab-users] How to execute MS-VECM code in Grocer toolbox of SciLab Message-ID: <1434951266908-4032483.post@n3.nabble.com> Hello, Kindly help me with the execution of MS-VECM code in Grocer toolbox of Scilab. I need to execute the MS-VECM code in Scilab. I am not able to find any code for MS-VECM in this. Kindly provide the help on this issue. Regards Ritika Jaiswal Doctoral Research Scholar National Institute of Technology Karnataka Surathkal, Mangalore- 575 025, Karnataka, India email: ritikastone at gmail.com ritika.hm13f07 at nitk.edu.in -- View this message in context: http://mailinglists.scilab.org/How-to-execute-MS-VECM-code-in-Grocer-toolbox-of-SciLab-tp4032483.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From ritikastone at gmail.com Mon Jun 22 13:14:02 2015 From: ritikastone at gmail.com (Ritika Jaiswal) Date: Mon, 22 Jun 2015 04:14:02 -0700 (MST) Subject: [Scilab-users] How to execute MS-VECM code in Grocer toolbox of SciLab In-Reply-To: <1434951266908-4032483.post@n3.nabble.com> References: <1434951266908-4032483.post@n3.nabble.com> Message-ID: <1434971642031-4032486.post@n3.nabble.com> Hello Kindly provide the help on the above issue. I need to implement MS-VECM code in Scilab, Grocer toolbox. Kindly advise the code for implementing MS-VECM code in Scilab Regards Ritika Jaiswal Doctoral Research Scholar National Institute of Technology Karnataka Surathkal, Mangalore- 575 025, Karnataka, India email: ritikastone at gmail.com ritika.hm13f07 at nitk.edu.in -- View this message in context: http://mailinglists.scilab.org/How-to-execute-MS-VECM-code-in-Grocer-toolbox-of-SciLab-tp4032483p4032486.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From ritikastone at gmail.com Tue Jun 23 06:40:58 2015 From: ritikastone at gmail.com (Ritika Jaiswal) Date: Mon, 22 Jun 2015 21:40:58 -0700 (MST) Subject: [Scilab-users] How to execute MS-VECM code in Grocer toolbox of SciLab In-Reply-To: <1434951266908-4032483.post@n3.nabble.com> References: <1434951266908-4032483.post@n3.nabble.com> Message-ID: <1435034458808-4032492.post@n3.nabble.com> Hello Kindly provide help on my query. Its actually very urgent. Regards Ritika Jaiswal Doctoral Research Scholar National Institute of Technology Karnataka Surathkal, Mangalore- 575 025, Karnataka, India email: ritikastone at gmail.com ritika.hm13f07 at nitk.edu.in -- View this message in context: http://mailinglists.scilab.org/How-to-execute-MS-VECM-code-in-Grocer-toolbox-of-SciLab-tp4032483p4032492.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From bmbouter at gmail.com Wed Jun 24 14:26:35 2015 From: bmbouter at gmail.com (Brian Bouterse) Date: Wed, 24 Jun 2015 08:26:35 -0400 Subject: [Scilab-users] How to execute MS-VECM code in Grocer toolbox of SciLab In-Reply-To: <1435034458808-4032492.post@n3.nabble.com> References: <1434951266908-4032483.post@n3.nabble.com> <1435034458808-4032492.post@n3.nabble.com> Message-ID: Ritika, I have no idea how to help you resolve the technical issue, but I do want to point out that the request contains almost no technical details. Search the internet for articles on how to write a good technical post or bug report and try to repost your request with more detail. Good luck, Brian On Tue, Jun 23, 2015 at 12:40 AM, Ritika Jaiswal wrote: > Hello > > > Kindly provide help on my query. Its actually very urgent. > > > Regards > Ritika Jaiswal > Doctoral Research Scholar > National Institute of Technology Karnataka > Surathkal, Mangalore- 575 025, > Karnataka, India > > email: ritikastone at gmail.com > ritika.hm13f07 at nitk.edu.in > > > > > -- > View this message in context: > http://mailinglists.scilab.org/How-to-execute-MS-VECM-code-in-Grocer-toolbox-of-SciLab-tp4032483p4032492.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 > -- Brian Bouterse -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Wed Jun 24 15:23:04 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Wed, 24 Jun 2015 06:23:04 -0700 (MST) Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix In-Reply-To: <558A69F8.9060608@free.fr> References: <1434982961621-4032488.post@n3.nabble.com> <55881F2C.8050408@scilab-enterprises.com> <1435117549376-4032494.post@n3.nabble.com> <558A69F8.9060608@free.fr> Message-ID: <1435152184501-4032497.post@n3.nabble.com> Hi,Samuel GOUGEON Thanks for reply. I suppose 'Yasp' is one which is in the following address. https://www.scilab.org/development/nightly_builds/yasp Am I right and can I install Yasp in the same PC where Scilab 5.5.2 has been already installed? Best Regards -- View this message in context: http://mailinglists.scilab.org/how-to-respond-the-shrinkage-of-hyper-matrix-tp4032488p4032497.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From sgougeon at free.fr Wed Jun 24 16:28:07 2015 From: sgougeon at free.fr (sgougeon at free.fr) Date: Wed, 24 Jun 2015 16:28:07 +0200 (CEST) Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix In-Reply-To: <1435152184501-4032497.post@n3.nabble.com> Message-ID: <739353198.130518001.1435156087534.JavaMail.root@zimbra75-e12.priv.proxad.net> >I suppose 'Yasp' is one which is in the following address. > >https://www.scilab.org/development/nightly_builds/yasp > >Am I right Yes you are >and can I install Yasp in the same PC where Scilab 5.5.2 has been already installed? Yes, you may have as many Scilab versions as needed on the same PC, and you may run simultaneously as many Scilab sessions in any available Scilab version, as needed. The only limit is just the required RAM ;) Best regards Samuel Gougeon From fujimoto2005 at gmail.com Wed Jun 24 17:03:08 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Wed, 24 Jun 2015 08:03:08 -0700 (MST) Subject: [Scilab-users] how to respond the shrinkage of hyper-matrix In-Reply-To: <739353198.130518001.1435156087534.JavaMail.root@zimbra75-e12.priv.proxad.net> References: <1434982961621-4032488.post@n3.nabble.com> <55881F2C.8050408@scilab-enterprises.com> <1435117549376-4032494.post@n3.nabble.com> <558A69F8.9060608@free.fr> <1435152184501-4032497.post@n3.nabble.com> <739353198.130518001.1435156087534.JavaMail.root@zimbra75-e12.priv.proxad.net> Message-ID: <1435158188589-4032499.post@n3.nabble.com> Dear Samuel Gougeon Thanks for your information. I could installed 'Yasp' and found the speed of calculation is much faster than 5.52. But I face one problem. I can't interrupt the run of the script. Neither 'ctrl x ' key nor the interrupt command work and the run does't stop. Are there any way to fix it? Best regards -- View this message in context: http://mailinglists.scilab.org/how-to-respond-the-shrinkage-of-hyper-matrix-tp4032488p4032499.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From tim at wescottdesign.com Wed Jun 24 18:22:22 2015 From: tim at wescottdesign.com (Tim Wescott) Date: Wed, 24 Jun 2015 09:22:22 -0700 Subject: [Scilab-users] Redefined function does not remain redefined in macro In-Reply-To: <1434965466660-4032484.post@n3.nabble.com> References: <1434965466660-4032484.post@n3.nabble.com> Message-ID: <1435162942.6474.20.camel@Servo> On Mon, 2015-06-22 at 02:31 -0700, DominicErnst wrote: > Dear all, > > I stumbled across a strange behavior I don't understand correctly. > Please also refer to the simplified code extract below. > > Apparently beta(x,y) refers to a built in SciLab function (beta function > (Euler integral of the first kind)). If I redefine the function by a simple > vector variable, I obviously get the warning that the function was redefined > (funcprot(1)). Which is what I would expect. > > Now, if I try to write-access a particular value in this now initialized > vector I get the Warning "Warning: obsolete use of '=' instead of '=='.". To > my understanding, this happens as scilab still interprets the function as > the default function beta and not the redefined vector. > > Is there any way to circumvent this? Or am I forced to choose another > variable name? > > Example code snippet: > > clc; > clear; > funcprot(1); > > b = 2; > > p= linspace(0,20,10)'; > > beta= zeros(1:10)'; > > for i = 1:length(p)-1 do > a(i)= b*cos(beta(i)/2); > beta(i+1)= p(i+1)*a(i); > end > > Any help is greatly appreciated > > Best, > > Dominic Hey Dominic: Yes, I think it's a bug. However, I'm of the opinion that redefining built-ins in your code is in Very Bad Style. You've had the suggestion to use Beta -- I'd use beta_, or b, or something else that avoids depending on capitalization to keep track of stuff (because I also think that is Very Bad Style). Do keep in mind that this is probably a quibble unless you're writing applications that need to be maintained by other people, or big applications that need to last over half a year (I find that after half a year of not looking at some piece of my own code I count as "other people" when I get back to it). But -- it's a good habit to get into. -- Tim Wescott www.wescottdesign.com Control & Communications systems, circuit & software design. Phone: 503.631.7815 Cell: 503.349.8432 From manasdas17 at gmail.com Thu Jun 25 03:46:36 2015 From: manasdas17 at gmail.com (manasdas17) Date: Wed, 24 Jun 2015 18:46:36 -0700 (MST) Subject: [Scilab-users] How to use writec_f block in xcos Message-ID: <1435196796002-4032501.post@n3.nabble.com> I am trying to use writec_f block as wfile_f block has been removed from scilab-5.5.2. But when i run the model it worked only for once in win 8 system and after that no file is being created. Same problem also occuring ubuntu system.Could you let me know the proper way of using it or may be suggest some different block with sam functionality. Regards, Manas -- View this message in context: http://mailinglists.scilab.org/How-to-use-writec-f-block-in-xcos-tp4032501.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From dominic.ernst at bystronic.com Thu Jun 25 09:22:30 2015 From: dominic.ernst at bystronic.com (DominicErnst) Date: Thu, 25 Jun 2015 00:22:30 -0700 (MST) Subject: [Scilab-users] Redefined function does not remain redefined in macro In-Reply-To: <1435162942.6474.20.camel@Servo> References: <1434965466660-4032484.post@n3.nabble.com> <1435162942.6474.20.camel@Servo> Message-ID: <1435216950891-4032502.post@n3.nabble.com> Hi Tim Thank you for the hint. I actually do totally agree to your suggestion that this is not very good style. I thought, that since I will most probably never use the function beta in this particular project I could risk this kind of usage. However, I will try to avoid redefining built-in functions in the future. But first I will have to define a naming convention if I come across such situations again. Luckily, there are not many functions that will cause this issue. In my opinion it is always preferable to have short but meaningful variable names, such that code remains readable, be it for me after a long time or for others. But with a good naming convention one can achieve this goal. Thank you very much for your input, I appreciate it very much. It gives the impression to be supported by the community, which is very valuable. Best, Dominic ----- -- Dominic Ernst CNC-PLC Software Developer @ Bystronic Laser AG -- View this message in context: http://mailinglists.scilab.org/Redefined-function-does-not-remain-redefined-in-macro-tp4032484p4032502.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From paul.carrico at esterline.com Thu Jun 25 17:17:12 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Thu, 25 Jun 2015 15:17:12 +0000 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> References: <3A6B7233274DB449A2A0053A47684F952B7B805B@BGS-EX01.auxitrol.ad> <3A6B7233274DB449A2A0053A47684F952B7B812A@BGS-EX01.auxitrol.ad> Message-ID: <3A6B7233274DB449A2A0053A47684F952B7CA26B@BGS-EX01.auxitrol.ad> Dear all I'm still working on my "eigs" issue topic and I'm still trying to understand what's going wrong; I run a test case : - same function is launched 10 times - each time the input data are recorded (K,M) - for each loop the eigenvalues u and the eigenvectors v are recorded Then the values of each loop are compared with the values of the loop 1 If K,M,u remains strictly identical, it is not the case for u (the eigenvectors) ... why ? I've ever check any initialization issue, but everything seems to be ok Paul PS : the results of this case Max delta v2 - v1 = 453.857 Max delta K2 - K1 = 0 Max delta M2 - M1 = 0 Max delta v3 - v1 = 549.214 Max delta K3 - K1 = 0 Max delta M3 - M1 = 0 Max delta v4 - v1 = 585.95 Max delta K4 - K1 = 0 Max delta M4 - M1 = 0 Max delta v5 - v1 = 379.702 Max delta K5 - K1 = 0 Max delta M5 - M1 = 0 Max delta v6 - v1 = 489.844 Max delta K6 - K1 = 0 Max delta M6 - M1 = 0 Max delta v7 - v1 = 439.221 Max delta K7 - K1 = 0 Max delta M7 - M1 = 0 Max delta v8 - v1 = 432.406 Max delta K8 - K1 = 0 Max delta M8 - M1 = 0 Max delta v9 - v1 = 351.752 Max delta K9 - K1 = 0 Max delta M9 - M1 = 0 Max delta v10 - v1 = 554.515 Max delta K10 - K1 = 0 Max delta M10 - M1 = 0 -----Message d'origine----- De?: Carrico, Paul Envoy??: mercredi 17 juin 2015 22:18 ??: International users mailing list for Scilab. Objet?: RE: [Scilab-users] eigs calculation Dear All Thanks for the answers. To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); Nota: ddl means dof I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... In any way, thanks for the debate Paul EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data From smp.meroni at gmail.com Thu Jun 25 22:07:00 2015 From: smp.meroni at gmail.com (Simone Meroni) Date: Thu, 25 Jun 2015 22:07:00 +0200 Subject: [Scilab-users] number od decimals with csvWrite Message-ID: Hi I want to write some data and export them with csvWrite I need to have no more than 2 decimals for aesthetic reasons. How can I set the number of decimals? Thanks, Simone Meroni -------------- next part -------------- An HTML attachment was scrubbed... URL: From Serge.Steer at inria.fr Fri Jun 26 09:09:13 2015 From: Serge.Steer at inria.fr (Serge Steer) Date: Fri, 26 Jun 2015 09:09:13 +0200 Subject: [Scilab-users] number od decimals with csvWrite In-Reply-To: References: Message-ID: <558CFA99.2020104@inria.fr> The simplest way is to use round A=rand(3,3); B=round(100*A)/100; Serge Steer Le 25/06/2015 22:07, Simone Meroni a ?crit : > Hi > I want to write some data and export them with > > csvWrite > > I need to have no more than 2 decimals for aesthetic reasons. > > How can I set the number of decimals? > > Thanks, > > Simone Meroni > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From amonmayr at laas.fr Fri Jun 26 10:50:13 2015 From: amonmayr at laas.fr (Antoine Monmayrant) Date: Fri, 26 Jun 2015 10:50:13 +0200 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7CA26B@BGS-EX01.auxitrol.ad> Message-ID: <418b-558d1280-5-69a6aa80@86134162> Hi Paul, I don't really like a function that gives different answers for the very same input. That sounds like a bug to me. That being said, from the data you showed, it is not clear that your eigenvector are really different. If what you show is just a difference in the norm of the difference between the eigenvalue at iteration 1 and N, that might be OK. Indeed, if v is an eigenvector, a.v with a non zero-scalar, is also an eigenvector. You should check whether v1 and vN are colinear: if they are, the results are not really different, they just differ by a scaling factor. But I would still call it a bug, as a function should always give the same answer when given the same input parameters. Cheers, Antoine Le Jeudi 25 Juin 2015 17:17 CEST, "Carrico, Paul" a ?crit: > Dear all > > I'm still working on my "eigs" issue topic and I'm still trying to understand what's going wrong; > > I run a test case : > - same function is launched 10 times > - each time the input data are recorded (K,M) > - for each loop the eigenvalues u and the eigenvectors v are recorded > > Then the values of each loop are compared with the values of the loop 1 > > If K,M,u remains strictly identical, it is not the case for u (the eigenvectors) ... why ? > > I've ever check any initialization issue, but everything seems to be ok > > Paul > > PS : the results of this case > > Max delta v2 - v1 = 453.857 > Max delta K2 - K1 = 0 > Max delta M2 - M1 = 0 > > Max delta v3 - v1 = 549.214 > Max delta K3 - K1 = 0 > Max delta M3 - M1 = 0 > > Max delta v4 - v1 = 585.95 > Max delta K4 - K1 = 0 > Max delta M4 - M1 = 0 > > Max delta v5 - v1 = 379.702 > Max delta K5 - K1 = 0 > Max delta M5 - M1 = 0 > > Max delta v6 - v1 = 489.844 > Max delta K6 - K1 = 0 > Max delta M6 - M1 = 0 > > Max delta v7 - v1 = 439.221 > Max delta K7 - K1 = 0 > Max delta M7 - M1 = 0 > > Max delta v8 - v1 = 432.406 > Max delta K8 - K1 = 0 > Max delta M8 - M1 = 0 > > Max delta v9 - v1 = 351.752 > Max delta K9 - K1 = 0 > Max delta M9 - M1 = 0 > > Max delta v10 - v1 = 554.515 > Max delta K10 - K1 = 0 > Max delta M10 - M1 = 0 > > -----Message d'origine----- > De?: Carrico, Paul > Envoy??: mercredi 17 juin 2015 22:18 > ??: International users mailing list for Scilab. > Objet?: RE: [Scilab-users] eigs calculation > > Dear All > > Thanks for the answers. > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > [u,v] = eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$),n,"SM"); > > Nota: ddl means dof > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... > > In any way, thanks for the debate > > Paul > > EXPORT CONTROL : > Cet email ne contient pas de donn?es techniques > This email does not contain technical data > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > From paul.carrico at esterline.com Fri Jun 26 11:16:15 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Fri, 26 Jun 2015 09:16:15 +0000 Subject: [Scilab-users] eigs calculation In-Reply-To: <418b-558d1280-5-69a6aa80@86134162> References: <3A6B7233274DB449A2A0053A47684F952B7CA26B@BGS-EX01.auxitrol.ad> <418b-558d1280-5-69a6aa80@86134162> Message-ID: <3A6B7233274DB449A2A0053A47684F952B7CA370@BGS-EX01.auxitrol.ad> Hi Antoine, Thanks for the comments; In attachment the results of the calculations (note that the modal effective depends on the eigenvectors) The calculation for the difference calculations are described herebellow .. Paul ###################################################################################" // prog principal nbre = 10; for k = 1 : nbre printf("\n*******\nIteration num %d\n",k); s1 = " ... [u" + string(k) + ",v" + string(k) + ",K" + string(k) + ",M" + string(k) + "] = calcul_v() , ... [nl,nc] = size(v" + string(k) + ") , ... save(''u" + string(k) + ".bin'',''u" + string(k) + "'') , ... save(''v" + string(k) + ".bin'',''v" + string(k) + "'') , ... clear nl , ... clear nc , ... save(''K" + string(k) + ".bin'',''K" + string(k) + "'') , ... save(''M" + string(k) + ".bin'',''M" + string(k) + "'') , ... "; execstr(s1) ; end printf("\n\n"); // difference v1 - v2 for k = 2 : nbre s2 = " ... printf(""(eigenvalues) Max delta u%d - u1 = %g\n"",k,abs(max(u" + string(k) + " - u1))) , ... printf(""(eigenvectors) Max delta v%d - v1 = %g\n"",k,abs(max(v" + string(k) + " - v1))) , ... printf(""(input matrix) Max delta K%d - K1 = %g\n"",k,abs(max(K" + string(k) + " - K1))) , ... printf(""(input matrix) Max delta M%d - M1 = %g\n"",k,abs(max(M" + string(k) + " - M1))) , ... "; execstr(s2) ; printf("\n"); end ######################################################################################" -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Antoine Monmayrant Envoy??: vendredi 26 juin 2015 10:50 ??: International users mailing list for Scilab. Objet?: Re: [Scilab-users] eigs calculation Hi Paul, I don't really like a function that gives different answers for the very same input. That sounds like a bug to me. That being said, from the data you showed, it is not clear that your eigenvector are really different. If what you show is just a difference in the norm of the difference between the eigenvalue at iteration 1 and N, that might be OK. Indeed, if v is an eigenvector, a.v with a non zero-scalar, is also an eigenvector. You should check whether v1 and vN are colinear: if they are, the results are not really different, they just differ by a scaling factor. But I would still call it a bug, as a function should always give the same answer when given the same input parameters. Cheers, Antoine Le Jeudi 25 Juin 2015 17:17 CEST, "Carrico, Paul" a ?crit: > Dear all > > I'm still working on my "eigs" issue topic and I'm still trying to > understand what's going wrong; > > I run a test case : > - same function is launched 10 times > - each time the input data are recorded (K,M) > - for each loop the eigenvalues u and the eigenvectors v are recorded > > Then the values of each loop are compared with the values of the loop > 1 > > If K,M,u remains strictly identical, it is not the case for u (the eigenvectors) ... why ? > > I've ever check any initialization issue, but everything seems to be > ok > > Paul > > PS : the results of this case > > Max delta v2 - v1 = 453.857 > Max delta K2 - K1 = 0 > Max delta M2 - M1 = 0 > > Max delta v3 - v1 = 549.214 > Max delta K3 - K1 = 0 > Max delta M3 - M1 = 0 > > Max delta v4 - v1 = 585.95 > Max delta K4 - K1 = 0 > Max delta M4 - M1 = 0 > > Max delta v5 - v1 = 379.702 > Max delta K5 - K1 = 0 > Max delta M5 - M1 = 0 > > Max delta v6 - v1 = 489.844 > Max delta K6 - K1 = 0 > Max delta M6 - M1 = 0 > > Max delta v7 - v1 = 439.221 > Max delta K7 - K1 = 0 > Max delta M7 - M1 = 0 > > Max delta v8 - v1 = 432.406 > Max delta K8 - K1 = 0 > Max delta M8 - M1 = 0 > > Max delta v9 - v1 = 351.752 > Max delta K9 - K1 = 0 > Max delta M9 - M1 = 0 > > Max delta v10 - v1 = 554.515 > Max delta K10 - K1 = 0 > Max delta M10 - M1 = 0 > > -----Message d'origine----- > De?: Carrico, Paul > Envoy??: mercredi 17 juin 2015 22:18 > ??: International users mailing list for Scilab. > Objet?: RE: [Scilab-users] eigs calculation > > Dear All > > Thanks for the answers. > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > [u,v] = > eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$) > ,n,"SM"); > > Nota: ddl means dof > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... > > In any way, thanks for the debate > > Paul > > EXPORT CONTROL : > Cet email ne contient pas de donn?es techniques This email does not > contain technical data _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_m > ailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJh > ZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=1q3UvWvNm8K7-IZGCUzIV3uKngKalGVN-2TH > izNu_Ek&s=Qd7_9c6GV5rmRHnsvDdDgvnnodut0J24q0KOnk1zqb8&e= > _______________________________________________ users mailing list users at lists.scilab.org https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=1q3UvWvNm8K7-IZGCUzIV3uKngKalGVN-2THizNu_Ek&s=Qd7_9c6GV5rmRHnsvDdDgvnnodut0J24q0KOnk1zqb8&e= EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data -------------- next part -------------- A non-text attachment was scrubbed... Name: diary.dat Type: application/octet-stream Size: 11492 bytes Desc: diary.dat URL: From amonmayr at laas.fr Fri Jun 26 12:26:26 2015 From: amonmayr at laas.fr (Antoine Monmayrant) Date: Fri, 26 Jun 2015 12:26:26 +0200 Subject: [Scilab-users] eigs calculation In-Reply-To: <3A6B7233274DB449A2A0053A47684F952B7CA370@BGS-EX01.auxitrol.ad> Message-ID: <55e0-558d2900-3-6755ca00@260716128> Le Vendredi 26 Juin 2015 11:16 CEST, "Carrico, Paul" a ?crit: > Hi Antoine, > > Thanks for the comments; In attachment the results of the calculations (note that the modal effective depends on the eigenvectors) > > The calculation for the difference calculations are described herebellow .. I'm confused, who's who: u = value or vector v = value or vector ? Antoine > > > Paul > > ###################################################################################" > // prog principal > nbre = 10; > for k = 1 : nbre > printf("\n*******\nIteration num %d\n",k); > s1 = " ... > [u" + string(k) + ",v" + string(k) + ",K" + string(k) + ",M" + string(k) + "] = calcul_v() , ... > [nl,nc] = size(v" + string(k) + ") , ... > save(''u" + string(k) + ".bin'',''u" + string(k) + "'') , ... > save(''v" + string(k) + ".bin'',''v" + string(k) + "'') , ... > clear nl , ... > clear nc , ... > save(''K" + string(k) + ".bin'',''K" + string(k) + "'') , ... > save(''M" + string(k) + ".bin'',''M" + string(k) + "'') , ... > "; > execstr(s1) ; > > end > > printf("\n\n"); > > // difference v1 - v2 > for k = 2 : nbre > s2 = " ... > printf(""(eigenvalues) Max delta u%d - u1 = %g\n"",k,abs(max(u" + string(k) + " - u1))) , ... > printf(""(eigenvectors) Max delta v%d - v1 = %g\n"",k,abs(max(v" + string(k) + " - v1))) , ... > printf(""(input matrix) Max delta K%d - K1 = %g\n"",k,abs(max(K" + string(k) + " - K1))) , ... > printf(""(input matrix) Max delta M%d - M1 = %g\n"",k,abs(max(M" + string(k) + " - M1))) , ... > "; > execstr(s2) ; > printf("\n"); > end > ######################################################################################" > > > > > > > > > > > > > > > > > > > -----Message d'origine----- > De?: users [mailto:users-bounces at lists.scilab.org] De la part de Antoine Monmayrant Envoy??: vendredi 26 juin 2015 10:50 ??: International users mailing list for Scilab. > Objet?: Re: [Scilab-users] eigs calculation > > Hi Paul, > > I don't really like a function that gives different answers for the very same input. > That sounds like a bug to me. > That being said, from the data you showed, it is not clear that your eigenvector are really different. > If what you show is just a difference in the norm of the difference between the eigenvalue at iteration 1 and N, that might be OK. > Indeed, if v is an eigenvector, a.v with a non zero-scalar, is also an eigenvector. > You should check whether v1 and vN are colinear: if they are, the results are not really different, they just differ by a scaling factor. > But I would still call it a bug, as a function should always give the same answer when given the same input parameters. > > Cheers, > > Antoine > > > Le Jeudi 25 Juin 2015 17:17 CEST, "Carrico, Paul" a ?crit: > > > Dear all > > > > I'm still working on my "eigs" issue topic and I'm still trying to > > understand what's going wrong; > > > > I run a test case : > > - same function is launched 10 times > > - each time the input data are recorded (K,M) > > - for each loop the eigenvalues u and the eigenvectors v are recorded > > > > Then the values of each loop are compared with the values of the loop > > 1 > > > > If K,M,u remains strictly identical, it is not the case for u (the eigenvectors) ... why ? > > > > I've ever check any initialization issue, but everything seems to be > > ok > > > > Paul > > > > PS : the results of this case > > > > Max delta v2 - v1 = 453.857 > > Max delta K2 - K1 = 0 > > Max delta M2 - M1 = 0 > > > > Max delta v3 - v1 = 549.214 > > Max delta K3 - K1 = 0 > > Max delta M3 - M1 = 0 > > > > Max delta v4 - v1 = 585.95 > > Max delta K4 - K1 = 0 > > Max delta M4 - M1 = 0 > > > > Max delta v5 - v1 = 379.702 > > Max delta K5 - K1 = 0 > > Max delta M5 - M1 = 0 > > > > Max delta v6 - v1 = 489.844 > > Max delta K6 - K1 = 0 > > Max delta M6 - M1 = 0 > > > > Max delta v7 - v1 = 439.221 > > Max delta K7 - K1 = 0 > > Max delta M7 - M1 = 0 > > > > Max delta v8 - v1 = 432.406 > > Max delta K8 - K1 = 0 > > Max delta M8 - M1 = 0 > > > > Max delta v9 - v1 = 351.752 > > Max delta K9 - K1 = 0 > > Max delta M9 - M1 = 0 > > > > Max delta v10 - v1 = 554.515 > > Max delta K10 - K1 = 0 > > Max delta M10 - M1 = 0 > > > > -----Message d'origine----- > > De?: Carrico, Paul > > Envoy??: mercredi 17 juin 2015 22:18 > > ??: International users mailing list for Scilab. > > Objet?: RE: [Scilab-users] eigs calculation > > > > Dear All > > > > Thanks for the answers. > > > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > > [u,v] = > > eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1):$) > > ,n,"SM"); > > > > Nota: ddl means dof > > > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... > > > > In any way, thanks for the debate > > > > Paul > > > > EXPORT CONTROL : > > Cet email ne contient pas de donn?es techniques This email does not > > contain technical data _______________________________________________ > > users mailing list > > users at lists.scilab.org > > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_m > > ailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJh > > ZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=1q3UvWvNm8K7-IZGCUzIV3uKngKalGVN-2TH > > izNu_Ek&s=Qd7_9c6GV5rmRHnsvDdDgvnnodut0J24q0KOnk1zqb8&e= > > > > > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=1q3UvWvNm8K7-IZGCUzIV3uKngKalGVN-2THizNu_Ek&s=Qd7_9c6GV5rmRHnsvDdDgvnnodut0J24q0KOnk1zqb8&e= > > EXPORT CONTROL : > Cet email ne contient pas de donn?es techniques > This email does not contain technical data From paul.carrico at esterline.com Fri Jun 26 13:28:25 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Fri, 26 Jun 2015 11:28:25 +0000 Subject: [Scilab-users] eigs calculation In-Reply-To: <55e0-558d2900-3-6755ca00@260716128> References: <3A6B7233274DB449A2A0053A47684F952B7CA370@BGS-EX01.auxitrol.ad> <55e0-558d2900-3-6755ca00@260716128> Message-ID: <3A6B7233274DB449A2A0053A47684F952B7CA3E3@BGS-EX01.auxitrol.ad> U = eigenvalues -> matrix (nx1) V = eigenvectors -> matrix (nxn) NB: n = 12 since I'm foccusing of the n first natural frequencies .... Paul -----Message d'origine----- De?: users [mailto:users-bounces at lists.scilab.org] De la part de Antoine Monmayrant Envoy??: vendredi 26 juin 2015 12:26 ??: International users mailing list for Scilab. Objet?: Re: [Scilab-users] eigs calculation Le Vendredi 26 Juin 2015 11:16 CEST, "Carrico, Paul" a ?crit: > Hi Antoine, > > Thanks for the comments; In attachment the results of the calculations > (note that the modal effective depends on the eigenvectors) > > The calculation for the difference calculations are described herebellow .. I'm confused, who's who: u = value or vector v = value or vector ? Antoine > > > Paul > > ###################################################################################" > // prog principal > nbre = 10; > for k = 1 : nbre > printf("\n*******\nIteration num %d\n",k); > s1 = " ... > [u" + string(k) + ",v" + string(k) + ",K" + string(k) + ",M" + string(k) + "] = calcul_v() , ... > [nl,nc] = size(v" + string(k) + ") , ... > save(''u" + string(k) + ".bin'',''u" + string(k) + "'') , ... > save(''v" + string(k) + ".bin'',''v" + string(k) + "'') , ... > clear nl , ... > clear nc , ... > save(''K" + string(k) + ".bin'',''K" + string(k) + "'') , ... > save(''M" + string(k) + ".bin'',''M" + string(k) + "'') , ... > "; > execstr(s1) ; > > end > > printf("\n\n"); > > // difference v1 - v2 > for k = 2 : nbre > s2 = " ... > printf(""(eigenvalues) Max delta u%d - u1 = %g\n"",k,abs(max(u" + string(k) + " - u1))) , ... > printf(""(eigenvectors) Max delta v%d - v1 = %g\n"",k,abs(max(v" + string(k) + " - v1))) , ... > printf(""(input matrix) Max delta K%d - K1 = %g\n"",k,abs(max(K" + string(k) + " - K1))) , ... > printf(""(input matrix) Max delta M%d - M1 = %g\n"",k,abs(max(M" + string(k) + " - M1))) , ... > "; > execstr(s2) ; > printf("\n"); > end > ######################################################################################" > > > > > > > > > > > > > > > > > > > -----Message d'origine----- > De?: users [mailto:users-bounces at lists.scilab.org] De la part de Antoine Monmayrant Envoy??: vendredi 26 juin 2015 10:50 ??: International users mailing list for Scilab. > Objet?: Re: [Scilab-users] eigs calculation > > Hi Paul, > > I don't really like a function that gives different answers for the very same input. > That sounds like a bug to me. > That being said, from the data you showed, it is not clear that your eigenvector are really different. > If what you show is just a difference in the norm of the difference between the eigenvalue at iteration 1 and N, that might be OK. > Indeed, if v is an eigenvector, a.v with a non zero-scalar, is also an eigenvector. > You should check whether v1 and vN are colinear: if they are, the results are not really different, they just differ by a scaling factor. > But I would still call it a bug, as a function should always give the same answer when given the same input parameters. > > Cheers, > > Antoine > > > Le Jeudi 25 Juin 2015 17:17 CEST, "Carrico, Paul" a ?crit: > > > Dear all > > > > I'm still working on my "eigs" issue topic and I'm still trying to > > understand what's going wrong; > > > > I run a test case : > > - same function is launched 10 times > > - each time the input data are recorded (K,M) > > - for each loop the eigenvalues u and the eigenvectors v are > > recorded > > > > Then the values of each loop are compared with the values of the > > loop > > 1 > > > > If K,M,u remains strictly identical, it is not the case for u (the eigenvectors) ... why ? > > > > I've ever check any initialization issue, but everything seems to be > > ok > > > > Paul > > > > PS : the results of this case > > > > Max delta v2 - v1 = 453.857 > > Max delta K2 - K1 = 0 > > Max delta M2 - M1 = 0 > > > > Max delta v3 - v1 = 549.214 > > Max delta K3 - K1 = 0 > > Max delta M3 - M1 = 0 > > > > Max delta v4 - v1 = 585.95 > > Max delta K4 - K1 = 0 > > Max delta M4 - M1 = 0 > > > > Max delta v5 - v1 = 379.702 > > Max delta K5 - K1 = 0 > > Max delta M5 - M1 = 0 > > > > Max delta v6 - v1 = 489.844 > > Max delta K6 - K1 = 0 > > Max delta M6 - M1 = 0 > > > > Max delta v7 - v1 = 439.221 > > Max delta K7 - K1 = 0 > > Max delta M7 - M1 = 0 > > > > Max delta v8 - v1 = 432.406 > > Max delta K8 - K1 = 0 > > Max delta M8 - M1 = 0 > > > > Max delta v9 - v1 = 351.752 > > Max delta K9 - K1 = 0 > > Max delta M9 - M1 = 0 > > > > Max delta v10 - v1 = 554.515 > > Max delta K10 - K1 = 0 > > Max delta M10 - M1 = 0 > > > > -----Message d'origine----- > > De?: Carrico, Paul > > Envoy??: mercredi 17 juin 2015 22:18 ??: International users mailing > > list for Scilab. > > Objet?: RE: [Scilab-users] eigs calculation > > > > Dear All > > > > Thanks for the answers. > > > > To give more information's on what I'm doing (That's quite new I confess), I'm performing a (basic) finite element calculation with beams using sparse matrix (stiffness matrix K and mass matrix M). > > [u,v] = > > eigs(K((ddl_elem+1):$,(ddl_elem+1):$),M((ddl_elem+1):$,(ddl_elem+1): > > $) > > ,n,"SM"); > > > > Nota: ddl means dof > > > > I'm calculated first the natural frequencies using (K - omega^2.M).x=0 ... the pulse (or circular frequencies) are no more and no less than the eigenvalues of the above system (u = omega^2). > > > > Just a "mechanical" remark: since the beam is clamped in one side (and free on the tip), it is absolutely normal that you find twice the same natural frequency (1rst mode in one direction, the second one in a new direction at 90?) .... I've been really surprised to find it, but happy at the same time ... > > > > The origin of my question was: since it obvious that the results are strongly sensitive to the "units" (i.e. the numbers), I'm wondering if there is a way to control the accuracy of the eigenvalues calculation using eigs keywords ... > > > > In any way, thanks for the debate > > > > Paul > > > > EXPORT CONTROL : > > Cet email ne contient pas de donn?es techniques This email does not > > contain technical data > > _______________________________________________ > > users mailing list > > users at lists.scilab.org > > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org > > _m > > ailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXf > > Jh > > ZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=1q3UvWvNm8K7-IZGCUzIV3uKngKalGVN-2 > > TH izNu_Ek&s=Qd7_9c6GV5rmRHnsvDdDgvnnodut0J24q0KOnk1zqb8&e= > > > > > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_m > ailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJh > ZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=1q3UvWvNm8K7-IZGCUzIV3uKngKalGVN-2TH > izNu_Ek&s=Qd7_9c6GV5rmRHnsvDdDgvnnodut0J24q0KOnk1zqb8&e= > > EXPORT CONTROL : > Cet email ne contient pas de donn?es techniques This email does not > contain technical data _______________________________________________ users mailing list users at lists.scilab.org https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=AwIGaQ&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=DtRDqCwpPaahLQM8E-MeOGMZ4CHFjo7bH4Zh_lxEzr0&s=72hXuchUSXoMknI3enblwt4qNp-LSwTkkeeIaXZI9wE&e= EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data From krotersv at gmail.com Fri Jun 26 14:38:27 2015 From: krotersv at gmail.com (=?UTF-8?B?0KHRgtCw0L3QuNGB0LvQsNCy?=) Date: Fri, 26 Jun 2015 17:38:27 +0500 Subject: [Scilab-users] number od decimals with csvWrite In-Reply-To: <558CFA99.2020104@inria.fr> References: <558CFA99.2020104@inria.fr> Message-ID: <558D47C3.9060201@gmail.com> Hi. I suppose no need to use the round function. Use csvWrite: A=rand(3,3); csvWrite(A, 'file_with_A_data.txt', ' ', '.', "%.2e") 7.26e-01 2.32e-01 8.83e-01 1.99e-01 2.31e-01 6.53e-01 5.44e-01 2.16e-01 3.08e-01 By the way, if you want to use a comma as a decimal mark the fourth parameter should be ',': csvWrite(A, 'file_with_A_data.txt', ' ', ',', "%.2e") Stanislav 26.06.2015 12:09, Serge Steer ?????: > The simplest way is to use round > > A=rand(3,3); > B=round(100*A)/100; > > Serge Steer > Le 25/06/2015 22:07, Simone Meroni a ?crit : >> Hi >> I want to write some data and export them with >> >> csvWrite >> >> I need to have no more than 2 decimals for aesthetic reasons. >> >> How can I set the number of decimals? >> >> Thanks, >> >> Simone Meroni >> >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From samuel.enibe at unn.edu.ng Sat Jun 27 23:20:24 2015 From: samuel.enibe at unn.edu.ng (Samuel Enibe) Date: Sat, 27 Jun 2015 22:20:24 +0100 Subject: [Scilab-users] Output Argument from the SCILAB ode function Message-ID: I would like to use the SCILAB ode function with the option of second or more output arguments activated. In the example shown below (also attached), I would like to pass the variable dd in line 3 of the function linear. When the ode function is called without the second output argument as in line 12, it works, but when the second output argument is included as in line 13, it does not work. Any help in correcting this will be greatly appreciated. /////////////////////////////////////////////////////////////////// function [xdot, dd] =linear(t, x, A, u, B, omega) xdot=A*x+B*u(t,omega); dd = cos(A); endfunction function ut=u(t, omega) ut=sin(omega*t) endfunction A=[1 1;0 2]; B=[1;1]; omega=5; y0=[1;0]; t0=0; t=[0.1,0.2,0.5,1]; [y] = ode(y0,t0,t,list(linear,A,u,B,omega)) [y1, dd] = ode(y0,t0,t,list(linear,A,u,B,omega)) /////////////////////////////////////////////////////////// -- Samuel Ogbonna Enibe BEng (Nig), MSc (Reading, England), PhD (Nig) Professor of Mechanical Engineering Dean, Faculty of Engineering University of Nigeria, Nsukka, Nigeria Tel: +2348063646798 Email: samuel.enibe at unn.edu.ng enibesam at yahoo.com -------------- next part -------------- A non-text attachment was scrubbed... Name: example-ode.sce Type: application/x-scilab-sce Size: 321 bytes Desc: not available URL: From sgougeon at free.fr Sun Jun 28 12:57:52 2015 From: sgougeon at free.fr (Samuel Gougeon) Date: Sun, 28 Jun 2015 12:57:52 +0200 Subject: [Scilab-users] Output Argument from the SCILAB ode function In-Reply-To: References: Message-ID: <558FD330.8090705@free.fr> Hello, Le 27/06/2015 23:20, Samuel Enibe a ?crit : > I would like to use the SCILAB ode function with the option of second > or more output arguments activated. In the example shown below (also > attached), I would like to pass the variable dd in line 3 of the > function linear. > > When the ode function is called without the second output argument as > in line 12, it works, but when the second output argument is included > as in line 13, it does not work. > . In normal mode (not for "root" or "discrete" ones), the second and third output arguments require to provide them also as (last) input arguments: [y, w, iw] = ode([type,] y0, t0, t [,rtol [,atol]], f [,jac] [,w, iw]) This "iterated" mode is not really documented. It lacks of examples. Anyway, it does not look to be what you want to do. In the case that you've posted, the need to define (and provide) u() is unclear. Indeed, omega and t are available in linear(), so sin(omega*) can be computed in it instead of an external u(). By the way, even keeping u() does not require to provide it in the list: linear() sees it anyway. Finally, A is not modified while ode() is run. Nor dd=cos(A). Therefore, cos(A) can be computed out of ode(), without being among output arguments. Hope this helps Samuel Gougeon > > Any help in correcting this will be greatly appreciated. > > /////////////////////////////////////////////////////////////////// > function [xdot, dd] =linear(t, x, A, u, B, omega) > xdot=A*x+B*u(t,omega); > dd = cos(A); > endfunction > function ut=u(t, omega) > ut=sin(omega*t) > endfunction > A=[1 1;0 2]; > B=[1;1]; > omega=5; > y0=[1;0]; > t0=0; > t=[0.1,0.2,0.5,1]; > [y] = ode(y0,t0,t,list(linear,A,u,B,omega)) > [y1, dd] = ode(y0,t0,t,list(linear,A,u,B,omega)) > /////////////////////////////////////////////////////////// > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From samuel.enibe at unn.edu.ng Mon Jun 29 07:24:46 2015 From: samuel.enibe at unn.edu.ng (Samuel Enibe) Date: Mon, 29 Jun 2015 06:24:46 +0100 Subject: [Scilab-users] Output Argument from the SCILAB ode function In-Reply-To: <558FD330.8090705@free.fr> References: <558FD330.8090705@free.fr> Message-ID: Thank you very much, *Samuel Gougeon* for your kind response. I really want to use the *ode* function in the "iterated" mode so that a matrix calculated in the function passed to *ode* function can be reused as input in an iteration until a solution can be found. I will appreciate any example of how to pass the arguments for the "iterated" mode. In the classic example [y, w, iw] = ode([type,] y0, t0, t [,rtol [,atol]], f [,jac] [,w, iw]) (1) should the 2nd and 3rd output arguments be of the same dimensions with y? (2) should these 2nd and 3rd output arguments not be passed via a *list* as in the examples given in SCILAB documentation? I will appreciate any example of use of the "iterated" mode. Samuel Ogbonna Enibe BEng (Nig), MSc (Reading, England), PhD (Nig) Professor of Mechanical Engineering Dean, Faculty of Engineering University of Nigeria, Nsukka, Nigeria Tel: +2348063646798 Email: samuel.enibe at unn.edu.ng enibesam at yahoo.com On Sun, Jun 28, 2015 at 11:57 AM, Samuel Gougeon wrote: > Hello, > > Le 27/06/2015 23:20, Samuel Enibe a ?crit : > > I would like to use the SCILAB ode function with the option of second > or more output arguments activated. In the example shown below (also > attached), I would like to pass the variable dd in line 3 of the > function linear. > > When the ode function is called without the second output argument as > in line 12, it works, but when the second output argument is included > as in line 13, it does not work. > . > > In normal mode (not for "root" or "discrete" ones), the second and third > output arguments require to provide them also as (last) input arguments: > [y, w, iw] = ode([type,] y0, t0, t [,rtol [,atol]], f [,jac] [,w, iw]) > > This "iterated" mode is not really documented. It lacks of examples. > Anyway, it does not look to be what you want to do. > > In the case that you've posted, the need to define (and provide) u() is > unclear. > Indeed, omega and t are available in linear(), so sin(omega*) can be > computed in it > instead of an external u(). > By the way, even keeping u() does not require to provide it in the list: > linear() sees it anyway. > Finally, A is not modified while ode() is run. Nor dd=cos(A). Therefore, > cos(A) can be > computed out of ode(), without being among output arguments. > > Hope this helps > Samuel Gougeon > > > Any help in correcting this will be greatly appreciated. > > /////////////////////////////////////////////////////////////////// > function [xdot, dd] =linear(t, x, A, u, B, omega) > xdot=A*x+B*u(t,omega); > dd = cos(A); > endfunction > function ut=u(t, omega) > ut=sin(omega*t) > endfunction > A=[1 1;0 2]; > B=[1;1]; > omega=5; > y0=[1;0]; > t0=0; > t=[0.1,0.2,0.5,1]; > [y] = ode(y0,t0,t,list(linear,A,u,B,omega)) > [y1, dd] = ode(y0,t0,t,list(linear,A,u,B,omega)) > /////////////////////////////////////////////////////////// > > > > _______________________________________________ > users mailing listusers at lists.scilab.orghttp://lists.scilab.org/mailman/listinfo/users > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Mon Jun 29 19:12:52 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Mon, 29 Jun 2015 10:12:52 -0700 (MST) Subject: [Scilab-users] setting the seed of grand Message-ID: <1435597972698-4032520.post@n3.nabble.com> I want to set the same seed for same k for every i,j. for i=1:100 for j=1:100 for k=1:100 x=grand(10,10,'norm',0,1) ????? end end end Pleas teach me how to set such seed. -- View this message in context: http://mailinglists.scilab.org/setting-the-seed-of-grand-tp4032520.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From grocer.toolbox at gmail.com Mon Jun 29 20:06:06 2015 From: grocer.toolbox at gmail.com (Eric Dubois) Date: Mon, 29 Jun 2015 20:06:06 +0200 Subject: [Scilab-users] setting the seed of grand In-Reply-To: <1435597972698-4032520.post@n3.nabble.com> References: <1435597972698-4032520.post@n3.nabble.com> Message-ID: Hi I do nbot understand why you should keep the same seed (it should give rise to 10000 times the same results), but teh way to do that is: S=grand("getsd") for i=1:100 for j=1:100 for k=1:100 grand("setsd",S") x=grand(10,10,'norm',0,1) ????? end end end ?ric. 2015-06-29 19:12 GMT+02:00 fujimoto2005 : > I want to set the same seed for same k for every i,j. > > for i=1:100 > for j=1:100 > for k=1:100 > x=grand(10,10,'norm',0,1) > ????? > end > end > end > Pleas teach me how to set such seed. > > > > > > -- > View this message in context: > http://mailinglists.scilab.org/setting-the-seed-of-grand-tp4032520.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive > at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From jrafaelbguerra at hotmail.com Mon Jun 29 23:55:32 2015 From: jrafaelbguerra at hotmail.com (Rafael Guera) Date: Mon, 29 Jun 2015 22:55:32 +0100 Subject: [Scilab-users] setting the seed of grand In-Reply-To: References: <1435597972698-4032520.post@n3.nabble.com> Message-ID: What about getting the k-loop outside? for k=1:100 x=grand(10,10,'nor',0,1) for i=1:100 for j=1:100 ????? end end end Not possible? Rafael From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Eric Dubois Sent: Monday, June 29, 2015 7:06 PM To: International users mailing list for Scilab. Subject: Re: [Scilab-users] setting the seed of grand Hi I do nbot understand why you should keep the same seed (it should give rise to 10000 times the same results), but teh way to do that is: S=grand("getsd") for i=1:100 for j=1:100 for k=1:100 grand("setsd",S") x=grand(10,10,'norm',0,1) ????? end end end ?ric. 2015-06-29 19:12 GMT+02:00 fujimoto2005 >: I want to set the same seed for same k for every i,j. for i=1:100 for j=1:100 for k=1:100 x=grand(10,10,'norm',0,1) ????? end end end Pleas teach me how to set such seed. -- View this message in context: http://mailinglists.scilab.org/setting-the-seed-of-grand-tp4032520.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Tue Jun 30 02:45:52 2015 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Mon, 29 Jun 2015 17:45:52 -0700 (MST) Subject: [Scilab-users] setting the seed of grand In-Reply-To: References: <1435597972698-4032520.post@n3.nabble.com> Message-ID: <1435625152862-4032523.post@n3.nabble.com> Dear, Erick and Rafael The actual script is like this. I want to exclude the effect of the different seed for the value of d(i,j,k) whit different i,j. for i=1:100 for j=1:100 for k=1:100 grand("setsd",S") x=grand(10,10,'norm',0,1) y(i,j,k)=sum(mean(f(p1(i),p2(j),k,X)),'r'),'c') ????? end end d(i,j,k)=(y(i,j,k)-y(i-1,j,k))/(p1(i)-p1(i-1)) end By the way, is the right position of 'the grand("setsd",S")' between 'for j=1:100' and 'for k=1:100'? for j=1:100 grand("setsd",S") for k=1:100 Before k-loop, there is a time consuming part independent of k-value . If k-loop is the most outside loop, time consuming part will be recalculated and I want to avoid that situation. Best regards -- View this message in context: http://mailinglists.scilab.org/setting-the-seed-of-grand-tp4032520p4032523.html Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. From grocer.toolbox at gmail.com Tue Jun 30 08:24:29 2015 From: grocer.toolbox at gmail.com (=?utf-8?Q?=C3=89ric_Dubois?=) Date: Tue, 30 Jun 2015 08:24:29 +0200 Subject: [Scilab-users] setting the seed of grand In-Reply-To: <1435625152862-4032523.post@n3.nabble.com> References: <1435597972698-4032520.post@n3.nabble.com> <1435625152862-4032523.post@n3.nabble.com> Message-ID: <8232A816-AC7A-46BF-B938-46E850FB43DC@gmail.com> If you want different draws for different k but the same i and j, yes the grand("setsd",S) should be outside the k loop. By the way, help grand should have given you the answer to your request. ?ric Envoy? de mon iPhone > Le 30 juin 2015 ? 02:45, fujimoto2005 a ?crit : > > Dear, Erick and Rafael > > The actual script is like this. > I want to exclude the effect of the different seed for the value of d(i,j,k) > whit different i,j. > > for i=1:100 > for j=1:100 > for k=1:100 > grand("setsd",S") > x=grand(10,10,'norm',0,1) > y(i,j,k)=sum(mean(f(p1(i),p2(j),k,X)),'r'),'c') > ????? > end > end > d(i,j,k)=(y(i,j,k)-y(i-1,j,k))/(p1(i)-p1(i-1)) > end > > By the way, is the right position of 'the grand("setsd",S")' between 'for > j=1:100' and 'for k=1:100'? > > for j=1:100 > grand("setsd",S") > for k=1:100 > > Before k-loop, there is a time consuming part independent of k-value . > If k-loop is the most outside loop, time consuming part will be recalculated > and I want to avoid that situation. > > Best regards > > > > -- > View this message in context: http://mailinglists.scilab.org/setting-the-seed-of-grand-tp4032520p4032523.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From grocer.toolbox at gmail.com Tue Jun 30 08:24:29 2015 From: grocer.toolbox at gmail.com (=?utf-8?Q?=C3=89ric_Dubois?=) Date: Tue, 30 Jun 2015 08:24:29 +0200 Subject: [Scilab-users] setting the seed of grand In-Reply-To: <1435625152862-4032523.post@n3.nabble.com> References: <1435597972698-4032520.post@n3.nabble.com> <1435625152862-4032523.post@n3.nabble.com> Message-ID: <8232A816-AC7A-46BF-B938-46E850FB43DC@gmail.com> If you want different draws for different k but the same i and j, yes the grand("setsd",S) should be outside the k loop. By the way, help grand should have given you the answer to your request. ?ric Envoy? de mon iPhone > Le 30 juin 2015 ? 02:45, fujimoto2005 a ?crit : > > Dear, Erick and Rafael > > The actual script is like this. > I want to exclude the effect of the different seed for the value of d(i,j,k) > whit different i,j. > > for i=1:100 > for j=1:100 > for k=1:100 > grand("setsd",S") > x=grand(10,10,'norm',0,1) > y(i,j,k)=sum(mean(f(p1(i),p2(j),k,X)),'r'),'c') > ????? > end > end > d(i,j,k)=(y(i,j,k)-y(i-1,j,k))/(p1(i)-p1(i-1)) > end > > By the way, is the right position of 'the grand("setsd",S")' between 'for > j=1:100' and 'for k=1:100'? > > for j=1:100 > grand("setsd",S") > for k=1:100 > > Before k-loop, there is a time consuming part independent of k-value . > If k-loop is the most outside loop, time consuming part will be recalculated > and I want to avoid that situation. > > Best regards > > > > -- > View this message in context: http://mailinglists.scilab.org/setting-the-seed-of-grand-tp4032520p4032523.html > Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com. > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users From pgsmntry at gmail.com Tue Jun 30 11:15:07 2015 From: pgsmntry at gmail.com (Parthageet Samantaray) Date: Tue, 30 Jun 2015 14:45:07 +0530 Subject: [Scilab-users] Mathematical Modelling using Xcos In-Reply-To: <1434644458.2689.94.camel@Servo> References: <1434611759685-4032469.post@n3.nabble.com> <1434644458.2689.94.camel@Servo> Message-ID: Hello, I am sorry for the delay in replying because of certain work. Actually I am very new to Xcos and Scilab. But as a part project I want to build up some palette of my own which I would be able to use for the mathematical modelling of PLL. If possible help me to build up my own blocks(that is as per my specs). Regards, Partha On Thu, Jun 18, 2015 at 9:50 PM, Tim Wescott wrote: > On Thu, 2015-06-18 at 00:15 -0700, Parthageet wrote: > > Dear All, > > > > I am doing a project and I am supposed to do the mathematical modelling > for > > a Phase Locked Loop. Can anyone please help me in this regard? Actually, > I > > am unable to figure out how to device a mathematical model for PLL using > > XCOS in Scilab. > > Are you having trouble with the mathematical modeling part, or the XCos > part? > > The biggest issue that I find in modeling phase-locked loops is first in > deciding what sort of model I want (i.e., do you want a full-on > continuous-time circuit model, a fully linearized model of the loop in > lock, or something in between), and then, if I'm modeling something more > abstract (e.g. a fully linearized model), how I want to represent > everything (this is particularly challenging if the fully linearized > model is in sampled time, and uses something non-traditional for the > phase detect, "NCO", or both). > > -- > > Tim Wescott > www.wescottdesign.com > Control & Communications systems, circuit & software design. > Phone: 503.631.7815 > Cell: 503.349.8432 > > > _______________________________________________ > 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.carrico at esterline.com Tue Jun 30 13:37:14 2015 From: paul.carrico at esterline.com (Carrico, Paul) Date: Tue, 30 Jun 2015 11:37:14 +0000 Subject: [Scilab-users] curiosity : topological optimization Message-ID: <3A6B7233274DB449A2A0053A47684F952B7CA960@BGS-EX01.auxitrol.ad> Dear All, I've been reading articles using Matlab for topological optimization, and naturally I'm wondering if there's somebody who's using successfully Scilab for such topic (for my curiosity) http://www.topopt.dtu.dk/?q=node/2 (99 lines release) http://www.topopt.dtu.dk/?q=node/751 (88 lines release) Regards Paul EXPORT CONTROL : Cet email ne contient pas de donn?es techniques This email does not contain technical data -------------- next part -------------- An HTML attachment was scrubbed... URL: