From sgougeon at free.fr Thu Mar 1 00:27:58 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Thu, 1 Mar 2018 00:27:58 +0100 Subject: [Scilab-users] int2d In-Reply-To: <1519856231321-0.post@n3.nabble.com> References: <1519856231321-0.post@n3.nabble.com> Message-ID: <132a8b8f-a4e3-46e3-1d99-f636eccaf374@free.fr> Hello Masahiro, > function rt=multNormDensFunc(x_V, m_V, cov_M); I am not sure about the semi-colon... ? > [I2,e]=int2d(X,Y,list(multNormDensFunc,[0,0],[1,0;0,1]),[0.0001, 1, 20*20, 4000, 0]) The function expects only one parameter after both 2D operands x_V, m_V, while the list provides 2 parameters. So either the function is missing a y_V, or the list has an extra parameter. I don't know how int2d() manages this. In practical, int2d() needs mesh2d() to build and feed its X and Y arguments. While mesh2d() is still an external function in the Metanet module, that is not Scilab 6 - ready: https://atoms.scilab.org/toolboxes/metanet ... However, mesh_2() from the CGLAB module https://atoms.scilab.org/toolboxes/cglab should also be able to do the job (hoping that the output is int2d() compatible). These tools should be inside Scilab, or int2d() outsite : http://bugzilla.scilab.org/8284 For int3d(), it is even worse. I am not sure there is any mesh3d() available at all anywhere for Scilab. Samuel Le 28/02/2018 ? 23:17, fujimoto2005 a ?crit : > In order to calculate the expected value of a certain two-variable function, > I have created a density function of two-dimensional normal distribution. > I created a code like the attached file to check the correctness of the > function, but it does not work. > Where is the code wrong? > > I have been confirmed int2D supports list format for a function call. > Furthermore running the code doesn't show "called multNormDensFunc1". > This means it fails to call "multNormDensFunc" function. > > Best regards. > > test_multNorm.sce > > > > > -- > Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html > _______________________________________________ > 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 Thu Mar 1 03:21:04 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Wed, 28 Feb 2018 19:21:04 -0700 (MST) Subject: [Scilab-users] int2d In-Reply-To: <132a8b8f-a4e3-46e3-1d99-f636eccaf374@free.fr> References: <1519856231321-0.post@n3.nabble.com> <132a8b8f-a4e3-46e3-1d99-f636eccaf374@free.fr> Message-ID: <1519870864015-0.post@n3.nabble.com> Hi Samuel I fixed it as you pointed out and it worked. There was another mistake. I misunderstood the function len () of another language as size () and used it. Was mesh2d internalized? I have not downloaded the Metanet module, but int2d worked fine. I have not used the ATOMS yet, but I will try using it little by little. Are there any guides that wrote the whole picture? Best regards test_biv.sce -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From stephane.mottelet at utc.fr Thu Mar 1 09:29:42 2018 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Thu, 1 Mar 2018 09:29:42 +0100 Subject: [Scilab-users] strange %h_set overloading Message-ID: Hello, I have a question concerning the overloading function %h_set and concerning overloading mechanism in general. Here is an example where the original %h_set is redefined in order to display "in %h_set" when it is called: function hdl=%h_set(varargin) ??? disp("in %h_set") ??? if (size(varargin) == 3) ??????? hdl = varargin(1) ??????? str = varargin(2) ??????? val = varargin(3) ??????? for i=1:size(hdl,"*") ??????????? set(hdl(i),str,val) ??????? end, ??? elseif (size(varargin) == 2) ??????? str = varargin(1) ??????? val = varargin(2) ??????? hdl=[] ??????? set(str,val) ??? end, endfunction set(gcf(),'tag','foo'); set([gcf(),gcf()],'tag','foo'); --> set(gcf(),'tag','foo'); --> set([gcf(),gcf()],'tag','foo'); ?in %h_set %h_set is not called wen a 1x1 handle is used, but only when a vector of more than 2 handles is used. I suppose this mechanism is implemented i the built-in "set" ? My other question/remark is about the naming convention of built-in functions like "set". When macros do exist wit the same name but in different libraries,? the library name used as a prefix plays the role of a namespace and allows to disambiguate calls: e.g."graphicslib.plot" vs "plolib.plot". However, there is no namespace for built-in functions. The reason for these questions is that I am currently re-engineering the plotlib and need a smooth mechanism to emulate new properties/fields for *native* Scilab handles. S. -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet From stephane.mottelet at utc.fr Thu Mar 1 19:47:15 2018 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Thu, 1 Mar 2018 19:47:15 +0100 Subject: [Scilab-users] strange %h_set overloading In-Reply-To: References: Message-ID: <2513320c-8812-20fc-cd27-9f64567b3277@utc.fr> OK, answering to myself... Le 01/03/2018 ? 09:29, St?phane Mottelet a ?crit?: > Hello, > > I have a question concerning the overloading function %h_set and > concerning overloading mechanism in general. > > Here is an example where the original %h_set is redefined in order to > display "in %h_set" when it is called: > > function hdl=%h_set(varargin) > ??? disp("in %h_set") > ??? if (size(varargin) == 3) > ??????? hdl = varargin(1) > ??????? str = varargin(2) > ??????? val = varargin(3) > ??????? for i=1:size(hdl,"*") > ??????????? set(hdl(i),str,val) > ??????? end, > ??? elseif (size(varargin) == 2) > ??????? str = varargin(1) > ??????? val = varargin(2) > ??????? hdl=[] > ??????? set(str,val) > ??? end, > endfunction > > set(gcf(),'tag','foo'); > set([gcf(),gcf()],'tag','foo'); > > --> set(gcf(),'tag','foo'); > > --> set([gcf(),gcf()],'tag','foo'); > > ?in %h_set > > %h_set is not called wen a 1x1 handle is used, but only when a vector > of more than 2 handles is used. I suppose this mechanism is > implemented i the built-in "set" ? > ok, seen in the source (modules/graphics/sci_gateway/c/sci_set.c) S. -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet From fujimoto2005 at gmail.com Fri Mar 2 01:41:18 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Thu, 1 Mar 2018 17:41:18 -0700 (MST) Subject: [Scilab-users] "optim" doesn't work Message-ID: <1519951278553-0.post@n3.nabble.com> I am trying to solve the minimization problem like the attached code. The problem of finding "shifted value" from [mLogY1Prod, mLogY2] which minimizes the objective function "ExSqShiftedVal". First I calculated the values of the objective function around [mLogY1Prod, mLogY2]. You can see that the objective function decreases as you shift from mLogY1Prod, mLogY2 to larger values. Analysis of the objective function shows the followings. 1. The value of the objective function is nonnegative. 2. When the shift becomes very large, the objective function value becomes larger than that in case of shift 0. Therefore, there is always an optimal solution. Next, I ran "optim" code, but it does not work. In the case of algorithm "qn", the optimum value does not change from the initial value. It might be because a numerical integration function "intg" is used in the objective function and the function value is not smooth and it fails. But I am not sure. In the case of algorithm "nd", Warning "Round-off error detected, the requested tolerance (or default) cannot be achieved. Try using bigger tolerances." appeared and "optim" stopped. Please let me know if you know how to manage this problem. Best regards test_optim.sce -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From fujimoto2005 at gmail.com Fri Mar 2 16:40:41 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Fri, 2 Mar 2018 08:40:41 -0700 (MST) Subject: [Scilab-users] mean of elements corresponding element is "true" Message-ID: <1520005241921-0.post@n3.nabble.com> In the attached code,"boolean_M" is a random Boolean matrix. I want to get the average of each column of A whose corresponding element of boolean_M has "true" as row vector x_V, but it does not work. Where am I wrong? Best regards test_meanWtboolean.sce -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From sgougeon at free.fr Fri Mar 2 16:54:10 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 2 Mar 2018 16:54:10 +0100 Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <1520005241921-0.post@n3.nabble.com> References: <1520005241921-0.post@n3.nabble.com> Message-ID: <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> Le 02/03/2018 ? 16:40, fujimoto2005 a ?crit : > In the attached code,"boolean_M" is a random Boolean matrix. > I want to get the average of each column of A whose corresponding element of > boolean_M has "true" as row vector x_V, but it does not work. > Where am I wrong? > n=10; > boolean_M=rand(n,n)>0.5; > disp(boolean_M) > A=rand(n,n); > disp(A); > x_V=mean(A(boolean_M),'r'); Instead, you may use n = 5; boolean_M = rand(n,n)>0.5 A = grand(n,n,"uin",0,9) A(boolean_M) = %nan nanmean(A,"r") --> n = 5; --> boolean_M = rand(n,n)>0.5 boolean_M = T F F F T T T T T T T T F T T F T F T F F F T T F --> A = grand(n,n,"uin",0,9) A = 6. 2. 5. 4. 6. 6. 9. 3. 6. 5. 3. 9. 3. 5. 4. 0. 8. 7. 0. 7. 1. 6. 8. 7. 2. --> A(boolean_M) = %nan A = Nan 2. 5. 4. Nan Nan Nan Nan Nan Nan Nan Nan 3. Nan Nan 0. Nan 7. Nan 7. 1. 6. Nan Nan 2. --> nanmean(A,"r") ans = 0.5 4. 5. 4. 4.5 Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Fri Mar 2 16:57:49 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 2 Mar 2018 16:57:49 +0100 Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> References: <1520005241921-0.post@n3.nabble.com> <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> Message-ID: <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> Le 02/03/2018 ? 16:54, Samuel Gougeon a ?crit : > Le 02/03/2018 ? 16:40, fujimoto2005 a ?crit : >> In the attached code,"boolean_M" is a random Boolean matrix. >> I want to get the average of each column of A whose corresponding element of >> boolean_M has "true" as row vector x_V, but it does not work. >> Where am I wrong? >> n=10; >> boolean_M=rand(n,n)>0.5; >> disp(boolean_M) >> A=rand(n,n); >> disp(A); >> x_V=mean(A(boolean_M),'r'); > > Instead, you may use > n = 5; > boolean_M = rand(n,n)>0.5 > A = grand(n,n,"uin",0,9) > A(boolean_M) = %nan > nanmean(A,"r") > --> n = 5; > --> boolean_M = rand(n,n)>0.5 > boolean_M = > T F F F T > T T T T T > T T F T T > F T F T F > F F T T F > > --> A = grand(n,n,"uin",0,9) > A = > 6. 2. 5. 4. 6. > 6. 9. 3. 6. 5. > 3. 9. 3. 5. 4. > 0. 8. 7. 0. 7. > 1. 6. 8. 7. 2. > > --> A(boolean_M) = %nan or --> A(~boolean_M) = %nan to stick to your requirement -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Fri Mar 2 17:03:56 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Fri, 2 Mar 2018 09:03:56 -0700 (MST) Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> References: <1520005241921-0.post@n3.nabble.com> <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> Message-ID: <1520006636913-0.post@n3.nabble.com> Dear Samuel Thank you for your answer. Best regards -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From jrafaelbguerra at hotmail.com Sat Mar 3 13:15:30 2018 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Sat, 3 Mar 2018 12:15:30 +0000 Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> References: <1520005241921-0.post@n3.nabble.com> <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> Message-ID: An alternative solution: n=10; M= bool2s(rand(n,n)>0.5); // Boolean to zero-ones matrix A= rand(n,n); x= sum(A.*M,'r') ./ sum(M,'r'); Regards, Rafael -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Sat Mar 3 13:53:50 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Sat, 3 Mar 2018 05:53:50 -0700 (MST) Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: References: <1520005241921-0.post@n3.nabble.com> <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> Message-ID: <1520081630656-0.post@n3.nabble.com> Dear Rafael, Thank you for another excellent solution. Best regards -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From sgougeon at free.fr Sat Mar 3 17:49:10 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Sat, 3 Mar 2018 17:49:10 +0100 Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: References: <1520005241921-0.post@n3.nabble.com> <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> Message-ID: <2022da0e-dc51-329b-aed1-502fae2731ca@free.fr> Le 03/03/2018 ? 13:15, Rafael Guerra a ?crit : > > An alternative solution: > > n=10; > > M=bool2s(rand(n,n)>0.5);/// Boolean to zero-ones matrix/ > > A=rand(n,n); > > x=sum(A.*M,'r')./sum(M,'r'); > Indeed. This is basically what nanmean() does, managing in addtition cases where no components remain on a range (all "%F"). -------------- next part -------------- An HTML attachment was scrubbed... URL: From herozigle at gmail.com Sun Mar 4 10:10:16 2018 From: herozigle at gmail.com (Hermes) Date: Sun, 4 Mar 2018 02:10:16 -0700 (MST) Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <1520005241921-0.post@n3.nabble.com> References: <1520005241921-0.post@n3.nabble.com> Message-ID: <1520154616236-0.post@n3.nabble.com> how about this way: n=10; boolean_M=rand(n,n)>0.5; disp(boolean_M) A=rand(n,n); disp(A); x_V=mean(boolean_M.*A,'r') disp(x_V) Hermes -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From jrafaelbguerra at hotmail.com Sun Mar 4 11:34:17 2018 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Sun, 4 Mar 2018 10:34:17 +0000 Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <1520154616236-0.post@n3.nabble.com> References: <1520005241921-0.post@n3.nabble.com> <1520154616236-0.post@n3.nabble.com> Message-ID: Hermes, With your solution all false elements (zeros) get averaged with the true ones. The averaging should be done only over the true-entries. Regards, Rafael -----Original Message----- From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Hermes Sent: Sunday, March 04, 2018 10:10 AM To: users at lists.scilab.org Subject: Re: [Scilab-users] mean of elements corresponding element is "true" how about this way: n=10; boolean_M=rand(n,n)>0.5; disp(boolean_M) A=rand(n,n); disp(A); x_V=mean(boolean_M.*A,'r') disp(x_V) Hermes -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users From jrafaelbguerra at hotmail.com Sun Mar 4 11:54:46 2018 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Sun, 4 Mar 2018 10:54:46 +0000 Subject: [Scilab-users] mean of elements corresponding element is "true" In-Reply-To: <2022da0e-dc51-329b-aed1-502fae2731ca@free.fr> References: <1520005241921-0.post@n3.nabble.com> <0d716791-670e-b81f-91d4-54077df90ac5@free.fr> <527376cc-deac-aa5e-191b-4a07340f50f1@free.fr> <2022da0e-dc51-329b-aed1-502fae2731ca@free.fr> Message-ID: "Indeed. This is basically what nanmean() does, managing in addition cases where no components remain on a range (all "%F")." If all components are %F then the results in Scilab 6 for both solutions are %nan, which is convenient. (Scilab 5 requires additional attention because of division by zero errors, though) -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Mon Mar 5 16:26:49 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 5 Mar 2018 16:26:49 +0100 Subject: [Scilab-users] Clone a function, continued In-Reply-To: <7eef23bf-9649-a37b-17d1-d0a31456b025@free.fr> References: <6613d71e-24e3-713b-a861-2c46802ff5ac@utc.fr> <7eef23bf-9649-a37b-17d1-d0a31456b025@free.fr> Message-ID: <8f1a8a43-38fa-e332-7341-bda09053a484@free.fr> Le 28/02/2018 ? 18:21, Samuel Gougeon a ?crit : > About > > http://mailinglists.scilab.org/Scilab-users-Scilab-and-hdf5-format-small-survey-tp4037711p4037718.html > : > Le 28/02/2018 ? 14:52, St?phane Mottelet a ?crit : >> .../... >> Doing things like this is not possible any more. It was was also to >> add persistent variables in libraries, with the same mechanism. > > External libraries are not persistent > . > They are also erased by the so-loved-by-matlabers > clear > instruction... My fault: External libraries that are autoloaded ARE protected and persistent. They are protected by the startup procedure. So actually functions belonging to them are persistent. Only the libraries that are loaded during the session may be cleared. Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From florent.halgatte at gmail.com Tue Mar 6 19:46:58 2018 From: florent.halgatte at gmail.com (Florent) Date: Tue, 6 Mar 2018 11:46:58 -0700 (MST) Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot Message-ID: <1520362018170-0.post@n3.nabble.com> Good morning, I'm attempting to make a Xcos version of a Simulink model of a loudspeaker that I found on Daniel San Aus?n's Thesis (page 8) . Please see the attached file for the signal flow chart. It's "quite simple", just additions and multiplications. The parameters Le, Re, Bl, Mms, Rms, Cms are the Thiele & Small parameters of a loudspeaker (it's just floating values). The input voltage can be a sinus, a noise or whatever. I'm getting a message error "Algebraic loop". Do you know how to solve that ? I think it can be a problem in the simulation setup as I'm new to Scilab/Xcos. Here is the Xcos diagram: linear_LS.zcos I'm also trying to display the Bode plot at different stages of the diagram that's why I'm using the CPGE toolbox. If you have an idea to solve the problem, feel free to contact me. You can answer in French. I'll also try a different approach with the state-space bloc. -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From cfuttrup at gmail.com Tue Mar 6 20:18:42 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Tue, 6 Mar 2018 20:18:42 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: <1520362018170-0.post@n3.nabble.com> References: <1520362018170-0.post@n3.nabble.com> Message-ID: Hi Florent I'm interested in the loudspeaker rmodel (I work in the loudspeaker industry). The Master Thesis you're referencing to, it uses Matlab. Did you get hold of the Matlab files? ... Are you simply taking some inspiration from this Thesis and experimenting with XCOS? I tried to download your zcos file, but it came as a ZIP file. When extracted it contains a few files, but no ZCOS file. Appologies if I'm a lost cause here, since I have no experience with XCOS, even basic stuff might throw me off. Please explain. Best regards, Claus -------------- next part -------------- An HTML attachment was scrubbed... URL: From florent.halgatte at gmail.com Tue Mar 6 20:36:02 2018 From: florent.halgatte at gmail.com (Florent) Date: Tue, 6 Mar 2018 12:36:02 -0700 (MST) Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: References: <1520362018170-0.post@n3.nabble.com> Message-ID: <1520364962127-0.post@n3.nabble.com> Hi Claus, No, unfortunately I didn't get the Simulink files. And we don't have Simulink where I'm working. Here is the zcos file linear_LS.zcos . I'm new to Xcos too, I started last week. Best regards, Florent -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From cfuttrup at gmail.com Tue Mar 6 20:40:00 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Tue, 6 Mar 2018 20:40:00 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: <1520364962127-0.post@n3.nabble.com> References: <1520362018170-0.post@n3.nabble.com> <1520364962127-0.post@n3.nabble.com> Message-ID: Hi Florian I see it's a zcos file, but GMAIL downloads it as a zip file. I just rename it from zip to zcos, and it loads into Scilab. I'll look at it tomorrow. :-) Best regards, Claus On Tue, Mar 6, 2018 at 8:36 PM, Florent wrote: > Hi Claus, > > No, unfortunately I didn't get the Simulink files. And we don't have > Simulink where I'm working. > > Here is the zcos file linear_LS.zcos > . > > I'm new to Xcos too, I started last week. > > Best regards, > Florent > > > > -- > Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing- > Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Tue Mar 6 20:42:37 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Tue, 6 Mar 2018 20:42:37 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: <1520364962127-0.post@n3.nabble.com> References: <1520362018170-0.post@n3.nabble.com> <1520364962127-0.post@n3.nabble.com> Message-ID: Hi Team Scilab I wish I could buy this book in English: https://www.d-booker.fr/scilab-2-xcos/112-scilab-xcos.html /Claus I'm new to Xcos too, I started last week. > > Best regards, > Florent > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Wed Mar 7 07:47:51 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Wed, 7 Mar 2018 07:47:51 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: References: <1520362018170-0.post@n3.nabble.com> <1520364962127-0.post@n3.nabble.com> Message-ID: Hi Florian I loaded the zcos file this morning and compared it to Fig. 2.5 in the Master Thesis. It looks like your Re-box is connected before the integrator, which means it becomes omega-dependent. This is not what we normally do. Did you intend it to be this way? I'm not sure this is actually your problem - I'm just trying to help - and clear the path. Best regards, Claus On Tue, Mar 6, 2018 at 8:40 PM, Claus Futtrup wrote: > Hi Florian > > I see it's a zcos file, but GMAIL downloads it as a zip file. I just > rename it from zip to zcos, and it loads into Scilab. I'll look at it > tomorrow. :-) > > Best regards, > Claus > > > On Tue, Mar 6, 2018 at 8:36 PM, Florent > wrote: > >> Hi Claus, >> >> No, unfortunately I didn't get the Simulink files. And we don't have >> Simulink where I'm working. >> >> Here is the zcos file linear_LS.zcos >> . >> >> I'm new to Xcos too, I started last week. >> >> Best regards, >> Florent >> >> >> >> -- >> Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists- >> Archives-f2602246.html >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users >> > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From florent.halgatte at gmail.com Wed Mar 7 17:56:31 2018 From: florent.halgatte at gmail.com (Florent) Date: Wed, 7 Mar 2018 09:56:31 -0700 (MST) Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: References: <1520362018170-0.post@n3.nabble.com> <1520364962127-0.post@n3.nabble.com> Message-ID: <1520441791794-0.post@n3.nabble.com> Hi Claus, I multiplied i(t) by Re in order to match with the equation: Or in an other form: What do you mean be omega-dependent? Is it because in the integration 1/s, s=jw? Best regards, Florent -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From cfuttrup at gmail.com Wed Mar 7 18:21:51 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Wed, 7 Mar 2018 18:21:51 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: <1520441791794-0.post@n3.nabble.com> References: <1520362018170-0.post@n3.nabble.com> <1520364962127-0.post@n3.nabble.com> <1520441791794-0.post@n3.nabble.com> Message-ID: Hi Florian Yes, because of the integration. I'm comparing to the source (the thesis), Fig. 2.5. Best regards, Claus On Wed, Mar 7, 2018 at 5:56 PM, Florent wrote: > Hi Claus, > > I multiplied i(t) by Re in order to match with the equation: > > Or in an other form: > > > What do you mean be omega-dependent? > Is it because in the integration 1/s, s=jw? > > Best regards, > Florent > > > > -- > Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing- > Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Wed Mar 7 18:24:00 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Wed, 7 Mar 2018 18:24:00 +0100 Subject: [Scilab-users] GUI programming In-Reply-To: References: Message-ID: Hi Philippe The book arrived in the mail today. Now - if someone would translate the book about XCOS ... then I'd buy that too! Cheers, Claus On Wed, Jan 31, 2018 at 11:04 PM, philippe wrote: > Hi > > Le 24/01/2018 ? 19:49, Claus Futtrup a ?crit : > > Hi Scilabers > > > > Is there a comprehensive manual or book (in English, or alt German) > > about programming GUI in Scilab? (or a comprehensive web-page) > > programming GUI in scilab relies on a good comprehension of graphics > handles. The 4th chapter of my book is dedicated to graphics and contain > few pages about GUI conception (part 23.4 pages 338 to 353) : > > https://www.d-booker.fr/scilab-book1/325-scilab-fundamentals.html > > You can find a simple example (with call-back inside) detailed : > > https://www.youtube.com/watch?v=zRXr22VZNJM > > It's sufficient to understand the basis of GUI. > > I think that what I need is not just examples. I need explanations, and > > preferably something coherent that I can combine and work with, into > > something that works. > > That was the aim of my book : giving simple examples with detailed > explanations, well ordered so that you can learn scilab with more > efficiency than jumping from one help page to another. > > Philippe > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From florent.halgatte at gmail.com Wed Mar 7 18:45:27 2018 From: florent.halgatte at gmail.com (Florent) Date: Wed, 7 Mar 2018 10:45:27 -0700 (MST) Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: References: <1520362018170-0.post@n3.nabble.com> <1520364962127-0.post@n3.nabble.com> <1520441791794-0.post@n3.nabble.com> Message-ID: <1520444727411-0.post@n3.nabble.com> Hoo, you're right I did a "copy/paste" mistake... Thanks for notice that. Now I don't get the Algebraic loop error anymore but i cannot display the Bode plot, I need to work on it. Best regards, Florent -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From lepezhr at yahoo.com.ar Fri Mar 9 11:57:48 2018 From: lepezhr at yahoo.com.ar (Hector R. Lepez) Date: Fri, 9 Mar 2018 10:57:48 +0000 (UTC) Subject: [Scilab-users] Help me plot y = x^(1/3); References: <1473516067.13167106.1520593068556.ref@mail.yahoo.com> Message-ID: <1473516067.13167106.1520593068556@mail.yahoo.com> I have problems to graph the function x=[-10:1:10] y = x^(1/3); plot(x,y); the same in plot of google gives a graph but scilab does anything with the sign. H?ctor L?pez -------------- next part -------------- An HTML attachment was scrubbed... URL: From stephane.mottelet at utc.fr Fri Mar 9 13:24:28 2018 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Fri, 9 Mar 2018 13:24:28 +0100 Subject: [Scilab-users] Help me plot y = x^(1/3); In-Reply-To: <1473516067.13167106.1520593068556@mail.yahoo.com> References: <1473516067.13167106.1520593068556.ref@mail.yahoo.com> <1473516067.13167106.1520593068556@mail.yahoo.com> Message-ID: Le 09/03/2018 ? 11:57, Hector R. Lepez a ?crit?: > I have problems to graph the function > x=[-10:1:10] y = x^(1/3); plot(x,y); the same in plot of google gives > a graph but scilab does anything with the sign. > > H?ctor L?pez > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users By default in scilab non-integer power elevation of negative number is made in the complex domain. The particular case of the real cubic root can be treated like this : y = sign(x).*abs(x)^(1/3); otherwise (-1)^(1/3) gives you one of the complex cubic root of -1, namely exp(%i*%pi*1/3) S. -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet -------------- next part -------------- An HTML attachment was scrubbed... URL: From denis.crete at thalesgroup.com Fri Mar 9 13:18:46 2018 From: denis.crete at thalesgroup.com (CRETE Denis) Date: Fri, 9 Mar 2018 13:18:46 +0100 Subject: [Scilab-users] Help me plot y = x^(1/3); In-Reply-To: <1473516067.13167106.1520593068556@mail.yahoo.com> References: <1473516067.13167106.1520593068556.ref@mail.yahoo.com> <1473516067.13167106.1520593068556@mail.yahoo.com> Message-ID: <8F1D40232A0E68409E3FC23A30C326620179DC771D06@THSONEA01CMS04P.one.grp> Hello, The operator ^ seems to return the real part of one of the complex roots of x^3=1. If x=[-10:1:10] and y = x^(1/3); Thus, for x<0, the values of y are half the opposite of values obtained with x>0. With y=x.^(1/3), the values of y are complex for x<0. HTH Denis De : users [mailto:users-bounces at lists.scilab.org] De la part de Hector R. Lepez Envoy? : vendredi 9 mars 2018 11:58 ? : users at lists.scilab.org Objet : [Scilab-users] Help me plot y = x^(1/3); I have problems to graph the function x=[-10:1:10] y = x^(1/3); plot(x,y); the same in plot of google gives a graph but scilab does anything with the sign. H?ctor L?pez -------------- next part -------------- An HTML attachment was scrubbed... URL: From lepezhr at yahoo.com.ar Fri Mar 9 13:58:15 2018 From: lepezhr at yahoo.com.ar (Hector R. Lepez) Date: Fri, 9 Mar 2018 12:58:15 +0000 (UTC) Subject: [Scilab-users] Help me plot y = x^(1/3); In-Reply-To: References: <1473516067.13167106.1520593068556.ref@mail.yahoo.com> <1473516067.13167106.1520593068556@mail.yahoo.com> Message-ID: <820186375.13246482.1520600295586@mail.yahoo.com> This solution worked, could you explain how you put it together? What does it mean to put a sign (x) * ?? Thank you H?ctor L?pez Cel: 154137337 El viernes, 9 de marzo de 2018 09:25:18 ART, St?phane Mottelet escribi?: Le 09/03/2018 ? 11:57, Hector R. Lepez a ?crit?: I have problems to graph the function x=[-10:1:10] y = x^(1/3); plot(x,y); the same in plot of google gives a graph but scilab does anything with the sign. H?ctor L?pez _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users By default in scilab non-integer power elevation of negative number is made in the complex domain. The particular case of the real cubic root can be treated like this : y = sign(x).*abs(x)^(1/3); otherwise (-1)^(1/3) gives you one of the complex cubic root of -1, namely exp(%i*%pi*1/3) S. -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet _______________________________________________ users mailing list users at lists.scilab.org http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Fri Mar 9 14:03:49 2018 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Fri, 9 Mar 2018 13:03:49 +0000 Subject: [Scilab-users] Help me plot y = x^(1/3); Message-ID: Hello Hector, > De : Hector R. Lepez > Envoy? : vendredi 9 mars 2018 13:58 > > What does it mean to put a sign (x) * ?? You'll find some answers here: https://help.scilab.org/docs/6.0.1/en_US/sign.html https://help.scilab.org/docs/6.0.1/en_US/hat.html https://help.scilab.org/docs/6.0.1/en_US/dot.html 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 florent.halgatte at gmail.com Fri Mar 9 22:29:01 2018 From: florent.halgatte at gmail.com (Florent) Date: Fri, 9 Mar 2018 14:29:01 -0700 (MST) Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: <1520362018170-0.post@n3.nabble.com> References: <1520362018170-0.post@n3.nabble.com> Message-ID: <1520630941799-0.post@n3.nabble.com> Good morning, Here is an update of my "problem". First of all I had a problem of compatibility between Scilab/Xcos and the toolbox CPGE. Before I was using Scilab 6.0.1 and CPGE Version: 1.6.0. Now, with Scilab 5.4.1 and CPGE Version: 1.5.2, the Bode plot works. But it doesn't give the result expected... For the Displacement vs Frequency, a Low Pass behavior is expected. For the Velocity vs Frequency, a Band Pass behavior is expected. For the Acceleration vs Frequency, a High Pass behavior is expected (similar to the pressure response). And here is what is get: -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From cfuttrup at gmail.com Sat Mar 10 09:51:11 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sat, 10 Mar 2018 09:51:11 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: <1520630941799-0.post@n3.nabble.com> References: <1520362018170-0.post@n3.nabble.com> <1520630941799-0.post@n3.nabble.com> Message-ID: Hi Florent Thank you for sharing your experiences with us. I don't speak but maybe 100 words in French and can't string together a sentence (but maybe the most elementary tourist-french). I looked up the CPGE in Atoms, but it looks French to me... https://atoms.scilab.org/toolboxes/CPGE/1.6.3/ You give the result but no details to replicate, etc. I'm afraid I can't help. Best regards, Claus On Fri, Mar 9, 2018 at 10:29 PM, Florent wrote: > Good morning, > > Here is an update of my "problem". > > First of all I had a problem of compatibility between Scilab/Xcos and the > toolbox CPGE. > > Before I was using Scilab 6.0.1 and CPGE Version: 1.6.0. > Now, with Scilab 5.4.1 and CPGE Version: 1.5.2, the Bode plot works. But it > doesn't give the result expected... > > For the Displacement vs Frequency, a Low Pass behavior is expected. > For the Velocity vs Frequency, a Band Pass behavior is expected. > For the Acceleration vs Frequency, a High Pass behavior is expected > (similar > to the pressure response). > > And here is what is get: > > > > > > -- > Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing- > Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From cfuttrup at gmail.com Sat Mar 10 20:09:42 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Sat, 10 Mar 2018 20:09:42 +0100 Subject: [Scilab-users] Xcos Loudspeaker simulation Algebraic loop and Bode plot In-Reply-To: References: <1520362018170-0.post@n3.nabble.com> <1520630941799-0.post@n3.nabble.com> Message-ID: Hi Florent (CC Scilab Users) I cannot run the zcos model with my current Scilab 6.0.1 setup (and I haven't installed CPGE 1.6.3 either - is there a way you could simulate a loudspeaker without the need for CPGE??). The bode plots you get is a function of the input + the loudspeaker model. We know that in principle the accelaration is equivalent to Loudspeaker Sound Pressure Level (SPL) response and it should show a flat area + roll off at both ends. I see that your plots don't do this. The velocity is like an integration and should "shift" (or "tilt") the response by 6 dB. The position is also like an integration and should further shift/tilt the reponse additionally 6 dB so that in the bandpass you get a 12 dB/octave slope. Your current loudspeaker model looks OK (I haven't dug into the details of each block), so my guess is - the output depends on the input. Now that you presumably have a working loudspeaker model - would it then maybe make sense to generate a very simple model - e.g. just a resistor (i.e. not a bandpass - but something that behaves from DC to lightning as a flat curve) - play around with the input and get it to behave correctly? Best regards, Claus On Sat, Mar 10, 2018 at 9:51 AM, Claus Futtrup wrote: > Hi Florent > > Thank you for sharing your experiences with us. I don't speak but maybe > 100 words in French and can't string together a sentence (but maybe the > most elementary tourist-french). I looked up the CPGE in Atoms, but it > looks French to me... https://atoms.scilab.org/toolboxes/CPGE/1.6.3/ > > You give the result but no details to replicate, etc. I'm afraid I can't > help. > > Best regards, > Claus > > > On Fri, Mar 9, 2018 at 10:29 PM, Florent > wrote: > >> Good morning, >> >> Here is an update of my "problem". >> >> First of all I had a problem of compatibility between Scilab/Xcos and the >> toolbox CPGE. >> >> Before I was using Scilab 6.0.1 and CPGE Version: 1.6.0. >> Now, with Scilab 5.4.1 and CPGE Version: 1.5.2, the Bode plot works. But >> it >> doesn't give the result expected... >> >> For the Displacement vs Frequency, a Low Pass behavior is expected. >> For the Velocity vs Frequency, a Band Pass behavior is expected. >> For the Acceleration vs Frequency, a High Pass behavior is expected >> (similar >> to the pressure response). >> >> And here is what is get: >> >> >> >> >> >> -- >> Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists- >> Archives-f2602246.html >> _______________________________________________ >> 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 Tue Mar 13 18:18:54 2018 From: samuel.enibe at unn.edu.ng (Samuel Enibe) Date: Tue, 13 Mar 2018 18:18:54 +0100 Subject: [Scilab-users] Function to read Excel xlsx file into SCILAB Message-ID: Dear SCILAB users, is there a direct function available in SCILAB to read an Excel .xlsx file (Excel 2007 or later) into SCILAB. I know I could open it in Excel and use the SaveAs menu to first convert it to .xls file, but a direct SCILAB function would be faster when handling many files. Thank you very much. God bless you. Samuel Enibe University of Nigeria, Nsukka, Nigeria -------------- next part -------------- An HTML attachment was scrubbed... URL: From jrafaelbguerra at hotmail.com Tue Mar 13 19:59:30 2018 From: jrafaelbguerra at hotmail.com (Rafael Guerra) Date: Tue, 13 Mar 2018 18:59:30 +0000 Subject: [Scilab-users] Function to read Excel xlsx file into SCILAB In-Reply-To: References: Message-ID: I am not aware of such direct function. A trick that I find quite useful is the possibility to copy data from Excel directly into Scilab via the Windows clipboard. See ones example below. // First copy an Excel data range to clipboard then run Scilab commands: str = clipboard("paste"); M = evstr(strsplit(str,ascii(10))); Regards, Rafael From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Enibe Sent: Tuesday, March 13, 2018 6:19 PM To: users at lists.scilab.org Subject: [Scilab-users] Function to read Excel xlsx file into SCILAB Dear SCILAB users, is there a direct function available in SCILAB to read an Excel .xlsx file (Excel 2007 or later) into SCILAB. I know I could open it in Excel and use the SaveAs menu to first convert it to .xls file, but a direct SCILAB function would be faster when handling many files. Thank you very much. God bless you. Samuel Enibe University of Nigeria, Nsukka, Nigeria -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Fri Mar 16 14:40:54 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 16 Mar 2018 14:40:54 +0100 Subject: [Scilab-users] Function to read Excel xlsx file into SCILAB In-Reply-To: References: Message-ID: Hello Samuel, https://atoms.scilab.org/toolboxes/xls_link is the module likely the most fitting to your need. There is no indication in the description about the support to .xlsx. But if you did not do it yet, you may try using it. Best regards Samuel Gougeon Le 13/03/2018 ? 18:18, Samuel Enibe a ?crit : > Dear SCILAB users, > > is there a direct function available in SCILAB to read an Excel .xlsx > file (Excel 2007 or later) into SCILAB. > > I know I could open it in Excel and use the SaveAs menu to first > convert it to .xls file, but a direct SCILAB function would be faster > when handling many files. > > Thank you very much. > > God bless you. > > Samuel Enibe > University of Nigeria, Nsukka, Nigeria > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users Hello Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Fri Mar 16 18:23:28 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 16 Mar 2018 18:23:28 +0100 Subject: [Scilab-users] Atoms login In-Reply-To: References: <04dc9350-1369-755e-6fd4-cc7ca60fd3ea@free.fr> <64253f89-c22f-358b-47f2-9bc0e61bc07c@free.fr> <6392edd5998b556a10872b645b32ee76@scilab-enterprises.com> <6022fd3b-9efa-193f-7f4e-49db84059c29@free.fr> <3c08316a-5a26-51cf-f5b8-3d254ddf4249@tritytech.com> <136360ed-18e9-b1f6-f88e-d56023bc678b@free.fr> Message-ID: Le 09/01/2018 ? 10:35, simon.marchetto at scilab-enterprises.com a ?crit : > Hello, > The ATOMS toolbox cache is once again automatically refreshed every 5 > *minutes*, as before. Your cron daemon is killed again, or intentionally stopped. I pushed and published 5 *days* ago a new version of a module with only macros for Scilab 5.5 and 6.0, packaged as usual. I am still expecting it in the ATOMS GUI after atomsSystemUpdate()... If ESI does not want to maintain this platform in a normal and workable way, i think that the best would be just to tell it. For my part, i do not believe too much on gremlins. Only for France, there are more than 100 universities all with Sciences and technologies departments, big public research institutions as CEA, CNES, CNES, INRA... So, it would be a pain not being able to set a new working platform out of ESI. We can call as well for proposals to implement and home a new platform in Japan, India, Brasil, or Russia. Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From skiba.g at gmail.com Fri Mar 16 20:12:44 2018 From: skiba.g at gmail.com (Grzegorz Skiba) Date: Fri, 16 Mar 2018 20:12:44 +0100 Subject: [Scilab-users] Atoms login In-Reply-To: References: <04dc9350-1369-755e-6fd4-cc7ca60fd3ea@free.fr> <64253f89-c22f-358b-47f2-9bc0e61bc07c@free.fr> <6392edd5998b556a10872b645b32ee76@scilab-enterprises.com> <6022fd3b-9efa-193f-7f4e-49db84059c29@free.fr> <3c08316a-5a26-51cf-f5b8-3d254ddf4249@tritytech.com> <136360ed-18e9-b1f6-f88e-d56023bc678b@free.fr> Message-ID: Same situation with our MicroDAQ 1.2.1v toolbox ( https://atoms.scilab.org/toolboxes/microdaq/1.2.1 ). Pushed even earlier on 05.03.2018 and still not available for Scilab and MicroDAQ users. Regards Grzegorz 2018-03-16 18:23 GMT+01:00 Samuel Gougeon : > Le 09/01/2018 ? 10:35, simon.marchetto at scilab-enterprises.com a ?crit : > > Hello, > The ATOMS toolbox cache is once again automatically refreshed every 5 > *minutes*, as before. > > > Your cron daemon is killed again, or intentionally stopped. > > I pushed and published 5 *days* ago a new version of a module with only > macros for Scilab 5.5 and 6.0, > packaged as usual. I am still expecting it in the ATOMS GUI after > atomsSystemUpdate()... > > If ESI does not want to maintain this platform in a normal and workable > way, > i think that the best would be just to tell it. > For my part, i do not believe too much on gremlins. > > Only for France, there are more than 100 universities all with Sciences > and technologies departments, > big public research institutions as CEA, CNES, CNES, INRA... > So, it would be a pain not being able to set a new working platform out of > ESI. > > We can call as well for proposals to implement and home a new platform in > Japan, India, Brasil, or Russia. > > Samuel > > > > > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Fri Mar 16 20:32:57 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Fri, 16 Mar 2018 20:32:57 +0100 Subject: [Scilab-users] Atoms login In-Reply-To: References: <04dc9350-1369-755e-6fd4-cc7ca60fd3ea@free.fr> <64253f89-c22f-358b-47f2-9bc0e61bc07c@free.fr> <6392edd5998b556a10872b645b32ee76@scilab-enterprises.com> <6022fd3b-9efa-193f-7f4e-49db84059c29@free.fr> <3c08316a-5a26-51cf-f5b8-3d254ddf4249@tritytech.com> <136360ed-18e9-b1f6-f88e-d56023bc678b@free.fr> Message-ID: Le 16/03/2018 ? 20:12, Grzegorz Skiba a ?crit : > Same situation with our MicroDAQ 1.2.1v toolbox > (https://atoms.scilab.org/toolboxes/microdaq/1.2.1 ). Pushed even > earlier on 05.03.2018 and still not available for Scilab and MicroDAQ > users. Thanks for your confirmation, Grzegorz! For sure, i was not claiming anything just for any small package of my own, but about the very dissuasive ATOMS for any contributor and finally user. And for years. Such a platform is a key component for users and contributors, for a huge amount of ressources. It is essential that it is reliable, and at least reachable, and easy to use. It is the aim of opensource. So we do not need to be 100 k.users among 1e6 agreeing about that, to change things to make them working. 1000 shall be enough. Best regards Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From grocer.toolbox at gmail.com Fri Mar 16 20:56:36 2018 From: grocer.toolbox at gmail.com (=?UTF-8?B?w4lyaWMgRHVib2lz?=) Date: Fri, 16 Mar 2018 20:56:36 +0100 Subject: [Scilab-users] Atoms login In-Reply-To: References: <04dc9350-1369-755e-6fd4-cc7ca60fd3ea@free.fr> <64253f89-c22f-358b-47f2-9bc0e61bc07c@free.fr> <6392edd5998b556a10872b645b32ee76@scilab-enterprises.com> <6022fd3b-9efa-193f-7f4e-49db84059c29@free.fr> <3c08316a-5a26-51cf-f5b8-3d254ddf4249@tritytech.com> <136360ed-18e9-b1f6-f88e-d56023bc678b@free.fr> Message-ID: Hello to all I can only share Samuel concerns. For toolbox developpers it can only be disheartening to see their toolbvoxes "being packed" for months and now even years. And even for potential users, the message cannot be conducive to use Scilab.... I do not understand why the situation, which was satisfying until 2 years ago, has been let deteriorated since then. Regards Eric 2018-03-16 20:12 GMT+01:00 Grzegorz Skiba : > Same situation with our MicroDAQ 1.2.1v toolbox ( > https://atoms.scilab.org/toolboxes/microdaq/1.2.1 ). Pushed even earlier > on 05.03.2018 and still not available for Scilab and MicroDAQ users. > > Regards Grzegorz > > > 2018-03-16 18:23 GMT+01:00 Samuel Gougeon : > >> Le 09/01/2018 ? 10:35, simon.marchetto at scilab-enterprises.com a ?crit : >> >> Hello, >> The ATOMS toolbox cache is once again automatically refreshed every 5 >> *minutes*, as before. >> >> >> Your cron daemon is killed again, or intentionally stopped. >> >> I pushed and published 5 *days* ago a new version of a module with only >> macros for Scilab 5.5 and 6.0, >> packaged as usual. I am still expecting it in the ATOMS GUI after >> atomsSystemUpdate()... >> >> If ESI does not want to maintain this platform in a normal and workable >> way, >> i think that the best would be just to tell it. >> For my part, i do not believe too much on gremlins. >> >> Only for France, there are more than 100 universities all with Sciences >> and technologies departments, >> big public research institutions as CEA, CNES, CNES, INRA... >> So, it would be a pain not being able to set a new working platform out >> of ESI. >> >> We can call as well for proposals to implement and home a new platform in >> Japan, India, Brasil, or Russia. >> >> Samuel >> >> >> >> >> >> >> >> _______________________________________________ >> 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 fujimoto2005 at gmail.com Sun Mar 18 17:31:42 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Sun, 18 Mar 2018 09:31:42 -0700 (MST) Subject: [Scilab-users] n dimensional normal distribution Message-ID: <1521390702445-0.post@n3.nabble.com> Let x1, x2, .., xn be an random variable of n dimensional normal distribution. Is there any function that gives the probability of {x1> = k1}&{x2> = k2}&...&{xn> = kn}? -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From Christophe.Dang at sidel.com Mon Mar 19 10:03:38 2018 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Mon, 19 Mar 2018 09:03:38 +0000 Subject: [Scilab-users] {EXT} n dimensional normal distribution In-Reply-To: <1521390702445-0.post@n3.nabble.com> References: <1521390702445-0.post@n3.nabble.com> Message-ID: Hello, > De : fujimoto2005 > Envoy? : dimanche 18 mars 2018 17:32 > > Let x1, x2, .., xn be an random variable of n dimensional normal distribution. > Is there any function that gives the probability of {x1> = k1}&{x2> = k2}&...&{xn> = kn}? I'm not sure I understand well your need. You can have the cumulative distribution function of the normal law with cdfnor() https://help.scilab.org/docs/6.0.1/en_US/cdfnor.html In your case, the syntax is [P,Q]=cdfnor("PQ",X,Mean,Std) So I guess the probability you're looking for would be the product of the Qs for all your single component. Additional functions can be found in the CASCI Atoms module: http://atoms.scilab.org/toolboxes/casci HTH -- Christophe Dang Ngoc Chan Mechanical calculation engineer This e-mail may contain confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error), please notify the sender immediately and destroy this e-mail. Any unauthorized copying, disclosure or distribution of the material in this e-mail is strictly forbidden. From paul.carrico at free.fr Mon Mar 19 11:35:54 2018 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Mon, 19 Mar 2018 11:35:54 +0100 Subject: [Scilab-users] unable to read a hdf5 file under windows Message-ID: <82650d3fc2f0f99114decc654ae5e5e3@free.fr> Hi All I'm still wondering why I cannot open a hdf5 file under windows, while it works find under linux (I'm using the latest scilab release in both cases - hdf5 file was made using h5py librarie) The path & file name are corect, so ????? Any trouble under windows? Paul #################### PATH = get_absolute_file_path("test_read_hdf5.sce"); fichier = "Results_gzip.h5"; resu = h5open(PATH + fichier); h5open : Impossible d'ajouter des donn?es au fichier (non HDF5) : D:\***\Results_gzip.h5. Description HDF5: unable to open file: name = 'D:\***\Results_gzip.h5', errno = 2, error message = 'No such file or directory', flags = 0, o_flags = 0. -------------- next part -------------- An HTML attachment was scrubbed... URL: From fujimoto2005 at gmail.com Mon Mar 19 13:15:14 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Mon, 19 Mar 2018 05:15:14 -0700 (MST) Subject: [Scilab-users] n dimensional normal distribution In-Reply-To: <1521390702445-0.post@n3.nabble.com> References: <1521390702445-0.post@n3.nabble.com> Message-ID: <1521461714506-0.post@n3.nabble.com> Dear Chan Thank for your reply. " the product of the Qs for all your single component. " is a right answer only when xis are independent. I am looking for the function when they are not independent. Best regards. -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From amonmayr at laas.fr Mon Mar 19 13:00:54 2018 From: amonmayr at laas.fr (amonmayr at laas.fr) Date: Mon, 19 Mar 2018 13:00:54 +0100 Subject: [Scilab-users] unable to read a hdf5 file under windows In-Reply-To: <82650d3fc2f0f99114decc654ae5e5e3@free.fr> References: <82650d3fc2f0f99114decc654ae5e5e3@free.fr> Message-ID: Le 19/03/2018 ? 11:35, paul.carrico at free.fr a ?crit?: > > Hi All > > I'm still wondering why I cannot open a hdf5 file under windows, while > it works find under linux > > (I'm using the latest scilab release in both cases - hdf5 file was > made using h5py librarie) > > The path & file name are corect, so ????? > > Any trouble under windows? > > Paul > > #################### > > PATH = get_absolute_file_path("test_read_hdf5.sce"); > fichier = "Results_gzip.h5"; > resu = h5open(PATH + fichier); > > h5open : Impossible d'ajouter des donn?es au fichier (non HDF5) : > D:\***\Results_gzip.h5. > Description HDF5: unable to open file: name = > 'D:\***\Results_gzip.h5', errno = 2, error message = 'No such file or > directory', flags = 0, o_flags = 0. > > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users That's weird. Can you share a small hdf5 file that exhibits this problem? I might try to read it from linux/windows. The only windows-specific hdf5 bug I am aware of is related to unknown types for some integers. Are you sure your problem is not something related to the difference in path definition & handling between linux and windows? Antoine -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.carrico at free.fr Mon Mar 19 13:42:48 2018 From: paul.carrico at free.fr (paul.carrico at free.fr) Date: Mon, 19 Mar 2018 13:42:48 +0100 Subject: [Scilab-users] unable to read a hdf5 file under windows In-Reply-To: References: <82650d3fc2f0f99114decc654ae5e5e3@free.fr> Message-ID: I finally found what went wrong? a latin letter in the path with '?' ? Now it works ... sorry Paul Le 2018-03-19 13:00, amonmayr at laas.fr a ?crit : > Le 19/03/2018 ? 11:35, paul.carrico at free.fr a ?crit : > >> Hi All >> >> I'm still wondering why I cannot open a hdf5 file under windows, while it works find under linux >> >> (I'm using the latest scilab release in both cases - hdf5 file was made using h5py librarie) >> >> The path & file name are corect, so ????? >> >> Any trouble under windows? >> >> Paul >> >> #################### >> >> PATH = get_absolute_file_path("test_read_hdf5.sce"); >> fichier = "Results_gzip.h5"; >> resu = h5open(PATH + fichier); >> >> h5open : Impossible d'ajouter des donn?es au fichier (non HDF5) : D:\***\Results_gzip.h5. >> Description HDF5: unable to open file: name = 'D:\***\Results_gzip.h5', errno = 2, error message = 'No such file or directory', flags = 0, o_flags = 0. >> >> _______________________________________________ >> users mailing list >> users at lists.scilab.org >> http://lists.scilab.org/mailman/listinfo/users > > That's weird. > Can you share a small hdf5 file that exhibits this problem? > I might try to read it from linux/windows. > The only windows-specific hdf5 bug I am aware of is related to unknown types for some integers. > Are you sure your problem is not something related to the difference in path definition & handling between linux and windows? > > Antoine > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From Christophe.Dang at sidel.com Mon Mar 19 14:12:23 2018 From: Christophe.Dang at sidel.com (Dang Ngoc Chan, Christophe) Date: Mon, 19 Mar 2018 13:12:23 +0000 Subject: [Scilab-users] {EXT} Re: n dimensional normal distribution In-Reply-To: <1521461714506-0.post@n3.nabble.com> References: <1521390702445-0.post@n3.nabble.com> <1521461714506-0.post@n3.nabble.com> Message-ID: Hello, > De : fujimoto2005 > Envoy? : lundi 19 mars 2018 13:15 > > " the product of the Qs for all your single component. " is a right answer only when xis are independent. > I am looking for the function when they are not independent. I'm afraid I can't help further; Maybe if you can share a small example? Otherwise, if you have strong needs I statistics, you may have a look at the R language https://www.r-project.org/ although the syntax is a bit weird when you come from Scilab. 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 stephane.mottelet at utc.fr Mon Mar 19 14:37:33 2018 From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=) Date: Mon, 19 Mar 2018 14:37:33 +0100 Subject: [Scilab-users] n dimensional normal distribution In-Reply-To: <1521461714506-0.post@n3.nabble.com> References: <1521390702445-0.post@n3.nabble.com> <1521461714506-0.post@n3.nabble.com> Message-ID: <452279b1-b1ed-0083-e85f-fe0885dd3d93@utc.fr> Le 19/03/2018 ? 13:15, fujimoto2005 a ?crit?: > Dear Chan > Thank for your reply. > > " the product of the Qs for all your single component. " is a right answer > only when xis are independent. > I am looking for the function when they are not independent. > > Best regards. > > > > -- > Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users Hello, R implements this in? the "TruncatedNormal" package (https://cran.r-project.org/web/packages/TruncatedNormal/) : The routines include (Quasi-) Monte Carlo estimator and a deterministic upper bound of the cumulative distribution function of the multivariate normal. hth S. -- St?phane Mottelet Ing?nieur de recherche EA 4297 Transformations Int?gr?es de la Mati?re Renouvelable D?partement G?nie des Proc?d?s Industriels Sorbonne Universit?s - Universit? de Technologie de Compi?gne CS 60319, 60203 Compi?gne cedex Tel : +33(0)344234688 http://www.utc.fr/~mottelet From sgougeon at free.fr Mon Mar 19 20:51:21 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Mon, 19 Mar 2018 20:51:21 +0100 Subject: [Scilab-users] n dimensional normal distribution In-Reply-To: <1521390702445-0.post@n3.nabble.com> References: <1521390702445-0.post@n3.nabble.com> Message-ID: <1ddb98fa-11cf-e90b-ce5a-25854c3cdbf1@free.fr> Hello Masahiro, Le 18/03/2018 ? 17:31, fujimoto2005 a ?crit : > Let x1, x2, .., xn be an random variable of n dimensional normal > distribution. Is there any function that gives the probability of {x1> = > k1}&{x2> = k2}&...&{xn> = kn}? I thougth at the first sight that it would be possible to get this probability with some lexicographic sorting, but it's not the case. I don't thing that we can avoid an explicit loop. Here is a possible direct calculation, from a list of actual samples : ns = 20; // number of samples nd = 3; // number of dimensions k = [5 3 7]; m = grand(ns,nd,"nor",8,3) for i = 1:nd [?,r] = gsort(m(:,i)); m = m(r,:); m = m(find(m(:,i)>=k(i)),:) end p = size(m,1)/ns // requested probability Example of run : m = 12.750144 7.1200267 5.900484 12.575508 5.410083 10.976422 12.399993 3.7475677 10.092495 12.220889 5.4940195 6.1479044 12.158487 5.2062957 11.651957 11.611694 6.6661928 5.9750098 10.718146 2.2454739 11.737011 10.343892 4.2714818 4.7199587 9.9907016 7.7903253 5.9802778 8.3135823 4.8094984 7.1769228 8.0939865 9.2484944 13.215993 8.0098652 3.0198012 7.5767533 7.5354006 10.715856 8.985266 7.4971339 17.821625 5.5456382 6.6502306 8.7791304 8.812858 6.5728805 12.299302 7.9823783 6.2940806 10.376389 7.8221558 5.2862072 6.0566186 11.102784 m = 7.4971339 17.821625 5.5456382 6.5728805 12.299302 7.9823783 7.5354006 10.715856 8.985266 6.2940806 10.376389 7.8221558 8.0939865 9.2484944 13.215993 6.6502306 8.7791304 8.812858 9.9907016 7.7903253 5.9802778 12.750144 7.1200267 5.900484 11.611694 6.6661928 5.9750098 5.2862072 6.0566186 11.102784 12.220889 5.4940195 6.1479044 12.575508 5.410083 10.976422 12.158487 5.2062957 11.651957 8.3135823 4.8094984 7.1769228 10.343892 4.2714818 4.7199587 m = 8.0939865 9.2484944 13.215993 12.158487 5.2062957 11.651957 5.2862072 6.0566186 11.102784 12.575508 5.410083 10.976422 7.5354006 10.715856 8.985266 6.6502306 8.7791304 8.812858 --> p = size(m,1)/ns // requested probability p = 0.3 Samuel -------------- next part -------------- An HTML attachment was scrubbed... URL: From lepezhr at yahoo.com.ar Wed Mar 21 13:18:35 2018 From: lepezhr at yahoo.com.ar (Hector R. Lepez) Date: Wed, 21 Mar 2018 12:18:35 +0000 (UTC) Subject: [Scilab-users] Problem with SCILAB 5 and 6 Version References: <1251921910.4639618.1521634715629.ref@mail.yahoo.com> Message-ID: <1251921910.4639618.1521634715629@mail.yahoo.com> In SCILAB 6 this code: function ingres() A=input ('Ingrese la matriz A '); disp(A); b=input ('Ingrese el vector independiente b ');disp(b); endfunction in console, with 6 version of SCILAB Input function faild, i don't know why?? because in older version works with no problems. --> ingres Ingrese la matriz A? ? [4,3,-2,1;3,2,1,5;-2,3,1,2;-5,0,1,1] ?? 4.?? 3.? -2.?? 1. ?? 3.?? 2.?? 1.?? 5. ? -2.?? 3.?? 1.?? 2. ? -5.?? 0.?? 1.?? 1. Ingrese el vector independiente b? ? --> [4;-8;-7;-8] Ingrese el vector independiente b? ? [4;-8;-7;-8] ? -8. H?ctor L?pez Cel: 154137337 -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgougeon at free.fr Wed Mar 21 14:02:42 2018 From: sgougeon at free.fr (Samuel Gougeon) Date: Wed, 21 Mar 2018 14:02:42 +0100 Subject: [Scilab-users] Problem with SCILAB 5 and 6 Version In-Reply-To: <1251921910.4639618.1521634715629@mail.yahoo.com> References: <1251921910.4639618.1521634715629.ref@mail.yahoo.com> <1251921910.4639618.1521634715629@mail.yahoo.com> Message-ID: <5ca96189-0e6f-5ed5-6164-51f45aa59b5d@free.fr> Hello H?ctor Are you using the latest Scilab 6.0.1 release? input() was bugged in previous Scilab 6.0.0+/- releases, but it is now fixed: --> function ingres() > > A=input ('Ingrese la matriz A '); > > disp(A); > > b=input ('Ingrese el vector independiente b '); > > disp(b); > > > endfunction --> ingres Ingrese la matriz A [4,3,-2,1;3,2,1,5;-2,3,1,2;-5,0,1,1] 4. 3. -2. 1. 3. 2. 1. 5. -2. 3. 1. 2. -5. 0. 1. 1. Ingrese el vector independiente b [4;-8;-7;-8] 4. -8. -7. -8. HTH Samuel Le 21/03/2018 ? 13:18, Hector R. Lepez a ?crit : > In SCILAB 6 this code: > > function ingres() > A=input ('Ingrese la matriz A '); > disp(A); > b=input ('Ingrese el vector independiente b '); > disp(b); > > endfunction > > > in console, with 6 version of SCILAB Input function faild, i don't > know why?? because in older version works with no problems. > > --> ingres > Ingrese la matriz A > [4,3,-2,1;3,2,1,5;-2,3,1,2;-5,0,1,1] > > 4. 3. -2. 1. > 3. 2. 1. 5. > -2. 3. 1. 2. > -5. 0. 1. 1. > Ingrese el vector independiente b > --> [4;-8;-7;-8] > Ingrese el vector independiente b > [4;-8;-7;-8] > > > -8. > > > > H?ctor L?pez Cel: 154137337 > > > _______________________________________________ > users mailing list > users at lists.scilab.org > http://lists.scilab.org/mailman/listinfo/users -------------- next part -------------- An HTML attachment was scrubbed... URL: From lepezhr at yahoo.com.ar Wed Mar 21 14:14:02 2018 From: lepezhr at yahoo.com.ar (Hector R. Lepez) Date: Wed, 21 Mar 2018 13:14:02 +0000 (UTC) Subject: [Scilab-users] Problem with SCILAB 5 and 6 Version In-Reply-To: <5ca96189-0e6f-5ed5-6164-51f45aa59b5d@free.fr> References: <1251921910.4639618.1521634715629.ref@mail.yahoo.com> <1251921910.4639618.1521634715629@mail.yahoo.com> <5ca96189-0e6f-5ed5-6164-51f45aa59b5d@free.fr> Message-ID: <832802185.4672995.1521638042616@mail.yahoo.com> YES SAMUEL!, thats the problem, i update SCILAB and fixed. Thanks!. H?ctor L?pez Cel: 154137337 El mi?rcoles, 21 de marzo de 2018 10:03:36 ART, Samuel Gougeon escribi?: Hello H?ctor Are you using the latest Scilab 6.0.1 release? input() was bugged in previous Scilab 6.0.0+/- releases, but it is now fixed: --> function ingres() ? > ? > A=input ('Ingrese la matriz A?? '); ? > ? > disp(A); ? > ? > b=input ('Ingrese el vector independiente b?? '); ? > ? > disp(b); ? > ? > ? > endfunction --> ingres Ingrese la matriz A?? [4,3,-2,1;3,2,1,5;-2,3,1,2;-5,0,1,1] ?? 4.?? 3.? -2.?? 1. ?? 3.?? 2.?? 1.?? 5. ? -2.?? 3.?? 1.?? 2. ? -5.?? 0.?? 1.?? 1. Ingrese el vector independiente b??? [4;-8;-7;-8] ?? 4. ? -8. ? -7. ? -8. HTH Samuel Le 21/03/2018 ? 13:18, Hector R. Lepez a ?crit?: In SCILAB 6 this code: function ingres() A=input ('Ingrese la matriz A '); disp(A); b=input ('Ingrese el vector independiente b '); disp(b); endfunction in console, with 6 version of SCILAB Input function faild, i don't know why?? because in older version works with no problems. --> ingres Ingrese la matriz A? ? [4,3,-2,1;3,2,1,5;-2,3,1,2;-5,0,1,1] ?? 4.?? 3.? -2.?? 1. ?? 3.?? 2.?? 1.?? 5. ? -2.?? 3.?? 1.?? 2. ? -5.?? 0.?? 1.?? 1. Ingrese el vector independiente b? ? --> [4;-8;-7;-8] Ingrese el vector independiente b? ? [4;-8;-7;-8] ? -8. H?ctor L?pez Cel: 154137337 _______________________________________________ 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 fujimoto2005 at gmail.com Tue Mar 27 00:05:51 2018 From: fujimoto2005 at gmail.com (fujimoto2005) Date: Mon, 26 Mar 2018 15:05:51 -0700 (MST) Subject: [Scilab-users] n dimensional normal distribution In-Reply-To: <1ddb98fa-11cf-e90b-ce5a-25854c3cdbf1@free.fr> References: <1521390702445-0.post@n3.nabble.com> <1ddb98fa-11cf-e90b-ce5a-25854c3cdbf1@free.fr> Message-ID: <1522101951559-0.post@n3.nabble.com> Dear all. Thank you for your answers. I am sorry that the reply is delayed. My original problem is v=??..?f(X1,X2,..,Xn,Y)n(X1,X2,..,Xn,Y)dX1dX2..dXndY where n(X1,X2,..,Xn,Y) is a n+1 dimensional normal density function. This is a problem of a mathematical finance and this integral is usually calculated by numerical integration using random numbers. In my problem f() has the following special form. f(X1,X2,..,Xn,Y)=g(Y)I(x1>=k1)...I(Xn>=Kn) where I() is an indication function. In this case, v can be decomposed as v=?g(Y){?..?n(X1,X2,..,Xn |Y)dX1dX2..dXn}n(Y)dY= ?g(Y)Prob[(X1>=K1)& ..?(Xn>=Kn)|Y}n(Y)dY where n(X1,X2,..,Xn |Y) is a normal density conditiona to Y, Prob(..|Y) is a probability conditional to Y and n(Y) is a one dimensional normal density function. This integration can be calculated firstly when n = 1 since Prob(x1>=K1|Y) has an analytical approximate function. When n = 2, it works still fine because Prob[(X1>=K1 &X2>=K2)|Y}=1-Prob[X1<=K1]-Prob[X2<=K2]+Prob[X1<=K1&X2<=K2] and I hear that Prob[X1<=K1&X2<=K2] have an analytical approximate functions. I posted the question because I can calculate it first if there is an analytical approximate function of Prob[(X1>=K1)& ..,(Xn>=Kn)|Y} for large n. However, as Samuel demonstrated explicitly and St?phane implicitly by referring to R, it seems that there is no choice but to calculate Prob[(X1>=K1)& ..,(Xn>=Kn)|Y} for large n by using random numbers. I understand my prospects doesn?t go well because there is no such function for large n. When using random numbers, it is the best way to calculate v in the original form. Best regards. -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html From cfuttrup at gmail.com Wed Mar 28 16:54:22 2018 From: cfuttrup at gmail.com (Claus Futtrup) Date: Wed, 28 Mar 2018 16:54:22 +0200 Subject: [Scilab-users] FEA and geometry Message-ID: Hi there It's Easter soon and I thought, maybe I could play around a bit with Scilab and figure out how to read a geometry and do some Finite Element Analysis? Maybe not finished during Easter ... first step would be to read a geometry from e.g. a CAD file and I wonder if Scilab supports this? I found this ATOMS module for reading STL files: https://fileexchange.scilab.org/toolboxes/490000 Are there any other options out there? ... or, how would you approach such a task? Second I found CGLab, which interfaces to CGAL (Computational Geometry Algorithms Library) which seems like it could analyze and define a mesh... but then again, an STL file in itself is like a mesh. I'm curious for input, to get started in the right direction. Best regards, Claus -------------- next part -------------- An HTML attachment was scrubbed... URL: