From delaleau at enib.fr  Tue May  1 00:50:22 2018
From: delaleau at enib.fr (Emmanuel Delaleau)
Date: Tue, 1 May 2018 00:50:22 +0200
Subject: [Scilab-users] Scliab 6.0.x compilation problem (jogl2)
In-Reply-To: <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
Message-ID: <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>

I have all jogl2 and gluegen2 packages included in Debian/stable 
installed and still the same error message.



Le 30/04/2018 ? 13:22, Masaru Nomiya a ?crit?:
> Hello,
>
> In the Message;
>
>    Subject    : Re: [Scilab-users] Scliab 6.0.x compilation problem (jogl2)
>    Message-ID : <b494dec9-123b-73d0-6e0d-22b3801c7f07 at enib.fr>
>    Date & Time: Mon, 30 Apr 2018 10:16:55 +0200
>
> [ED] == Emmanuel Delaleau <delaleau at enib.fr> has written:
>
> ED>  Dear Masaru Nomiya, thank you for your response. See result below (I am under
> ED>  Debian and I use dpkg instead of rpm) :
>
> ED>  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
>
> ED>  delaleau at ed-lapt:~$ dpkg -L libjogl2-java
> ED>  /.
> ED>  /usr
> ED>  /usr/share
> ED>  /usr/share/doc
> ED>  /usr/share/doc/libjogl2-java
> ED>  /usr/share/doc/libjogl2-java/README.Debian
> ED>  /usr/share/doc/libjogl2-java/changelog.Debian.gz
> ED>  /usr/share/doc/libjogl2-java/changelog.gz
> ED>  /usr/share/doc/libjogl2-java/copyright
> ED>  /usr/share/java
> ED>  /usr/share/java/jogl2-2.3.2.jar
> ED>  /usr/share/maven-repo
> ED>  /usr/share/maven-repo/org
> ED>  /usr/share/maven-repo/org/jogamp
> ED>  /usr/share/maven-repo/org/jogamp/jogl
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all/2.3.2
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all/2.3.2/jogl-all-2.3.2.pom
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all/debian
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all/debian/jogl-all-debian.pom
> ED>  /usr/share/java/jogl2.jar
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all/2.3.2/jogl-all-2.3.2.jar
> ED>  /usr/share/maven-repo/org/jogamp/jogl/jogl-all/debian/jogl-all-debian.jar
>
> ED>  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
>
> ED>  delaleau at ed-lapt:~$ dpkg -L libgluegen2-rt-java
> ED>  /.
> ED>  /usr
> ED>  /usr/share
> ED>  /usr/share/doc
> ED>  /usr/share/doc/libgluegen2-rt-java
> ED>  /usr/share/doc/libgluegen2-rt-java/changelog.Debian.gz
> ED>  /usr/share/doc/libgluegen2-rt-java/copyright
> ED>  /usr/share/java
> ED>  /usr/share/java/gluegen2-rt-2.3.2.jar
> ED>  /usr/share/maven-repo
> ED>  /usr/share/maven-repo/org
> ED>  /usr/share/maven-repo/org/jogamp
> ED>  /usr/share/maven-repo/org/jogamp/gluegen
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt/2.3.2
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt/2.3.2/gluegen-rt-2.3.2.pom
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt/debian
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt/debian/gluegen-rt-debian.pom
> ED>  /usr/share/java/gluegen2-rt.jar
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt/2.3.2/gluegen-rt-2.3.2.jar
> ED>  /usr/share/maven-repo/org/jogamp/gluegen/gluegen-rt/debian/gluegen-rt-debian.jar
>
> As you might know, gulegen2 and jogl2 ara key factors for building
> scilab 6.0.x.
>
> You don't have gluegen2-rt.so, do you?
>
> This is strange. Without gluegen2-rt.so, you can't compile scilab.
>
> Another problem.?Do as follows;
>
> # mkdir /usr/lib64/java
> # ln -s /usr/share/java/jogl2.jar /usr/lib64/java
>
> Regards,
>
> ---
> ????? Masaru Nomiya             mail-to: nomiya @ galaxy.dti.ne.jp
> ????
> ???? "Bill! You married with Computers.
>            Not with Me!"
>           "No..., with money."
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users



From heinznabielek at me.com  Thu May  3 10:57:18 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Thu, 03 May 2018 10:57:18 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull random
 variables with prior given autocorrelation
In-Reply-To: <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
Message-ID: <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>

Friends:

urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.

Uncorrelated is easy with 
vc*((-log(grand(1,5,"def"))).^(1/m)) producing
0.9666848   1.4646143   10.402793   5.9513675   3.3241823
with characteristic velocity vc=7 and modulus m=2.

But I need Weibull deviates with given autocorrelation, say rho=0.80

Literature has dozen of references on the subject which either do not work or I do not know how to code.

Any help here?
Heinz

From heinznabielek at me.com  Thu May  3 11:15:29 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Thu, 03 May 2018 11:15:29 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull
 random variables with prior given autocorrelation
In-Reply-To: <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
 <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
Message-ID: <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>

Uniform U(0,1) random deviates with given autocorrelation would be good enough?...


> On 03.05.2018, at 10:57, Heinz Nabielek <heinznabielek at me.com> wrote:
> 
> Friends:
> 
> urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.
> 
> Uncorrelated is easy with 
> vc*((-log(grand(1,5,"def"))).^(1/m)) producing
> 0.9666848   1.4646143   10.402793   5.9513675   3.3241823
> with characteristic velocity vc=7 and modulus m=2.
> 
> But I need Weibull deviates with given autocorrelation, say rho=0.80
> 
> Literature has dozen of references on the subject which either do not work or I do not know how to code.
> 
> Any help here?
> Heinz



From stephane.mottelet at utc.fr  Thu May  3 11:37:38 2018
From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=)
Date: Thu, 3 May 2018 11:37:38 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull
 random variables with prior given autocorrelation
In-Reply-To: <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
 <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
 <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>
Message-ID: <61156ad9-6c46-19fb-3460-1d624e87e586@utc.fr>

Why don't start with a simple autoregressive process ? See e.g.

https://en.wikipedia.org/wiki/Autoregressive_model


Le 03/05/2018 ? 11:15, Heinz Nabielek a ?crit?:
> Uniform U(0,1) random deviates with given autocorrelation would be good enough?...
>
>
>> On 03.05.2018, at 10:57, Heinz Nabielek <heinznabielek at me.com> wrote:
>>
>> Friends:
>>
>> urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.
>>
>> Uncorrelated is easy with
>> vc*((-log(grand(1,5,"def"))).^(1/m)) producing
>> 0.9666848   1.4646143   10.402793   5.9513675   3.3241823
>> with characteristic velocity vc=7 and modulus m=2.
>>
>> But I need Weibull deviates with given autocorrelation, say rho=0.80
>>
>> Literature has dozen of references on the subject which either do not work or I do not know how to code.
>>
>> Any help here?
>> Heinz
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users



From heinznabielek at me.com  Thu May  3 15:49:04 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Thu, 03 May 2018 15:49:04 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull
 random variables with prior given autocorrelation
In-Reply-To: <61156ad9-6c46-19fb-3460-1d624e87e586@utc.fr>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
 <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
 <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>
 <61156ad9-6c46-19fb-3460-1d624e87e586@utc.fr>
Message-ID: <B6BA05AB-D19E-446B-9DAF-DCA184A6C1BB@me.com>

Many thanks. <https://en.wikipedia.org/wiki/Autoregressive_model#Example:_An_AR(1)_process> works fine on an originating normal distribution [where i would not need it, because we have grand(n, "mn", Mean, Cov) in SciLab].

With an initiating U(-0.5, 0,5), however, I get a triangle function, but not a autocorrelated uniform random distribution.

What to do?
Heinz




> On 03.05.2018, at 11:37, St?phane Mottelet <stephane.mottelet at utc.fr> wrote:
> 
> Why don't start with a simple autoregressive process ? See e.g.
> 
> https://en.wikipedia.org/wiki/Autoregressive_model
> 
> 
> Le 03/05/2018 ? 11:15, Heinz Nabielek a ?crit :
>> Uniform U(0,1) random deviates with given autocorrelation would be good enough?...
>> 
>> 
>>> On 03.05.2018, at 10:57, Heinz Nabielek <heinznabielek at me.com> wrote:
>>> 
>>> Friends:
>>> 
>>> urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.
>>> 
>>> Uncorrelated is easy with
>>> vc*((-log(grand(1,5,"def"))).^(1/m)) producing
>>> 0.9666848   1.4646143   10.402793   5.9513675   3.3241823
>>> with characteristic velocity vc=7 and modulus m=2.
>>> 
>>> But I need Weibull deviates with given autocorrelation, say rho=0.80
>>> 
>>> Literature has dozen of references on the subject which either do not work or I do not know how to code.
>>> 
>>> Any help here?
>>> Heinz



From stephane.mottelet at utc.fr  Thu May  3 17:40:39 2018
From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=)
Date: Thu, 3 May 2018 17:40:39 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull
 random variables with prior given autocorrelation
In-Reply-To: <B6BA05AB-D19E-446B-9DAF-DCA184A6C1BB@me.com>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
 <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
 <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>
 <61156ad9-6c46-19fb-3460-1d624e87e586@utc.fr>
 <B6BA05AB-D19E-446B-9DAF-DCA184A6C1BB@me.com>
Message-ID: <2bb062ce-3002-2581-10d9-d0172adfafc7@utc.fr>

Le 03/05/2018 ? 15:49, Heinz Nabielek a ?crit?:
> Many thanks. <https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/en.wikipedia.org/wiki/Autoregressive_model#Example:_An_AR(1)_process> works fine on an originating normal distribution [where i would not need it, because we have grand(n, "mn", Mean, Cov) in SciLab].
grand(n, "mn", Mean, Cov) allows to generate realizations of a 
multi-normal variable, i.e. covariance matrix describes the covariance 
between each component? of the (vector) variable.

AR processes are different stuff (random process <> random variable). 
For an AR process the non-idependance is w.r.t. to time (discrete or 
continuous)

Your wind data is a random process.
> With an initiating U(-0.5, 0,5), however, I get a triangle function, but not a autocorrelated uniform random distribution.
what do you mean by " autocorrelated uniform random distribution" ? Here 
your triangle is just the distribution of a sum of two uniform vraiables.

S.
> What to do?
> Heinz
>
>
>
>
>> On 03.05.2018, at 11:37, St?phane Mottelet <stephane.mottelet at utc.fr> wrote:
>>
>> Why don't start with a simple autoregressive process ? See e.g.
>>
>> https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/en.wikipedia.org/wiki/Autoregressive_model
>>
>>
>> Le 03/05/2018 ? 11:15, Heinz Nabielek a ?crit :
>>> Uniform U(0,1) random deviates with given autocorrelation would be good enough?...
>>>
>>>
>>>> On 03.05.2018, at 10:57, Heinz Nabielek <heinznabielek at me.com> wrote:
>>>>
>>>> Friends:
>>>>
>>>> urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.
>>>>
>>>> Uncorrelated is easy with
>>>> vc*((-log(grand(1,5,"def"))).^(1/m)) producing
>>>> 0.9666848   1.4646143   10.402793   5.9513675   3.3241823
>>>> with characteristic velocity vc=7 and modulus m=2.
>>>>
>>>> But I need Weibull deviates with given autocorrelation, say rho=0.80
>>>>
>>>> Literature has dozen of references on the subject which either do not work or I do not know how to code.
>>>>
>>>> Any help here?
>>>> Heinz
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users



From stephane.mottelet at utc.fr  Thu May  3 18:29:42 2018
From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=)
Date: Thu, 3 May 2018 18:29:42 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull
 random variables with prior given autocorrelation
In-Reply-To: <B6BA05AB-D19E-446B-9DAF-DCA184A6C1BB@me.com>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
 <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
 <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>
 <61156ad9-6c46-19fb-3460-1d624e87e586@utc.fr>
 <B6BA05AB-D19E-446B-9DAF-DCA184A6C1BB@me.com>
Message-ID: <9d5d4a61-c10b-9f45-8511-72fcad501ad1@utc.fr>

Le 03/05/2018 ? 15:49, Heinz Nabielek a ?crit?:
> Many thanks. <https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/en.wikipedia.org/wiki/Autoregressive_model#Example:_An_AR(1)_process> works fine on an originating normal distribution [where i would not need it, because we have grand(n, "mn", Mean, Cov) in SciLab].
grand(n, "mn", Mean, Cov) allows to generate realizations of a 
multi-normal variable, i.e. covariance matrix describes the covariance 
between each component? of the (vector) variable.

AR processes are different stuff (random process <> random variable). 
For an AR process the non-idependance is w.r.t. to time (discrete or 
continuous)

Your wind data is a random process.
> With an initiating U(-0.5, 0,5), however, I get a triangle function, but not a autocorrelated uniform random distribution.
what do you mean by " autocorrelated uniform random distribution" ? Here 
your triangle is just the distribution of a sum of two uniform variables.

S.
> What to do?
> Heinz
>
>
>
>
>> On 03.05.2018, at 11:37, St?phane Mottelet <stephane.mottelet at utc.fr> wrote:
>>
>> Why don't start with a simple autoregressive process ? See e.g.
>>
>> https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/en.wikipedia.org/wiki/Autoregressive_model
>>
>>
>> Le 03/05/2018 ? 11:15, Heinz Nabielek a ?crit :
>>> Uniform U(0,1) random deviates with given autocorrelation would be good enough?...
>>>
>>>
>>>> On 03.05.2018, at 10:57, Heinz Nabielek <heinznabielek at me.com> wrote:
>>>>
>>>> Friends:
>>>>
>>>> urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.
>>>>
>>>> Uncorrelated is easy with
>>>> vc*((-log(grand(1,5,"def"))).^(1/m)) producing
>>>> 0.9666848   1.4646143   10.402793   5.9513675   3.3241823
>>>> with characteristic velocity vc=7 and modulus m=2.
>>>>
>>>> But I need Weibull deviates with given autocorrelation, say rho=0.80
>>>>
>>>> Literature has dozen of references on the subject which either do not work or I do not know how to code.
>>>>
>>>> Any help here?
>>>> Heinz
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users



From heinznabielek at me.com  Thu May  3 18:30:51 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Thu, 03 May 2018 18:30:51 +0200
Subject: [Scilab-users] needed: SciLab routine for generating Weibull
 random variables with prior given autocorrelation
In-Reply-To: <2bb062ce-3002-2581-10d9-d0172adfafc7@utc.fr>
References: <6b2a8014-6244-9c3f-3c74-1b25d79cbc64@enib.fr>
 <8736zd8llt.wl-nomiya@galaxy.dti.ne.jp>
 <b494dec9-123b-73d0-6e0d-22b3801c7f07@enib.fr>
 <87lgd5djln.wl-nomiya@galaxy.dti.ne.jp>
 <65f1016b-cc20-4d0d-1668-d6156c8e590b@enib.fr>
 <A7E8D8DA-98D9-444B-98DD-A03526BD8F05@me.com>
 <FF78A374-CBB0-43BF-AE7F-059BD622A542@me.com>
 <61156ad9-6c46-19fb-3460-1d624e87e586@utc.fr>
 <B6BA05AB-D19E-446B-9DAF-DCA184A6C1BB@me.com>
 <2bb062ce-3002-2581-10d9-d0172adfafc7@utc.fr>
Message-ID: <66650484-E817-4E27-9BB2-0DA6E2304DB8@me.com>

thanks. 

but still: how construct autocorrelated uniform randoms?

Sent from Heinz Nabielek

