Re: [Scilab-users] Compatibility With

2012-10-21 Thread Stéphane Mottelet
Hello, Seems to me that this problem is the same as the one reported here : http://bugzilla.scilab.org/show_bug.cgi?id=11994 Le 22 oct. 2012 à 05:39, jtibau a écrit : > whew... finally figured it out. > > The problem is actually scinote's, not scilab's. I don't have a Lion machine > to check

Re: [Scilab-users] root calculation

2012-11-05 Thread Stéphane Mottelet
Le 05/11/12 14:53, Adrien Vogt-Schilb a écrit : On 05/11/2012 14:45, Paul Carrico wrote: All, Thanks for the feedback ... ... a basic plot shows the curve is monotonic (see attachment) ; I should plot a 3D curve in order to study the influence of the exponents for example to verify there's on

Re: [Scilab-users] How to plot 3d implicit surface

2012-11-23 Thread Stéphane Mottelet
Le 23/11/12 09:14, Orbeman a écrit : Hello, I would like to plot a 3d implicit surface defines with an equation like f(x, y, z)=0. There's no parametric transformation and is not possible to express for example z=g(x, y). Does anyone have an idea ? Thank a lot. -- View this message in conte

Re: [Scilab-users] u .^ k ?

2012-12-19 Thread Stéphane Mottelet
Le 19/12/12 11:17, Lamy Alain a écrit : Hi, Does anyone know how to perform efficiently an operation like : u .^k with: u = [u1, u2, ..., un] // ui : real numbers and: k = [k1; k2; kp] // ki : positive integers The result being : [ u1^k1, u2^k1, ..., un^k1; u1^k2, u2^k2, ..., un^k2; ...

Re: [Scilab-users] Problem with ATOMS

2013-01-28 Thread Stéphane Mottelet
Hello, this problem is fixed in the current developpement version of the plotlib, a new version will be released this week. S. Le 28/01/13 08:52, Vincent COUVERT a écrit : Hello, It seems that this error is due to plotlib since the error message indicates that the error occurred at line 46

Re: [Scilab-users] Scilab is not a pipe

2013-02-27 Thread Stéphane Mottelet
Le 26/02/13 17:45, Adrien Vogt-Schilb a écrit : hi is there any reasons why I can't type from a console: scilab myfile.sce and expect scilab to execute myfile.sce and redirect output to the standard output? I can almost to that with scilab-cli -f and scilab-cli -e but i still have to exit

Re: [Scilab-users] Surface smoothing in Scilab, immune to outliers

2013-03-04 Thread Stéphane Mottelet
Hello, Replacing the squared L2 norm by the L1 norm in the linear regression gives a good robustness to outliers (cf. Donoho and al. papers). The problem is then non differentiable but you can implement it by iteratively reweighting the classical L2 method (IRLS method), or by writing an equi

Re: [Scilab-users] Surface smoothing in Scilab, immune to outliers

2013-03-04 Thread Stéphane Mottelet
onvolution does not seem, a priori, to be the most effective method to remove them. Regards, Rafael -Original Message- From: users-boun...@lists.scilab.org [mailto:users-boun...@lists.scilab.org] On Behalf Of Stéphane Mottelet Sent: Monday, March 04, 2013 10:14 AM To: users@lists.s

[Scilab-users] parallel_run under OSX

2013-03-06 Thread Stéphane Mottelet
Hello, is parallel_run supposed to work under MacOS X (tested here with Scilab 5.4.0) ? -->function a=g(arg1) --> a=arg1*arg1 -->endfunction --> -->res=parallel_run(1:10, g); A previous error has been detected while loading libsciparallel.dylib: !--error 999 Imposs

Re: [Scilab-users] parallel_run under OSX

2013-03-06 Thread Stéphane Mottelet
Le 06/03/13 15:30, Stéphane Mottelet a écrit : Hello, is parallel_run supposed to work under MacOS X (tested here with Scilab 5.4.0) ? -->function a=g(arg1) --> a=arg1*arg1 -->endfunction --> -->res=parallel_run(1:10, g); A previous error has been detecte

Re: [Scilab-users] parallel_run under OSX

2013-03-06 Thread Stéphane Mottelet
Le 06/03/13 15:30, Stéphane Mottelet a écrit : Hello, is parallel_run supposed to work under MacOS X (tested here with Scilab 5.4.0) ? -->function a=g(arg1) --> a=arg1*arg1 -->endfunction --> -->res=parallel_run(1:10, g); A previous error has been detected while loading libsc

Re: [Scilab-users] parallel_run under OSX

2013-03-08 Thread Stéphane Mottelet
Le 06/03/13 18:27, Sylvestre Ledru a écrit : Le 06/03/2013 15:43, Stéphane Mottelet a écrit : Le 06/03/13 15:30, Stéphane Mottelet a écrit : Hello, is parallel_run supposed to work under MacOS X (tested here with Scilab 5.4.0) ? There seems to be ideas about a fix for Darwin : http

[Scilab-users] declaration of uint8 matrix