> On 03 May 2018, at 17:40, St?phane Mottelet <stephane.mottelet at utc.fr> wrote:
> 
>> Le 03/05/2018 ? 15:49, Heinz Nabielek a ?crit :
>> Many thanks. <https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/en.wikipedia.org/wiki/Autoregressive_model#Example:_An_AR(1)_process> works fine on an originating normal distribution [where i would not need it, because we have grand(n, "mn", Mean, Cov) in SciLab].
> grand(n, "mn", Mean, Cov) allows to generate realizations of a multi-normal variable, i.e. covariance matrix describes the covariance between each component  of the (vector) variable.
> 
> AR processes are different stuff (random process <> random variable). For an AR process the non-idependance is w.r.t. to time (discrete or continuous)
> 
> Your wind data is a random process.
>> With an initiating U(-0.5, 0,5), however, I get a triangle function, but not a autocorrelated uniform random distribution.
> what do you mean by " autocorrelated uniform random distribution" ? Here your triangle is just the distribution of a sum of two uniform vraiables.
> 
> S.
>> What to do?
>> Heinz
>> 
>> 
>> 
>> 
>>> On 03.05.2018, at 11:37, St?phane Mottelet <stephane.mottelet at utc.fr> wrote:
>>> 
>>> Why don't start with a simple autoregressive process ? See e.g.
>>> 
>>> https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/en.wikipedia.org/wiki/Autoregressive_model
>>> 
>>> 
>>>> Le 03/05/2018 ? 11:15, Heinz Nabielek a ?crit :
>>>> Uniform U(0,1) random deviates with given autocorrelation would be good enough?...
>>>> 
>>>> 
>>>>> On 03.05.2018, at 10:57, Heinz Nabielek <heinznabielek at me.com> wrote:
>>>>> 
>>>>> Friends:
>>>>> 
>>>>> urgently needed: a SciLab routine for generating Weibull random variables with prior given autocorrelation.
>>>>> 
>>>>> Uncorrelated is easy with
>>>>> vc*((-log(grand(1,5,"def"))).^(1/m)) producing
>>>>> 0.9666848   1.4646143   10.402793   5.9513675   3.3241823
>>>>> with characteristic velocity vc=7 and modulus m=2.
>>>>> 
>>>>> But I need Weibull deviates with given autocorrelation, say rho=0.80
>>>>> 
>>>>> Literature has dozen of references on the subject which either do not work or I do not know how to code.
>>>>> 
>>>>> Any help here?
>>>>> Heinz
>> _______________________________________________
>> users mailing list
>> users at lists.scilab.org
>> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users
> 
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users


From sgougeon at free.fr  Fri May  4 15:31:43 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Fri, 4 May 2018 15:31:43 +0200
Subject: [Scilab-users] [SEP] factorial() extension on [171, 10^14]
In-Reply-To: <192586ac-7c55-8c6d-80a0-daba885f21f7@free.fr>
References: <192586ac-7c55-8c6d-80a0-daba885f21f7@free.fr>
Message-ID: <8749f83a-67e8-44af-f607-1c54008f6556@free.fr>

Le 14/04/2018 ? 15:44, Samuel Gougeon a ?crit :
> While fixing the bug 7562 about a bad memory usage by factorial(),
>
> i was wondering about the poor factorial() result, since currently
> it can process only integers in [0, 170] before saturating and 
> returning %inf.
>
> This Scilab Enhancement Proposed (SEP) aims to extend factorial()
> to all integers in [0, 10^4].
>

The proposed extension goes up to n=10^14 (not 10^4)

Samuel



From jrafaelbguerra at hotmail.com  Fri May  4 16:13:55 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Fri, 4 May 2018 14:13:55 +0000
Subject: [Scilab-users] [SEP] factorial() extension on [171, 10^14]
In-Reply-To: <8749f83a-67e8-44af-f607-1c54008f6556@free.fr>
References: <192586ac-7c55-8c6d-80a0-daba885f21f7@free.fr>
 <8749f83a-67e8-44af-f607-1c54008f6556@free.fr>
Message-ID: <VI1PR04MB41095223E03BE09765444386CC860@VI1PR04MB4109.eurprd04.prod.outlook.com>

Hi Samuel,

How are in the proposed extension the resulting very large numbers represented?
In other words, how are the limitations of double precision overcome and can we do math with such very large integers?

PS: In Scilab the gammaln() function allows computing the logarithm of the factorial for large integers.

Regards,
Rafael

-----Original Message-----
From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Gougeon
Sent: Friday, May 04, 2018 3:32 PM
To: International users mailing list for Scilab. <users at lists.scilab.org>
Subject: Re: [Scilab-users] [SEP] factorial() extension on [171, 10^14]

Le 14/04/2018 ? 15:44, Samuel Gougeon a ?crit :
> While fixing the bug 7562 about a bad memory usage by factorial(),
>
> i was wondering about the poor factorial() result, since currently it 
> can process only integers in [0, 170] before saturating and returning 
> %inf.
>
> This Scilab Enhancement Proposed (SEP) aims to extend factorial() to 
> all integers in [0, 10^4].
>

The proposed extension goes up to n=10^14 (not 10^4)

Samuel

_______________________________________________
users mailing list
users at lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users


From sgougeon at free.fr  Fri May  4 16:23:43 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Fri, 4 May 2018 16:23:43 +0200
Subject: [Scilab-users] [SEP] factorial() extension on [171, 10^14]
In-Reply-To: <VI1PR04MB41095223E03BE09765444386CC860@VI1PR04MB4109.eurprd04.prod.outlook.com>
References: <192586ac-7c55-8c6d-80a0-daba885f21f7@free.fr>
 <8749f83a-67e8-44af-f607-1c54008f6556@free.fr>
 <VI1PR04MB41095223E03BE09765444386CC860@VI1PR04MB4109.eurprd04.prod.outlook.com>
Message-ID: <057f169c-99c6-b4a5-e227-ac7d6fc6ac5a@free.fr>

Le 04/05/2018 ? 16:13, Rafael Guerra a ?crit :
> Hi Samuel,
>
> How are in the proposed extension the resulting very large numbers represented?
> In other words, how are the limitations of double precision overcome and can we do math with such very large integers?
>
> PS: In Scilab the gammaln() function allows computing the logarithm of the factorial for large integers.

Hello Rafael,

Please see the link to the help page posted in my first message. It is 
documented in the page.
The algorithm can be seen on CodeReview (ref in the message/bugzilla 
entry). It uses a careful
processing with log(). It is more accurate than gammaln() by a factor of 
up to 100 (it is also
documented in the test).

Regards
Samuel



From sgougeon at free.fr  Mon May  7 14:36:45 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Mon, 7 May 2018 14:36:45 +0200
Subject: [Scilab-users] gcf().anti_aliasing : Confirmation of (un)effect?
Message-ID: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>

Hello,

The gcf().anti_aliasing figure property was implemented many years ago 
(in the early Scilab 5.x versions).

It aims to reduce the "staircase effect" for lines drawn very close to 
vertical or horizontal directions.

Yet, has anyone ever seen an actual effect on graphic rendering, as 
expected and described in the related SEP#16 
<http://gitweb.scilab.org/?p=scilab.git;a=blob;f=SEP/SEP_016_anti_aliasing.odt;h=e267088c7ed7e5a4cd0ae29b744616332423224e;hb=HEAD> 
?
I never did, whatever were the computer and its graphic card i used.
Still now with Scilab 6.0, i don't see any effect.

The staircase effect can easily be seen with the following code:

x = [0 1];
y = (0:0.02:1)';
y = [y y+0.02]';
clf
plot2d(x,y)
f = gcf();
f.anti_aliasing  // should display "off" (by default)

f = scf();
f.anti_aliasing = "8x";
plot2d(x,y)

Then compare both figures.
I get the following, whatever is the .anti_aliasing value.


Does anyone get anything else, closer to the expected "dithered" rendering?
Either it will be possible to better document the effect of this property,
as requested in the bug report 6777 
<http://bugzilla.scilab.org/show_bug.cgi?id=6777>, or to reopen the bug 
11311 <http://bugzilla.scilab.org/show_bug.cgi?id=11311> reporting
that it is out of work.

Thanks
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180507/6cd008d9/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dbhpbgjbciobikhj.png
Type: image/png
Size: 7020 bytes
Desc: not available
URL: <https://lists.scilab.org/pipermail/users/attachments/20180507/6cd008d9/attachment.png>

From amonmayr at laas.fr  Mon May  7 15:09:24 2018
From: amonmayr at laas.fr (amonmayr at laas.fr)
Date: Mon, 7 May 2018 15:09:24 +0200
Subject: [Scilab-users] gcf().anti_aliasing : Confirmation of (un)effect?
In-Reply-To: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>
References: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>
Message-ID: <d70ea5af-05d8-41db-5c97-d9318ef47f5a@laas.fr>

Hello Samuel,

With your example, I don't see any difference on the various computers I 
have.
I think I used anti_aliasing in the past but I cannot find back a script 
where I used it.
I think I used it on 3D plots (surf?), maybe for the export to graphic 
files but I am not sure.
If I can find something I'll let you know.

Antoine


Le 07/05/2018 ? 14:36, Samuel Gougeon a ?crit?:
>
> Hello,
>
> The gcf().anti_aliasing figure property was implemented many years ago 
> (in the early Scilab 5.x versions).
>
> It aims to reduce the "staircase effect" for lines drawn very close to 
> vertical or horizontal directions.
>
> Yet, has anyone ever seen an actual effect on graphic rendering, as 
> expected and described in the related SEP#16 
> <http://gitweb.scilab.org/?p=scilab.git;a=blob;f=SEP/SEP_016_anti_aliasing.odt;h=e267088c7ed7e5a4cd0ae29b744616332423224e;hb=HEAD> 
> ?
> I never did, whatever were the computer and its graphic card i used.
> Still now with Scilab 6.0, i don't see any effect.
>
> The staircase effect can easily be seen with the following code:
>
> x = [0 1];
> y = (0:0.02:1)';
> y = [y y+0.02]';
> clf
> plot2d(x,y)
> f = gcf();
> f.anti_aliasing  // should display "off" (by default)
>
> f = scf();
> f.anti_aliasing = "8x";
> plot2d(x,y)
> Then compare both figures.
> I get the following, whatever is the .anti_aliasing value.
>
>
> Does anyone get anything else, closer to the expected "dithered" 
> rendering?
> Either it will be possible to better document the effect of this property,
> as requested in the bug report 6777 
> <http://bugzilla.scilab.org/show_bug.cgi?id=6777>, or to reopen the 
> bug 11311 <http://bugzilla.scilab.org/show_bug.cgi?id=11311> reporting
> that it is out of work.
>
> Thanks
> Samuel
>
>
>
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users


-- 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++

  Antoine Monmayrant LAAS - CNRS
  7 avenue du Colonel Roche
  BP 54200
  31031 TOULOUSE Cedex 4
  FRANCE

  Tel:+33 5 61 33 64 59
  
  email : antoine.monmayrant at laas.fr
  permanent email : antoine.monmayrant at polytechnique.org

+++++++++++++++++++++++++++++++++++++++++++++++++++++++

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180507/08060f07/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dbhpbgjbciobikhj.png
Type: image/png
Size: 7020 bytes
Desc: not available
URL: <https://lists.scilab.org/pipermail/users/attachments/20180507/08060f07/attachment.png>

From jrafaelbguerra at hotmail.com  Mon May  7 16:23:42 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Mon, 7 May 2018 07:23:42 -0700 (MST)
Subject: [Scilab-users] gcf().anti_aliasing : Confirmation of (un)effect?
In-Reply-To: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>
References: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>
Message-ID: <1525703022019-0.post@n3.nabble.com>

Hi Samuel,

Could not see any anti-aliasing being performed, either in Scilab 5 or
Scilab 6.
A quick trial on customizing the anti-aliasing settings in NVIDIA Quadro
K4100M for the Scilab program did not help either:
<http://mailinglists.scilab.org/file/t495698/NVIDIA_antialiasing.gif> 

It may need more tweaking.

Regards,
Rafael




--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From cfuttrup at gmail.com  Mon May  7 20:45:11 2018
From: cfuttrup at gmail.com (Claus Futtrup)
Date: Mon, 07 May 2018 18:45:11 +0000
Subject: [Scilab-users] GSOC 2018
Message-ID: <CA+Q5PwtZ4FLSeQJFAnRbs20PEjvTnihNqm2QxudcqoPUJi2i6Q@mail.gmail.com>

Hi there

Can Scilab confirm that there are four GSOC projects this year? I read from:

https://summerofcode.withgoogle.com/organizations/4726304075677696/#!

1. About Machine Learning
2. MacOS signed app
3. Code Quality
4. Working demo for Matlab MEX modules

... and will anyone post any news about it? Maybe here:
https://www.scilab.org/en/projects/gsoc

Best regards,
Claus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180507/9c105bdf/attachment.htm>

From sgougeon at free.fr  Mon May  7 21:01:04 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Mon, 7 May 2018 21:01:04 +0200
Subject: [Scilab-users] GSOC 2018
In-Reply-To: <CA+Q5PwtZ4FLSeQJFAnRbs20PEjvTnihNqm2QxudcqoPUJi2i6Q@mail.gmail.com>
References: <CA+Q5PwtZ4FLSeQJFAnRbs20PEjvTnihNqm2QxudcqoPUJi2i6Q@mail.gmail.com>
Message-ID: <16a3bb71-9377-89e8-17a0-3735eded3e88@free.fr>

Hello Claus,

AFAIK, there is this portal on the wiki : 
https://wiki.scilab.org/GSoC_project_proposal
and this dedicated and active mailing list : 
http://mailinglists.scilab.org/Scilab-GSOC-Mailing-Lists-Archives-f2646148.html

But you are right : it would be worthwhile having official news here (on 
this list) about the selected projects and about their outcomes, status, 
and chances and schedule to see them available as ATOMS modules, or even 
afterwards  being included in Scilab for some of them.

Regards
Samuel

Le 07/05/2018 ? 20:45, Claus Futtrup a ?crit :
> Hi there
>
> Can Scilab confirm that there are four GSOC projects this year? I read 
> from:
>
> https://summerofcode.withgoogle.com/organizations/4726304075677696/#!
>
> 1. About Machine Learning
> 2. MacOS signed app
> 3. Code Quality
> 4. Working demo for Matlab MEX modules
>
> ... and will anyone post any news about it? Maybe here: 
> https://www.scilab.org/en/projects/gsoc
>
> Best regards,
> Claus
>
>
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180507/9de9a03b/attachment.htm>

From delaleau at enib.fr  Fri May 11 09:42:27 2018
From: delaleau at enib.fr (Emmanuel Delaleau)
Date: Fri, 11 May 2018 09:42:27 +0200
Subject: [Scilab-users] How to close/kill "messagebox"
Message-ID: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>

Hi,


Does it exists a way to close/kill windows created by "messagebox".

"close" or "delete" do not work as messageboxes do not have a handle.

Thanks

Emmanuel



From amonmayr at laas.fr  Fri May 11 10:26:14 2018
From: amonmayr at laas.fr (amonmayr at laas.fr)
Date: Fri, 11 May 2018 10:26:14 +0200
Subject: [Scilab-users] How to close/kill "messagebox"
In-Reply-To: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
Message-ID: <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>


Hello,


Indeed, messagebox() seems to live outside the realm of figures and 
other graphic entities.
It's a bit weird in my opinion.
Maybe a workaround for you would be to emulate what messagebox() is 
doing by using a figure and a bunch of uicontrols and callbacks.
That way you'll be able to treat this emulated version like any othe 
figure and close it using delete(winsid_of_said_figure).

Hope it helps,

Antoine



Le 11/05/2018 ? 09:42, Emmanuel Delaleau a ?crit?:



> Hi,
>
>
> Does it exists a way to close/kill windows created by "messagebox".
>
> "close" or "delete" do not work as messageboxes do not have a handle.
>
> Thanks
>
> Emmanuel
>
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
>

-- 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++

  Antoine Monmayrant LAAS - CNRS
  7 avenue du Colonel Roche
  BP 54200
  31031 TOULOUSE Cedex 4
  FRANCE

  Tel:+33 5 61 33 64 59
  
  email : antoine.monmayrant at laas.fr
  permanent email : antoine.monmayrant at polytechnique.org

+++++++++++++++++++++++++++++++++++++++++++++++++++++++



From sgougeon at free.fr  Fri May 11 11:30:05 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Fri, 11 May 2018 11:30:05 +0200
Subject: [Scilab-users] gcf().anti_aliasing : Confirmation of (un)effect?
In-Reply-To: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>
References: <3e7e61f7-ff70-62d5-9a15-ab4a772e3ec0@free.fr>
Message-ID: <25f10094-ccdf-9790-b03f-6ee496ca53ed@free.fr>

Le 07/05/2018 ? 14:36, Samuel Gougeon a ?crit :
>
> Hello,
>
> The gcf().anti_aliasing figure property was implemented many years ago 
> (in the early Scilab 5.x versions).
>
> It aims to reduce the "staircase effect" for lines drawn very close to 
> vertical or horizontal directions.
>
> Yet, has anyone ever seen an actual effect on graphic rendering, as 
> expected and described in the related SEP#16 
> <http://gitweb.scilab.org/?p=scilab.git;a=blob;f=SEP/SEP_016_anti_aliasing.odt;h=e267088c7ed7e5a4cd0ae29b744616332423224e;hb=HEAD> 
> ?
> I never did, whatever were the computer and its graphic card i used.
> Still now with Scilab 6.0, i don't see any effect.
>
> The staircase effect can easily be seen with the following code:
>
> x = [0 1];
> y = (0:0.02:1)';
> y = [y y+0.02]';
> clf
> plot2d(x,y)
> f = gcf();
> f.anti_aliasing  // should display "off" (by default)
>
> f = scf();
> f.anti_aliasing = "8x";
> plot2d(x,y)
> Then compare both figures.
> I get the following, whatever is the .anti_aliasing value.
>
>
> Does anyone get anything else, closer to the expected "dithered" 
> rendering?
> Either it will be possible to better document the effect of this property,
> as requested in the bug report 6777 
> <http://bugzilla.scilab.org/show_bug.cgi?id=6777>, or to reopen the 
> bug 11311 <http://bugzilla.scilab.org/show_bug.cgi?id=11311> reporting
> that it is out of work.

Thanks Rafael & Antoine.
I have finally got it. It works

  * only with usecanvas(%t). The bug 11311
    <http://bugzilla.scilab.org/show_bug.cgi?id=11311> claimed that it
    was killed in usecanvas(%t) mode. Maybe fixing it has killed it in
    usecanvas(%f) mode...

  * only up to Scilab 5.4.1 included. Scilab 5.5.0 has killed any
    graphic rendering in usecanvas(%t) mode (and has changed the
    usecanvas() type of output from boolean to integer).

With 5.4.1, usecanvas(%t) displays:


which is the expected rendering.

I am going to see what's on bugzilla about usecanvas()...

Samuel



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180511/2d498d27/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: jpndpeifadohoabc.png
Type: image/png
Size: 10769 bytes
Desc: not available
URL: <https://lists.scilab.org/pipermail/users/attachments/20180511/2d498d27/attachment.png>

From Clement.David at esi-group.com  Tue May 15 17:36:38 2018
From: Clement.David at esi-group.com (=?utf-8?B?Q2zDqW1lbnQgRGF2aWQ=?=)
Date: Tue, 15 May 2018 15:36:38 +0000
Subject: [Scilab-users] How to close/kill "messagebox"
In-Reply-To: <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
 <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
Message-ID: <1526398598.2108.4.camel@esi-group.com>

Hello,

Exactly messagebox() use a semi-hardcoded Java-based GUI that totally differs from the uicontrol
code. If you want to mix messagebox() and regular figure(), you should (as Antoine suggested)
implement a popup-like by yourself (should be ~20 lines of code).

Thanks,

--
Cl?ment

Le vendredi 11 mai 2018 ? 10:26 +0200, amonmayr at laas.fr a ?crit :
> Hello,
> 
> 
> Indeed, messagebox() seems to live outside the realm of figures and 
> other graphic entities.
> It's a bit weird in my opinion.
> Maybe a workaround for you would be to emulate what messagebox() is 
> doing by using a figure and a bunch of uicontrols and callbacks.
> That way you'll be able to treat this emulated version like any othe 
> figure and close it using delete(winsid_of_said_figure).
> 
> Hope it helps,
> 
> Antoine
> 
> 
> 
> Le 11/05/2018 ? 09:42, Emmanuel Delaleau a ?crit :
> 
> 
> 
> > Hi,
> > 
> > 
> > Does it exists a way to close/kill windows created by "messagebox".
> > 
> > "close" or "delete" do not work as messageboxes do not have a handle.
> > 
> > Thanks
> > 
> > Emmanuel
> > 
> > _______________________________________________
> > users mailing list
> > users at lists.scilab.org
> > http://lists.scilab.org/mailman/listinfo/users
> > 
> 
> 

From heinznabielek at me.com  Wed May 16 00:20:16 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Wed, 16 May 2018 00:20:16 +0200
Subject: [Scilab-users] any suggestion to make my Scilab code faster would
	be welcome:
In-Reply-To: <1526398598.2108.4.camel@esi-group.com>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
 <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
 <1526398598.2108.4.camel@esi-group.com>
Message-ID: <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>

Friends and colleagues:

I have already a pretty efficient Scilab code, but any suggestion to make it faster would be welcome:

I have position measurements of 20,000 objects in 3d space stored in column vectors X, Y, Z and want to explore the distribution of nearest neighbour distances. This works well like below, but an acceleration would help me a lot.

Heinz

n=length(X);
d=[];
for i=1:n; 
x=X(i);y=Y(i);z=Z(i);
SqDist=(X-x)^2 + (Y-y)^2 + (Z-z)^2;
SqDist(i)=[];
d=[d; sqrt(min(SqDist))];
end;


.. and 'd' gives me the desired distribution in a column vector of n elements

From kopac.jakub at gmail.com  Wed May 16 08:40:13 2018
From: kopac.jakub at gmail.com (kjubo)
Date: Tue, 15 May 2018 23:40:13 -0700 (MST)
Subject: [Scilab-users] any suggestion to make my Scilab code faster
 would be welcome:
In-Reply-To: <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
 <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
 <1526398598.2108.4.camel@esi-group.com>
 <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>
Message-ID: <1526452813137-0.post@n3.nabble.com>

Hello,

d = zeros(n,1);
for ii=1:n;
//    SqDist   = (X-X(i))^2 + (Y-Y(i))^2 + (Z-Z(i))^2;
    SqDist   = (X-X(ii)).*(X-X(ii)) + (Y-Y(ii)).*(Y-Y(ii)) +
(Z-Z(ii)).*(Z-Z(ii));
    SqDist(ii)= [];
    d(ii)    = sqrt(min(SqDist));
end; 

some comments:
* preallocate matrix
* avoid unnecessary assign of variables (this slow down scilab more than
expected)
* as seen in previous meesages here, use X*X instead of X^2

for me speed up factor is ~10

Regards.



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From wozai at protonmail.ch  Wed May 16 10:53:47 2018
From: wozai at protonmail.ch (Wozai)
Date: Wed, 16 May 2018 01:53:47 -0700 (MST)
Subject: [Scilab-users] any suggestion to make my Scilab code faster
 would be welcome:
In-Reply-To: <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
 <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
 <1526398598.2108.4.camel@esi-group.com>
 <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>
Message-ID: <1526460827855-0.post@n3.nabble.com>

Hi,

A very similar problem was discussed
here
<http://mailinglists.scilab.org/Scilab-users-need-a-more-efficient-and-faster-code-suggestions-welcome-tc4037482.html#a4037515>  
in length. You'll find there numerous hints how to make your code faster.

Cheers

Wozai



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From heinznabielek at me.com  Wed May 16 19:44:18 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Wed, 16 May 2018 19:44:18 +0200
Subject: [Scilab-users] any suggestion to make my Scilab code faster
 would be welcome:
In-Reply-To: <1526452813137-0.post@n3.nabble.com>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
 <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
 <1526398598.2108.4.camel@esi-group.com>
 <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>
 <1526452813137-0.post@n3.nabble.com>
Message-ID: <1B2CD69D-63BD-41E0-81A0-E3B6BF07A13B@me.com>

Yes- factor 10. Great many thanks for all the help.
Heinz

> On 16.05.2018, at 08:40, kjubo <kopac.jakub at gmail.com> wrote:
> 
> Hello,
> 
> d = zeros(n,1);
> for ii=1:n;
> //    SqDist   = (X-X(i))^2 + (Y-Y(i))^2 + (Z-Z(i))^2;
>    SqDist   = (X-X(ii)).*(X-X(ii)) + (Y-Y(ii)).*(Y-Y(ii)) +
> (Z-Z(ii)).*(Z-Z(ii));
>    SqDist(ii)= [];
>    d(ii)    = sqrt(min(SqDist));
> end; 
> 
> some comments:
> * preallocate matrix
> * avoid unnecessary assign of variables (this slow down scilab more than
> expected)
> * as seen in previous meesages here, use X*X instead of X^2
> 
> for me speed up factor is ~10
> 
> Regards.



From heinznabielek at me.com  Wed May 16 19:58:03 2018
From: heinznabielek at me.com (Heinz Nabielek)
Date: Wed, 16 May 2018 19:58:03 +0200
Subject: [Scilab-users] any suggestion to make my Scilab code faster
 would be welcome:
In-Reply-To: <1B2CD69D-63BD-41E0-81A0-E3B6BF07A13B@me.com>
References: <a203dcdd-f8bb-b379-fea3-7efb19f126a2@enib.fr>
 <d1e10304-42f4-e371-2518-57be454eac77@laas.fr>
 <1526398598.2108.4.camel@esi-group.com>
 <789033FA-02C6-4CAC-81EB-F6600C151760@me.com>
 <1526452813137-0.post@n3.nabble.com>
 <1B2CD69D-63BD-41E0-81A0-E3B6BF07A13B@me.com>
Message-ID: <9428E462-D562-4455-B87A-F48207E867C8@me.com>

Histogram shows distribution of d and comparison with my model prediction.
Heinz


> On 16.05.2018, at 19:44, Heinz Nabielek <heinznabielek at me.com> wrote:
> 
> Yes- factor 10. Great many thanks for all the help.
> Heinz
> 
>> On 16.05.2018, at 08:40, kjubo <kopac.jakub at gmail.com> wrote:
>> 
>> Hello,
>> 
>> d = zeros(n,1);
>> for ii=1:n;
>> //    SqDist   = (X-X(i))^2 + (Y-Y(i))^2 + (Z-Z(i))^2;
>>   SqDist   = (X-X(ii)).*(X-X(ii)) + (Y-Y(ii)).*(Y-Y(ii)) +
>> (Z-Z(ii)).*(Z-Z(ii));
>>   SqDist(ii)= [];
>>   d(ii)    = sqrt(min(SqDist));
>> end; 
>> 
>> some comments:
>> * preallocate matrix
>> * avoid unnecessary assign of variables (this slow down scilab more than
>> expected)
>> * as seen in previous meesages here, use X*X instead of X^2
>> 
>> for me speed up factor is ~10
>> 
>> Regards.
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180516/02bc19e7/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-1.tiff
Type: image/tiff
Size: 51494 bytes
Desc: not available
URL: <https://lists.scilab.org/pipermail/users/attachments/20180516/02bc19e7/attachment.tiff>

From j.s.strom at hslmg.de  Sat May 19 12:35:48 2018
From: j.s.strom at hslmg.de (Jens Simon Strom)
Date: Sat, 19 May 2018 12:35:48 +0200
Subject: [Scilab-users] Exclamation mark in variable names
Message-ID: <5AFFFE04.4030207@hslmg.de>

Hallo Scilab intimates,
Is *!a*     a valid variable name in the latest version of Scilab?  
Vers. 6.0.1 sees an*Error: Unexpected token '!'*

Kind regards
Jens
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180519/216db738/attachment.htm>

From sgougeon at free.fr  Sat May 19 13:19:21 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sat, 19 May 2018 13:19:21 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <5AFFFE04.4030207@hslmg.de>
References: <5AFFFE04.4030207@hslmg.de>
Message-ID: <b2466d98-ba7f-f649-8c9b-fd6ff941d48b@free.fr>

Hello Jens,

Le 19/05/2018 ? 12:35, Jens Simon Strom a ?crit :
> Hallo Scilab intimates,
> Is *!a*     a valid variable name in the latest version of Scilab?  
> Vers. 6.0.1 sees an*Error: Unexpected token '!'*

Yes it is. You have reported it as the bug 15126 
<http://bugzilla.scilab.org/show_bug.cgi?id=15126>, and Clement has 
confirmed to you that this naming should still be supported.

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180519/9ad458df/attachment.htm>

From sgougeon at free.fr  Mon May 21 14:10:19 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Mon, 21 May 2018 14:10:19 +0200
Subject: [Scilab-users] det([]) changed in Scilab 6.0
Message-ID: <c4b3cebc-d7fb-9458-03a4-8abe1b866210@free.fr>

Hello,

While updating "help empty" after changes in Scilab 6 about inequality 
comparisons with [],
i have found that det([]) is changed:

  * 5.5: det([]) => 1
  * 6.0: det([]) => 0

I don't know if this modification is intentional.
I am wondering about the added value of this change:

  * If we consider the property for square matrices A & B:
            det(A*B) = det(A)*det(B)    (1)
    setting det([])=0 better complies with it. Indeed, since []*B => [],
    we have now det([]*B) = det([])*det(B) = 0*k = 0 = det([])  : OK

  * In another hand, we should also have :
    det(1/A) = 1/det(A)     (2)
    Since inv([]) => [], setting det([]) => 1  complied with (2),
    while det([]) => 0  no longer does it.

Settingdet([]) => [] would comply with both (1) and (2).

Let's note that Octave returns 1.

So, shall we restore det([]) ==1, or keep 0, or set it to [] ?

What's the improvement in changing det([]) to 0 ?

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180521/1207601e/attachment.htm>

From sgougeon at free.fr  Mon May 21 15:46:21 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Mon, 21 May 2018 15:46:21 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
Message-ID: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>

Hello,

Another case with []:
The bug 15579 <http://bugzilla.scilab.org/show_bug.cgi?id=15579> reports 
that *cond([])* returns *1*, while *cond([],2)* -- that computes the 
same 2-norm conditional number for any other matrix than [] -- returns *0*.

Actually, *cond([], p)* returns 0 for any accepted p value.

At least since Scilab 4, and likely since the beginning of Scilab, 
*cond([])* is set to 1 as a specific case.
I am wondering why.

In my opinion, the mismatch between *cond([])* and the equivalent 
*cond([],2)* must be fixed.
There are at least two way to do it :

  * either keep *cond([])* to *1* and set all *cond([], p)* to 1 instead
    of 0
  * or set *cond([])* to *0*.

I don't see any clear reason enforcing a choice rather than the other.
Do you?

Let's note that

  * Octave returns 0 in both cond([]) and cond([],p) cases
  * Julia, Python, and R (with kappa()) yield an error

Regards
Samuel


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180521/3de4e049/attachment.htm>