2013-03-13 Thread Stéphane Mottelet
Hello, Is it possible to declare a uint8 matrix without first declaring it as a double ? Having to write m=uint8(zeros(n,n)) makes it impossible to declare a matrix that would fit in the stack as a uint8 but not as a double (8 times bigger). S. _

Re: [Scilab-users] declaration of uint8 matrix

2013-03-13 Thread Stéphane Mottelet
Le 13/03/13 18:57, sgoug...@free.fr a écrit : Hello Stéphane, De: "Stéphane Mottelet" Is it possible to declare a uint8 matrix without first declaring it as a double ? ../.. Yes, you can do: -->ui = resize_matrix(uint8(0),5,3) ui = 0 0 0 0 0 0 0 0 0 0

Re: [Scilab-users] Data Extraction Issue in Scilab

2013-03-14 Thread Stéphane Mottelet
Hello, this small script solves your problem : fd=mopen('engine.in','r'); s=mgetl(fd); mclose(fd); data1=[]; for i=1:size(s,1) w=strsplit(s(i),'/\s+/'); for j=1:size(w,1) data1(i,j)=w(j); end end [r,c]=find(data=='bore'); bore=eval(da

Re: [Scilab-users] Data Extraction Issue in Scilab

2013-03-14 Thread Stéphane Mottelet
Sorry for the small typing errors ("data" instead of "data1") : fd=mopen('engine.in','r'); s=mgetl(fd); mclose(fd); data1=[]; for i=1:size(s,1) w=strsplit(s(i),'/\s+/'); for j=1:size(w,1) data1(i,j)=w(j); end end [r,c]=find(data1=='bore'); bore=eva

Re: [Scilab-users] Data Extraction Issue in Scilab

2013-03-14 Thread Stéphane Mottelet
Sorry for the small typing errors ("data" instead of "data1") : fd=mopen('engine.in','r'); s=mgetl(fd); mclose(fd); data1=[]; for i=1:size(s,1) w=strsplit(s(i),'/\s+/'); for j=1:size(w,1) data1(i,j)=w(j); end end [r,c]=find(data1=='bore'); bore=eva

Re: [Scilab-users] Data Extraction Issue in Scilab

2013-03-14 Thread Stéphane Mottelet
ta.%s=%s',w(2),w(1))); end -->data data = "" 0.08 stroke: 0.08 squish: 0.001671 conrod: 0.14835 wrist: 0 rpm: 4000 swirl: 1.7 swirl_profile: 3.11 zhead: 0 pist_id: 1 liner_id: 4 head_id: 5 crevice_flag: 0 S.

Re: [Scilab-users] Subplot question

2013-05-13 Thread Stéphane Mottelet
Le 13/05/13 16:38, Adrien Vogt-Schilb a écrit : Hi It is possible, just use subplot(321) subplot(322) subplot(323) subplot(324) and then subplot(313) On 13/05/2013 16:30, SCHULZ Wolfgang wrote: Hello, I want plot data and divided the plot region with the subplot command in 6 regions (subpl

Re: [Scilab-users] Subplot question

2013-05-13 Thread Stéphane Mottelet
Le 13/05/13 16:38, Adrien Vogt-Schilb a écrit : Hi It is possible, just use subplot(321) subplot(322) subplot(323) subplot(324) and then subplot(313) On 13/05/2013 16:30, SCHULZ Wolfgang wrote: Hello, I want plot data and divided the plot region with the subplot command in 6 regions (subplo

Re: [Scilab-users] Subplot question

2013-05-16 Thread Stéphane Mottelet
Hello, Le 13/05/13 16:38, Adrien Vogt-Schilb a écrit : Hi It is possible, just use subplot(321) subplot(322) subplot(323) subplot(324) and then subplot(313) then you get margins that are not homogeneous (this is ugly), as you can see it in the attached file scilab_subplot.png (although plot

Re: [Scilab-users] Name of function's input arguments

2013-05-29 Thread Stéphane Mottelet
You can do it with the macrovar primitive. S. Le 29/05/13 22:56, Adrien Vogt-Schilb a écrit : you cannot. But you can provide the function directly with the name of your variable, and transform te name of the variable into the variable by evaluating it... function f(nameofvar1,nameofvar2)

Re: [Scilab-users] 4D gridding and interpolation of scattered data points

2013-06-03 Thread Stéphane Mottelet
Hello, try splin3d or linear_interpn. S. Le 03/06/13 16:21, Ariel a écrit : > Do you mean that he V values you have are not placed on a 3d grid, > but on randomly scattered (x, y, z) points? Exactly.I have randomly scattered (x,y,z) points with corresponding value V and I want to map those

Re: [Scilab-users] Sequential Quadratic Programming in Scilab

2013-08-20 Thread Stéphane Mottelet
Hello, I have updated the fsqp toolbox for Scilab 5.4 an can provide the package (without the cfsqp.c file which has to be obtained from the authors). S. 2013/8/19 Michael BAUDIN > Hi, > > What kind of optimization problem are you willing to solve ? SQP is > typical for nonlinear programming

Re: [Scilab-users] how to get list of object that are present in scilab script file

2013-09-03 Thread Stéphane Mottelet
Le 03/09/13 17:27, Serge Steer a écrit : Le 03/09/2013 12:20, hilife5 a écrit : I want to know is there exist any way to get the list of objects like variable names, matrices, plots that are present in scilab script file. for script file, no but it is possible for functions using the macrovar f

Re: [Scilab-users] how to get list of object that are present in scilab script file

2013-09-03 Thread Stéphane Mottelet
Le 03/09/13 17:27, Serge Steer a écrit : Le 03/09/2013 12:20, hilife5 a écrit : I want to know is there exist any way to get the list of objects like variable names, matrices, plots that are present in scilab script file. for script file, no but it is possible for functions using the macrovar f

Re: [Scilab-users] [Users-fr] plot avec scilab 5.4.1 ss ubuntu 13.04

2013-09-16 Thread Stéphane Mottelet
Le 16 sept. 2013 à 20:50, Antoine Monmayrant a écrit : > Le 16/09/13 17:57, Calixte a écrit : >> On 16/09/2013 15:18, Antoine Monmayrant wrote: >>> On 09/16/2013 02:05 PM, quantpa...@numericable.fr wrote: j'ai executé x=[0:0.1:2*%pi];plot(x,sin(x)); sur 5.4.1 cela se graphe

[Scilab-users] Scilab 5.5.0 beta1 does not start under OSX 10.8

2013-10-02 Thread Stéphane Mottelet
Hello all, The 5.5 beta seems to be broken : # open /Applications/scilab-5.5.0-beta-1.app/ LSOpenURLsWithRole() failed with error -10810 for the file /Applications/scilab-5.5.0-beta-1.app. Trying to lauch the binary gives : # /Applications/scilab-5.5.0-beta-1.app/Contents/MacOS/bin/scilab dy

Re: [Scilab-users] Scilab 5.5.0 beta1 does not start under OSX 10.8

2013-10-04 Thread Stéphane Mottelet
nts/MacOS/bin/scilab: line 911: 14256 Trace/BPT trap: 5 "$SCILABBIN" "$@" S. Le 04/10/13 17:46, Bruno JOFRET a écrit : Hi, Can you try this nightly [1] please ? [1] http://downloadarea.scilab.org/download/2013-10-04/scilab-master-1380893237-x86_64.dmg On 10/02/

Re: [Scilab-users] Scilab 5.5.0 beta1 does not start under OSX 10.8

2013-10-07 Thread Stéphane Mottelet
you try something with "umfpack" ? (It should complain about dlopen) Thx On 10/04/2013 06:01 PM, Stéphane Mottelet wrote: Same behaviour : macmottelet-gi-0:Applications mottelet$ scilab-master-1380893237.app/Contents/MacOS/bin/scilab dyld: Library not loaded: @executable_path/../lib/

Re: [Scilab-users] Solving a ecuation

2013-11-08 Thread Stéphane Mottelet
Hello, spec(Kx,M) will give you all the solutions. S. > Le 8 nov. 2013 à 20:30, Chuox a écrit : > > I need to solve the following equation > > det(Kx-w*M)=0, Were Kx, and M are know matrix, in my texas vollage 200 (R) i > solve it like this: > > Kx= > > >457458.17 - 187960.850.

[Scilab-users] XMLlab 1.7.6 is out !

2014-01-24 Thread Stéphane Mottelet
Dear Scilab users, A new version of XMLlab (1.7.6) has just been released ! For those who never heard about XMLlab, the official site http://xmllab.org which is up again after a fatal server crash just before Christmas, will give you some piece of information about it. Basically, XMLlab can

Re: [Scilab-users] 3d curve plot

2014-02-13 Thread Stéphane Mottelet
Le 13/02/2014 10:51, Serge Steer a écrit : Le 12/02/2014 12:49, Jeibros a écrit : Hi all, I'm totally new to Scilab (I used to work with Matlab), and I have a rather simple question. I have two arrays (x,y). Next, I make a conversion so that to get a new variable u out of x and y. And finall

Re: [Scilab-users] evaluate error on each parameter calculated with leastsq

2014-02-21 Thread Stéphane Mottelet
Here is a simple example where the derivative is approximated with the "derivative" macro of Scilab. S. 2014-02-20 9:43 GMT+01:00 Stéphane Mottelet : > Hello, > > If you can make the hypothesis that your data is corrupted by gaussian > noise, then you can approximate th

Re: [Scilab-users] evaluate error on each parameter calculated with leastsq

2014-02-21 Thread Stéphane Mottelet
Hello, If you can make the hypothesis that your data is corrupted by gaussian noise, then you can approximate the covariance matrix of your estimated parameters. Let p be the vector of parameters and r(p) the residual vector given by r(p)=sigma^(-1)*(y-Y(p)) where y is your measuement vector, Y(

Re: [Scilab-users] evaluate error on each parameter calculated with leastsq

2014-02-21 Thread Stéphane Mottelet
Hello, If you can make the hypothesis that your data is corrupted by gaussian noise, then you can approximate the covariance matrix of your estimated parameters. Let p be the vector of parameters and r(p) the residual vector given by r(p)=sigma^(-1)*(y-Y(p)) where y is your measuement vector, Y(

[Scilab-users] Test

2014-02-21 Thread Stéphane Mottelet
This is just a test, sorry for the disturbance. S. ___ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] Test

2014-02-21 Thread Stéphane Mottelet
This a test, sorry. S. ___ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] Test

2014-02-26 Thread Stéphane Mottelet
Sorry, this is a test... S. ___ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] fmincon won't work

2014-03-14 Thread Stéphane Mottelet
nk you for your attention, Natasha On 26/02/2014 20:25, "Stéphane Mottelet (UTC)" wrote: Hello, Currently, sci_ipopt seems to be broken on Linux and MacOSX when using the atoms package for Scilab 5.4.1. Since fmincon relies on sci_ipot, fmicon is equally broken... BTW, users needing

Re: [Scilab-users] fmincon won't work

2014-04-24 Thread Stéphane Mottelet
Le 23/04/2014 19:48, Gate007 a écrit : Dear users, I have the same trouble in Scilab 5.3.3 Ubuntu 12.04. I have installed IPOPT from repository, then compiled ipopt from sci-ipopt archive, maked it and maked istall. Loading sci-ipopt works fine, without any errors or warnings. Loading fmincon al

[Scilab-users] problem/bug with xmlNs

2014-07-08 Thread Stéphane Mottelet
Hello, xmlNS is not able to create an XML namespace without a prefix, although this is completely correct in XML, as it allows inheritance of it. Example: http://name.space";> 1 here the "b" element is in the "http://name.space"; namespace. The only workaround I have found is the followi

[Scilab-users] Problème avec xmltoformat

2014-07-11 Thread Stéphane Mottelet
Bonjour, Depuis la 5.5.0 j'ai le problème suivant lors de la construction d'un module atoms : Création des fichiers d'aide... Génération du document maître : /Users/mottelet/svn/plotlib/trunk/plotlib/help/en_US Génération du fichier d'aide [javaHelp] in /Users/mottelet/svn/plotlib/trunk/p

[Scilab-users] xsl processor

2014-07-11 Thread Stéphane Mottelet
Hello, Is there any low level access to saxon in Scilab ? Functions of the XML module are fine, but it will be great to be also able to make xsl transformations. S. ___ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listin

Re: [Scilab-users] xsl processor

2014-07-11 Thread Stéphane Mottelet
calls whereas xsltproc doesn't. For example, I have an xsl stylesheet which computes the transitive closure of a graph which runs with Saxon but fails with xsltproc. S. -- Clément Le vendredi 11 juillet 2014 à 17:11 +0200, Stéphane Mottelet a écrit : Hello, Is there any low level acce

Re: [Scilab-users] Scilab leastsq exponential fitting

2014-07-23 Thread Stéphane Mottelet
Le 23/07/2014 11:45, chloe.kykam a écrit : k=6.63e-34*3e8/1.38e-23 x=[1;2;3;4;5;6;7;8;9;10] y=[280;320;369.22772;391.25743;414.74257;439.75248;466.06931;493.60396;523.87129;530] w=[0;0;1;1;1;1;1;1;1;0] function y=yth(x,a) y=a(1)*exp(-k/x/a(2)) endfunction a0=[1.0;1.0] function e=myfun(a,x

Re: [Scilab-users] How to solve sparse nonlinear system by scilabe

2014-08-11 Thread Stéphane Mottelet
Hello, fsolve cannot handle a sparse jacobian, but if you use the Newton's method (less than 5 lines of Scilab code), the sparsity of the jacobian will be taken into account in the linear system resolution step. S. Le 11/08/2014 14:29, Michael J McCann a écrit : Since the network is passive

Re: [Scilab-users] Optimization question

2014-08-13 Thread Stéphane Mottelet
Le 13/08/2014 03:11, ebmb.sci a écrit : Hello All, please, could someone guide me on how to proceed to minimize the following polynomial using functions optim and numdiff. Since the parameters A and B have different constraints (i.e.: A = [0 1] and B = [1000 2000]), how must I sign the values

Re: [Scilab-users] Recognition of cursor key as input

2014-08-13 Thread Stéphane Mottelet
Hello, I presume you are not using scilab 5.5.0, change the function as follows : function key=waitforkeypress() global _keypress set(gcf(),'event_handler','handler'); set(gcf(),'event_handler_enable','on') _keypress=0; while ~_keypress sleep(1) end set(gcf()

Re: [Scilab-users] Recognition of cursor key as input

2014-08-13 Thread Stéphane Mottelet
Le 13/08/2014 15:53, Jens a écrit : It works now on my installation, probably the way you designed it. In my case the graphic window does not open after starting the script by F5 really strange... I think that you mean that the windows, if already open (or minimized, or docked), does not show up

Re: [Scilab-users] Fwd: Help modeling ODE system

2014-09-15 Thread Stéphane Mottelet
Hello, Why use a cannon to kill a fly ? Maybe this would be easier like this : function dvdt=f(t,v) dvdt=[0.01 -0.01;-0.01 0.01]*v; endfunction t=0:0.1:10; v=ode([1;5000],0,t,f); plot(t,v) S. Le 15/09/2014 13:05, Berns Buenaobra a écrit : You could simulate those in XCOS like an anal

Re: [Scilab-users] printf bug?

2014-10-20 Thread Stéphane Mottelet
Hello, I think that the decimal separator is controlled by the LC_NUMERIC or the LANG environment variable. With my local install (LANG=fr_FR.UTF-8) I get -->msprintf("%0.1f",1.2345) ans = 1.2 S. Le 20/10/2014 17:48, Stefan Du Rietz a écrit : Hello, look here: -->x = 1.2345 x =

Re: [Scilab-users] printf bug?

2014-10-20 Thread Stéphane Mottelet
z a écrit : But why the difference between msprintf() and string()? Stefan On 2014-10-20 18:19, Stéphane Mottelet wrote: Hello, I think that the decimal separator is controlled by the LC_NUMERIC or the LANG environment variable. With my local install (LANG=fr_FR.UTF-8) I get -->msprintf(&q

Re: [Scilab-users] 3d interpolation griddata

2014-11-13 Thread Stéphane Mottelet
Hello, the cshep2d macro should fit your needs. S. Le 12/11/2014 22:40, simond a écrit : Hi, I am trying to convert a matlab code using griddata. Unfortunately, I havent succeed to solve the problem with the available scilal functions (linear_interpn, interpd..). I have a data set of disorde

Re: [Scilab-users] Plotlib

2014-11-28 Thread Stéphane Mottelet
Le 28/11/2014 15:45, Collewet Guylaine a écrit : Hello, I would like to plot an image using a colormap with plotlib (pcolor function) In order to "zoom" on the values between A and B, I use the function caxis([A B]) But I would like also the pixels with values lower than A or greater tha

Re: [Scilab-users] Parallel execution on Mac

2015-02-15 Thread Stéphane Mottelet
Le 15/02/2015 12:01, Arvid Rosén a écrit : Anyone out there who is using any successful ways of running heavy scilab jobs in parallel across several cores on Mac OS? 8 cores, with 1 in use is, is a bit annoying when you are sitting there waiting for the processing to finish. Cheers, Arvid ___

[Scilab-users] bug #13158 /parallel_run on MacOSX

2015-04-07 Thread Stéphane Mottelet
Hello, Is there any plan to solve bug #13158 ? http://bugzilla.scilab.org/show_bug.cgi?id=13158 There is no assignee since the bug has been signaled (in 2014)... As far as I am concerned, I have access to a Linux machine where parallel_run is completely functional, but my usual machine is a M

Re: [Scilab-users] scilab 5.5.2 Ubuntu 14.10 plot error

2015-04-08 Thread Stéphane Mottelet
Hello, You are using the plotlib, please file a ticket at the plotlib Atoms page: https://atoms.scilab.org/toolboxes/plotlib S. Le 08/04/2015 08:53, Vincent COUVERT a écrit : Hello, The plot function you use is not Scilab's one (_update_leaf is not defined in Scilab). Did you install some

Re: [Scilab-users] bug #13158 /parallel_run on MacOSX

2015-04-14 Thread Stéphane Mottelet
Le 07/04/2015 15:39, Stéphane Mottelet a écrit : Hello, Is there any plan to solve bug #13158 ? http://bugzilla.scilab.org/show_bug.cgi?id=13158 There is no assignee since the bug has been signaled (in 2014)... As far as I am concerned, I have access to a Linux machine where parallel_run is

[Scilab-users] Ways to speed up simple things in Scilab ?

2015-04-23 Thread Stéphane Mottelet
Hello, I am currently working on a project where Scilab code is automatically generated, and after many code optimization, the remaining bottleneck is the time that Scilab spends to execute simple code like this (full script (where the vector has 839 lines) with timings is attached) : M1_v=[

Re: [Scilab-users] lsqrsolve

2015-04-23 Thread Stéphane Mottelet
Le 24/04/2015 01:56, fujimoto2005 a écrit : function ret=err_test(param,argVector1,argVector2) err(1)=param(1)^2; err(2)=param(2)^2; err(3)=param(3)^2; err(4)=param(1)*param(2); ret=err; endfunction para=[1,2,3]; argVector1=zeros(1,5); argVector2=zeros(1,3) lsqrsolve(par

Re: [Scilab-users] Ways to speed up simple things in Scilab ?

2015-04-24 Thread Stéphane Mottelet
diff, UK, CF15 7AB._ __www.hidglobal.com_ <http://www.hidglobal.com/> From: Stéphane Mottelet To: "International users mailing list for Scilab." Date: 23/04/2015 22:52 Subject:[Scilab-users] Wa

Re: [Scilab-users] Ways to speed up simple things in Scilab ?

2015-04-24 Thread Stéphane Mottelet
Hello, this is not trivial indexing, in fact some terms are linear combination of v's components M1_v=[v(17) v(104) v(149) -(v(18)+v(63)+v(103)) -(v(18)+v(63)+v(103)) v(17) ... v(104) v(149) ] How do you take this into account in your proposed method ? These combinations are sums of influxes

Re: [Scilab-users] Ways to speed up simple things in Scilab ?

2015-04-24 Thread Stéphane Mottelet
t in Scilab. > > Antoine > > > Le 04/23/2015 11:51 PM, Stéphane Mottelet a écrit : >> Hello, >> >> I am currently working on a project where Scilab code is automatically >> generated, and after many code optimization, the remaining bottleneck is the >&g

Re: [Scilab-users] Ways to speed up simple things in Scilab ?

2015-04-24 Thread Stéphane Mottelet
Le 24/04/2015 15:35, Samuel Gougeon a écrit : Le 24/04/2015 14:50, Stéphane Mottelet a écrit : Hello Antoine I think your idea is excellent in my context, I Will try it and keep you informed of the actual improvement, S. Le 24 avr. 2015 à 13:31, Antoine Monmayrant mailto:antoine.monmayr

Re: [Scilab-users] Ways to speed up simple things in Scilab ?

2015-04-24 Thread Stéphane Mottelet
ainly will result in correctness problems. It's likely better to make your improvements in the software that is generating the scilab code. I think the main feature there is to have it generate a vectorized implementation instead of a procedural one (as you have now). -Brian On Fri, Apr 24,

Re: [Scilab-users] Accidentally displaying huge matrices

2015-05-21 Thread Stéphane Mottelet
Le 22/05/2015 02:42, Antoine Monmayrant a écrit : Le Jeudi 21 Mai 2015 17:48 CEST, Tim Wescott a écrit: If you're working with something humongous and you do some operation that displays, it takes forever before your mistake is done displaying. E.g., if x is 10 elements long and you t

Re: [Scilab-users] Accidentally displaying huge matrices

2015-05-26 Thread Stéphane Mottelet
btw how can you overload the display of double matrices ? S. Le 26/05/2015 11:34, Serge Steer a écrit : Le 25/05/2015 15:29, Antoine Monmayrant a écrit : Le Samedi 23 Mai 2015 00:05 CEST, Samuel Gougeon a écrit: Hello Tim, Le 21/05/2015 17:48, Tim Wescott a écrit : .../... First, is t

Re: [Scilab-users] How to speed up making big structure?

2015-07-03 Thread Stéphane Mottelet
Le 03/07/2015 11:01, Samuel Gougeon a écrit : Le 03/07/2015 09:42, jaipur a écrit : [...] *Conclusion*: for structures arrays, - Scilab 6 (aka YaSp) is ~12 times faster than the best of Scilab 5.5.2 - initializing the array size still makes things faster. It decreases by a factor ~2 the time

[Scilab-users] nightly builds of scilab 6 ?

2015-09-16 Thread Stéphane Mottelet
Hello, When will be recent nighty builds of Scilab 6 available ? Some important bugs have been fixed since the alpha version so it would be nice to benefit from these fixes. S. -- Département de Génie Informatique EA 4297 Transformations Intégrées de la Matière Renouvelable Université de Tec

Re: [Scilab-users] Fwd: How to get fjac from lsqrsolve

2015-09-17 Thread Stéphane Mottelet
Hello, You can use the numderivative macro, with your residual function as first argument. hth S. Le 17/09/2015 16:18, Pablo Caron a écrit : Dear user group, I need to extract the Jacobian from the lsqrsolve function. The fortran routine lmdif (http://www.netlib.org/minpack/lmdif.f) prov

Re: [Scilab-users] Function

2015-10-01 Thread Stéphane Mottelet
Le 01/10/2015 12:36, Rafael Guerra a écrit : I can only think of Da Vinci's wise words: "simplicity is the ultimate sophistication" +1, definitively > To: users@lists.scilab.org > From: cfutt...@gmail.com > Date: Wed, 30 Sep 2015 23:51:11 +0200 > Subject: Re: [Scilab-users] Function > > Hi Chr

Re: [Scilab-users] Backslash operator with square singular matrix

2015-10-27 Thread Stéphane Mottelet
Hello, Le 27/10/2015 13:27, Pierre Vuillemin a écrit : A = [2 4;1 2]; disp(det(A)) b = [1;1]; x = A \ b; -->A = [2 4;1 2]; -->disp(det(A)) 0. -->b = [1;1]; -->x = A \ b; Attention : La matrice est presque singulière ou mal conditionnée. rcond = 0.D+00 Calcul de la solution des mo

Re: [Scilab-users] Backslash operator with square singular matrix

2015-10-27 Thread Stéphane Mottelet
I forgot to mention that I tested this with 5.4.1 S. > Le 27 oct. 2015 à 15:14, JLan a écrit : > > Different in 6.0.0-alpha: > --> A = [2 4;1 2]; > --> disp(det(A)) > det: LAPACK error n°2. > --> b = [1;1]; > --> x = A \ b; > --> disp(norm(A*x-b)); >0.4472136 > > J > > > > -- > View th

Re: [Scilab-users] Problem installing modules manually (atomsInstall)

2016-01-15 Thread Stéphane Mottelet
Le 15/01/2016 14:25, Lester Anderson a écrit : Hello Antoine, How does one get a "cached" version of the repository list and where is it placed? -->atomsRepositoryList() ans = !http://atoms.scilab.org/5.5 official ! -->atomsVersion() ans = 5.5 Is there an option to disable fetchin

Re: [Scilab-users] Conditional random variable

2016-01-23 Thread Stéphane Mottelet
Hello, You will find the Matlab code for the "trandn" function there : http://www.mathworks.com/matlabcentral/fileexchange/53180-truncated-normal-generator/content/trandn.m hth S. Le 23/01/2016 13:25, fujimoto2005 a écrit : Is there any function generating random variable X with a standard

Re: [Scilab-users] Conditional random variable

2016-01-23 Thread Stéphane Mottelet
You also have this distribution (and many more) in the atoms package https://forge.scilab.org/index.php/p/distfun/ S. Le 23/01/2016 13:25, fujimoto2005 a écrit : Is there any function generating random variable X with a standard normal density conditional on a range (a,b] where 0http://mailing

Re: [Scilab-users] optim

2016-01-25 Thread Stéphane Mottelet
Hello, your costfunction has to provide f and its gradient : function [f, g, ind]=costfunction(x, ind, para1, para2) f=x(1)^2+x(2)^2; g=2*x; endfunction x0=[100,100]'; para1=[1,2]'; para2=[2,3]'; costf=list(costfunction,para1,para2) [fopt,xopt]=optim(costf,x0) S. Le 25/01/2016 10:06

Re: [Scilab-users] optim

2016-01-25 Thread Stéphane Mottelet
btw the error message is incorrect as you had an incorrect number of *output* arguments (Scilab 6.0 bug ?) S. Le 25/01/2016 10:06, fujimoto2005 a écrit : I'm useing scilab 6.0. I want to use 'optim' function . But it fail to work in my code with an error message 'costfunction: Wrong number of

Re: [Scilab-users] 'fsolver'

2016-01-27 Thread Stéphane Mottelet
Le 27/01/2016 15:24, fujimoto2005 a écrit : if f(x)=0 has multiple solutions and I want to get the smallest solution, is there any way to get a such solution by 'fsolver'? In my actual problem, 'fsolver' give the largest solution. is the solution set countable or not ? I mean, is the problem u

Re: [Scilab-users] 'fsolver'

2016-01-27 Thread Stéphane Mottelet
Le 27/01/2016 16:54, fujimoto2005 a écrit : The feature of my f(x) defined for x>0 is as follows. f(x)<0 for x0 for x1=x2 'fsolver' gives some x such as x>x2 as a solution. I want to get x1 as a solution. -- View this message in context: http://mailinglists.scilab.org/fsolver-tp4033340p4033

Re: [Scilab-users] 'fsolver'

2016-01-29 Thread Stéphane Mottelet
Hello, Can you give us the actual definition of the function f(x) (i suppose you are trying to find x such that f(x)=0 ) ? S. Le 29/01/2016 08:38, fujimoto2005 a écrit : Hi, motterlet and Steer Thanks a lot your helps. Unfortunately 'lsqrsolve' did't give the smallest initial point for my

Re: [Scilab-users] Logical zero

2016-02-01 Thread Stéphane Mottelet
Le 01/02/2016 09:10, shamikam a écrit : Is there a Scilab equivalent for the false function in Matlab? I want to create a matrix of logical zeros. Is there any other way to do it. Shamika -- View this message in context: http://mailinglists.scilab.org/Logical-zero-tp4033361.html Sent from

Re: [Scilab-users] Logical zero

2016-02-01 Thread Stéphane Mottelet
BTY, I am worried about the Scilab 6.0 timing (3.3228213 seconds) of a(1:n,1:m)=%f; i.e. more than ten times slower than Scilab 5.5.1 (0.2964019 seconds). Does somebody have an explanation ? S. Le 01/02/2016 14:40, Jan-Åge Langeland a écrit : I ran a little speed test with the different alt

[Scilab-users] Scilab 6 bad performances

2016-02-19 Thread Stéphane Mottelet
Hello, I am currently porting a Scilab package (Sysmetab, on the Scilab forge but not yet Atoms packaged) and I have noticed that Scilab 6 is two times slower than Scilab 5.5.2. The involved computations are very diverse but are vectorized products of submatrices, sparse matrix algebra and sp

Re: [Scilab-users] Scilab 6 bad performances

2016-02-19 Thread Stéphane Mottelet
Le 19/02/2016 11:17, Samuel Gougeon a écrit : Hello, Le 19/02/2016 11:05, Stéphane Mottelet a écrit : Hello, I am currently porting a Scilab package (Sysmetab, on the Scilab forge but not yet Atoms packaged) and I have noticed that Scilab 6 is two times slower than Scilab 5.5.2. The

Re: [Scilab-users] Scilab 6 bad performances

2016-02-19 Thread Stéphane Mottelet
Le 19/02/2016 11:41, Samuel Gougeon a écrit : Le 19/02/2016 11:31, Stéphane Mottelet a écrit : Le 19/02/2016 11:17, Samuel Gougeon a écrit : .../... bench_run() is designed for that: https://help.scilab.org/docs/5.5.2/en_US/bench_run.html On the same computer, you may run it on 5.5.2 for

Re: [Scilab-users] leastsq question : what is 'gopt' useful for?

2016-03-22 Thread Stéphane Mottelet
Hello, Le 22/03/2016 10:41, antoine.monmayr...@laas.fr a écrit : Hi everyone, I have a very general and naive question concerning leastsq: what am I to do with "gopt", the "gradient of f at xopt"? Is there a way to link it to the confidence interval for each parameter of my fit? Not really,

Re: [Scilab-users] "Smoothing" very localised discontinuities in curves.

2016-04-03 Thread Stéphane Mottelet
Why not use median filtering with a small window ? > Le 4 avr. 2016 à 07:57, t...@wescottdesign.com a écrit : > > As written (with or without the commenting out) that code won't produce phase > shifts. In technical terms it's a FIR filter that is symmetrical around zero > delay, and such filt

Re: [Scilab-users] "Smoothing" very localised discontinuities in (scilab: to exclusive) (scilab: to exclusive) curves.

2016-04-04 Thread Stéphane Mottelet
Hello, The last time I had used a median filter, it was to locate peaks in a frequency DSP. The idea was to substract the median-filtered spectrum to the original one, and treshold the difference. This was enough to locate the peaks. Maybe you could use the same idea there. S. Le 04/04/2016

Re: [Scilab-users] Function int() strange behavior

2016-04-17 Thread Stéphane Mottelet
It seems to be the normal behavior, as int() and floor() behave the same for positive numbers (and 1.45 cannot be stored exactly with IEEE 754) S. Le 18/04/2016 08:17, ol.bond a écrit : Thank you Samuel, I got it also 2016-04-18 15:13 GMT+09:00 Samuel GOUGEON [via Scilab / Xcos - Mailing

Re: [Scilab-users] create random values between [-1,1]

2016-06-10 Thread Stéphane Mottelet
Le 10/06/2016 à 16:19, Stéphane Mottelet a écrit : Le 10/06/2016 à 15:59, Jens Simon Strom a écrit : rand(n,1)–rand(n,1) is definitely NOT uniform. by convolution the density of such a variable seems triangle shaped : f(x)=1/2-abs(x)/4, if x \in [-2,2] =0, elsewhere oops, rather f(x

Re: [Scilab-users] create random values between [-1,1]

2016-06-10 Thread Stéphane Mottelet
Le 10/06/2016 à 15:59, Jens Simon Strom a écrit : rand(n,1)–rand(n,1) is definitely NOT uniform. by convolution the density of such a variable seems triangle shaped : f(x)=1/2-abs(x)/4, if x \in [-2,2] =0, elsewhere S. Jens ---

Re: [Scilab-users] issue with parallel_run

2016-06-14 Thread Stéphane Mottelet
Hello, 1-parallel_run is designed to cope with input-output data organized in columns 2-you have to precise the size of the output of yout fct() function. Here is your corrected script (with comments) : mode(0); clear; function random_matrix=fct(n) M = grand(n,12,'unf',-1,1); M =

Re: [Scilab-users] parallel_run mixed with vectorization

2016-06-14 Thread Stéphane Mottelet
Hello, IMHO parallel_run is very unstable and should be considered as an experimental feature. I never managed to make it work with multiple arguments. The workaround is the following : function U=fct(i) U1 = R1(i) * i1(i); U2 = R2(i) * i2(i); U = [U1 U2]'; endfunction n = 1000;

Re: [Scilab-users] Scilab vs matlab

2016-09-26 Thread Stéphane Mottelet
The price :-D Le 26/09/2016 à 16:20, Chijioke Kenechukwu a écrit : What is the advantage of scilab over matlab ___ users mailing list users@lists.scilab.org http://lists.scilab.org/mailman/listinfo/users _

Re: [Scilab-users] issue with "sum"

2016-09-27 Thread Stéphane Mottelet
paul, remember the transpose is *conjugate transpose*, use .' instead. S. Le 27/09/2016 à 15:04, paul.carr...@free.fr a écrit : Hi again In the following example, the sum lead to an opposite sign for the complex part … Am I doing something wrong ? thanks for any highlight Paul ###

Re: [Scilab-users] scilab and "parallelization" objectives

2016-10-22 Thread Stéphane Mottelet
Hello, As parallel_run is using child Scilab processes, the std output of all children is the terminal where you started Scilab and not the Scilab console. S. > Le 22 oct. 2016 à 15:11, paul.carr...@free.fr a écrit : > > Thanks claus > > I completely forgot this features (that I've ever test

Re: [Scilab-users] scilab and "parallelization" objectives

2016-10-23 Thread Stéphane Mottelet
gt; What OS are you running? parallel_run does not work on macOS. It does work on > Linux though. > > Cheers, > Arvid > > Get Outlook for iOS > > > From: users on behalf of Stéphane Mottelet > > Sent: Saturday, October 22, 2016 4:33:22 PM > To: Users maili

Re: [Scilab-users] The sound of mathematical formulas

2016-11-05 Thread Stéphane Mottelet
Hello, Cascaded frequency modulation gives nices sounds. For example signals like sin((w0+a1*sin(w1*t))*t) w1 itself can also be modulated, and so on. This kind of synthesis was used by the famous Yamaha DX-7 and gave typical sounds (bass guitar, bells, ...) S. Le 05/11/2016 à 13:43, Rafa

  1   2   3   4   5   6   7   >