From jrafaelbguerra at hotmail.com  Mon May 21 20:25:20 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Mon, 21 May 2018 18:25:20 +0000
Subject: [Scilab-users] det([]) changed in Scilab 6.0
In-Reply-To: <c4b3cebc-d7fb-9458-03a4-8abe1b866210@free.fr>
References: <c4b3cebc-d7fb-9458-03a4-8abe1b866210@free.fr>
Message-ID: <VI1PR04MB4109EDC13C0AC3FD1200E90ACC950@VI1PR04MB4109.eurprd04.prod.outlook.com>

Hi Samuel,

I am not a mathematician but I suspect that determinant rules you have enunciated do not hold for empty matrices.

Fyi, the following Wikipedia?s entry states that det([]) should be 1 (see ?Empty matrices? paragraph):
https://en.wikipedia.org/wiki/Matrix_(mathematics)#Empty_matrices

Regards,
Rafael

From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Gougeon
Sent: Monday, May 21, 2018 2:10 PM
To: International users mailing list for Scilab. <users at lists.scilab.org>
Subject: [Scilab-users] det([]) changed in Scilab 6.0

Hello,

While updating "help empty" after changes in Scilab 6 about inequality comparisons with [],
i have found that det([]) is changed:

  *   5.5: det([]) => 1
  *   6.0: det([]) => 0

I don't know if this modification is intentional.
I am wondering about the added value of this change:

  *   If we consider the property for square matrices A & B:
       det(A*B) = det(A)*det(B)    (1)
setting det([])=0 better complies with it. Indeed, since []*B => [],
we have now det([]*B) = det([])*det(B) = 0*k = 0 = det([])  : OK
  *   In another hand, we should also have :
            det(1/A) = 1/det(A)     (2)
Since inv([]) => [], setting det([]) => 1  complied with (2),
while det([]) => 0  no longer does it.

Setting det([]) => [] would comply with both (1) and (2).

Let's note that Octave returns 1.

So, shall we restore det([]) ==1, or keep 0, or set it to [] ?

What's the improvement in changing det([]) to 0 ?

Regards
Samuel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180521/84fc1f5f/attachment.htm>

From jrafaelbguerra at hotmail.com  Mon May 21 20:58:51 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Mon, 21 May 2018 18:58:51 +0000
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
Message-ID: <VI1PR04MB410964CF69AFD2028482DCCBCC950@VI1PR04MB4109.eurprd04.prod.outlook.com>

Hi Samuel,

Please read the following article on this interesting topic of the algebra of empty matrices:
https://pdfs.semanticscholar.org/0f3b/c36f19d5c6a761c19fbc3c4ebde2f31b0a10.pdf

It states that the condition number of the empty matrix should be 0.

Regards,
Rafael

From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Gougeon
Sent: Monday, May 21, 2018 3:46 PM
To: International users mailing list for Scilab. <users at lists.scilab.org>
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)

Hello,

Another case with []:
The bug 15579<http://bugzilla.scilab.org/show_bug.cgi?id=15579> reports that cond([]) returns 1, while cond([],2) -- that computes the same 2-norm conditional number for any other matrix than [] -- returns 0.

Actually, cond([], p) returns 0 for any accepted p value.

At least since Scilab 4, and likely since the beginning of Scilab,  cond([]) is set to 1 as a specific case.
I am wondering why.

In my opinion, the mismatch between cond([]) and the equivalent cond([],2) must be fixed.
There are at least two way to do it :

  *   either keep cond([]) to 1 and set all cond([], p) to 1 instead of 0
  *   or set cond([]) to 0.

I don't see any clear reason enforcing a choice rather than the other.
Do you?

Let's note that

  *   Octave returns 0 in both cond([]) and cond([],p) cases
  *   Julia, Python, and R (with kappa()) yield an error
Regards
Samuel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180521/636559fb/attachment.htm>

From sgougeon at free.fr  Mon May 21 21:25:50 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Mon, 21 May 2018 21:25:50 +0200
Subject: [Scilab-users] det([]) changed in Scilab 6.0
In-Reply-To: <VI1PR04MB4109EDC13C0AC3FD1200E90ACC950@VI1PR04MB4109.eurprd04.prod.outlook.com>
References: <c4b3cebc-d7fb-9458-03a4-8abe1b866210@free.fr>
 <VI1PR04MB4109EDC13C0AC3FD1200E90ACC950@VI1PR04MB4109.eurprd04.prod.outlook.com>
Message-ID: <e4c86243-d556-314a-b08e-f0d6215ee4a8@free.fr>

Le 21/05/2018 ? 20:25, Rafael Guerra a ?crit :
>
> Hi Samuel,
>
> I am not a mathematician but I suspect that determinant rules you have 
> enunciated do not hold for empty matrices.
>
> Fyi, the following Wikipedia?s entry states that det([]) should be 1 
> (see ?Empty matrices? paragraph):
>
> https://en.wikipedia.org/wiki/Matrix_(mathematics)#Empty_matrices 
> <https://en.wikipedia.org/wiki/Matrix_%28mathematics%29#Empty_matrices>
>

Thanks Rafael!

This bug is then reported as the bug 15580 
<http://bugzilla.scilab.org/show_bug.cgi?id=15580>.

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180521/99ca43d1/attachment.htm>

From j_holland at msn.com  Tue May 22 12:28:40 2018
From: j_holland at msn.com (James Holland)
Date: Tue, 22 May 2018 03:28:40 -0700 (MST)
Subject: [Scilab-users] Plot2d and legends
Message-ID: <1526984920350-0.post@n3.nabble.com>

I am plotting three columns of data using plot2d and I'm adding a legend box
using legends. I have declared the line colours for each plot to make the
graph easier to view. The problem that I have is that the legends box shows
all three of the lines in the same colour (blue). Is this an issue with the
way my code is written?
TIA James

<code>
// plot data 

// the third column: -->m(:,3)
x1 = M1(:,1);   // frequency
y1 = M1(:,2);   // M1 Average
y2 = M2(:,2);   // M2 Average
x2 = M3(:,1);   // Limits frequency
y3 = M3(:,2);   // Limits - average


plot2d(x1, y1, style= color('orange'), logflag = "ln");
plot2d(x1, y2, style= color('blue'), logflag = "ln");
plot2d(x2, y3, style= color('red'), logflag = "ln");
f = gca();                     // get the current axes
bg_color = color('lightgray');              
f.background= bg_color;        //RGB
xtitle("EN55022 Conducted Emissions", "Frequency (MHz)", "dBuV");
legends(['uA7805';'7805SR'; 'EN55022B'],[2,2,2],opt="ur");
xgrid;
</code>



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From stephane.mottelet at utc.fr  Tue May 22 12:44:01 2018
From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=)
Date: Tue, 22 May 2018 12:44:01 +0200
Subject: [Scilab-users] Plot2d and legends
In-Reply-To: <1526984920350-0.post@n3.nabble.com>
References: <1526984920350-0.post@n3.nabble.com>
Message-ID: <a6c8ef16-1f8f-a81b-602d-b0ed1f05363b@utc.fr>

Hello,

Le 22/05/2018 ? 12:28, James Holland a ?crit?:
> I am plotting three columns of data using plot2d and I'm adding a legend box
> using legends. I have declared the line colours for each plot to make the
> graph easier to view. The problem that I have is that the legends box shows
> all three of the lines in the same colour (blue). Is this an issue with the
> way my code is written?
> TIA James
>
> <code>
> // plot data
>
> // the third column: -->m(:,3)
> x1 = M1(:,1);   // frequency
> y1 = M1(:,2);   // M1 Average
> y2 = M2(:,2);   // M2 Average
> x2 = M3(:,1);   // Limits frequency
> y3 = M3(:,2);   // Limits - average
>
>
> plot2d(x1, y1, style= color('orange'), logflag = "ln");
> plot2d(x1, y2, style= color('blue'), logflag = "ln");
> plot2d(x2, y3, style= color('red'), logflag = "ln");
> f = gca();                     // get the current axes
> bg_color = color('lightgray');
> f.background= bg_color;        //RGB
> xtitle("EN55022 Conducted Emissions", "Frequency (MHz)", "dBuV");
> legends(['uA7805';'7805SR'; 'EN55022B'],[2,2,2],opt="ur");

here [2 2 2] overwrites the colors you specified in the plot2d calls. I 
would suggest:

legend('uA7805','7805SR','EN55022B','in_upper_right')

S.
> xgrid;
> </code>
>
>
>
> --
> Sent from: https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users


-- 
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 j_holland at msn.com  Thu May 24 15:40:43 2018
From: j_holland at msn.com (James Holland)
Date: Thu, 24 May 2018 06:40:43 -0700 (MST)
Subject: [Scilab-users] Sockets in scilab 6.0
Message-ID: <1527169243922-0.post@n3.nabble.com>

I want to communicate with an oscilloscope through an Ethernet connection. I
have downloaded the files linked from the Atoms page.
If I try to run the loader then I get the message 'lib: Old binary lib
detected. Please recompile it for Scilab 6.' I don't see any files that
aren't Scilab so I ran 'builder.sce'. I get the following message:
'--> exec('C:\Users\James.Holland\Documents\Scilab\socket
src\builder.sce',-1)
Building macros...
-- Creation of [socket_toolboxlib] (Macros) --
Building help...

Building the master document:
	C:\Users\James.Holland\Documents\Scilab\socket src\help\en_US

Building the manual file [javaHelp] in
C:\Users\James.Holland\Documents\Scilab\socket src\help\en_US.

Building the master document:
	C:\Users\James.Holland\Documents\Scilab\socket src\help\fr_FR

Building the manual file [javaHelp] in
C:\Users\James.Holland\Documents\Scilab\socket src\help\fr_FR.
WARNING: tbx_build_loader: tbx_build_loader(name, path) is obsolete. Please
use tbx_build_loader(path) instead.
Generating loader.sce...
Generating unloader.sce...
WARNING: tbx_build_cleaner: tbx_build_cleaner(name, path) is obsolete.
Please use tbx_build_cleaner(path) instead.
Generating cleaner.sce...
'

After that I ran the loader which returned:
--> exec('C:\Users\James.Holland\Documents\Scilab\socket src\loader.sce',-1)
Start Socket Toolbox
	Load macros
	Load help
	Load demos
I then get the following error:

<code>
--> exec('C:\Users\James.Holland\Documents\Scilab\socket
src\demos\ex_1.sce',-1)

    []

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat

Undefined variable: mat
mat=["<!doctype html><html itemscope=""""
itemtype=""http://schema.org/WebPage"" lang=""en-GB""><head><meta
content=""text/html; charset=UTF-8"" http-equiv=""Content-Type""><meta
content=""/images/branding/googleg/1x/googleg_standard_color_128dp.png""
itemprop=""image""><title>Google</title><style>#gb{font:13px/27px
Arial,sans-serif;height:30px}#gbz,#gbg{position:absolute;white-space:nowrap;top:0;height:30px;z-index:1000}#gbz{left:0;padding-left:4px}#gbg{right:0;padding-right:5px}#gbs{background:transparent;position:absolute;top:-999px;visibility:hidden;z-index:998;right:0}.gbto
#gbs{background:#fff}#gbx3,#gbx4{background-color:#2d2d2d;background-image:none;_background-image:none;background-position:0
-138px;background-repeat:repeat-x;border-bottom:1px solid
#000;font-size:24px;height:29px;_height:30px;opacity:1;filter:alpha(opacity=100);position:absolute;top:0;width:100%;z-index:990}#gbx3{left:0}#gbx4{right:0}#gbb{position:relative}#gbbw{left:0;position:absolute;top:30px;width:100%}.gbtcb{position:absolute;visibility:hidden}#gbz
.gbtcb{right:0}#gbg .gbtcb{left:0}.gbxx{display:none
!important}.gbxo{opacity:0 !important;filter:alpha(opacity=0)
!important}.gbm{position:absolute;z-index:999;top:-999px;visibility:hidden;text-align:left;border:1px
soli
                                                                                                                                                                                                                                                                                                                                                                             
^^
Error: Heterogeneous string detected, starting with " and ending with '.
</ code>



Has anyone successfully used socket_toolbox_2.0.1_5.5 within Scilab 6.0?
How do I get the toolbox to load on start up like other Atoms modules?

TIA James






--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From j_holland at msn.com  Thu May 24 17:18:45 2018
From: j_holland at msn.com (James Holland)
Date: Thu, 24 May 2018 08:18:45 -0700 (MST)
Subject: [Scilab-users] Sockets in scilab 6.0
In-Reply-To: <1527169243922-0.post@n3.nabble.com>
References: <1527169243922-0.post@n3.nabble.com>
Message-ID: <1527175125191-0.post@n3.nabble.com>

I've played about with this some more. The toolbox did install and is
accessed from the 'toolboxes' tab at the top - I haven't used that before.
The code seems to work except for the syntax problem ' and '' and the
undefined variable that keeps coming up - maybe a caused by there not being
a delay between the write and the read.
If I change the example code from google.fr to another website it doesn't
work, no idea why, I will try it with an oscilloscope next.



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From rouxph.22 at gmail.com  Thu May 24 23:20:55 2018
From: rouxph.22 at gmail.com (philippe)
Date: Thu, 24 May 2018 23:20:55 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
Message-ID: <pe7a7j$et5$1@blaine.gmane.org>

Le 21/05/2018 ? 15:46, Samuel Gougeon a ?crit?:
> 
> There are at least two way to do it :
> 
>   * either keep *cond([])* to *1* and set all *cond([], p)* to 1 instead
>     of 0
>   * or set *cond([])* to *0*.
> 
> I don't see any clear reason enforcing a choice rather than the other.
> Do you?

from a mathematical point of view  the Condition number of a matrix A is
defined by

cond(A)=||A|| . ||A^(-1)||
its interest is that when solving A*x=y rounding errors on y (eps) are
amplified to be cond(A)*eps. The optimal value of cond(A) is 1 (for
identity matrix) so for me it looks natural that cond([])=1 .

Philippe



From wozai at protonmail.ch  Fri May 25 10:17:25 2018
From: wozai at protonmail.ch (Wozai)
Date: Fri, 25 May 2018 01:17:25 -0700 (MST)
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <5AFFFE04.4030207@hslmg.de>
References: <5AFFFE04.4030207@hslmg.de>
Message-ID: <1527236245681-0.post@n3.nabble.com>

Hi,

I have had the same problem, and found a solution (at least if you build
scilab by yourself), apply the following changes to the source and compile
Scilab.

1.  Changes Hotey "!" to ":" for Recalling History in Java Console
<https://github.com/rdbyk/balisc/pull/320>  

2.  Fixes Parser to Accept "!" in Variable Names
<https://github.com/rdbyk/balisc/pull/323>  

Of course the downside of this fixing approach would be, that your code only
works in "your" Scilab.
 
Cheers,

Wozai




--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From Clement.David at esi-group.com  Fri May 25 10:45:25 2018
From: Clement.David at esi-group.com (=?utf-8?B?Q2zDqW1lbnQgRGF2aWQ=?=)
Date: Fri, 25 May 2018 08:45:25 +0000
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <1527236245681-0.post@n3.nabble.com>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
Message-ID: <1527237924.2237.52.camel@esi-group.com>

Hi,

Thanks for your experimental work on that topic. I guess we could merge the parser fix easily
however the `!` usage to recall a previous command will remain problematic. This usage does not seem
to be documented somewhere in the help page, maybe we could discard it. What's your opinion on that
?

Thanks,

--
Cl?ment

Le vendredi 25 mai 2018 ? 01:17 -0700, Wozai a ?crit :
> Hi,
> 
> I have had the same problem, and found a solution (at least if you build
> scilab by yourself), apply the following changes to the source and compile
> Scilab.
> 
> 1.  Changes Hotey "!" to ":" for Recalling History in Java Console
> <https://github.com/rdbyk/balisc/pull/320>? 
> 
> 2.  Fixes Parser to Accept "!" in Variable Names
> <https://github.com/rdbyk/balisc/pull/323>? 
> 
> Of course the downside of this fixing approach would be, that your code only
> works in "your" Scilab.
>  
> Cheers,
> 
> Wozai
> 
> 
> 
> 
> --
> 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 sgougeon at free.fr  Fri May 25 11:02:20 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Fri, 25 May 2018 11:02:20 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <1527237924.2237.52.camel@esi-group.com>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
Message-ID: <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>

Hello,

Le 25/05/2018 ? 10:45, Cl?ment David a ?crit :
> Hi,
>
> Thanks for your experimental work on that topic. I guess we could merge the parser fix easily
> however the `!` usage to recall a previous command will remain problematic. This usage does not seem
> to be documented somewhere in the help page, maybe we could discard it. What's your opinion on that
> ?

It is an historical and standard feature, used not only in Scilab language.
Please see the entry http://bugzilla.scilab.org/10776

The only issue is that it is still undocumented. It could be described 
in the Console's page.

Samuel



From wozai at protonmail.ch  Fri May 25 11:08:25 2018
From: wozai at protonmail.ch (Wozai)
Date: Fri, 25 May 2018 02:08:25 -0700 (MST)
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
Message-ID: <1527239305614-0.post@n3.nabble.com>

Hi,

I think it should be possible to execute things like "!A=123" in the Java
Console.

Cheers,

Wozai



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From sgougeon at free.fr  Fri May 25 11:16:44 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Fri, 25 May 2018 11:16:44 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <1527239305614-0.post@n3.nabble.com>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <1527239305614-0.post@n3.nabble.com>
Message-ID: <e3084182-fcf3-5692-6a4f-a679fe527396@free.fr>

Le 25/05/2018 ? 11:08, Wozai a ?crit :
> Hi,
>
> I think it should be possible to execute things like "!A=123" in the Java
> Console.

Just prepend a space, and it works.

Samuel



From cfuttrup at gmail.com  Fri May 25 19:17:20 2018
From: cfuttrup at gmail.com (Claus Futtrup)
Date: Fri, 25 May 2018 19:17:20 +0200
Subject: [Scilab-users] Simulating phase modulation
Message-ID: <CA+Q5PwtAYBv9zeq07O6x9VkULVHCczP9G-jXPdQ1tb6iDskkvA@mail.gmail.com>

Hi there

In a loudspeaker the driver can move several millimeter in an attempt to
reproduce a low-frequency note. If the speaker also at the same time
produce a higher tone, this second tone is phase modulated by the first
one. This is a distortion of the original signal which I'd like to simulate
/ illustrate with some simple Scilab code, if possible.

In Matlab this can be simulated with pmmod.
https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/

Is there a similar function in Scilab? (name - please ?)

Here's the code I have written so far - this is the part that shows the
input signal (the un-distorted signal):

sample_rate=20000;t = 0:1/sample_rate:0.6;N=size(t,'*'); //number of
samplesy1 = sin(2*%pi*50*t);y2 = 0.5*sin(2*%pi*500*t);// y2 =
0.5*sin(2*%pi*500*t+%pi/4);s=y1+y2+grand(1,N,'nor',0,1);
// Plot time-domainendplot = round(N/15);twoplots = scf(); // Set
Current Figure (Graphics
Window)subplot(211);plot(t(1:endplot),y1(1:endplot),t(1:endplot),y2(1:endplot));subplot(212);plot(t(1:endplot),y1(1:endplot)+y2(1:endplot));
y=fft(s);ymax = max(abs(y));y = y ./ ymax; // Normalize
// s is real so the fft response is conjugate symmetric// and we
retain only the first N/2 pointsf=sample_rate*(0:(N/2))/N;
//associated frequency vectorn=size(f,'*');fftplots =
scf();plot(f(2:$),abs(y(2:n))); // drop first datapoint, f = 0 (it
prevents log-plot)a = gca();a.log_flags = "lnn";


Best regards,
Claus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180525/fed687ec/attachment.htm>

From jrafaelbguerra at hotmail.com  Fri May 25 20:42:23 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Fri, 25 May 2018 18:42:23 +0000
Subject: [Scilab-users] Simulating phase modulation
In-Reply-To: <CA+Q5PwtAYBv9zeq07O6x9VkULVHCczP9G-jXPdQ1tb6iDskkvA@mail.gmail.com>
References: <CA+Q5PwtAYBv9zeq07O6x9VkULVHCczP9G-jXPdQ1tb6iDskkvA@mail.gmail.com>
Message-ID: <VI1PR04MB4109CEA1679F0CE45050B684CC690@VI1PR04MB4109.eurprd04.prod.outlook.com>

Hi Claus,

I am not aware of such function. However, you can find simple code here below for both phase modulation and demodulation, which is straightforward to translate in Scilab:
https://www.gaussianwaves.com/2017/06/phase-demodulation-using-hilbert-transform-application-of-analytic-signal/

Note that the phase modulation is coded differently from you snippet below.

Regards,
Rafael

From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Claus Futtrup
Sent: Friday, May 25, 2018 7:17 PM
To: International users mailing list for Scilab. <users at lists.scilab.org>
Subject: [Scilab-users] Simulating phase modulation

Hi there

In a loudspeaker the driver can move several millimeter in an attempt to reproduce a low-frequency note. If the speaker also at the same time produce a higher tone, this second tone is phase modulated by the first one. This is a distortion of the original signal which I'd like to simulate / illustrate with some simple Scilab code, if possible.

In Matlab this can be simulated with pmmod.
https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/

Is there a similar function in Scilab? (name - please ?)

Here's the code I have written so far - this is the part that shows the input signal (the un-distorted signal):


sample_rate=20000;

t = 0:1/sample_rate:0.6;

N=size(t,'*'); //number of samples

y1 = sin(2*%pi*50*t);

y2 = 0.5*sin(2*%pi*500*t);

// y2 = 0.5*sin(2*%pi*500*t+%pi/4);

s=y1+y2+grand(1,N,'nor',0,1);



// Plot time-domain

endplot = round(N/15);

twoplots = scf(); // Set Current Figure (Graphics Window)

subplot(211);

plot(t(1:endplot),y1(1:endplot),t(1:endplot),y2(1:endplot));

subplot(212);

plot(t(1:endplot),y1(1:endplot)+y2(1:endplot));



y=fft(s);

ymax = max(abs(y));

y = y ./ ymax; // Normalize



// s is real so the fft response is conjugate symmetric

// and we retain only the first N/2 points

f=sample_rate*(0:(N/2))/N; //associated frequency vector

n=size(f,'*');

fftplots = scf();

plot(f(2:$),abs(y(2:n))); // drop first datapoint, f = 0 (it prevents log-plot)

a = gca();

a.log_flags = "lnn";

Best regards,
Claus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180525/952c8c1d/attachment.htm>

From sgougeon at free.fr  Sat May 26 14:10:55 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sat, 26 May 2018 14:10:55 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
Message-ID: <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr>

Le 25/05/2018 ? 11:02, Samuel Gougeon a ?crit :
> Hello,
>
> Le 25/05/2018 ? 10:45, Cl?ment David a ?crit :
>> Hi,
>>
>> Thanks for your experimental work on that topic. I guess we could 
>> merge the parser fix easily
>> however the `!` usage to recall a previous command will remain 
>> problematic. This usage does not seem
>> to be documented somewhere in the help page, maybe we could discard 
>> it. What's your opinion on that
>> ?
>
> It is an historical and standard feature, used not only in Scilab 
> language.
> Please see the entry http://bugzilla.scilab.org/10776
>
> The only issue is that it is still undocumented. It could be described 
> in the Console's page.

"The only issue" : hmm, not really. Even if the missing documentation is 
the main one, IMO this (already) very handy feature could be improved in 
at least two ways:

  * when no match is found in the history, nothing happens. For people
    who does not know the feature (or who has forgotten that it is
    available), this could be easily puzzling and felt as a bug.
    => when no match is found, Scilab shall print a message "No match
    found in history"

  * A listbox listing matches could be displayed, to select the only
    entry searched for.
    This might be done in two ways :

      o either typing "!" + one char as the very first characters would
        enter a completion mode based on the history, with matches
        listed in a listbox as described hereabove

      o or this could be done only after validating with <enter>, or
        CTRL+space as for the feature completion (but here targeted to
        the history with the leading "!")

If a completion mode is implemented, then the message "no match found in 
history" could be displayed in the listbox instead of in the console 
(note that the general feature completion triggered with CTRL+space 
displays nothing when no match is found).

By the way, as for this "!" history recall, the existing interactive 
usage of the completion is not documented in the Console page. Let's add 
both in the page, as they presently are.

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180526/d44a62cc/attachment.htm>

From cfuttrup at gmail.com  Sat May 26 16:06:48 2018
From: cfuttrup at gmail.com (Claus Futtrup)
Date: Sat, 26 May 2018 16:06:48 +0200
Subject: [Scilab-users] Scilab code published in scientific article
Message-ID: <CA+Q5PwuEQh5uhM1=5FE+p0HKPWyt50b_Fcoap3m58EwFmnM8_Q@mail.gmail.com>

Dear Scilabers

Just FYI, I had a scientific article published this weekend. You can see
it  here:
http://www.aes.org/e-lib/browse.cfm?elib=19566

I'm not allowed to freely copy the article, but I hope nobody minds that I
show the Scilab code included in the article. Here it is:

Qts = 0.5;alpha = 0; // alpha = Vas/Vb
           // alpha = 0, open baffleh = 0; // bass reflex system tuning ratio
       // h = 0 for closed box systemsfunction y=resp(s)y =
s^2/(s^2+s/Qts+1+alpha);endfunction // Eq (15)
N0    = 12; // resolutionsteps = 8;  // steps per unit time,t_max = 5;
 // time (dimensionless)t_step = 1/steps;t = t_step:t_step:t_max;
mu_c = max([1 h]);t_c = %pi*N0/(12*mu_c);nt = length(t);
mprintf("   t      -  Contour");mprintf(" method  -  Exact\n");for i=1:nt do
    ti = t(i);
    // Selection defined in
    // Eqs (31) and (32)
    if ti < t_c then
        mu = %pi*N0/(12*ti);
        N = N0;
        d = 3/N;
    else
        mu = mu_c;
        N = ceil(N0*ti/t_c);
        d = 3/N;
    end

    // Perform sum in Eq (24)
    xsk = 0.0;
    for k=-N:N do
        u = k*d;
        s = mu*(%i*u+1)^2;
        dsdu = 2.0*mu*(%i-u);
        xsk = xsk+imag(exp(s*ti)* ..
            resp(s)/s*dsdu);
    end
    xsk = xsk*d/(2*%pi);
    xs(i) = xsk; // Collect x-values
    exact(i) = exp(-ti)*(1-ti);
    // Collect exact result, compare
    mprintf("%5.4e ",ti);
    mprintf("%+12.11e ",xsk);
    mprintf("%+12.11e\n",exact(i));end
g_tr = scf(); // Graph Time Responseplot([0 t],[1;
xs],"-xr");xgrid(color("grey70"));// Plot contour method numerical
errorerr = abs(xs - exact);g_e = scf(); // Graph Errorplot([0 t],[0;
err],"-r");xgrid(color("grey70"));



Best regards,
Claus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180526/efb1e923/attachment.htm>

From cfuttrup at gmail.com  Sat May 26 16:40:19 2018
From: cfuttrup at gmail.com (Claus Futtrup)
Date: Sat, 26 May 2018 16:40:19 +0200
Subject: [Scilab-users] Simulating phase modulation
In-Reply-To: <VI1PR04MB4109CEA1679F0CE45050B684CC690@VI1PR04MB4109.eurprd04.prod.outlook.com>
References: <CA+Q5PwtAYBv9zeq07O6x9VkULVHCczP9G-jXPdQ1tb6iDskkvA@mail.gmail.com>
 <VI1PR04MB4109CEA1679F0CE45050B684CC690@VI1PR04MB4109.eurprd04.prod.outlook.com>
Message-ID: <CA+Q5PwtdVW2+zeh8gwhMu47nRhXG5P9NU2q+2i1_mtDXi2yJMw@mail.gmail.com>

Hi Rafael

Thank you, I shall print and study. :-)

Cheers,
Claus

On Fri, May 25, 2018 at 8:42 PM Rafael Guerra <jrafaelbguerra at hotmail.com>
wrote:

> Hi Claus,
>
>
>
> I am not aware of such function. However, you can find simple code here
> below for both phase modulation and demodulation, which is straightforward
> to translate in Scilab:
>
>
> https://www.gaussianwaves.com/2017/06/phase-demodulation-using-hilbert-transform-application-of-analytic-signal/
>
>
>
> Note that the phase modulation is coded differently from you snippet below.
>
>
>
> Regards,
>
> Rafael
>
>
>
> *From:* users [mailto:users-bounces at lists.scilab.org] *On Behalf Of *Claus
> Futtrup
> *Sent:* Friday, May 25, 2018 7:17 PM
> *To:* International users mailing list for Scilab. <users at lists.scilab.org
> >
> *Subject:* [Scilab-users] Simulating phase modulation
>
>
>
> Hi there
>
>
>
> In a loudspeaker the driver can move several millimeter in an attempt to
> reproduce a low-frequency note. If the speaker also at the same time
> produce a higher tone, this second tone is phase modulated by the first
> one. This is a distortion of the original signal which I'd like to simulate
> / illustrate with some simple Scilab code, if possible.
>
>
>
> In Matlab this can be simulated with pmmod.
>
> https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/
>
>
>
> Is there a similar function in Scilab? (name - please ?)
>
>
>
> Here's the code I have written so far - this is the part that shows the
> input signal (the un-distorted signal):
>
>
>
> sample_rate=20000;
>
> t = 0:1/sample_rate:0.6;
>
> N=size(t,'*'); *//number of samples*
>
> y1 = sin(2*%pi*50*t);
>
> y2 = 0.5*sin(2*%pi*500*t);
>
> *// y2 = 0.5*sin(2*%pi*500*t+%pi/4);*
>
> s=y1+y2+grand(1,N,'nor',0,1);
>
>
>
> *// Plot time-domain*
>
> endplot = round(N/15);
>
> twoplots = scf(); *// Set Current Figure (Graphics Window)*
>
> subplot(211);
>
> plot(t(1:endplot),y1(1:endplot),t(1:endplot),y2(1:endplot));
>
> subplot(212);
>
> plot(t(1:endplot),y1(1:endplot)+y2(1:endplot));
>
>
>
> y=fft(s);
>
> ymax = max(abs(y));
>
> y = y ./ ymax; *// Normalize*
>
>
>
> *// s is real so the fft response is conjugate symmetric*
>
> *// and we retain only the first N/2 points*
>
> f=sample_rate*(0:(N/2))/N; *//associated frequency vector*
>
> n=size(f,'*');
>
> fftplots = scf();
>
> plot(f(2:$),abs(y(2:n))); *// drop first datapoint, f = 0 (it prevents log-plot)*
>
> a = gca();
>
> a.log_flags = "lnn";
>
>
>
> Best regards,
>
> Claus
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180526/c1e878bc/attachment.htm>

From porpoiseseeker at gmail.com  Sat May 26 22:54:45 2018
From: porpoiseseeker at gmail.com (Gary Nelson)
Date: Sat, 26 May 2018 13:54:45 -0700
Subject: [Scilab-users] Simulating phase modulation
In-Reply-To: <CA+Q5PwtdVW2+zeh8gwhMu47nRhXG5P9NU2q+2i1_mtDXi2yJMw@mail.gmail.com>
References: <CA+Q5PwtAYBv9zeq07O6x9VkULVHCczP9G-jXPdQ1tb6iDskkvA@mail.gmail.com>
 <VI1PR04MB4109CEA1679F0CE45050B684CC690@VI1PR04MB4109.eurprd04.prod.outlook.com>
 <CA+Q5PwtdVW2+zeh8gwhMu47nRhXG5P9NU2q+2i1_mtDXi2yJMw@mail.gmail.com>
Message-ID: <5b09c999.1c69fb81.23238.a038@mx.google.com>

Claus,
Back in the late 60s, I did my PhD research using analytic signal. Yes,it works, and the implementation details are more complex that we find on the link. A few years ago, I implemented a system for analyzing harbor porpoise vocalizations using scilab. For example, instantaneous phase(t) wants to be differentiable if instantaneous frequency is to be positive. I found that octave bands are the widest that preserves this feature. Yes, you can calculate phase(t) = atan(imaginary/real), but you have to add 2PI when atan wraps around in order to make phase monotonically increase.

Also, I am not convinced that a loudspeaker does phase modulation.  Certainly, if the speaker is linear, then superposition applies. If phase modulation occurs, it is a non-linear effect. Perhaps that effect is real, but we need to see a model to show how it comes to be. 

I might be able to help you write analytic signal  code. 

Good wishes
Gary Nelson 


Sent from my Windows 10 phone

From: Claus Futtrup
Sent: Saturday, May 26, 2018 7:41 AM
To: International users mailing list for Scilab.
Subject: Re: [Scilab-users] Simulating phase modulation

Hi Rafael

Thank you, I shall print and study. :-)

Cheers,
Claus

On Fri, May 25, 2018 at 8:42 PM Rafael Guerra <jrafaelbguerra at hotmail.com> wrote:
Hi Claus,
?
I am not aware of such function. However, you can find simple code here below for both phase modulation and demodulation, which is straightforward to translate in Scilab:
https://www.gaussianwaves.com/2017/06/phase-demodulation-using-hilbert-transform-application-of-analytic-signal/
?
Note that the phase modulation is coded differently from you snippet below.
?
Regards,
Rafael
?
From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Claus Futtrup
Sent: Friday, May 25, 2018 7:17 PM
To: International users mailing list for Scilab. <users at lists.scilab.org>
Subject: [Scilab-users] Simulating phase modulation
?
Hi there
?
In a loudspeaker the driver can move several millimeter in an attempt to reproduce a low-frequency?note. If the speaker also at the same time produce a higher tone, this second tone is phase modulated by the first one. This is a distortion of the original signal which I'd like to simulate / illustrate with some simple Scilab code, if possible.
?
In Matlab this can be simulated with pmmod.
https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/
?
Is there a similar function in Scilab? (name - please ?)
?
Here's the code I have written so far - this is the part that shows the input signal (the un-distorted signal):
?
sample_rate=20000;
t = 0:1/sample_rate:0.6;
N=size(t,'*'); //number of samples
y1 = sin(2*%pi*50*t);
y2 = 0.5*sin(2*%pi*500*t);
// y2 = 0.5*sin(2*%pi*500*t+%pi/4);
s=y1+y2+grand(1,N,'nor',0,1);
?
// Plot time-domain
endplot = round(N/15);
twoplots = scf(); // Set Current Figure (Graphics Window)
subplot(211);
plot(t(1:endplot),y1(1:endplot),t(1:endplot),y2(1:endplot));
subplot(212);
plot(t(1:endplot),y1(1:endplot)+y2(1:endplot));
?
y=fft(s);
ymax = max(abs(y));
y = y ./ ymax; // Normalize
?
// s is real so the fft response is conjugate symmetric
// and we retain only the first N/2 points
f=sample_rate*(0:(N/2))/N; //associated frequency vector
n=size(f,'*');
fftplots = scf();
plot(f(2:$),abs(y(2:n))); // drop first datapoint, f = 0 (it prevents log-plot)
a = gca();
a.log_flags = "lnn";
?
Best regards,
Claus
_______________________________________________
users mailing list
users at lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180526/e46d898f/attachment.htm>

From j.s.strom at hslmg.de  Sun May 27 10:52:32 2018
From: j.s.strom at hslmg.de (Jens Simon Strom)
Date: Sun, 27 May 2018 10:52:32 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr>
Message-ID: <5B0A71D0.3000304@hslmg.de>

Hi,
The `!` usage to recall a previous command is discussed here in detail. 
I don't miss it.  My question was different.

How about implementing again the "!" as a valid character in variable 
names as documented? How long is the reasonable response time for fixing 
the reported bug?

Regards Jens

--------------------------------------------------------------




Am 26.05.2018 14:10, schrieb Samuel Gougeon:
> Le 25/05/2018 ? 11:02, Samuel Gougeon a ?crit :
>> Hello,
>>
>> Le 25/05/2018 ? 10:45, Cl?ment David a ?crit :
>>> Hi,
>>>
>>> Thanks for your experimental work on that topic. I guess we could 
>>> merge the parser fix easily
>>> however the `!` usage to recall a previous command will remain 
>>> problematic. This usage does not seem
>>> to be documented somewhere in the help page, maybe we could discard 
>>> it. What's your opinion on that
>>> ?
>>
>> It is an historical and standard feature, used not only in Scilab 
>> language.
>> Please see the entry http://bugzilla.scilab.org/10776
>>
>> The only issue is that it is still undocumented. It could be 
>> described in the Console's page.
>
> "The only issue" : hmm, not really. Even if the missing documentation 
> is the main one, IMO this (already) very handy feature could be 
> improved in at least two ways:
>
>   * when no match is found in the history, nothing happens. For people
>     who does not know the feature (or who has forgotten that it is
>     available), this could be easily puzzling and felt as a bug.
>     => when no match is found, Scilab shall print a message "No match
>     found in history"
>
>   * A listbox listing matches could be displayed, to select the only
>     entry searched for.
>     This might be done in two ways :
>
>       o either typing "!" + one char as the very first characters
>         would enter a completion mode based on the history, with
>         matches listed in a listbox as described hereabove
>
>       o or this could be done only after validating with <enter>, or
>         CTRL+space as for the feature completion (but here targeted to
>         the history with the leading "!")
>
> If a completion mode is implemented, then the message "no match found 
> in history" could be displayed in the listbox instead of in the 
> console (note that the general feature completion triggered with 
> CTRL+space displays nothing when no match is found).
>
> By the way, as for this "!" history recall, the existing interactive 
> usage of the completion is not documented in the Console page. Let's 
> add both in the page, as they presently are.
>
> Regards
> Samuel
>
>
>
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/b170bf39/attachment.htm>

From sgougeon at free.fr  Sun May 27 12:20:21 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sun, 27 May 2018 12:20:21 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <5B0A71D0.3000304@hslmg.de>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
Message-ID: <12893834-717e-9958-3bd3-19805729acf8@free.fr>

Le 27/05/2018 ? 10:52, Jens Simon Strom a ?crit :
> Hi,
> The `!` usage to recall a previous command is discussed here in 
> detail. I don't miss it.  My question was different.

I was answering to the "What's your opinion?" question by Cl?ment.

Yes, i definitely agree on that both questions are completely 
independent, and that since the fix for variables names is deemed 
/easy/, it should be applied ASAP and fastly made available in the 
continuous built.

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/bc507cf6/attachment.htm>

From wozai at protonmail.ch  Sun May 27 13:45:42 2018
From: wozai at protonmail.ch (Wozai)
Date: Sun, 27 May 2018 04:45:42 -0700 (MST)
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <12893834-717e-9958-3bd3-19805729acf8@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
 <12893834-717e-9958-3bd3-19805729acf8@free.fr>
Message-ID: <1527421542771-0.post@n3.nabble.com>

Hi,

In theory they are "completely independent", but currently as a matter of
fact there is a clash wrt. to the
semantics of the "!" character.

Does is it really make sense allowing "!" in identifier names, while keeping
the the current history
recall feature, and assuming, that the user prepends a space in front of
expressions like "!A=1" to
get them executed in the Java console?

My conclusion: if both features are going to be kept, then the clash wrt "!"
should be resolved ASAP,
anyhow.

Cheers,

Wozai




--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From cfuttrup at gmail.com  Sun May 27 15:15:58 2018
From: cfuttrup at gmail.com (Claus Futtrup)
Date: Sun, 27 May 2018 15:15:58 +0200
Subject: [Scilab-users] Simulating phase modulation
In-Reply-To: <5b09c999.1c69fb81.23238.a038@mx.google.com>
References: <CA+Q5PwtAYBv9zeq07O6x9VkULVHCczP9G-jXPdQ1tb6iDskkvA@mail.gmail.com>
 <VI1PR04MB4109CEA1679F0CE45050B684CC690@VI1PR04MB4109.eurprd04.prod.outlook.com>
 <CA+Q5PwtdVW2+zeh8gwhMu47nRhXG5P9NU2q+2i1_mtDXi2yJMw@mail.gmail.com>
 <5b09c999.1c69fb81.23238.a038@mx.google.com>
Message-ID: <CA+Q5PwstedkKqwM2J3aakA_bTzfLaR4snDyU76SY_3xepcjXdg@mail.gmail.com>

Hi Gary

Thanks for your input. I went the analytic route for now, and simplified a
bit (as you often do with analytic solutions). You're right that from an
electro-mechanical point of view, superposition applies, but as the waves
transmit in acoustics in air - not, I'm afraid, this is where we see the
distortion, simply because the speaker cone moves several millimeter.

Best regards,
Claus

On Sat, May 26, 2018 at 10:55 PM Gary Nelson <porpoiseseeker at gmail.com>
wrote:

> Claus,
>
> Back in the late 60s, I did my PhD research using analytic signal. Yes,it
> works, and the implementation details are more complex that we find on the
> link. A few years ago, I implemented a system for analyzing harbor porpoise
> vocalizations using scilab. For example, instantaneous phase(t) wants to be
> differentiable if instantaneous frequency is to be positive. I found that
> octave bands are the widest that preserves this feature. Yes, you can
> calculate phase(t) = atan(imaginary/real), but you have to add 2PI when
> atan wraps around in order to make phase monotonically increase.
>
>
>
> Also, I am not convinced that a loudspeaker does phase modulation.
> Certainly, if the speaker is linear, then superposition applies. If phase
> modulation occurs, it is a non-linear effect. Perhaps that effect is real,
> but we need to see a model to show how it comes to be.
>
>
>
> I might be able to help you write analytic signal  code.
>
>
>
> Good wishes
>
> Gary Nelson
>
>
>
>
>
> Sent from my Windows 10 phone
>
>
>
> *From: *Claus Futtrup <cfuttrup at gmail.com>
> *Sent: *Saturday, May 26, 2018 7:41 AM
> *To: *International users mailing list for Scilab.
> <users at lists.scilab.org>
> *Subject: *Re: [Scilab-users] Simulating phase modulation
>
>
>
> Hi Rafael
>
>
>
> Thank you, I shall print and study. :-)
>
>
>
> Cheers,
>
> Claus
>
>
>
> On Fri, May 25, 2018 at 8:42 PM Rafael Guerra <jrafaelbguerra at hotmail.com>
> wrote:
>
> Hi Claus,
>
>
>
> I am not aware of such function. However, you can find simple code here
> below for both phase modulation and demodulation, which is straightforward
> to translate in Scilab:
>
>
> https://www.gaussianwaves.com/2017/06/phase-demodulation-using-hilbert-transform-application-of-analytic-signal/
>
>
>
> Note that the phase modulation is coded differently from you snippet below.
>
>
>
> Regards,
>
> Rafael
>
>
>
> *From:* users [mailto:users-bounces at lists.scilab.org] *On Behalf Of *Claus
> Futtrup
> *Sent:* Friday, May 25, 2018 7:17 PM
> *To:* International users mailing list for Scilab. <users at lists.scilab.org
> >
> *Subject:* [Scilab-users] Simulating phase modulation
>
>
>
> Hi there
>
>
>
> In a loudspeaker the driver can move several millimeter in an attempt to
> reproduce a low-frequency note. If the speaker also at the same time
> produce a higher tone, this second tone is phase modulated by the first
> one. This is a distortion of the original signal which I'd like to simulate
> / illustrate with some simple Scilab code, if possible.
>
>
>
> In Matlab this can be simulated with pmmod.
>
> https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/
>
>
>
> Is there a similar function in Scilab? (name - please ?)
>
>
>
> Here's the code I have written so far - this is the part that shows the
> input signal (the un-distorted signal):
>
>
>
> sample_rate=20000;
>
> t = 0:1/sample_rate:0.6;
>
> N=size(t,'*'); *//number of samples*
>
> y1 = sin(2*%pi*50*t);
>
> y2 = 0.5*sin(2*%pi*500*t);
>
> *// y2 = 0.5*sin(2*%pi*500*t+%pi/4);*
>
> s=y1+y2+grand(1,N,'nor',0,1);
>
>
>
> *// Plot time-domain*
>
> endplot = round(N/15);
>
> twoplots = scf(); *// Set Current Figure (Graphics Window)*
>
> subplot(211);
>
> plot(t(1:endplot),y1(1:endplot),t(1:endplot),y2(1:endplot));
>
> subplot(212);
>
> plot(t(1:endplot),y1(1:endplot)+y2(1:endplot));
>
>
>
> y=fft(s);
>
> ymax = max(abs(y));
>
> y = y ./ ymax; *// Normalize*
>
>
>
> *// s is real so the fft response is conjugate symmetric*
>
> *// and we retain only the first N/2 points*
>
> f=sample_rate*(0:(N/2))/N; *//associated frequency vector*
>
> n=size(f,'*');
>
> fftplots = scf();
>
> plot(f(2:$),abs(y(2:n))); *// drop first datapoint, f = 0 (it prevents log-plot)*
>
> a = gca();
>
> a.log_flags = "lnn";
>
>
>
> Best regards,
>
> Claus
>
> _______________________________________________
> 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: <https://lists.scilab.org/pipermail/users/attachments/20180527/4d715dee/attachment.htm>

From sgougeon at free.fr  Sun May 27 15:47:29 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sun, 27 May 2018 15:47:29 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <1527421542771-0.post@n3.nabble.com>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
 <12893834-717e-9958-3bd3-19805729acf8@free.fr>
 <1527421542771-0.post@n3.nabble.com>
Message-ID: <83c7f360-afc8-eedf-197f-b2a0b1320a43@free.fr>

Le 27/05/2018 ? 13:45, Wozai a ?crit :
> Hi,
>
> In theory they are "completely independent", but currently as a matter of
> fact there is a clash wrt. to the semantics of the "!" character.

Wozai, i am afraid i don't catch your point.
Your point is not specifically about Scilab 6, is it?
You already disagreed with Scilab 5 behavior about using "!" to recall 
the history, while "!" was actually accepted in identifiers, didn't you?

If so, then, your point is clearly distinct from Jens (and any user 
using variable names including "!") urgent claim,
that is specific to a Scilab 5 => 6 change.

If "!" is reintroduced among allowed characters in identifiers -- 
including as the first one --,
and "!" is not changed into ":" to recall the history, do we get with 
the fixed Scilab 6 anything else than the Scilab < 6 behavior?

Samuel



From wozai at protonmail.ch  Sun May 27 16:30:09 2018
From: wozai at protonmail.ch (Wozai)
Date: Sun, 27 May 2018 07:30:09 -0700 (MST)
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <83c7f360-afc8-eedf-197f-b2a0b1320a43@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
 <12893834-717e-9958-3bd3-19805729acf8@free.fr>
 <1527421542771-0.post@n3.nabble.com>
 <83c7f360-afc8-eedf-197f-b2a0b1320a43@free.fr>
Message-ID: <1527431409578-0.post@n3.nabble.com>

Dear Samuel,

my point is, it wasn't possible to execute "!A=1" in Scilab 5.x (without
prepending
a space) and IMHO one should really *improve* the situation in Scilab 6,
rather than
just restoring the status quo.

Cheers,

Wozai




--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html


From sgougeon at free.fr  Sun May 27 17:14:18 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sun, 27 May 2018 17:14:18 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <1527431409578-0.post@n3.nabble.com>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
 <12893834-717e-9958-3bd3-19805729acf8@free.fr>
 <1527421542771-0.post@n3.nabble.com>
 <83c7f360-afc8-eedf-197f-b2a0b1320a43@free.fr>
 <1527431409578-0.post@n3.nabble.com>
Message-ID: <8a92bc2d-200e-069b-e403-2758e152b4a5@free.fr>

Le 27/05/2018 ? 16:30, Wozai a ?crit :
> Dear Samuel,
>
> my point is, it wasn't possible to execute "!A=1" in Scilab 5.x (without prepending
> a space) and IMHO one should really *improve* the situation in Scilab 6,
> rather than just restoring the status quo.

Thank you Wozai for having confirmed my reading of your messages.

What Jens and i see is that Scilab 6 has not kept a status quo: by 
removing "!" from allowed characters, and by not fixing this, Scilab 6 
is downgraded now for almost 3 years it is released. Unfortunately, this 
is not a status quo.

I definitively agree that improving Scilab is necessary. This is why i 
am involved in it. But i definitely disagree about the fact that what is 
considered as an improvement or a downgrade is the decision of one 
person or a very few people just declaring that "This is better than 
existing features. It must be changed!" (with a frequent trailing 
exclamation mark).

Other users than you and me, more numerous than you alone, could easily 
consider that keeping the usage of "!" as the history call is 
preferable, because they use to use it in other environments, and then 
would consider this change as a *regression*, not an improvement.

You seem to consider that using very occasionally a leading space to 
address a variable is really a big deal. But how many coders name 
preferably their variables with a leading "!" such that having to 
prepend a space actually becomes a big deal? I am sorry, but in all the 
external codes i had to work on for 20 years, variables with a leading 
"!" were rare enough that i don't remember having met a single one.

The usual way to propose an improvement is to post a wish on Bugzilla. 
Yours is welcome!

I let Clement pushing the fix for the parser, and i hope it will be 
straightforwardly merged.

I will wait 2-3 weeks before documenting the history recall and 
CTRL+Space completion
features in the Console page. If within this time it is decided to 
change for ":", the page won't
have to be possibly updated afterwards.

Samuel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/f49f0caa/attachment.htm>

From jrafaelbguerra at hotmail.com  Sun May 27 17:56:24 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Sun, 27 May 2018 15:56:24 +0000
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <8a92bc2d-200e-069b-e403-2758e152b4a5@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
 <12893834-717e-9958-3bd3-19805729acf8@free.fr>
 <1527421542771-0.post@n3.nabble.com>
 <83c7f360-afc8-eedf-197f-b2a0b1320a43@free.fr>
 <1527431409578-0.post@n3.nabble.com>
 <8a92bc2d-200e-069b-e403-2758e152b4a5@free.fr>
Message-ID: <VI1PR04MB4109FAB4DDF988D38AB95BFCCC6F0@VI1PR04MB4109.eurprd04.prod.outlook.com>

Another two cents.
Allowing a leading exclamation mark in variable names seems to be, IMHO, a bad idea because (1) most computer languages do not seem to allow it and (2) some when encountering it, will interpret as logical negation.

Regards,
Rafael

From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Gougeon
Sent: Sunday, May 27, 2018 5:14 PM
To: Users mailing list for Scilab <users at lists.scilab.org>
Subject: Re: [Scilab-users] Exclamation mark in variable names

Le 27/05/2018 ? 16:30, Wozai a ?crit :

Dear Samuel,



my point is, it wasn't possible to execute "!A=1" in Scilab 5.x (without prepending

a space) and IMHO one should really *improve* the situation in Scilab 6,

rather than just restoring the status quo.

Thank you Wozai for having confirmed my reading of your messages.

What Jens and i see is that Scilab 6 has not kept a status quo: by removing "!" from allowed characters, and by not fixing this, Scilab 6 is downgraded now for almost 3 years it is released. Unfortunately, this is not a status quo.

I definitively agree that improving Scilab is necessary. This is why i am involved in it. But i definitely disagree about the fact that what is considered as an improvement or a downgrade is the decision of one person or a very few people just declaring that "This is better than existing features. It must be changed!" (with a frequent trailing exclamation mark).

Other users than you and me, more numerous than you alone, could easily consider that keeping the usage of "!" as the history call is preferable, because they use to use it in other environments, and then would consider this change as a regression, not an improvement.

You seem to consider that using very occasionally a leading space to address a variable is really a big deal. But how many coders name preferably their variables with a leading "!" such that having to prepend a space actually becomes a big deal? I am sorry, but in all the external codes i had to work on for 20 years, variables with a leading "!" were rare enough that i don't remember having met a single one.

The usual way to propose an improvement is to post a wish on Bugzilla. Yours is welcome!

I let Clement pushing the fix for the parser, and i hope it will be straightforwardly merged.

I will wait 2-3 weeks before documenting the history recall and CTRL+Space completion
features in the Console page. If within this time it is decided to change for ":", the page won't
have to be possibly updated afterwards.

Samuel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/ccd347ed/attachment.htm>

From sgougeon at free.fr  Sun May 27 19:23:39 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sun, 27 May 2018 19:23:39 +0200
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr>
Message-ID: <cc45b486-3d33-f6f5-95a8-849dad0cb73d@free.fr>

Le 26/05/2018 ? 14:10, Samuel Gougeon a ?crit :
> Le 25/05/2018 ? 11:02, Samuel Gougeon a ?crit :
>> Hello,
>>
>> Le 25/05/2018 ? 10:45, Cl?ment David a ?crit :
>>> Hi,
>>>
>>> Thanks for your experimental work on that topic. I guess we could 
>>> merge the parser fix easily
>>> however the `!` usage to recall a previous command will remain 
>>> problematic. This usage does not seem
>>> to be documented somewhere in the help page, maybe we could discard 
>>> it. What's your opinion on that
>>> ?
>>
>> It is an historical and standard feature, used not only in Scilab 
>> language.
>> Please see the entry http://bugzilla.scilab.org/10776
>>
>> The only issue is that it is still undocumented. It could be 
>> described in the Console's page.
>
> "The only issue" : hmm, not really. Even if the missing documentation 
> is the main one, IMO this (already) very handy feature could be 
> improved in at least two ways:
>
>   * when no match is found in the history, nothing happens. For people
>     who does not know the feature (or who has forgotten that it is
>     available), this could be easily puzzling and felt as a bug.
>     => when no match is found, Scilab shall print a message "No match
>     found in history"
>
>   * A listbox listing matches could be displayed, to select the only
>     entry searched for.
>     This might be done in two ways :
>
>       o either typing "!" + one char as the very first characters
>         would enter a completion mode based on the history, with
>         matches listed in a listbox as described hereabove
>
>       o or this could be done only after validating with <enter>, or
>         CTRL+space as for the feature completion (but here targeted to
>         the history with the leading "!")
>
> If a completion mode is implemented, then the message "no match found 
> in history" could be displayed in the listbox instead of in the 
> console (note that the general feature completion triggered with 
> CTRL+space displays nothing when no match is found).

This request for improvements is now reported here: 
http://bugzilla.scilab.org/15592

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/1505163c/attachment.htm>

From sgougeon at free.fr  Sun May 27 19:51:36 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Sun, 27 May 2018 19:51:36 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <pe7a7j$et5$1@blaine.gmane.org>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
Message-ID: <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>

Hello Philippe,

Le 24/05/2018 ? 23:20, philippe a ?crit :
> Le 21/05/2018 ? 15:46, Samuel Gougeon a ?crit :
>> There are at least two way to do it :
>>
>>    * either keep *cond([])* to *1* and set all *cond([], p)* to 1 instead
>>      of 0
>>    * or set *cond([])* to *0*.
>>
>> I don't see any clear reason enforcing a choice rather than the other.
>> Do you?
> from a mathematical point of view  the Condition number of a matrix A is
> defined by
>
> cond(A)=||A|| . ||A^(-1)||
> its interest is that when solving A*x=y rounding errors on y (eps) are
> amplified to be cond(A)*eps. The optimal value of cond(A) is 1 (for
> identity matrix) so for me it looks natural that cond([])=1 .

I must confess that, unless claiming that *norm([])***is NOT 0,  i do 
not catch clearly the logical of your conclusion

  * /"The optimal value of cond(A) is 1 (for identity matrix) so for me
    it looks natural that cond([])=1"/

Indeed, cond([]) = norm([])*norm(inv([])) = 0*norm([]) = 0*0 = 0

By the way, according to the clear explanation you give about the 
meaning of the condition number, the value 0 is even more optimal: it 
says that small variations are not getting amplified, but killed.

Please let us know more about your proposal for 1 instead of 0.

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/401b0ce7/attachment.htm>

From rouxph.22 at gmail.com  Sun May 27 23:23:47 2018
From: rouxph.22 at gmail.com (philippe)
Date: Sun, 27 May 2018 23:23:47 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
 <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
Message-ID: <pef7h0$prd$1@blaine.gmane.org>

Le 27/05/2018 ? 19:51, Samuel Gougeon a ?crit?:
> [...]
> 
> Please let us know more about your proposal for 1 instead of 0.

I can't say it now, the best value also depends on the compatibility of
this value with others assigned to matrix functions with argument [].  I
would like to check it before, but with the example you mention 0 looks
a better choice.

Philippe



From jrafaelbguerra at hotmail.com  Sun May 27 23:49:34 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Sun, 27 May 2018 21:49:34 +0000
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <pef7h0$prd$1@blaine.gmane.org>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
 <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
 <pef7h0$prd$1@blaine.gmane.org>
Message-ID: <VI1PR04MB4109140FFC4173DD40B969ADCC6F0@VI1PR04MB4109.eurprd04.prod.outlook.com>

You guys may have missed the paper sent before on the algebra of empty matrices by Carl de Boor (Emeritus Professor in Mathematics and Computer Science), so I am sending it again:

https://pdfs.semanticscholar.org/0f3b/c36f19d5c6a761c19fbc3c4ebde2f31b0a10.pdf



He argues that the condition number of the square empty matrix [] = 0 and its det([]) = 1



Regards,

Rafael



-----Original Message-----
From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of philippe
Sent: Sunday, May 27, 2018 11:24 PM
To: users at lists.scilab.org
Subject: Re: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)



Le 27/05/2018 ? 19:51, Samuel Gougeon a ?crit :

> [...]

>

> Please let us know more about your proposal for 1 instead of 0.



I can't say it now, the best value also depends on the compatibility of this value with others assigned to matrix functions with argument [].  I would like to check it before, but with the example you mention 0 looks a better choice.



Philippe








-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180527/6ee33da6/attachment.htm>

From sgougeon at free.fr  Mon May 28 00:09:10 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Mon, 28 May 2018 00:09:10 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <VI1PR04MB4109140FFC4173DD40B969ADCC6F0@VI1PR04MB4109.eurprd04.prod.outlook.com>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
 <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
 <pef7h0$prd$1@blaine.gmane.org>
 <VI1PR04MB4109140FFC4173DD40B969ADCC6F0@VI1PR04MB4109.eurprd04.prod.outlook.com>
Message-ID: <285c8bd8-1275-4afd-04a3-538d6f8d643d@free.fr>

Le 27/05/2018 ? 23:49, Rafael Guerra a ?crit :
>
> You guys may have missed the paper sent before on the algebra of empty 
> matrices by Carl de Boor (Emeritus Professor in Mathematics and 
> Computer Science), so I am sending it again:
>
> https://pdfs.semanticscholar.org/0f3b/c36f19d5c6a761c19fbc3c4ebde2f31b0a10.pdf
>
> He argues that the condition number of the square empty matrix [] = 0 
> and its det([]) = 1
>

No no, i did not miss it. His argument is that norm([]) is 0:

/Any norm of an empty matrix is zero, as the supremum of the empty set 
of nonnegative numbers. This implies that the condition number of the 
square empty matrix [] is 0./

This does not prevent asking to Philippe why he proposes 1, noticeably 
if he also considers that norm([])==0, or not.

The De Boor argument "/as the supremum of the empty set of nonnegative 
numbers."/ is not clear to me./

/Is it clear for you? Would you have a second convergent independent 
reference?

Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180528/1046b123/attachment.htm>

From jrafaelbguerra at hotmail.com  Mon May 28 10:58:09 2018
From: jrafaelbguerra at hotmail.com (Rafael Guerra)
Date: Mon, 28 May 2018 08:58:09 +0000
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <285c8bd8-1275-4afd-04a3-538d6f8d643d@free.fr>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
 <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
 <pef7h0$prd$1@blaine.gmane.org>
 <VI1PR04MB4109140FFC4173DD40B969ADCC6F0@VI1PR04MB4109.eurprd04.prod.outlook.com>
 <285c8bd8-1275-4afd-04a3-538d6f8d643d@free.fr>
Message-ID: <VI1PR04MB4109565E527F915B95D56B3FCC6E0@VI1PR04MB4109.eurprd04.prod.outlook.com>

No, it is not clear but then, I am not an Emeritus Professor of Mathematics.
Sorry but I have no other references on this empty subject.

Regards,
Rafael

From: users [mailto:users-bounces at lists.scilab.org] On Behalf Of Samuel Gougeon
Sent: Monday, May 28, 2018 12:09 AM
To: Users mailing list for Scilab <users at lists.scilab.org>
Subject: Re: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)

Le 27/05/2018 ? 23:49, Rafael Guerra a ?crit :

You guys may have missed the paper sent before on the algebra of empty matrices by Carl de Boor (Emeritus Professor in Mathematics and Computer Science), so I am sending it again:

https://pdfs.semanticscholar.org/0f3b/c36f19d5c6a761c19fbc3c4ebde2f31b0a10.pdf



He argues that the condition number of the square empty matrix [] = 0 and its det([]) = 1

No no, i did not miss it. His argument is that norm([]) is 0:

Any norm of an empty matrix is zero, as the supremum of the empty set of nonnegative numbers. This implies that the condition number of the square empty matrix [] is 0.

This does not prevent asking to Philippe why he proposes 1, noticeably if he also considers that norm([])==0, or not.

The De Boor argument "as the supremum of the empty set of nonnegative numbers." is not clear to me.

Is it clear for you? Would you have a second convergent independent reference?

Samuel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180528/28c5df2e/attachment.htm>

From stephane.mottelet at utc.fr  Mon May 28 11:26:47 2018
From: stephane.mottelet at utc.fr (=?UTF-8?Q?St=c3=a9phane_Mottelet?=)
Date: Mon, 28 May 2018 11:26:47 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <VI1PR04MB4109565E527F915B95D56B3FCC6E0@VI1PR04MB4109.eurprd04.prod.outlook.com>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
 <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
 <pef7h0$prd$1@blaine.gmane.org>
 <VI1PR04MB4109140FFC4173DD40B969ADCC6F0@VI1PR04MB4109.eurprd04.prod.outlook.com>
 <285c8bd8-1275-4afd-04a3-538d6f8d643d@free.fr>
 <VI1PR04MB4109565E527F915B95D56B3FCC6E0@VI1PR04MB4109.eurprd04.prod.outlook.com>
Message-ID: <b9603fd0-3ff9-950d-3efe-df32cc244e70@utc.fr>

Le 28/05/2018 ? 10:58, Rafael Guerra a ?crit?:
>
> No, it is not clear but then, I am not an Emeritus Professor of 
> Mathematics.
>
> Sorry but I have no other references on this empty subject.
>
> Regards,
>
> Rafael
>
> *From:*users [mailto:users-bounces at lists.scilab.org] *On Behalf Of 
> *Samuel Gougeon
> *Sent:* Monday, May 28, 2018 12:09 AM
> *To:* Users mailing list for Scilab <users at lists.scilab.org>
> *Subject:* Re: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
>
> Le 27/05/2018 ? 23:49, Rafael Guerra a ?crit?:
>
>     You guys may have missed the paper sent before on the algebra of
>     empty matrices by Carl de Boor (Emeritus Professor in Mathematics
>     and Computer Science), so I am sending it again:
>
>     https://pdfs.semanticscholar.org/0f3b/c36f19d5c6a761c19fbc3c4ebde2f31b0a10.pdf
>     <https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/pdfs.semanticscholar.org/0f3b/c36f19d5c6a761c19fbc3c4ebde2f31b0a10.pdf>
>
>     He argues that the condition number of the square empty matrix []
>     = 0 and its det([]) = 1
>
>
> No no, i did not miss it. His argument is that norm([]) is 0:
>
> /Any norm of an empty matrix is zero, as the supremum of the empty set 
> of nonnegative numbers. This implies that the condition number of the 
> square empty matrix [] is 0./
>
> This does not prevent asking to Philippe why he proposes 1, noticeably 
> if he also considers that norm([])==0, or not.
>
> The De Boor argument "/as the supremum of the empty set of nonnegative 
> numbers."/ is not clear to me./
> /
>
Every property (even trivially false) is true for all elements of the 
empty set. Hence, any positive number can be taken as the supremum of 
the empty set. However, I don't see in which sense this could justify 
the discusses convention...

S.
>
> /
> /Is it clear for you? Would you have a second convergent independent 
> reference?
>
> Samuel
>
>
>
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users


-- 
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: <https://lists.scilab.org/pipermail/users/attachments/20180528/befec8f8/attachment.htm>

From Clement.David at esi-group.com  Mon May 28 12:09:35 2018
From: Clement.David at esi-group.com (=?utf-8?B?Q2zDqW1lbnQgRGF2aWQ=?=)
Date: Mon, 28 May 2018 10:09:35 +0000
Subject: [Scilab-users] Exclamation mark in variable names
In-Reply-To: <12893834-717e-9958-3bd3-19805729acf8@free.fr>
References: <5AFFFE04.4030207@hslmg.de> <1527236245681-0.post@n3.nabble.com>
 <1527237924.2237.52.camel@esi-group.com>
 <63e19c5f-9d19-0f2b-5378-8eb5aca5fcf6@free.fr>
 <7b963591-5c30-934a-c6e1-568cbb226f2a@free.fr> <5B0A71D0.3000304@hslmg.de>
 <12893834-717e-9958-3bd3-19805729acf8@free.fr>
Message-ID: <1527502174.3835.15.camel@esi-group.com>

Hello,

Le dimanche 27 mai 2018 ? 12:20 +0200, Samuel Gougeon a ?crit :
> Le 27/05/2018 ? 10:52, Jens Simon Strom a ?crit :
> > Hi,
> > The `!` usage to recall a previous command is discussed here in detail. I don't miss it.  My
> > question was different.
>  
> I was answering to the "What's your opinion?" question by Cl?ment.
> 
> Yes, i definitely agree on that both questions are completely independent, and that since the fix
> for variables names is deemed easy, it should be applied ASAP and fastly made available in the
> continuous built.

"Fixes Parser to Accept "!" in Variable Names
<https://github.com/rdbyk/balisc/pull/323>?" has been merged through
https://codereview.scilab.org/#/c/20069/ .

Thanks for your contribution !

--
Cl?ment

From sgougeon at free.fr  Wed May 30 07:16:30 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Wed, 30 May 2018 07:16:30 +0200
Subject: [Scilab-users] A way to manage and use homonymous builtin functions?
Message-ID: <4666d3d3-7991-3d8a-d1d3-b6fb2394266a@free.fr>

Hello,

Would anyone know a way to call a builtin function in a specific gateway?

I mean: for a function in Scilab language, if several homonymous 
functions -- say thisfun() -- are available in distinct available 
libraries -- say scilib and mylib --, it is possible to call either the 
scilib one with scilib.thisfun(), or the mylib one with mylib.thisfun().

Is there a similar mechanism for builtin functions?

The use case is the following: the external scicv module redefines the 
write() Scilab builtin function as a new builtin. But the scicv instance 
is not compatible with the Scilab write(). And is called first (i guess 
because it has been loaded /after/ the io gateway of Scilab's write()). 
This makes Scilab some programs failing when scicv is loaded. This is 
reported on the scicv forge/ATOMS page.

So, i am looking for a way to point to the Scilab's write() in any case.

Thanks for any hints

Regards

Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180530/e7f6f6ae/attachment.htm>

From sgougeon at free.fr  Wed May 30 19:44:56 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Wed, 30 May 2018 19:44:56 +0200
Subject: [Scilab-users] Is cond([]) 0 or 1 ? (bug 15579)
In-Reply-To: <pef7h0$prd$1@blaine.gmane.org>
References: <c4af5d46-ab50-d6e8-31a9-7ca118a695d2@free.fr>
 <pe7a7j$et5$1@blaine.gmane.org>
 <72a2c293-eb5f-2bcf-58b6-33a45d7c3324@free.fr>
 <pef7h0$prd$1@blaine.gmane.org>
Message-ID: <cb1a3011-77db-3b04-696f-9f4a67717971@free.fr>

Le 27/05/2018 ? 23:23, philippe a ?crit :
> Le 27/05/2018 ? 19:51, Samuel Gougeon a ?crit :
>> [...]
>>
>> Please let us know more about your proposal for 1 instead of 0.
> I can't say it now, the best value also depends on the compatibility of
> this value with others assigned to matrix functions with argument [].
> I would like to check it before, but with the example you mention 0 looks
> a better choice.

As far as i understand, rcond(X) should return the inverse of cond(X,1).

In Scilab, rcond([]) returns []. This looks pretty inconsistent with any 
Scilab cond([],p) result,
and does not help much for this longly empty cond([]) discussion.

In Octave, rcond([]) returns %inf, that's consistent with cond([])==0.

However, still, a guy whose det() is 1 and norm() is 0 is really a 
special guy.

If we have no evidence/reference on the fact norm([]) is 1 rather than 0,
i guess that we will change Scilab's cond([]) to 0.
May the Scilab creator having initially set 1 read us and correct us if 
it's needed.

Cheers
Samuel



From sgougeon at free.fr  Wed May 30 21:35:11 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Wed, 30 May 2018 21:35:11 +0200
Subject: [Scilab-users] surf & isoview: bug or unexpected "feature"
In-Reply-To: <5500d27f-12bc-c505-2a09-7b27ec28307c@utc.fr>
References: <b59644b5-5522-b0a5-1b4f-e3d0e7cfaea5@laas.fr>
 <9d780bce-0cc1-229d-ffa5-71357640b029@free.fr>
 <5500d27f-12bc-c505-2a09-7b27ec28307c@utc.fr>
Message-ID: <04f239bd-ec9d-7046-5ec4-a71d878e21e3@free.fr>

Le 21/02/2018 ? 13:47, St?phane Mottelet a ?crit :
> Great ! But in this case I would add also
>
> isoview('off') missed switching .cube_scaling='on'

Done. The fixed behaviour of isoview() is now available in Scilab 6.1.0- 
(by now, continuous built of the master)

Samuel



From sgougeon at free.fr  Wed May 30 22:02:57 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Wed, 30 May 2018 22:02:57 +0200
Subject: [Scilab-users] .isoview (and .cube_scaling) design <= Re: surf &
 isoview: bug or unexpected "feature"
In-Reply-To: <98256e3d-9110-e818-7456-d7dc1d601c77@utc.fr>
References: <b59644b5-5522-b0a5-1b4f-e3d0e7cfaea5@laas.fr>
 <91aa43f0-f41c-5754-83dc-6fad5050af07@laas.fr>
 <98256e3d-9110-e818-7456-d7dc1d601c77@utc.fr>
Message-ID: <27a588a5-92a7-b557-6ac8-1f2b57a51c00@free.fr>

Le 21/02/2018 ? 12:05, St?phane Mottelet a ?crit :
> Hello,
>
> The cube_scaling appeared with version 5.2 or 5.3 (I don't remember) 
> because the previous default behavior, which corresponds to
>
> cube_scaling="off"
>
> could produce surfaces which were very hard to visualize with default 
> angles and hard to homogeneously "rotate" in the elevation angle, when 
> the surface is very flat. See this example:
>
> clf
> [x,y]=meshgrid(-1:0.1:1,-2:0.1:2);
> surf(x,y,x.*y*1e-5)
> gca().cube_scaling="off";
>
> The default combination cube_scaling="on" and isoview="off" 
> corresponds to the default behavior in Matlab. Introducing the 
> cube_scaling property and its default value helps "migrating users" as 
> they obtain the same thing in both softwares.
>
> Anyway, I don't think that the combination cube_scaling="off" and 
> isoview="off" has still an interest as it can produce the annoying 
> behavior of the above example. The default combination 
> cube_scaling="on" and isoview="off" allows the 3D object to fill the 
> whole 3D box whatever the proportions of the enclosing 2D canvas, 
> which is the expected default behavior.
>
> Antoine's remark points out that cube_scaling="on" and isoview="on" 
> produces a plot whose proportions are completely wrong with respect of 
> the value of isoview="on". Adding in the help page of 
> "axes_properties/isoview " a sentence such as
>
> If you want *real* isometric scales on all axes in 3D you have to set 
> cube_scaling property to 'off'
>
> is not admissible.Setting the value of cube_scaling  to the opposite 
> of the isoview value when the latter is changed could be a solution. 
> At least, when the high level isoview *command* is used.

The implementation of this *gca().cube_scaling* property is 
unfortunately an example of the (too) raw and fast importation of 
a feature answering in a truly awkward way to a true need.

Extending some isometrical tunings in 3D by /adding/ /an axes property 
/instead of by /adding some new values to the existing .isoview property 
/IMO was a costy error. Why ?

  * adding a property whose job is not orthogonal to the /.isoview/ one
    makes the whole thing puzzling, and i guess the coding as well. In
    3D, what does this property, and what does the other one in a
    dependent way? This is really annoying.

  * Due to a missing or poor initial analysis, the result is rather
    poor. The point is that we may need isometrical modes coupling any
    possible pair of directions, or even all of them. This is why, from
    the "on"|"off" values available in 2D, new .isoview values like
    "xy", "xz", "yz", "xyz" -- without creating any new property --
    would have been more handy and powerful. And likely clearer and
    easier to code.
    Then, to actually get a cubic shape whatever are the x,y and z
    scales, a last value "cube" could be possible, still without any new
    property.

Analyzing a feature before implementing it has certainly a (low) cost. 
The use of some poorly built features wastes a lot more time for 
everyone and for a long time.
It's a pity.

Regards
Samuel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180530/9c11d1ce/attachment.htm>

From sgougeon at free.fr  Thu May 31 11:39:08 2018
From: sgougeon at free.fr (Samuel Gougeon)
Date: Thu, 31 May 2018 11:39:08 +0200
Subject: [Scilab-users] circshift() : Scilab Enhancement Proposal
In-Reply-To: <VI1PR04MB3229F14F74882E9FB67B4FF8CC100@VI1PR04MB3229.eurprd04.prod.outlook.com>
References: <1b44e4b4-a924-06f2-b4ee-c1c416008947@free.fr>
 <e48726de-36fc-8c15-22f4-dc01ea65dded@laas.fr>
 <VI1PR04MB3229E6B86A7EFD519D276A97CC1E0@VI1PR04MB3229.eurprd04.prod.outlook.com>
 <8f751180-20e9-3c28-83d2-f6801271caf6@free.fr>
 <VI1PR04MB3229F14F74882E9FB67B4FF8CC100@VI1PR04MB3229.eurprd04.prod.outlook.com>
Message-ID: <d8d18c74-9823-b7d7-c5e8-57a2911f9931@free.fr>

Hello Rafael,

Le 10/01/2018 ? 00:01, Rafael Guerra a ?crit :
>
> Hi Samuel,
>
> .../...
>
> Regarding extending circshift() to non-integers I had in mind 
> n-dimensional matrices and the discrete Fourier transform.
>
> These matrices could represent regularly sampled 1D time signals, or 
> 2D/3D spatial functions, etc.
>
> In this context the shift theorem gives a specific meaning to the 
> circular shifts of the input signals (corresponding to linear phase 
> shifts in their spectrum).
>

Using fractional shifts means that we somewhat address some frequencies
higher than the sampling one. Then, as expected, if such a feature is
implemented through the DFT, and unless the input signal is
* a truly periodic one
* sampled over an integer number of its period,
which are very restrictive conditions, we will face issues such as the Gibb
artefact, deserving or even demanding some damping window options, etc.

This is far beyond the main purpose of circshift() whose elementary job 
misses for years now.
We would go then into signal processing, with its usual complications,
only applicable to arrays of real or complex numbers, while circshift() is
basically applicable to all arrays of any datatype.

To me, your proposal to extend circshift() to fractional shifts needs some
further analysis, to compare several possible implementations.
It could make an additional SEP with its own new options, after merging
the basic but very expected circshift().

Using the DFT to interpolate input data has the main advantage to not be
a local method. As a trade-off, it yields some artefacts, at least with 
a simple
implementation as provided herebelow.
Please do not hesitate to go on with it. I have not tried to recover the 
undamped
shifted signal after using some damping window to get ride of Gibb.

Using classical local interpolations -- linear or splined ones -- is 
also meaningful
and applicable only to real or complex numerical data.

A benchmark could be set to compare performances of DFT versus local 
interpolations,
iterating small fractional shifts summed to an integer, and measuring 
the distance
of the result to the true expected one (after an integer shift).
I let you do that and build your extension proposal.

Cheers
Samuel


// Fractional circshift: tests
n  =  100;
sig  =  sin(linspace(0,  2.5,  n));  // Built discontinuity from 2.5=>0
clf
subplot(1,  3,  1)
plot(sig)
title("Input 1D signal",  "fontsize",3);

subplot(1,  3,  2)
v  =  frac_shift(sig,  35);
plot(real(v))
title("Integer shift = 35",  "fontsize",3);

subplot(1,  3,  3)
v  =  frac_shift(sig,  35.4);   // fractional shift = 35.4
plot(real(v))
title("Fractional shift = 35.4",  "fontsize",3);
gca().tight_limits="on";
// As expected, we get some Gibb artefact oscillations near the
// discontinuity. Avoiding them required to damp wings of the input
// signal in order to damp the discontinuity, requiring to use a
// damping window and profile.

function  v= frac_shift(x, n)
     // x: 1D signal
     // n: shift, in decimal number of samples
     // v: x shifted by i samples
     N  =  length(x)
     sp  =  fft(x,  -1).';
     frequencies  =  (floor(-(N-1)/2):floor((N-1)/2))/N
     c  =  exp(2*%pi*%i*n*frequencies).';
     shifted_sp  =  ifftshift(fftshift(sp).*c);
     v  =  fft(shifted_sp,1);
     mprintf("%5.2f max(abs(imag(shifted sig)))= %f\n",  n,  max(abs(imag(v))));
endfunction

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180531/3ee75b8b/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: lhameelmlamodjop.png
Type: image/png
Size: 10079 bytes
Desc: not available
URL: <https://lists.scilab.org/pipermail/users/attachments/20180531/3ee75b8b/attachment.png>