Re: [Scilab-users] I need a numerial procedure to generate Exponentially Modified Gaussian deviates

[Heinz Nabielek]

Very useful. I had done a similar thing with the EMG, but much more clumsy
Heinz

> On 29.04.2024, at 16:41, Federico Miyara  wrote:
>
>
> Heinz,
>
> I don't know if this might be useful. The function I'm attaching allows to 
> generate random numbers according to any distribution, either empirical or 
> analytical. The results may be scaled if necessary.
>
> Regards,
>
> Federico Miyara
>
>
> On 29/4/2024 00:11, Heinz Nabielek wrote:
>> Colleagues:
>>
>> bird flight altitude probabilities are given by the exponentially modified 
>> Gaussian distribution EMG f=f(h), in my case
>>
>> f = (a/2)*exp((a/2)*((a*σ^2)-2*h)) .* erfc((a*σ^2-h)/σ/sqrt(2))
>>
>> with h=hmeasured-85 in meters, a=1/60m, σ=30m. See: 
>> <https://en.wikipedia.org/wiki/Exponentially_modified_Gaussian_distribution>.
>>
>> For bird flight Monte-Carlo simulations, I need random deviates representing 
>> this distribution.
>> Normally, I generate the cumulative function F of f. Then, I project uniform 
>> random numbers U(0,1) on F^-1, the inverse of the cumulative.
>>
>> With the EMG, closed solutions are not available and I need a numerial 
>> procedure to generate EMG deviates.
>> What is the best way to do this?
>> Best greetings
>> Heinz
>> This email and any attachments are intended solely for the use of the 
>> individual or entity to whom it is addressed and may be confidential and/or 
>> privileged.
>>
>> If you are not one of the named recipients or have received this email in 
>> error,
>>
>> (i) you should not read, disclose, or copy it,
>>
>> (ii) please notify sender of your receipt by reply email and delete this 
>> email and all attachments,
>>
>> (iii) Dassault Systèmes does not accept or assume any liability or 
>> responsibility for any use of or reliance on this email.
>>
>>
>> Please be informed that your personal data are processed according to our 
>> data privacy policy as described on our website. Should you have any 
>> questions related to personal data protection, please contact 3DS Data 
>> Protection Officer https://www.3ds.com/privacy-policy/contact/
>>
>>
>>
>>
>>
>> ___
>> users mailing list - users@lists.scilab.org
>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>> https://lists.scilab.org/mailman/listinfo/users
>
>
> --
> Este correo electrónico ha sido analizado en busca de virus por el software 
> antivirus de Avast.
> www.avast.com
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

function x = rand_weighted(w, N)
   // This function generates N random numbers in the unit 
   // interval according to an empirical statistical 
   // distribution given by vector w
   //
   // Usage:
   //x = rand_weighted(w, N)
   //
   // where
   //x: vector containing random numbers
   //w: vector specifying the statistical distribution 
   //   as the probabilities in equally-wide bins along
   //   the unit interval
   //N: number of random numbers to be created  
   //
   // NOTES: 1) Vector w is a probability function, so sum(w) 
   //   must be equal to 1 
   //2) A simple linear transformation such as 
   //   y = a + x*(b - a)
   //   may be used to get a distribution over the 
   //   interval from a to b
   //3) The number M of bins is arbitrary, but as
   //   an interpolation is performed, as M gets
   //   larger the accusacy improves. Usually M 
   //   depends on the size of the sample data used
   //   to estimate the empirical distribution 
   //
   // 
   // Author: Federico Miyara
   // Date:   2017-02-06
   // 2020-03-12
   // 2023-11-05

   M = length(w);

   // Cumulative probability including 0
   W = cumsum([0,w]);

   // The inverse fu

Re: [Scilab-users] I need a numerial procedure to generate Exponentially Modified Gaussian deviates

[Heinz Nabielek]

On 29.04.2024, at 09:24, Stéphane Mottelet  wrote:
>
> Sorry,
>
> This is just a manifestation of another occurence of the xy problem. You
> told about y (drawing from a given PDF) but the original problem was x :
> drawing a random variable which is the sum of two random variables for
> which we already have  generators. In fact an EMG random variable Z may
> be expressed as Z = X + Y, where X and Y are independent, X is Gaussian
> with mean μ and variance σ2, and Y is exponential of rate λ
> (https://en.wikipedia.org/wiki/Exponentially_modified_Gaussian_distribution).
>
> So you just have to use grand
> (https://help.scilab.org/docs/2024.0.0/en_US/grand.html) to make
> independant draws of X and Y then make the sum afterwards


Marvellous. I should have thought of that myself.
Heinz

Shall we try to repeat it in the DISCOURSE?
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] I need a numerial procedure to generate Exponentially Modified Gaussian deviates

[Heinz Nabielek]

Colleagues:

bird flight altitude probabilities are given by the exponentially modified 
Gaussian distribution EMG f=f(h), in my case

f = (a/2)*exp((a/2)*((a*σ^2)-2*h)) .* erfc((a*σ^2-h)/σ/sqrt(2))

with h=hmeasured-85 in meters, a=1/60m, σ=30m. See: 
<https://en.wikipedia.org/wiki/Exponentially_modified_Gaussian_distribution>.

For bird flight Monte-Carlo simulations, I need random deviates representing 
this distribution.
Normally, I generate the cumulative function F of f. Then, I project uniform 
random numbers U(0,1) on F^-1, the inverse of the cumulative.

With the EMG, closed solutions are not available and I need a numerial 
procedure to generate EMG deviates.
What is the best way to do this?
Best greetings
Heinz
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



Exponentially Modified Gaussian.pdf
Description: Adobe PDF document



___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

[Scilab-users] How to compute the non-central t-distribution

[Heinz Nabielek]

https://scilab.discourse.group/t/non-central-t-distribution/446?u=heinznabielek

MATLAB has the function ‘nctpdf’ function to compute the non-central 
t-distribution. What would be the SciLab equivalent? Or how would you write 
Scilab coding to arrive at the non-central t-distribution? Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] macOS Sonoma 14.4 SciLab still fails to start

[Heinz Nabielek]

macOS Sonoma 14.4 full version now, but SciLab still fails to start. What can 
we do?
Heinz


heinznabielek@AppleSilicon ~ % 
/Applications/scilab-2024.0.0.app/Contents/MacOS/scilab
Warning: Localization issue. Does not support the locale ‘’
Returned: NULL
Current system locale: C
Did you install the system locales?
Picked up _JAVA_OPTIONS: -Djava.awt.headless=false
2024-03-13 14:37:29.878 scilab-bin[16883:143280] WARNING: Secure coding is not 
enabled for restorable state! Enable secure coding by implementing 
NSApplicationDelegate.applicationSupportsSecureRestorableState: and returning 
YES.
/Applications/scilab-2024.0.0.app/Contents/MacOS/scilab: line 928: 16883 
Trace/BPT trap: 5 arch -arch arm64 “SCILABBIN""@”
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scilab crashes every time

[Heinz Nabielek]

On 04.02.2024, at 19:04, Stéphane Mottelet  wrote:
>
> Hello,
>
> Why are you using beta versions of macOS ?


Sorry, I was not even aware of it. How can I backstep to normal macOS?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] What is the simplest way to plot a diagram with two x-axes?

[Heinz Nabielek]

What is the simplest way to plot a diagram with two x-axes using using 
"newaxes()"?
Heinz

[cid:E15015B9-E43B-432E-B482-3BB22849F10F@home]

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab crashes every time

[Heinz Nabielek]

 0x19d93150c 
NSPerformVisuallyAtomicChange + 108
6   AppKit 0x19e3a7690 -[NSWindow 
_reallyDoOrderWindowOutRelativeTo:] + 448
7   AppKit 0x19e3a7a54 -[NSWindow 
_reallyDoOrderWindow:] + 80
8   AppKit 0x19e3a7ca4 -[NSWindow 
_doOrderWindow:] + 264
9   AppKit 0x19e3a30b4 -[NSWindow 
_finishClosingWindow] + 472
10  AppKit 0x19db76a0c -[NSWindow _close] + 
408
11  libnativewindow_macosx.dylib   0x109a48ac4 
Java_jogamp_nativewindow_macosx_OSXUtil_DestroyNSWindow0 + 88
12  ???0x175d118ac ???
13  ???0x175d0dfc8 ???

and so on…...

> On 30.01.2024, at 08:03, Stéphane Mottelet  wrote:
> 
> Hello,
> 
> What did you change ? OS version, startup files, toolboxes installed ? First 
> step as usual, launch Scilab from a terminal to see something more.
> 
> S.
> 
>> Le 30 janv. 2024 à 04:55, Heinz Nabielek  a écrit :
>> 
>> Suddenly, Scilab crashes every time I try to start. What happened? Equally 
>> so after reinstallation.
>> Heinz
>> 
>> Scilab 2024.0.0 - macOS 64 bits (ARM) (dmg)
>> Apple M1 macOS Sonoma 14.4 Beta
>> ___
>> users mailing list - users@lists.scilab.org
>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>> https://antiphishing.vadesecure.com/v4?f=RnlZMDJMbzlZZ1hsTmtFVkDccrhAoM7SzRmBq20IJL9IM1Dfz-BjDqrY7ltRDXXI=anNVMGFsTWdEV0k3d09ZUybl5aVhDiFd5rBCIEVNi9U=AkQR=U1U2bFJicVpYWHFiUTEzTbUwYCs7bp9qVR7fTqMpOXVr0W44Q2tUuCYKa-ZRQSGt=3777ffbb2d16b9e5b77d440234152d85ea3be33ee6846548c550f3c07fc90edd=https%3A%2F%2Flists.scilab.org%2Fmailman%2Flistinfo%2Fusers
>> This email and any attachments are intended solely for the use of the 
>> individual or entity to whom it is addressed and may be confidential and/or 
>> privileged.
>> 
>> If you are not one of the named recipients or have received this email in 
>> error,
>> 
>> (i) you should not read, disclose, or copy it,
>> 
>> (ii) please notify sender of your receipt by reply email and delete this 
>> email and all attachments,
>> 
>> (iii) Dassault Systèmes does not accept or assume any liability or 
>> responsibility for any use of or reliance on this email.
>> 
>> 
>> Please be informed that your personal data are processed according to our 
>> data privacy policy as described on our website. Should you have any 
>> questions related to personal data protection, please contact 3DS Data 
>> Protection Officer 
>> https://antiphishing.vadesecure.com/v4?f=RnlZMDJMbzlZZ1hsTmtFVkDccrhAoM7SzRmBq20IJL9IM1Dfz-BjDqrY7ltRDXXI=anNVMGFsTWdEV0k3d09ZUybl5aVhDiFd5rBCIEVNi9U=AkQR=U1U2bFJicVpYWHFiUTEzTbUwYCs7bp9qVR7fTqMpOXVr0W44Q2tUuCYKa-ZRQSGt=88aa59a4ef9dc5e12ffc7e7833ddead939352f23c13d533ebab1478db8fc9eed=https%3A%2F%2Fwww.3ds.com%2Fprivacy-policy%2Fcontact%2F
>> 
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
> 
> If you are not one of the named recipients or have received this email in 
> error,
> 
> (i) you should not read, disclose, or copy it,
> 
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
> 
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
> 
> 
> Please be informed that your personal data are processed according to our 
> data privacy policy as described on our website. Should you have any 
> questions related to personal data protection, please contact 3DS Data 
> Protection Officer https://www.3ds.com/privacy-policy/contact/
> 

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab crashes every time

[Heinz Nabielek]

There was a recent minor OS Update. Νοw
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



Error messages:
/Applications/scilab-2024.0.0.app/Contents/bin/scilab
Warning: Localization issue. Does not support the locale ''
Returned: NULL
Current system locale: C
Did you install the system locales?
Picked up _JAVA_OPTIONS: -Djava.awt.headless=false
/Applications/scilab-2024.0.0.app/Contents/bin/scilab: line 928: 47472 
Trace/BPT trap: 5   arch -arch arm64 "$SCILABBIN" "$@“

Error window:
’Scilab-bin quit unexpectedly’
Process:   scilab-bin [47396]
Path:  /Applications/scilab-2024.0.0.app/Contents/bin/scilab-bin
Identifier:scilab-bin
Version:   ???
Code Type: ARM-64 (Native)
Parent Process:Exited process [47341]
Responsible:   Terminal [3476]
User ID:   501

Date/Time: 2024-01-30 10:23:37.3711 +0100
OS Version:macOS 14.4 (23E5180j)
Report Version:12
Anonymous UUID:38D39807-3D05-7015-C951-3FAC62AAEA16

Sleep/Wake UUID:   3F154BC0-177A-4E5D-A1CE-FA16346EB121


On 30.01.2024, at 08:03, Stéphane Mottelet  wrote:
> 
> Hello,
> 
> What did you change ? OS version, startup files, toolboxes installed ? First 
> step as usual, launch Scilab from a terminal to see something more.
> 
> S.
> 
>> Le 30 janv. 2024 à 04:55, Heinz Nabielek  a écrit :
>> 
>> Suddenly, Scilab crashes every time I try to start. What happened? Equally 
>> so after reinstallation.
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

[Scilab-users] Scilab crashes every time

[Heinz Nabielek]

Suddenly, Scilab crashes every time I try to start. What happened? Equally so 
after reinstallation.
Heinz

Scilab 2024.0.0 - macOS 64 bits (ARM) (dmg)
Apple M1 macOS Sonoma 14.4 Beta
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Scilab keeps freezing when trying to save diagram as *.pdf

[Heinz Nabielek]

sw_vers

ProductName: macOS

ProductVersion: 14.3

BuildVersion: 23D5033f

  "Scilab Version: ""2024.0.0.1698152278"
  "Operating System: "  "Mac OS X 14.3"
  "Java version: "  "17.0.8.1"
  "Java runtime information: "  "OpenJDK Runtime Environment (build 
17.0.8.1+1)"
  "Java Virtual Machine information: "  "OpenJDK 64-Bit Server VM (build 
17.0.8.1+1, mixed mode)"
  "Vendor specification: "  "Oracle Corporation"

histplot Export pdf —>. SCILAB FROZEN
Only since 3 days, and I had not changed anything on the M1 Mac. I have no clue 
what could have happened, everything else the same during last 5 weeks.
Heinz



sw_vers
ProductName: Mac OS X
ProductVersion: 10.15.7
BuildVersion: 19H2026

 "Scilab Version: ""2024.0.0.1698152278"
 "Operating System: "  "Mac OS X 10.15.7"
 "Java version: "  "17.0.8.1"
 "Java runtime information: "  "OpenJDK Runtime Environment (build 
17.0.8.1+1)"
 "Java Virtual Machine information: "  "OpenJDK 64-Bit Server VM (build 
17.0.8.1+1, mixed mode, sharing)"
 "Vendor specification: "  "Oracle Corporation"

histplot Export pdf NEVER A PROBLEM


OK. You should give more context and steps, because I am not able to
reproduce the freeze with default choices, i.e. export with default file
name and default path. What did you change in your computer yesterday ?

S.

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab keeps freezing when trying to save diagram as *.pdf

[Heinz Nabielek]

On 18.12.2023, at 08:41, Stéphane Mottelet  wrote:
>
> Hello,
>
> Do you have a simple example, with a precise description of the
> different steps to reproduce the problem ?

histplot
Export to pdf. —> FREEZE Scilab (not responding)

Now, I have used this export to pdf command successfully for years on a variety 
of Macs and PCS. The FREEZE is new since yesterday.
With save as - same problem.
Heinz




sw_vers
ProductName: macOS
ProductVersion: 14.3
BuildVersion: 23D5033f

uname -a
Darwin AppleSilicon.local 23.3.0 Darwin Kernel Version 23.3.0: Fri Dec  1 
03:16:57 PST 2023; root:xnu-10002.80.11~58/RELEASE_ARM64_T8103 arm64

--> ver
  "Scilab Version: ""2024.0.0.1698152278"
  "Operating System: "  "Mac OS X 14.3"
  "Java version: "  "17.0.8.1"
  "Java runtime information: "  "OpenJDK Runtime Environment (build 
17.0.8.1+1)"
  "Java Virtual Machine information: "  "OpenJDK 64-Bit Server VM (build 
17.0.8.1+1, mixed mode)"
  "Vendor specification: "  "Oracle Corporation"

MESSAGE AFTER FORCED QUIT
Date/Time:2023-12-18 16:23:29.418 +0100
End time: 2023-12-18 16:24:55.663 +0100
OS Version:   macOS 14.3 (Build 23D5033f)
Architecture: arm64e
Report Version:   44
Incident Identifier: 31710C48-FF69-4BBB-9559-1BADEE0DB9AE

Data Source:  Stackshots
Shared Cache: 01283273-9123-3D4C-8494-E5FA1711DEC5 slid base address 
0x182bcc000, slide 0x2bcc000 (System Primary)
Shared Cache: 2612A53C-5071-3956-BFD7-6B75EF2C2607 slid base address 
0x1a4abc000, slide 0x24abc000 (DriverKit)
Shared Cache: F228BFFA-E0C6-3B72-B877-1994D9CBA12F slid base address 
0x7ff817774000, slide 0x17774000 (Rosetta)

Command:  scilab-bin
Path: /Applications/scilab-2024.0.0.app/Contents/bin/scilab-bin
Codesigning ID:   scilab-bin
Team ID:  A2A9VZ4WJT
Architecture: arm64
Parent:   bash [29588] [unique pid 128961]
Responsible:  bash [29588] [unique pid 128961]
PID:  29643
Time Since Fork:  117s

Event:hang
Duration: 86.24s
Duration Sampled: 1.10s (process was unresponsive for 85 seconds before 
sampling)
Steps:11 (100ms sampling interval)

Hardware model:   iMac21,1
Active cpus:  8
HW page size: 16384
VM page size: 16384
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] On computational speed

[Heinz Nabielek]

Thanks. In the 1960s FORTRAN, we would always write A2=A*A rather than A2=A**2.

Below the Federico examples on my M1 Apple silicon iMac:
Heinz

--> tic, u = rand(1,1e7); toc
 ans  =   0.111413

--> tic, v = sqrt(u); toc
 ans  =   0.0384130

-->  tic, v = u.*u; toc
 ans  =   0.023148

--> tic, v = u.*%pi; toc
 ans  =   0.017872

--> tic, v = u.^2; toc
 ans  =   0.156129

--> tic, v = u.^%pi; toc
 ans  =   0.1537830

--> tic, v = log(u); toc
 ans  =   0.093586

--> tic, v = exp(u); toc
 ans  =   0.0821170

--> tic, v = exp(2*log(u)); toc
 ans  =   0.139802





> Am 30.11.2023 um 00:46 schrieb Federico Miyara :
>
> Dear All,
>
> I was trying to make some simple experiments regarding computational speed of 
> several functions and operations and some results puzzled me. To perform the 
> tests I used massive data to minimize the effect of idle time due to priority 
> handling by the operating system or whatever.
>
> --> tic, u = rand(1,1e7); toc
>  ans  =
>0.4948326
>
> This creates ten million random numbers in half a second. Very decent 
> performance. Then
>
> --> tic, v = sqrt(u); toc
>  ans  =
>0.282802
>
> The algorithm for square root seems quite fast. Multiplication
>
> --> tic, v = u.*u; toc
>  ans  =
>0.1218028
>
> --> tic, v = u.*%pi; toc
>  ans  =
>0.119245
>
> is really fast. But when it comes to power or exponentiation:
>
> --> tic, v = u.^2; toc
>  ans  =
>1.9633435
>
> --> tic, v = u.^%pi; toc
>  ans  =
>1.958321
>
> it is quite slow. Logarithm is much faster:
>
> --> tic, v = log(u); toc
>  ans  =
>0.5166959
>
> Even the exponential function with base e is much faster:
>
> --> tic, v = exp(u); toc
>  ans  =
>0.5451482
>
> I wonder why special cases such as small integer powers are not dealt with by 
> simple multiplication. Particularly the square is very frequently used on 
> large vectors representing discrete signals (for instance to get energy 
> measures) and no warning is issued in the documentation.
>
> I also notice that computing the square using this formula (a cannon to kill 
> a mosquito!)
>
> --> tic, v = exp(2*log(u)); toc
>  ans  =
>1.2789763
>
> is faster than using the power operator ^.
>
> This seems to be an inconsistent approach, since many special functions are 
> painstakingly tailored to be very efficient, yet one of the simplest of all 
> operations, raising to power 2, is sluggishly slow...
>
> Regards,
>
> Federico Miyara
>
>
> Libre de virus.www.avast.com This email and any attachments are intended 
> solely for the use of the individual or entity to whom it is addressed and 
> may be confidential and/or privileged.
> If you are not one of the named recipients or have received this email in 
> error,
> (i) you should not read, disclose, or copy it,
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
>
> Please be informed that your personal data are processed according to our 
> data privacy policy as described on our website. Should you have any 
> questions related to personal data protection, please contact 3DS Data 
> Protection Officer https://www.3ds.com/privacy-policy/contact/
>
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] my history file from the previous Scilab version has lost much of its first part

[Heinz Nabielek]

Great explanation, thank you. Would I have known
Heinz

What exactly should I do here? Sorry for being dumb
[cid:F854B4EB-87E6-4A23-814A-333E7B1D59EA@home]

On 11.11.2023, at 02:34, Federico Miyara 
mailto:fmiy...@fceia.unr.edu.ar>> wrote:


Heinz,

The history file has a maximum number of lines that can be set by the
user from the Scilab preferences, Command history. I've seen that you
can set it at a rather high value. I have it at 10. But eventually
it will reach that maximum and I guess it will be cleared following a
first-in first-out protocol so that the oldest information is
periodically cleared.

I'm afraid that if you reached the limit the information is definitely
and irreversibly lost. The only way it would be somehow retrievable
would be if you had noticed such situation right away when it started to
happen and that nothing more has been done so that the probability that
the information of the old version of the file, even if not accessible
from the file system, is still there and can be located by a specialized
recovery application such as Recuva (see
https://www.ccleaner.com/es-es/recuva/download, there is a free
version). But I guess one discovers that when it is too late.

Anyway, it isn't a good policy to trust important information to a
feature that is not intended for that purpose. What I do is to have a
file named drafts.sce where I paste ideas, fragments of scripts,
miscellaneous data and so on.

Regards,

Federico Miyara

On 10/11/2023 21:35, Heinz Nabielek wrote:
Scilab does a good job writing all commands into a history file specific for 
every Scilab version.
Now my history file from the previous Scilab version has lost much of its first 
part. Could that have been because of truncation because too big?
Is any of this retrievable? Contains important data!
Heinz
___
users mailing list - users@lists.scilab.org<mailto:users@lists.scilab.org>
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/




--
Este correo electrónico ha sido analizado en busca de virus por el software 
antivirus de Avast.
www.avast.com<http://www.avast.com>
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-un

[Scilab-users] my history file from the previous Scilab version has lost much of its first part

[Heinz Nabielek]

Scilab does a good job writing all commands into a history file specific for 
every Scilab version.
Now my history file from the previous Scilab version has lost much of its first 
part. Could that have been because of truncation because too big?
Is any of this retrievable? Contains important data!
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] 2024 Scilab versions do not allow copy and paste into filenames

[Heinz Nabielek]

On 09.11.2023, at 00:27, Stéphane Mottelet  wrote:
>
> Hello Heinz,
>
> You may have noticed that save/load dialogs are now native macOS
> dialogs, this was a long awaited feature. There may be some features
> that have not been tested and it seems that the classical copy/paste
> shortcuts do not work between strings in these dialogs, but you should
> be able to copy/paste by using right-click on the string to display a
> popup menu.

PERFECT!
Thanks
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] 2024 Scilab versions do not allow copy and paste into filenames

[Heinz Nabielek]

The 2024 Scilab versions do not allow copy and paste into filenames when 
storing a diagram. Both Intel Macs and Apple silicon Macs.

This was no problem in the last 10 years: I always store a graphic file with 
the name from the title of the diagram and this does not work any more.
What happened?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Plot sequence: Problem only, one colour is white. What is the easiest way to change the colour sequence?

[Heinz Nabielek]

So sorry, my fault. From 9 data series, I saw only 8 on my plot. But, when 
checking, this was because one was on top of another and not because one was 
white.

Thanks for the patient explanation.
Heinz

> On 30.10.2023, at 08:30, Federico Miyara  wrote:
>
>
> Heinz,
>
> I think you asked this very question in the past. I don't experience your 
> problem.
>
> According to https://help.scilab.org/LineSpec if you don't specify the color 
> it cycles through the table
>
> R G   B
> 0.0.  1.
> 0.0.5 0.
> 1.0.  0.
> 0.0.750.75
> 0.75  0.  0.75
> 0.75  0.750.
> 0.25  0.250.25
>
> None of these colors is white (which would be RGB = [1 1 1]). This code
>
> scf(1);
> clf(1);
> // Create sample x and y
> x = 0:0.01:1;
> n = 9;
> y = sin(x'*(1:n));
> plot(x',y)
>
> generates this 9 curve plot
> 
> which cycles from blue, green, ... to dark gray and starts over, blue, green. 
> Only 7 colors, two of them repeated.
>
> So I guess you have somehow assigned a different color map. In my example, 
> entering
>
> gce().children.foreground
>
> yields
>
>  ans  =
>33.
>2.
>37.
>36.
>35.
>34.
>5.
>33.
>2.
>
> If now you execute
>
> getcolor
>
> you get the following interactive color chart,
>
> 
>
> which contains a sample of each color in the color map, which are indexed 
> from top to bottom and left to right. Clicking on any one you get at the 
> bottom the index and the RGB formula. Referring to the table above, the first 
> curve (starting from the bottom) has the index 2, which correspondes to blue, 
> the second has the index 33, which you can see it is the 3rd color of the 6th 
> column, i.e., green; then comes 5, the 5th of the 1st column, i.e., red, so 
> the indices are consistent with the colors of the curves and with the colors 
> in the color map.
>
> WARNING: Don't forget to close the interactive chart or else Scilab will seem 
> to freeze!
>
> Now you can force the cycling to another order, for instance instead of 2, 
> 33, 5, ..., 37 you could use 1, 2, 3, ..., 37. I guess you did this which by 
> the way can be done this way:
>
> gce().children.foreground = (1:n)';
>
> This assigns the indices 1 to 9 to the foreground color or line color of each 
> curve. The problem is that, indeed, index 8 corresponds to white. One easy 
> way to get rid of the white is
>
> gce().children.foreground = [1:7, 9:10]';
>
> an you'll be done!
>
> However, the problem of further customizing the color map is an interesting 
> one if you want to have complete control of your colors. You can do so by 
> setting a color map from scratch using  gcf().color_map. You just equate it 
> to any mx3 matrix, where m is the number of colors, one per row, where each 
> row is the RGB formula of the color given by three numbers from 0 (darkest) 
> to 1 (most intense).
>
> However, if you don't have some criterion your colors may result caotic. One 
> way to prevent that is to use one of the color maps offered by Scilab. For 
> instance I like a lot this one:
>
> gcf().color_map = jetcolormap(n);
>
> They are called swatches, and you can find the available ones at 
> https://help.scilab.org/colormap
>
> Regards,
>
> Federico Miyara
>
>
>
> On 28/10/2023 22:45, Heinz Nabielek wrote:
>> I have a vector x with  --> size(x) = 13.   1.
>> and an array y with --> size(y) = 13.   9.
>>
>> Scilab has that beautiful power that I can plot them all at once with 9 
>> different colours
>> --> plot(x,y,’o’);
>>
>> Problem only, one colour is white.
>> What is the easiest way to change the colour sequence?
>> Greetings
>> Heinz

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Plot sequence: Problem only, one colour is white. What is the easiest way to change the colour sequence?

[Heinz Nabielek]

I have a vector x with  --> size(x) = 13.   1.
and an array y with --> size(y) = 13.   9.

Scilab has that beautiful power that I can plot them all at once with 9 
different colours
--> plot(x,y,’o’);

Problem only, one colour is white.
What is the easiest way to change the colour sequence?
Greetings
Heinz



___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] SCILAB Apple Silicon iMac Factor 4 faster

[Heinz Nabielek]

a=[];timer();for j=1:3;y=zeros(1,12);i=0;while 
~and(y);z=grand(1,1,'uin',1,12);y(z)=y(z)+1;i=i+1;end;a=[a,i];end; timer()
 ans  =   26.895136 INTEL iMac
 ans  =6.4323840Apple Silicon iMac

Factor 4 faster

Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] ARM version works great and fast….

[Heinz Nabielek]

ARM version works great and fast….
Heinz

--> timer();A=rand(1,1);timer()
 ans  =0.8124390

--> ver
 ans  =
 column 1
  "Scilab Version: "
  "Operating System: "
  "Java version: "
  "Java runtime information: "
  "Java Virtual Machine information: "
  "Vendor specification: "

 column 2
  "2024.0.0.1698152278"
  "Mac OS X 14.1"
  "17.0.8.1"
  "OpenJDK Runtime Environment (build 17.0.8.1+1)"
  "OpenJDK 64-Bit Server VM (build 17.0.8.1+1, mixed mode)"
  "Oracle Corporation"
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scilab 2024.0.0 refuses to start up on my Intel iMac

[Heinz Nabielek]

I am not certain I did it right.
Heinz

(base) heinz@Heinzs-iMac ~ % 
/Applications/scilab-2024.0.0.app/Contents/bin/scilab ; exit;
dyld: Library not loaded: @rpath/libklu.1.dylib
  Referenced from: /Applications/scilab-2024.0.0.app/Contents/bin/scilab-bin
  Reason: no suitable image found.  Did find:

/Applications/scilab-2024.0.0.app/Contents/bin/../lib/thirdparty/libklu.1.dylib:
 code signing blocked mmap() of 
'/Applications/scilab-2024.0.0.app/Contents/bin/../lib/thirdparty/libklu.1.dylib'

/Applications/scilab-2024.0.0.app/Contents/lib/thirdparty/libklu.1.dylib: code 
signing blocked mmap() of 
'/Applications/scilab-2024.0.0.app/Contents/lib/thirdparty/libklu.1.dylib'
/Applications/scilab-2024.0.0.app/Contents/bin/scilab: line 928:  2287 Abort 
trap: 6   "$SCILABBIN" "$@"

> On 25.10.2023, at 18:53, Stéphane Mottelet  wrote:
>
> Hello Heinz,
>
> Can you start Scilab from the command line and tell me what you see ?
>
> S.
>
> On 25/10/2023 18:47, Heinz Nabielek wrote:
>> I have downloaded Scilab 2024.0.0 - macOS 64 bits (Intel) (dmg).
>> But Scilab 2024.0.0 refuses to start up on my Intel iMac macOS 10.15.7 
>> (19H2026) Darwin 19.6.0
>>
>> What can I do?
>> Heinz
>> ___
>> users mailing list - users@lists.scilab.org
>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>> https://antiphishing.vadesecure.com/v4?f=bHdDQW5tZDVCemI1ZVczSVYUQ55Sw8XPuQVe9hm6kqysUQCM2UGjdKIXACCktR0y=UERGdHg5cm1GRGl1YjhpePkNkpYwiHjXsr1SPKiLzhc=19x6=ZEtPTklHeGR1a0VPT25scWufIKyLsieczpW0GQvrmzYhJoXZ-Jppo4DD58kxr3mQ=90d0593ce5c0ca95e0f26ec20c4d05a234e96994f7fdcc71b8ac1bdbb000f9b3=https%3A%2F%2Flists.scilab.org%2Fmailman%2Flistinfo%2Fusers
>> This email and any attachments are intended solely for the use of the 
>> individual or entity to whom it is addressed and may be confidential and/or 
>> privileged.
>>
>> If you are not one of the named recipients or have received this email in 
>> error,
>>
>> (i) you should not read, disclose, or copy it,
>>
>> (ii) please notify sender of your receipt by reply email and delete this 
>> email and all attachments,
>>
>> (iii) Dassault Systèmes does not accept or assume any liability or 
>> responsibility for any use of or reliance on this email.
>>
>>
>> Please be informed that your personal data are processed according to our 
>> data privacy policy as described on our website. Should you have any 
>> questions related to personal data protection, please contact 3DS Data 
>> Protection Officer 
>> https://antiphishing.vadesecure.com/v4?f=bHdDQW5tZDVCemI1ZVczSVYUQ55Sw8XPuQVe9hm6kqysUQCM2UGjdKIXACCktR0y=UERGdHg5cm1GRGl1YjhpePkNkpYwiHjXsr1SPKiLzhc=19x6=ZEtPTklHeGR1a0VPT25scWufIKyLsieczpW0GQvrmzYhJoXZ-Jppo4DD58kxr3mQ=ea4e5e0d14dbdc22b8975351f3b2de774acd096f08f838269726d8d7b511bbba=https%3A%2F%2Fwww.3ds.com%2Fprivacy-policy%2Fcontact%2F
>>
> --
> Stéphane Mottelet
> Ingénieur de recherche
> EA 4297 Transformations Intégrées de la Matière Renouvelable
> Département Génie des Procédés Industriels
> Sorbonne Universités - Université de Technologie de Compiègne
> CS 60319, 60203 Compiègne cedex
> Tel : +33(0)344234688
> http://www.utc.fr/~mottelet
>
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
>
> If you are not one of the named recipients or have received this email in 
> error,
>
> (i) you should not read, disclose, or copy it,
>
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
>
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
>
>
> Please be informed that your personal data are processed according to our 
> data privacy policy as described on our website. Should you have any 
> questions related to personal data protection, please contact 3DS Data 
> Protection Officer https://www.3ds.com/privacy-policy/contact/
>

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privil

[Scilab-users] Scilab 2024.0.0 refuses to start up on my Intel iMac

[Heinz Nabielek]

I have downloaded Scilab 2024.0.0 - macOS 64 bits (Intel) (dmg).
But Scilab 2024.0.0 refuses to start up on my Intel iMac macOS 10.15.7 
(19H2026) Darwin 19.6.0

What can I do?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Apple silicon speed improvements

[Heinz Nabielek]

Friends:

below code to simulate the Sophienalpe (I can explain if interested) problem:

timer();a=[];for j=1:3;y=zeros(1,12);i=0;while 
~and(y);z=grand(1,1,'uin',1,12);y(z)=y(z)+1;i=i+1;end;a=[a,i];end; timer()

On my new Apple silicon M1 iMac it takes 6.7 s;
on my old i7 iMac 3.2 GHz it always took 27.5 s.

This speed improvement is very useful for my Monte-Carlo simulations.
And thanks for providing an Apple silicon SciLab version.
Greetings
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Apple silicon speed improvements

[Heinz Nabielek]

Friends:

below code to simulate the Sophienalpe (I can explain if interested) problem:

timer();a=[];for j=1:3;y=zeros(1,12);i=0;while 
~and(y);z=grand(1,1,'uin',1,12);y(z)=y(z)+1;i=i+1;end;a=[a,i];end; timer()

On my new Apple silicon M1 iMac it takes 6.7 s;
on my old i7 iMac 3.2 GHz it always took 27.5 s.

This speed improvement is very useful for my Monte-Carlo simulations.
And thanks for providing an Apple silicon SciLab version.
Greetings
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] [Scilab community] [Specific domains/Graphics and plotting] How can I make errbar thicker, colours, provide endcap?

[Heinz Nabielek]

Great many thanks.
Heinz

> On 21.08.2023, at 15:21, philipp via Scilab community 
>  wrote:
>
>   philipp
> August 21
> What about:
>
> errbar(); // NOTE: add vertical error bars on a 2D plot
>
> a = gca();
> e = a.children(1);
> e.thickness = 2; // changing error bar thickness
> e.segs_color = 2; // changing error bar colors
>
> /*
> for colors, individual error bars can be reached like this:
> e.segs_color(i) = …
> with i = 1:nrOfErrorBars
> */
>
> // Some kind of endcaps:
>
> e.mark_mode = ‘on’
> e.mark_style = 2;
> e.mark_size = 0;
> e.mark_foreground = 2;
>
> If you like to have horizontal lines as endcaps, a nasty option would be to 
> plot these line separatly since there is no “-” in LineSpec → marker type. 
> (Scilab 6.1.0)
>
> Besides: I guess errbar() plots the vertikal lines all as single lines…see: 
> e.data
>
> BR
> Philipp
>
> Visit Topic or reply to this email to respond.
>
> To unsubscribe from these emails, click here.
>

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] linear least-squares with errors

[Heinz Nabielek]

On 15.08.2023, at 12:00, users-requ...@lists.scilab.org wrote:
>   1.  linear least-squares with errors (Heinz Nabielek)
> --
>
>
> I am doing my linear least-squares fits without any subroutine, just by 
> inversion of a rectangular array. How do I handle this, when every y-value 
> has its individual error range?

Dead simple: divide your measurement vector and you rectangular matrix by the 
error vector. How could I forget?
Heinz

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



Weighted Least Squares.pdf
Description: Adobe PDF document
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

[Scilab-users] linear least-squares with errors

[Heinz Nabielek]

I am doing my linear least-squares fits without any subroutine, just by 
inversion of a rectangular array. How do I handle this, when every y-value has 
its individual error range?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] how do I plot 9 curves with automatic colour sequencing?

[Heinz Nabielek]

Colleagues:

[sorry, has been asked before but forgot]: how do I plot 9 curves with 
automatic colour sequencing?

o   Problem #1: colour 8 is white and I cannot see my curve
o   Problem #2: colours >= 9 become rather non distinct

Is there an easy solution without complex SciLab acrobatics?
Heinz

Specifier   Color
#1  r   Red
#2  g   Green
#3  b   Blue
#4  c   Cyan
#5  m   Magenta
#6  y   Yellow
#7  k   Black
#8  w   White
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scilab on the M1/M2 Mac

[Heinz Nabielek]

I had tested both native builds and I had reported that both crash immediately 
with a call to cdfnor
Heinz



>> scilab-6.1.1-accelerate-arm64.dmg installed and started well
>>  "6.1.1.988271013"
>>  "Mac OS X 13.3.1"
>>  "1.8.0_292"
>>  "OpenJDK Runtime Environment (build 1.8.0_292-b10)"
>>  "OpenJDK 64-Bit Server VM (build 25.292-b10, mixed mode)"
>>  "Oracle Corporation“
>>
>> When I attempted execution speed comparisons, it failed quickly…
>> NNN=300;
>> WT=12; hub=150; radius=110; r2=radius^2;
>> n=60;r=0.9868;a=ones(1:n);m=2;vm=7; vc=vm/(sqrt(%pi)/2);
>> Z=grand(1,n,'nor',0,1);V=Z;for i=2:n;V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);end;
>> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);v=vc*((-log(VZ))^(1/m))’;
>>
>> UNEXPECTED TERMINATION….
>> Process:   scilab-bin [715]
>> Path:  
>> /Applications/scilab-6.1.1.app/Contents/MacOS/scilab-bin
>> Identifier:scilab-bin
>> Version:   ???
>> Code Type: ARM-64 (Native)
>> Parent Process:Exited process [655]
>> User ID:   501
>>
>> Date/Time: 2023-05-10 01:36:45.6480 +0200
>> OS Version:macOS 13.3.1 (22E772610a)
>> Report Version:12
>> Anonymous UUID:38D39807-3D05-7015-C951-3FAC62AAEA16
>>
>> What happened?
>> Heinz
>>
>> heinznabielek@iMac-von-Heinz ~ % sw_vers
>> ProductName: macOS
>> ProductVersion: 13.3.1
>> ProductVersionExtra: (a)
>> BuildVersion: 22E772610a
>> heinznabielek@iMac-von-Heinz ~ % uname -a
>> Darwin iMac-von-Heinz.local 22.4.0 Darwin Kernel Version 22.4.0:
>> Mon Mar  6 21:00:41 PST 2023; root:xnu-8796.101.5~3/RELEASE_ARM64_T8103 arm64

> On 11.05.2023, at 18:15, Stéphane Mottelet  wrote:
>
> Can you also test the two native builds ?
>
> S.
>
> Le 11/05/2023 à 05:27, Heinz Nabielek a écrit :
>> scilab-2023.0.0-x86_64.dmg is running perfectly well on my new M1 iMac, but 
>> my complex wind-bird swarm-avian mortality simulation code is 30% slower 
>> than on the old i7 iMac!
>> Heinz
>>
>>
>>
>>
>>
>>> Am 10/05/2023 um 02:56 schrieb Heinz Nabielek :
>>>
>>> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);
>>> is the command that kills scilab-6.1.1-accelerate-arm64
>>> Heinz
>>>
>>>
>>>
>>>
>>>
>>>
>>> NNN=300;
>>> WT=12;
>>> hub=150;
>>> radius=110;
>>> r2=radius^2;
>>> n=60;
>>> r=0.9868;
>>> a=ones(1:n);
>>> m=2;
>>> vm=7;
>>> vc=vm/(sqrt(%pi)/2);
>>> Z=grand(1,n,'nor',0,1);
>>> V=Z;
>>> for i=2:n;
>>> V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);
>>> end;
>>> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);
>>>
>>>> Am 10/05/2023 um 01:50 schrieb Heinz Nabielek :
>>>>
>>>> scilab-6.1.1-accelerate-arm64.dmg installed and started well
>>>> "6.1.1.988271013"
>>>> "Mac OS X 13.3.1"
>>>> "1.8.0_292"
>>>> "OpenJDK Runtime Environment (build 1.8.0_292-b10)"
>>>> "OpenJDK 64-Bit Server VM (build 25.292-b10, mixed mode)"
>>>> "Oracle Corporation“
>>>>
>>>> When I attempted execution speed comparisons, it failed quickly…
>>>> NNN=300;
>>>> WT=12; hub=150; radius=110; r2=radius^2;
>>>> n=60;r=0.9868;a=ones(1:n);m=2;vm=7; vc=vm/(sqrt(%pi)/2);
>>>> Z=grand(1,n,'nor',0,1);V=Z;for i=2:n;V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);end;
>>>> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);v=vc*((-log(VZ))^(1/m))’;
>>>>
>>>> UNEXPECTED TERMINATION….
>>>> Process:   scilab-bin [715]
>>>> Path:  
>>>> /Applications/scilab-6.1.1.app/Contents/MacOS/scilab-bin
>>>> Identifier:scilab-bin
>>>> Version:   ???
>>>> Code Type: ARM-64 (Native)
>>>> Parent Process:Exited process [655]
>>>> User ID:   501
>>>>
>>>> Date/Time: 2023-05-10 01:36:45.6480 +0200
>>>> OS Version:macOS 13.3.1 (22E772610a)
>>>> Report Version:12
>>>> Anonymous UUID:38D39807-3D05-7015-C951-3FAC62AAEA16
>>>>
>>>> What happened?
>>>> Heinz
>>>>
>>>> heinznabielek@iMac-von-Heinz ~ % sw_vers
>>

Re: [Scilab-users] Scilab on the M1/M2 Mac

[Heinz Nabielek]

scilab-2023.0.0-x86_64.dmg is running perfectly well on my new M1 iMac, but my 
complex wind-bird swarm-avian mortality simulation code is 30% slower than on 
the old i7 iMac!
Heinz





> Am 10/05/2023 um 02:56 schrieb Heinz Nabielek :
>
> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);
> is the command that kills scilab-6.1.1-accelerate-arm64
> Heinz
>
>
>
>
>
>
> NNN=300;
> WT=12;
> hub=150;
> radius=110;
> r2=radius^2;
> n=60;
> r=0.9868;
> a=ones(1:n);
> m=2;
> vm=7;
> vc=vm/(sqrt(%pi)/2);
> Z=grand(1,n,'nor',0,1);
> V=Z;
> for i=2:n;
> V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);
> end;
> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);
>
>> Am 10/05/2023 um 01:50 schrieb Heinz Nabielek :
>>
>> scilab-6.1.1-accelerate-arm64.dmg installed and started well
>> "6.1.1.988271013"
>> "Mac OS X 13.3.1"
>> "1.8.0_292"
>> "OpenJDK Runtime Environment (build 1.8.0_292-b10)"
>> "OpenJDK 64-Bit Server VM (build 25.292-b10, mixed mode)"
>> "Oracle Corporation“
>>
>> When I attempted execution speed comparisons, it failed quickly…
>> NNN=300;
>> WT=12; hub=150; radius=110; r2=radius^2;
>> n=60;r=0.9868;a=ones(1:n);m=2;vm=7; vc=vm/(sqrt(%pi)/2);
>> Z=grand(1,n,'nor',0,1);V=Z;for i=2:n;V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);end;
>> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);v=vc*((-log(VZ))^(1/m))’;
>>
>> UNEXPECTED TERMINATION….
>> Process:   scilab-bin [715]
>> Path:  
>> /Applications/scilab-6.1.1.app/Contents/MacOS/scilab-bin
>> Identifier:scilab-bin
>> Version:   ???
>> Code Type: ARM-64 (Native)
>> Parent Process:    Exited process [655]
>> User ID:   501
>>
>> Date/Time: 2023-05-10 01:36:45.6480 +0200
>> OS Version:macOS 13.3.1 (22E772610a)
>> Report Version:12
>> Anonymous UUID:38D39807-3D05-7015-C951-3FAC62AAEA16
>>
>> What happened?
>> Heinz
>>
>> heinznabielek@iMac-von-Heinz ~ % sw_vers
>> ProductName: macOS
>> ProductVersion: 13.3.1
>> ProductVersionExtra: (a)
>> BuildVersion: 22E772610a
>> heinznabielek@iMac-von-Heinz ~ % uname -a
>> Darwin iMac-von-Heinz.local 22.4.0 Darwin Kernel Version 22.4.0:
>> Mon Mar  6 21:00:41 PST 2023; root:xnu-8796.101.5~3/RELEASE_ARM64_T8103 arm64
>>
>>
>>
>>> Am 09/05/2023 um 20:59 schrieb Stéphane Mottelet :
>>>
>>> Hello,
>>>
>>> You can already test the Intel build and the older M1 builds of the 6.1.1 
>>> version, which will already run quite fast (in emulation mode for the Intel 
>>> build). I am not sure we will have an official M1/M2 release before  the 
>>> 2024.0.0, but it is likely that we will build it only with refBlas.
>>>
>>> S.
>>>
>>>> Le 8 mai 2023 à 01:15, Heinz Nabielek  a écrit :
>>>>
>>>> Dear Stéphane,
>>>>
>>>> I have ordered a 24" M1 iMac (not yet arrived).
>>>>
>>>> From the M1 options described in 
>>>> <https://www.utc.fr/~mottelet/scilab_for_macOS.html>, which is most 
>>>> straightforward for me?
>>>> -   I need fast vector operations
>>>> -   I do not use atoms
>>>> -   I use complex Scilab statistical functions
>>>>
>>>> Heinz
>

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scilab on the M1/M2 Mac

[Heinz Nabielek]

[VZ VZ2]=cdfnor("PQ",V,0*a,1*a);
is the command that kills scilab-6.1.1-accelerate-arm64
Heinz






NNN=300;
WT=12;
hub=150;
radius=110;
r2=radius^2;
n=60;
r=0.9868;
a=ones(1:n);
m=2;
vm=7;
vc=vm/(sqrt(%pi)/2);
Z=grand(1,n,'nor',0,1);
V=Z;
for i=2:n;
V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);
end;
[VZ VZ2]=cdfnor("PQ",V,0*a,1*a);

> Am 10/05/2023 um 01:50 schrieb Heinz Nabielek :
>
> scilab-6.1.1-accelerate-arm64.dmg installed and started well
>  "6.1.1.988271013"
>  "Mac OS X 13.3.1"
>  "1.8.0_292"
>  "OpenJDK Runtime Environment (build 1.8.0_292-b10)"
>  "OpenJDK 64-Bit Server VM (build 25.292-b10, mixed mode)"
>  "Oracle Corporation“
>
> When I attempted execution speed comparisons, it failed quickly…
> NNN=300;
> WT=12; hub=150; radius=110; r2=radius^2;
> n=60;r=0.9868;a=ones(1:n);m=2;vm=7; vc=vm/(sqrt(%pi)/2);
> Z=grand(1,n,'nor',0,1);V=Z;for i=2:n;V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);end;
> [VZ VZ2]=cdfnor("PQ",V,0*a,1*a);v=vc*((-log(VZ))^(1/m))’;
>
> UNEXPECTED TERMINATION….
> Process:   scilab-bin [715]
> Path:  
> /Applications/scilab-6.1.1.app/Contents/MacOS/scilab-bin
> Identifier:scilab-bin
> Version:   ???
> Code Type: ARM-64 (Native)
> Parent Process:Exited process [655]
> User ID:   501
>
> Date/Time: 2023-05-10 01:36:45.6480 +0200
> OS Version:    macOS 13.3.1 (22E772610a)
> Report Version:12
> Anonymous UUID:38D39807-3D05-7015-C951-3FAC62AAEA16
>
> What happened?
> Heinz
>
> heinznabielek@iMac-von-Heinz ~ % sw_vers
> ProductName: macOS
> ProductVersion: 13.3.1
> ProductVersionExtra: (a)
> BuildVersion: 22E772610a
> heinznabielek@iMac-von-Heinz ~ % uname -a
> Darwin iMac-von-Heinz.local 22.4.0 Darwin Kernel Version 22.4.0:
> Mon Mar  6 21:00:41 PST 2023; root:xnu-8796.101.5~3/RELEASE_ARM64_T8103 arm64
>
>
>
>> Am 09/05/2023 um 20:59 schrieb Stéphane Mottelet :
>>
>> Hello,
>>
>> You can already test the Intel build and the older M1 builds of the 6.1.1 
>> version, which will already run quite fast (in emulation mode for the Intel 
>> build). I am not sure we will have an official M1/M2 release before  the 
>> 2024.0.0, but it is likely that we will build it only with refBlas.
>>
>> S.
>>
>>> Le 8 mai 2023 à 01:15, Heinz Nabielek  a écrit :
>>>
>>> Dear Stéphane,
>>>
>>> I have ordered a 24" M1 iMac (not yet arrived).
>>>
>>> From the M1 options described in 
>>> <https://www.utc.fr/~mottelet/scilab_for_macOS.html>, which is most 
>>> straightforward for me?
>>> -   I need fast vector operations
>>> -   I do not use atoms
>>> -   I use complex Scilab statistical functions
>>>
>>> Heinz

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scilab on the M1/M2 Mac

[Heinz Nabielek]

scilab-6.1.1-accelerate-arm64.dmg installed and started well
  "6.1.1.988271013"
  "Mac OS X 13.3.1"
  "1.8.0_292"
  "OpenJDK Runtime Environment (build 1.8.0_292-b10)"
  "OpenJDK 64-Bit Server VM (build 25.292-b10, mixed mode)"
  "Oracle Corporation“

When I attempted execution speed comparisons, it failed quickly…
NNN=300;
WT=12; hub=150; radius=110; r2=radius^2;
n=60;r=0.9868;a=ones(1:n);m=2;vm=7; vc=vm/(sqrt(%pi)/2);
Z=grand(1,n,'nor',0,1);V=Z;for i=2:n;V(i)=r*V(i-1)+sqrt(1-r^2)*Z(i);end;
[VZ VZ2]=cdfnor("PQ",V,0*a,1*a);v=vc*((-log(VZ))^(1/m))’;

UNEXPECTED TERMINATION….
Process:   scilab-bin [715]
Path:  /Applications/scilab-6.1.1.app/Contents/MacOS/scilab-bin
Identifier:scilab-bin
Version:   ???
Code Type: ARM-64 (Native)
Parent Process:Exited process [655]
User ID:   501

Date/Time: 2023-05-10 01:36:45.6480 +0200
OS Version:macOS 13.3.1 (22E772610a)
Report Version:12
Anonymous UUID:    38D39807-3D05-7015-C951-3FAC62AAEA16

What happened?
Heinz

heinznabielek@iMac-von-Heinz ~ % sw_vers
ProductName: macOS
ProductVersion: 13.3.1
ProductVersionExtra: (a)
BuildVersion: 22E772610a
heinznabielek@iMac-von-Heinz ~ % uname -a
Darwin iMac-von-Heinz.local 22.4.0 Darwin Kernel Version 22.4.0:
Mon Mar  6 21:00:41 PST 2023; root:xnu-8796.101.5~3/RELEASE_ARM64_T8103 arm64



> Am 09/05/2023 um 20:59 schrieb Stéphane Mottelet :
>
> Hello,
>
> You can already test the Intel build and the older M1 builds of the 6.1.1 
> version, which will already run quite fast (in emulation mode for the Intel 
> build). I am not sure we will have an official M1/M2 release before  the 
> 2024.0.0, but it is likely that we will build it only with refBlas.
>
> S.
>
>> Le 8 mai 2023 à 01:15, Heinz Nabielek  a écrit :
>>
>> Dear Stéphane,
>>
>> I have ordered a 24" M1 iMac (not yet arrived).
>>
>> From the M1 options described in 
>> <https://www.utc.fr/~mottelet/scilab_for_macOS.html>, which is most 
>> straightforward for me?
>> -   I need fast vector operations
>> -   I do not use atoms
>> -   I use complex Scilab statistical functions
>>
>> Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Scilab on the M1/M2 Mac

[Heinz Nabielek]

Dear Stéphane,

I have ordered a 24" M1 iMac (not yet arrived).

From the M1 options described in 
<https://www.utc.fr/~mottelet/scilab_for_macOS.html>, which is most 
straightforward for me?
-   I need fast vector operations
-   I do not use atoms
-   I use complex Scilab statistical functions

Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Plotting on left and right y axes

[Heinz Nabielek]

Friends:


here is one working example
Heinz



B=read('155508.txt',155508,2);
v=B(:,2);
plot(35,.2,'k+');
histplot(-0.5:30.5,v,style=2);
xgrid();
h1=gca();h1.font_size=4;
xlabel('wind speed (m/s)','fontsize',4);
ylabel('wind speed probability ()','fontsize',4);
title('Power Generation is given by wind speed probability X turbine performance','fontsize',4);
// Axis y2
h2=newaxes();
c=color("red");
h2.font_color=c;
h2.font_size=4;
plot(35,2,'k+');
function [y]=eff(x)
y=min(((x/10.7).^3),1) .* (x>3) .*(x<=25)
endfunction
xx=0:.1:31;EFF=eff(xx);
plot(xx,EFF,'r-','Linewidth',5);
h2.filled="off";
h2.axes_visible(1)="off";
h2.y_location="right";
h2.children(1).children(1).thickness=2;
ylabel('relative turbine performance','fontsize',4);




This email and any attachments are intended solely for the use of the individual or entity to whom it is addressed and may be confidential and/or privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data privacy policy as described on our website. Should you have any questions related to personal data protection, please contact 3DS Data Protection Officer
https://www.3ds.com/privacy-policy/contact/







Power Generation is given by wind speed probability X turbine performance.pdf
Description: Adobe PDF document
Date: Fri, 28 Apr 2023 17:32:02 +0100From: Samuel Enibe <samuel.en...@unn.edu.ng>To: Users mailing list for Scilab <users@lists.scilab.org>Subject: [Scilab-users] Plotting on left and right y axesMessage-ID:       <CAJKcLoW1r0dOZTUcOf=dj2vgjp_+nbudajdqjk5dtoe0xyh...@mail.gmail.com>Content-Type: text/plain; charset="utf-8"How can I plot two graphs on the same window,  with graph 1 on the left yaxis,  and graph 2 on the right y axis.For example,x = 1:4;y1 = x.^2;y2 = x.^3;plot(x, y1);// y1 on left y axisplot (x. y2);// y2 on right y axisYour assistance will be greatly appreciated.Samuel Enibe___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] animaGIF does not work. Why?

[Heinz Nabielek]

In the Internet, https://convertio.co/ transfers my Scilab series of gif 
figures into a .mov movie. Other format options are also available.
Heinz



> On 14.04.2023, at 21:03, Samuel Gougeon  wrote:
>
> Le 14/04/2023 à 01:00, Heinz Nabielek a écrit :
>> On 13.04.2023, at 20:35, Samuel Gougeon  wrote:
>>> .../...
>>> Installing is not loading (that's intentional).
>>> Just after installing, without restarting Scilab, you need to load the 
>>> module, either with atomsLoad(..), or interactively thanks to the Toolboxes 
>>> console's menu.
>> Silly me. Thanks. Works, of course.
>> In my present coding, the movie is stored as a *.gif file. Graphics programs 
>> see this as stationary figure,
>
>
> It depends on the software. On Windows, for instance the popular
> Irfanview can display the animation. I don't know the MacOS softwares.

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] animaGIF does not work. Why?

[Heinz Nabielek]

On 13.04.2023, at 20:35, Samuel Gougeon  wrote:
>
> Le 13/04/2023 à 04:59, Heinz Nabielek a écrit :
>> Colleagues:
>>
>> animaGIF is loaded
>>
>>> --> atomsInstall('animaGIF');
>>>   animaGIF (1.0) is already installed in the 'allusers' section and 
>>> up-to-date
>>>
>> but does not work
>>
>>> --> idGif = animaGIF(gcf(), outgif, 1,2);
>>> Undefined variable: animaGIF
>>>
>> Why?
>> Heinz
>>
>
>
> Installing is not loading (that's intentional).
> Just after installing, without restarting Scilab, you need to load the 
> module, either with atomsLoad(..), or interactively thanks to the Toolboxes 
> console's menu.

Silly me. Thanks. Works, of course.
In my present coding, the movie is stored as a *.gif file. Graphics programs 
see this as stationary figure, with web browsers I can see the full movie.

What would be the most obious method, to convert it to a *.mov, *.mp4 file (for 
the Mac)?
Heinz


___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] animaGIF does not work. Why?

[Heinz Nabielek]

Colleagues:

animaGIF is loaded
> --> atomsInstall('animaGIF');
>   animaGIF (1.0) is already installed in the 'allusers' section and 
> up-to-date

but does not work
> --> idGif = animaGIF(gcf(), outgif, 1,2);
> Undefined variable: animaGIF

Why?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Generic plot command draws automatically _which_ coloured curves

[Heinz Nabielek]

On 21.03.2023, at 17:17, Jean-Yves Baudais  
wrote:
>
> Hello,
>
> Le 20/03/2023 à 14:40, Samuel Gougeon a écrit :
>> getcolor() tells it to you, in the infobar for the selected color:
>
> > plot((1:10)',rand(10,6),"-*");legend("1","2","3","4","5","6")


Sorry, colleagues, I am still at a loss. At first, I am plotting 7 arrays in a 
sequence of pre-defined colours (that I find nowhwere in the Scilab help 
functions)

A=rand(10,7);
plot((1:10)',A,"-*");
xgrid(color('grey70'));
legend("1","2","3","4","5","6","7");

Now, I would like to redraw array #6 with a thicker line of the same colour as 
in the original colour sequence, thus:
plot((1:10)',A(:,6),"-",'linewidth',3);

How can I specify the 6th colour in the plot command?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Generic plot command draws automatically _which_ coloured curves

[Heinz Nabielek]

 size(E4)=   81.   1.
size(D) =81.   7.
plot (E4,D,'-')

This generic plot command draws automatically coloured curves in blue, green, 
red, etcetera

What are the following colours? I know that 8 is white

Heinz

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scatter plot with fixed marker dimension and overlap analysis

[Heinz Nabielek]

Very interesting display and very interesting piece of coding. I think I will 
be able to use it in one of my display of spherical objects.
Heinz

> On 18.03.2023, at 06:15, Federico Miyara  wrote:
>
> // Image size
> n = 500;
> m = 800;
> // Impact area radius
> R = 20;
> // Number of impacts
> N = 200;
> // Generate test random impact points
> IPx = 1 + floor(m*rand(1,N));
> IPy = 1 + floor(n*rand(1,N));
> X=IPx; [min(X) mean(X) max(X) stdev(X)]
> X=IPy; [min(X) mean(X) max(X) stdev(X)]
> // Plot impacts
> scf(1);
> clf(1);
> plot(IPx, IPy,"o");
> xgrid(color('grey70'));
> // Initialize matrix to contain image
> A = zeros(n,m);
> // Navigate impacts
> for k=1:N
>// Horizontal range for kth impact
>x1 = max(IPx(k) - R, 1);
>x2 = min(IPx(k) + R, m);
>for i=x1:x2
>   // Vertical range for kth impact and ith horizontal pixel
>   y1 = IPy(k) - sqrt(R^2 - (i - IPx(k))^2);
>   y2 = IPy(k) + sqrt(R^2 - (i - IPx(k))^2);
>   y1 = max(y1, 1);
>   y2 = min(y2, n);
>   // Overlap
>   A(y1:y2, i) = A(y1:y2, i) + 1;
>end
> end
> scf(2);
> clf(2);
> gcf().color_map = jetcolormap(64);
> Sgrayplot([1:m],[1:n],A')
> isoview on
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] diagrams that should have identical x-y real picture plotting frames

[Heinz Nabielek]

Thanks a lot.
Heinz



> On 15.03.2023, at 06:55, Federico Miyara  wrote:
>
>
> You could use after each plot the following instruction:
>
> gca().data_bounds = [xmin ymin; xmax, ymax]
>
> where gca() gets the current axis handle and .data_bounds invokes the data 
> bounds field of the handle. xmin ymin are the coordinates of the left bottom 
> corner of the diagram and xmax ymax ar the corresponding ones to the right 
> top corner.
>
> Regards,
>
> Federico Miyara
>
>
> On 14/03/2023 19:04, Heinz Nabielek wrote:
>> Collegues,
>>
>> I am generating dozens of diagrams that should have identical x-y real 
>> picture plotting frames. How to make certain that they all have exactly same 
>> dimension in real space?
>> Heinz
>>
>>
>> ___
>> users mailing list -
>> users@lists.scilab.org
>>
>> Click here to unsubscribe:
>> <mailto:users-unsubscr...@lists.scilab.org>
>> https://lists.scilab.org/mailman/listinfo/users
>>
>> This email and any attachments are intended solely for the use of the 
>> individual or entity to whom it is addressed and may be confidential and/or 
>> privileged.
>>
>> If you are not one of the named recipients or have received this email in 
>> error,
>>
>> (i) you should not read, disclose, or copy it,
>>
>> (ii) please notify sender of your receipt by reply email and delete this 
>> email and all attachments,
>>
>> (iii) Dassault Systèmes does not accept or assume any liability or 
>> responsibility for any use of or reliance on this email.
>>
>>
>> Please be informed that your personal data are processed according to our 
>> data privacy policy as described on our website. Should you have any 
>> questions related to personal data protection, please contact 3DS Data 
>> Protection Officer
>> https://www.3ds.com/privacy-policy/contact/
>>
>>
>>
>>
>>
>
>
>
>
> El software de antivirus Avast ha analizado este correo electrónico en busca 
> de virus.
> www.avast.com
>
>
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
> If you are not one of the named recipients or have received this email in 
> error,
> (i) you should not read, disclose, or copy it,
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
>
> Please be informed that your personal data are processed according to our 
> data privacy policy as described on our website. Should you have any 
> questions related to personal data protection, please contact 3DS Data 
> Protection Officer https://www.3ds.com/privacy-policy/contact/
>
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] diagrams that should have identical x-y real picture plotting frames

[Heinz Nabielek]

Collegues,

I am generating dozens of diagrams that should have identical x-y real picture 
plotting frames. How to make certain that they all have exactly same dimension 
in real space?
Heinz


___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Scilab 2023.0.0 release

[Heinz Nabielek]

Thanks for the advice.
Heinz



> On 10.03.2023, at 23:41, Stéphane Mottelet  wrote:
>
> Hello Heinz,
>
> Yes please use the latest official  release. The version you have was just a 
> continuous build artefact which is now outdated !
>
> S.
>
>> Le 10 mars 2023 à 23:25, Heinz Nabielek  a écrit :
>>
>> Friends,
>>
>> presently very happy with scilab-mr3-190bf2ee.app on my Mac.
>> "6.1.2.1676533517"
>> "Mac OS X 10.15.7"
>>
>> Should there be a good reason to change to Scilab 2023.0.0 - macOS 64 bits 
>> (dmg)?
>> Heinz
>>
>>
>>
>>> On 10.03.2023, at 18:50, COUVERT Vincent  wrote:
>>>
>>> Dear Scilab users,
>>>
>>> We are pleased to announce the release of Scilab 2023.0.0 as a joint effort 
>>> between Scilab contributors and the Scilab team at Dassault Systèmes.
>>>
>>> Scilab 2023.0.0 is the first version released following the new numbering 
>>> and schedule.
>>> It fixes more than a hundred bugs, improve stability and contains new 
>>> features.
>>>
>>> For the complete list of changes and bugs fixed, please take a look 
>>> athttps://help.scilab.org/docs/2023.0.0/en_US/CHANGES.html.
>>>
>>> Download this brand new version at https://www.scilab.org/download/2023.0.0.
>>>
>>> Best regards,
>>> Vincent
>>>
>>>
>>> This email and any attachments are intended solely for the use of the 
>>> individual or entity to whom it is addressed and may be confidential and/or 
>>> privileged.
>>> If you are not one of the named recipients or have received this email in 
>>> error,
>>> (i) you should not read, disclose, or copy it,
>>> (ii) please notify sender of your receipt by reply email and delete this 
>>> email and all attachments,
>>> (iii) Dassault Systèmes does not accept or assume any liability or 
>>> responsibility for any use of or reliance on this email.
>>>
>>> Please be informed that your personal data are processed according to our 
>>> data privacy policy as described on our website. Should you have any 
>>> questions related to personal data protection, please contact 3DS Data 
>>> Protection Officer https://www.3ds.com/privacy-policy/contact/
>>>
>>> ___
>>> users mailing list - users@lists.scilab.org
>>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>>> https://lists.scilab.org/mailman/listinfo/users
>>
>> ___
>> users mailing list - users@lists.scilab.org
>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>> https://lists.scilab.org/mailman/listinfo/users
>> This email and any attachments are intended solely for the use of the 
>> individual or entity to whom it is addressed and may be confidential and/or 
>> privileged.
>>
>> If you are not one of the named recipients or have received this email in 
>> error,
>>
>> (i) you should not read, disclose, or copy it,
>>
>> (ii) please notify sender of your receipt by reply email and delete this 
>> email and all attachments,
>>
>> (iii) Dassault Systèmes does not accept or assume any liability or 
>> responsibility for any use of or reliance on this email.
>>
>>
>> Please be informed that your personal data are processed according to our 
>> data privacy policy as described on our website. Should you have any 
>> questions related to personal data protection, please contact 3DS Data 
>> Protection Officer https://www.3ds.com/privacy-policy/contact/
>>
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
>
> If you are not one of the named recipients or have received this email in 
> error,
>
> (i) you should not read, disclose, or copy it,
>
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
>
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
>
>
> Please be informed that your personal data are processed according to our 
> data pr

Re: [Scilab-users] Scilab 2023.0.0 release

[Heinz Nabielek]

Friends,

presently very happy with scilab-mr3-190bf2ee.app on my Mac.
  "6.1.2.1676533517"
  "Mac OS X 10.15.7"

Should there be a good reason to change to Scilab 2023.0.0 - macOS 64 bits 
(dmg)?
Heinz



> On 10.03.2023, at 18:50, COUVERT Vincent  wrote:
>
> Dear Scilab users,
>
> We are pleased to announce the release of Scilab 2023.0.0 as a joint effort 
> between Scilab contributors and the Scilab team at Dassault Systèmes.
>
> Scilab 2023.0.0 is the first version released following the new numbering and 
> schedule.
> It fixes more than a hundred bugs, improve stability and contains new 
> features.
>
> For the complete list of changes and bugs fixed, please take a look 
> athttps://help.scilab.org/docs/2023.0.0/en_US/CHANGES.html.
>
> Download this brand new version at https://www.scilab.org/download/2023.0.0.
>
> Best regards,
> Vincent
>
>
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
> If you are not one of the named recipients or have received this email in 
> error,
> (i) you should not read, disclose, or copy it,
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
>
> Please be informed that your personal data are processed according to our 
> data privacy policy as described on our website. Should you have any 
> questions related to personal data protection, please contact 3DS Data 
> Protection Officer https://www.3ds.com/privacy-policy/contact/
>
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Arrow on line from (x1, y1) to (x2, y2)

[Heinz Nabielek]

A simple demo:
yy=0.01;
x1=[580 56]; x2=[yy yy];
xarrows(x1,x2,555,6);
Heinz

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



xarrows.pdf
Description: Adobe PDF document


> On 06.03.2023, at 18:53, Samuel Enibe  wrote:
> 
> I need to attach an arrow to a line pointing from (x1, y1) to (x2, y2).
> 
> How can this be done? 
> 
> Your assistance will be greatly appreciated. 
> 
> Samuel Enibe 
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
> If you are not one of the named recipients or have received this email in 
> error,
> (i) you should not read, disclose, or copy it,
> (ii) please notify sender of your receipt by reply email and delete this 
> email and all attachments,
> (iii) Dassault Systèmes does not accept or assume any liability or 
> responsibility for any use of or reliance on this email.
> 
> Please be informed that your personal data are processed according to our 
> data privacy policy as described on our website. Should you have any 
> questions related to personal data protection, please contact 3DS Data 
> Protection Officer https://www.3ds.com/privacy-policy/contact/
> 
> ___
> users mailing list - users@lists.scilab.org
> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
> https://lists.scilab.org/mailman/listinfo/users

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] error in starting scilab from terminal

[Heinz Nabielek]

heinz@Heinzs-iMac ~ % java -version

openjdk version "11.0.15" 2022-04-19
OpenJDK Runtime Environment (build 11.0.15+10)
OpenJDK 64-Bit Server VM (build 11.0.15+10, mixed mode)

  System Version:   macOS 10.15.7 (19H2026)
  Kernel Version:   Darwin 19.6.0




> On 02.03.2023, at 09:57, Stéphane Mottelet  wrote:
>
> What is the output of
>
> java -version
>
> on the command line ?
>
> S.
>
>
> Le 02/03/2023 à 09:45, Wolfgang Engelmann a écrit :
>>
>>
>> Am 02.03.23 um 08:36 schrieb Stéphane Mottelet:
>>> Hello,
>>>
>>> What is your Linux distribution (I guess Ubuntu or Debian ?) and which
>>> Scilab version did you install ?
>>
>> Debian and Scilab 6.1.0, sorry, forgot to mention it
>>
>> Wolfgang
>>
>>>
>>> S.
>>>
>>> Le 02/03/2023 à 07:27, Wolfgang Engelmann a écrit :
>>>> If I try to start scilab from a terminal, it won't start and I get this
>>>> error
>>>>
>>>> .
>>>> g.scilab.modules.types.jar:/usr/share/scilab/modules/xcos/jar/org.scilab.modules.xcos.jar:/usr/share/scilab/modules/localization/jar/org.scilab.modules.localization.jar:/usr/share/scilab/modules/scinotes/jar/org.scilab.modules.scinotes.jar:/usr/share/scilab/modules/gui/jar/org.scilab.modules.gui.jar:/usr/share/scilab/modules/core/jar/org.scilab.modules.core.jar:
>>>>
>>>>
>>>> WARNING: An illegal reflective access operation has occurred
>>>> WARNING: Illegal reflective access by
>>>> org.scilab.modules.jvm.LibraryPath
>>>> (file:/usr/share/scilab/modules/jvm/jar/org.scilab.modules.jvm.jar) to
>>>> method java.lang.ClassLoader.initLibraryPaths()
>>>> WARNING: Please consider reporting this to the maintainers of
>>>> org.scilab.modules.jvm.LibraryPath
>>>> WARNING: Use --illegal-access=warn to enable warnings of further
>>>> illegal
>>>> reflective access operations
>>>> WARNING: All illegal access operations will be denied in a future
>>>> release
>>>>
>>>> What is wrong?
>>>>
>>>> Wolfgang
>>>> ___
>>>> users mailing list - users@lists.scilab.org
>>>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>>>> https://lists.scilab.org/mailman/listinfo/users
>>>> This email and any attachments are intended solely for the use of the
>>>> individual or entity to whom it is addressed and may be confidential
>>>> and/or privileged.
>>>>
>>>> If you are not one of the named recipients or have received this email
>>>> in error,
>>>>
>>>> (i) you should not read, disclose, or copy it,
>>>>
>>>> (ii) please notify sender of your receipt by reply email and delete
>>>> this email and all attachments,
>>>>
>>>> (iii) Dassault Systèmes does not accept or assume any liability or
>>>> responsibility for any use of or reliance on this email.
>>>>
>>>>
>>>> Please be informed that your personal data are processed according to
>>>> our data privacy policy as described on our website. Should you have
>>>> any questions related to personal data protection, please contact 3DS
>>>> Data Protection Officer https://www.3ds.com/privacy-policy/contact/
>>>>
>>> --
>>> Stéphane Mottelet
>>> Ingénieur de recherche
>>> EA 4297 Transformations Intégrées de la Matière Renouvelable
>>> Département Génie des Procédés Industriels
>>> Sorbonne Universités - Université de Technologie de Compiègne
>>> CS 60319, 60203 Compiègne cedex
>>> Tel : +33(0)344234688
>>> http://www.utc.fr/~mottelet
>>>
>>> ___
>>> users mailing list - users@lists.scilab.org
>>> Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
>>> https://lists.scilab.org/mailman/listinfo/users
>>> This email and any attachments are intended solely for the use of the
>>> individual or entity to whom it is addressed and may be confidential
>>> and/or privileged.
>>>
>>> If you are not one of the named recipients or have received this email
>>> in error,
>>>
>>> (i) you should not read, disclose, or copy it,
>>>
>>> (ii) please notify sender of your receipt by reply email and delete this
>>> email and all attachments,
>

Re: [Scilab-users] Can I program a table into a diagram?

[Heinz Nabielek]

Hi Samuel:

mytab = ["Turbine Power" "Turbine hub height" "Turbine radius";
"WT=12 MW" "150 m" "110 m";
"WT= 7 MW" "100 m" "82 m";
"WT= 2 MW" "70 m" "40 m"];
Lat = prettyprint(mytab,,"");   // "" is to remove overall parentheses
Lat = strsubst(Lat, "${", "$\large{");
xstring(0.0015, 12, Lat)

works just fine with me. Absolutely great. And great many thanks.
Heinz

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



Bird survival depends on avoidance, turbine power and wind speed.pdf
Description: Adobe PDF document


> On 27.02.2023, at 20:26, Samuel Gougeon  wrote:
> 
> Hello Heinz,
> 
> Le 25/02/2023 à 00:13, Heinz Nabielek a écrit :
>> Colleagues:
>> 
>> can I program a table like below into a Scilab generated diagram?
> 
> Sending some HTML-ized code to xstring() would have been the more versatile 
> and tunable way to proceed, to get a styled table, with colored backgrounds, 
> etc. Unfortunately, unlike for the GUI, xstring() does not interpret HTML and 
> renders it as is.
> 
> The first actual way might be to use an editable uitable (uicontrol 
> style="table") in a figure.
> Then you will need gui2bitmap to export the figure as an image, including the 
> uitable, and cells won't be filled with individual background colors.
> 
> The second way might be to use a LaTeX table,
> that can be prepared with prettyprint(), then be a bit customized, and 
> finally be sent to xstring():
> 
> mytab = ["Turbine Power" "Turbine hub height" "Turbine radius"
> ;
> 
> "7 MW" "100 m" "82 m";
> "12 MW" "150 m" "110 m"
> ;
> 
> "2 MW" "70 m" "40 m"];
> Lat = prettyprint(mytab,,"");   // ""
>  is to remove overall parentheses
> 
> Lat = strsubst(Lat, "${", "$\large{");
> clf, xstring(0.1, 0.1, Lat
> )
> 
> 
> 
> 
> You can tune the size of the text rendered in the table by changing \large 
> into \Large or \LARGE etc.
> You can add some row and column  frames by switching to a {tabular} object 
> instead of a {matrix} one:
> 
> mytab = ["Turbine Power" "Turbine hub height" "Turbine radius"
> "7 MW" "100 m" "82 m"
> "12 MW" "150 m" "110 m"
> "2 MW" "70 m" "40 m"];
> Lat = prettyprint(mytab,,"");   // "" is to remove overall parentheses
> Lat = strsubst(Lat, "${", "$\large{");
> Lat = strsubst(Lat, "begin{matrix}", "begin{tabular}{|c|c|c|}\hline");
> Lat = strsubst(Lat, "end{matrix}", "end{tabular}");
> Lat = strsubst(Lat, "\cr", "\cr\hline");
> I have not found any trivial direct way to split titles on 2 rows in cells:
> since strsubst() does not handle groups replacement,  the Lat string can be 
> processed by hand:
> In it, you may replace "Turbine Power" with "\substack{Turbine\\Power}", etc. 
> Then you get a string rendered as
> 
> 
> 
> ... that you can pass to xstring().
> 
> Regards
> Samuel
> 
> PS : Since in the figure you get the displayed legend is quite tall, it might 
> be worth to try the mc_legends() macro available in the FileExchange, to get 
> a 2-column legend.
> 
> PS2:  IMHO, prettyprint() would deserve getting more LaTeX-related options, 
> to output a matrix through the {tabular} environment instead of the {matrix} 
> one, to pass some {|c|l|c...} justification and border options, \hline 
> options, process in-cell line breaks, etc.
> 
> This email and any attachments are intended solely for the use of the 
> individual or entity to whom it is addressed and may be confidential and/or 
> privileged.
> If you are not one of the named recipients or have 

[Scilab-users] Fwd: LATEX in Scilab

[Heinz Nabielek]

[cid:B2ACD25C-2E57-4D8E-8F73-7BBC685F550B@home]

Terrific. Thanks. Heinz

Begin forwarded message:

From: Samuel Gougeon mailto:sgoug...@free.fr>>
Subject: Re: [Scilab-users] LATEX in Scilab
Date: 19. February 2023 at 21:13:04 CET
To: users@lists.scilab.org<mailto:users@lists.scilab.org>
Reply-To: sgoug...@free.fr, Users mailing list for Scilab 


Le 19/02/2023 à 20:49, Heinz Nabielek a écrit :
[cid:part1.hrz1sFOc.8Gza9hCA@free.fr]

In a sophisticated Scilab diagram, I would like to position this formula in the 
middle.
I can manage xstring just fine, but can someone help me with LATEX?


Something like this might help:

"$\Phi = 1 - \exp \left(-\left[ \int_{0}^{t}\!\!k\,dt^\prime\,\right]^m\right)$
"


[cid:part2.8q5f60jy.0Hum4RkM@free.fr]

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.
If you are not one of the named recipients or have received this email in error,
(i) you should not read, disclose, or copy it,
(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,
(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.

Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users


This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

[Scilab-users] LATEX in Scilab

[Heinz Nabielek]

[cid:9524BF53-5A8A-41C4-97EF-DC22C96DF05B@home]

In a sophisticated Scilab diagram, I would like to position this formula in the 
middle.
I can manage xstring just fine, but can someone help me with LATEX?
Heinz

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] macOS Scilab users, help needed ! (Stéphane Mottelet)

[Heinz Nabielek]

-->  jre_path
  "/Applications/scilab-mr3-190bf2ee.app/Contents/thirdparty/jre"
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: 
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] macOS Scilab users, help needed !

[Heinz Nabielek]

Works!
Any special properties?
Heinz


iMac 21.5-inch, Late 2013 3,1 GHz Quad-Core Intel Core i7
Darwin Heinzs-iMac.local 19.6.0 Darwin Kernel Version 19.6.0: Tue Jun 21 
21:18:39 PDT 2022; root:xnu-6153.141.66~1/RELEASE_X86_64 x86_64
System Software Overview:
  System Version: macOS 10.15.7 (19H2026)
  Kernel Version: Darwin 19.6.0
  Boot Volume: Macintosh HD
  Boot Mode: Normal
  Computer Name: Heinz’s iMac



[cid:DD5A3167-D1FF-46E1-81C1-5CF2F11F5D2E@home]

--> r = ver()
 column 1

  "Scilab Version: "
  "Operating System: "
  "Java version: "
  "Java runtime information: "
  "Java Virtual Machine information: "
  "Vendor specification: "

 column 2

  "6.1.2.1676533517"
  "Mac OS X 10.15.7"
  "1.8.0_332"
  "OpenJDK Runtime Environment (build 1.8.0_332-b09)"

  "OpenJDK 64-Bit Server VM (build 25.332-b09, mixed mode)"

  "Oracle Corporation"

On 16.02.2023, at 08:59, Stéphane Mottelet 
mailto:stephane.motte...@utc.fr>> wrote:

Hello macOS Scilab users,

Can you test the latest CI (continuous integration) buid of macOS
version of Scilab ? The dmg archive is there :

https://gitlab.com/mottelet/scilab/-/jobs/3780913723/artifacts/file/scilab-mr3-190bf2ee-unsigned-x86_64.dmg

As the archive is not signed you will have to allow the execution of
Scilab in System Preferences/Security and Privacy. Please report if
Scilab opens or not and mention your macOS version.

Thanks in advance

--
Stéphane Mottelet
Ingénieur de recherche
EA 4297 Transformations Intégrées de la Matière Renouvelable
Département Génie des Procédés Industriels
Sorbonne Universités - Université de Technologie de Compiègne
CS 60319, 60203 Compiègne cedex
Tel : +33(0)344234688
http://www.utc.fr/~mottelet

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/



This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

[Scilab-users] Poisson confidence range of small numbers

[Heinz Nabielek]

Friends:

no problem this time, but the solution with thanks to the SciLab community:

Statistical small numbers are Poisson distributed, but the inverse of the 
cumulative Poisson distribution is not an easy affair and we have to remember 
the Chi-Square distribution.

Example: the 95% confidence range of numbers 1 to 5 is given by:

conf=0.95;
n=(1:5)';
P = (1 - (1-conf)/2)*ones(n); Q=1-P;
// 95% confidence range of small numbers
[cdfchi("X",2*n,Q,P)/2  n  cdfchi("X",2*(n+1),P,Q)/2 ]

   0.0253178   1.   5.5716434
   0.2422093   2.   7.2246877
   0.6186721   3.   8.7672731
   1.0898654   4.   10.241589
   1.6234864   5.   11.668332

Thanks for all the help
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] 1. Re: cdfchi=cumulative distribution function chi-square distribution (Dang Ngoc Chan, Christophe)

[Heinz Nabielek]

On 10.02.2023, at 09:22, Dang Ngoc Chan, Christophe  
wrote:
>
> Hello Heinz,
>
>> De : Heinz Nabielek 
>> Envoyé : jeudi 9 février 2023 13:36
>>
>> Result is then
>> [...]
>>
>> However, EXCEL conf=0.95 gives me:
>> [...]
>
> Sorry, I don't know particularly the chi-square distribution, I suppose you 
> don't ask the same question to both softwares.
>
> If you use for XL
>
> Xmax=CHIINV((1-conf)/2,2*(n+1))/2
>
> Then for Scilab, you should use
>
> N = 1:5
>
> P = (1 - 0.5*0.05)*ones(n)
>
> Q = 1 - P
>
> [Xmax] = 0.5*cdfchi("X", 2*(n+1), P, Q)
>
> Which the gives you the same result

n = 1:5
P = (1 - 0.5*0.05)*ones(n)
Q = 1 - P
Upper95=0.5*cdfchi("X", 2*(n+1), P, Q)'
   5.5716434
   7.2246877
   8.7672731
   10.241589
   11.668332

Perfect agreement! I should have thought of it myself..
But: why does it have to be so difficult for an experimental physicist when 
trying to assess the uncertainties of low numbers?
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] 1. Re: cdfchi=cumulative distribution function chi-square distribution (Dang Ngoc Chan, Christophe)

[Heinz Nabielek]

> Hello Heinz,
>
>> De : Heinz Nabielek
>> Envoy? : mardi 7 f?vrier 2023 13:07
>>
>> cdfchi = cumulative distribution function chi-square distribution.
>> I want to generate the 95% confidence limit of small numbers like below
>
> According to https://help.scilab.org/doc/6.0.0/en_US/cdfchi.html
> I think the following does what you want:
>
> n = 1:5
> P = 0.95*ones(n)
> Q = 1 - P
> [Xmax] = cdfchi("X", n, P, Q)
> [Xmin] = cdfchi("X", n, Q, P)

Great many thanks. Result is then

--> [Xmin' n' Xmax']
   0.0039321   1.   3.8414588
   0.1025866   2.   5.9914645
   0.3518463   3.   7.8147279
   0.7107234.   9.487729
   1.1454762   5.   11.070498

However, EXCEL conf=0.95 gives me:
0.02531781  1   5.57164339
0.24220928  2   7.22468767
0.61867212  3   8.76727307
1.08986537  4   10.2415887
1.62348639  5   11.6683321

for Xmin=CHIINV(1-(1-conf)/2,2*n)/2
for Xmax=CHIINV((1-conf)/2,2*(n+1))/2

Also, for large numbers, upper 95% limit should be approx. n+1.96√n
n   n+1.96√n  cdfchi
50  63.86   67.50
100 119.60  124.34
150 174.00  179.58

Where do I go wrong with the  chi-square distribution?
Heinz
  __
Dr Heinz Nabielek
Schüttelstrasse 77A/11
1020 Wien, Österreich
Tel +43 1 276 56 13
cell +43 677 616 349 22


___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] cdfchi=cumulative distribution function chi-square distribution

[Heinz Nabielek]

Dear colleagues:

I have problems computing the inverse-chi-squared distribution in SciLab with

cdfchi = cumulative distribution function chi-square distribution.

I want to generate the 95% confidence limit of small numbers like below


number nlower limit upper limit
0   0.003.69
1   0.035.57
2   0.247.22
3   0.628.77
4   1.0910.24
5   1.6211.67

and this is easily calculated in EXCEL with conf=0.95 and

=CHIINV(1-(1-conf)/2,2*n)/2 for the lower limit; and

=CHIINV((1-conf)/2,2*(n+1))/2 .

But how do I do this in SciLab?
Heinz

This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Problems with Scilab routine "conv"

[Heinz Nabielek]

This is my latest code version: the 'convoluted secondary failure fraction' is 
nicely below primary failure,  but seems too low.
Heinz


m=2; // Weibull modulus in mechanism  #1
k=1E-7;  // corrosion rate(s-1) in mechanism  #1
n=500;   // number of time steps in hourly intervals
t=(1:n)';// timesteps in hours
ts= 3600*t;  // timesteps in seconds
PHI= 1 - exp(-((k*ts)^m)); // failure fraction in mechanism #1
deriv=m*k*((k*ts)^(m-1)) .* exp(-((k*ts)^m)); //derivative of PHI
phi=3600*deriv;  //derivative of PHI integrated over an hour
f= 1 - exp(-ts/3E5); // "pure" failure fraction in mechanism #2
F=(convol(phi,f)/n)(1:n);   //"convoluted" failure fraction mechanism #2
plot("nl",t, [PHI phi f F'],'-');xgrid;
legend('PHI primary failure','phi primary failure rate','pure secondary failure 
fraction','convoluted secondary failure fraction',4);
xlabel('time (hours)');
ylabel('failure fraction/ failure rate');



> On 03.02.2023, at 11:24, Heinz Nabielek  wrote:
>
> On 03.02.2023, at 11:13, Stéphane Mottelet  wrote:
>>
>> Thanks for the code.
>>
>> Just a remark on the notations, you should write :
>>
>> F(T)  = Int_{0}^{T}  PHI(t) . f(T-t) . dt
>>
>> i.e. not F(t) since t is mute.
>>
>> However, you should pay attention to the delay notion associated with 
>> convolution and the relationships between discrete convolution and 
>> continuous convolution. I am not sure that the output of conv matches with a 
>> given discretization of the integral above. Maybe rectangle method, but I am 
>> not sure at all. Anyway, you should have F(0)=0 which does not seem to be 
>> the case in your graph.
>
> F(t), of course.
>
> But the formula is fundamentally wrong and one should have seen it already 
> from a dimensional consideration.
>
> Correct should be F(t)  = Int_{0}^{T}  (d PHI(t)/dt) . f(T-t) . dt
>
> Hope you can help between conv and convol and possible something else.
> Heinz
>
>
>
>
>
>
>> Dear SciLab Friends:
>>
>> I have an object consisting of many (~10,000) small components that can fail 
>> in a statistical way during long-term operation at extreme conditions.
>>
>> My primary failure model is described by PHI(t) going monotonically from 
>> zero to one at times from t=0 to T. In the computer, this is realized in n 
>> timesteps.
>>
>> Mechanism 2: There is a secondary failure mechanism that starts only when a 
>> component has previously failed by mechanism 1 and occurs after some delay. 
>> The delay function is f(t) and final expression for the evolution of the 
>> mechanism 2 failure is given by:
>>
>> F(t)  = Int_{0}^{T}  PHI(t) . f(T-t) . dt
>>
>> a classical convolution Faltungsintegral.
>>
>> In Scilab, I write
>>
>> F= conv(PHI, f, "same")/ n;
>>
>> and the resulting function F looks reasonable for most of the time. Under 
>> some conditions, however, I can have
>>
>> F>PHI
>>
>> at early stages, but this is physically impossible. Because it comes later, 
>> F must always be below PHI.
>> What do I do wrong?
>> Heinz
>>
>> _
>> PS: massively simplified code below..
>>
>> k=1E-7;  // corrosion rate(s-1) in mechanism #1
>> n=100;   // number of timesteps
>> t=(1:n)'; // time in days
>> ts= 3600*t;  // time in seconds
>> PHI= 1 - exp(-((k*ts)^2));   // failure fraction in mechanism #1
>> plot("nl", t, PHI, 'r--');
>> f= 1 - exp(-ts/3E5); // "pure" failure fraction in mechanism #2
>> F=conv(PHI,f,"same")/length(t); //"convoluted" failure fraction mechanism #2
>> plot(t,F,'b--');
>> legend('failure fraction #1','subsequent failure fraction #2',4);
>> xlabel('time (hours)');
>> ylabel('failure fraction');
>> title('Component failure evolution #1 is followed by #2');

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] Problems with Scilab routine "conv"

[Heinz Nabielek]

On 03.02.2023, at 11:13, Stéphane Mottelet  wrote:
>
> Thanks for the code.
>
> Just a remark on the notations, you should write :
>
> F(T)  = Int_{0}^{T}  PHI(t) . f(T-t) . dt
>
> i.e. not F(t) since t is mute.
>
> However, you should pay attention to the delay notion associated with 
> convolution and the relationships between discrete convolution and continuous 
> convolution. I am not sure that the output of conv matches with a given 
> discretization of the integral above. Maybe rectangle method, but I am not 
> sure at all. Anyway, you should have F(0)=0 which does not seem to be the 
> case in your graph.

F(t), of course.

But the formula is fundamentally wrong and one should have seen it already from 
a dimensional consideration.

Correct should be F(t)  = Int_{0}^{T}  (d PHI(t)/dt) . f(T-t) . dt

Hope you can help between conv and convol and possible something else.
Heinz






> Dear SciLab Friends:
>
> I have an object consisting of many (~10,000) small components that can fail 
> in a statistical way during long-term operation at extreme conditions.
>
> My primary failure model is described by PHI(t) going monotonically from zero 
> to one at times from t=0 to T. In the computer, this is realized in n 
> timesteps.
>
> Mechanism 2: There is a secondary failure mechanism that starts only when a 
> component has previously failed by mechanism 1 and occurs after some delay. 
> The delay function is f(t) and final expression for the evolution of the 
> mechanism 2 failure is given by:
>
> F(t)  = Int_{0}^{T}  PHI(t) . f(T-t) . dt
>
> a classical convolution Faltungsintegral.
>
> In Scilab, I write
>
> F= conv(PHI, f, "same")/ n;
>
> and the resulting function F looks reasonable for most of the time. Under 
> some conditions, however, I can have
>
> F>PHI
>
> at early stages, but this is physically impossible. Because it comes later, F 
> must always be below PHI.
> What do I do wrong?
> Heinz
>
> _
> PS: massively simplified code below..
>
> k=1E-7;  // corrosion rate(s-1) in mechanism #1
> n=100;   // number of timesteps
> t=(1:n)'; // time in days
> ts= 3600*t;  // time in seconds
> PHI= 1 - exp(-((k*ts)^2));   // failure fraction in mechanism #1
> plot("nl", t, PHI, 'r--');
> f= 1 - exp(-ts/3E5); // "pure" failure fraction in mechanism #2
> F=conv(PHI,f,"same")/length(t); //"convoluted" failure fraction mechanism #2
> plot(t,F,'b--');
> legend('failure fraction #1','subsequent failure fraction #2',4);
> xlabel('time (hours)');
> ylabel('failure fraction');
> title('Component failure evolution #1 is followed by #2');
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] Problems with Scilab routine "conv"

[Heinz Nabielek]

Dear SciLab Friends:

I have an object consisting of many (~10,000) small components that can fail in 
a statistical way during long-term operation at extreme conditions.

My primary failure model is described by PHI(t) going monotonically from zero 
to one at times from t=0 to T. In the computer, this is realized in n timesteps.

Mechanism 2: There is a secondary failure mechanism that starts only when a 
component has previously failed by mechanism 1 and occurs after some delay. The 
delay function is f(t) and final expression for the evolution of the mechanism 
2 failure is given by:

F(t)  = Int_{0}^{T}  PHI(t) . f(T-t) . dt

a classical convolution Faltungsintegral.

In Scilab, I write

F= conv(PHI, f, "same")/ n;

and the resulting function F looks reasonable for most of the time. Under some 
conditions, however, I can have

F>PHI

at early stages, but this is physically impossible. Because it comes later, F 
must always be below PHI.
What do I do wrong?
Heinz

_
PS: massively simplified code below..

k=1E-7;  // corrosion rate(s-1) in mechanism #1
n=100;   // number of timesteps
t=(1:n)'; // time in days
ts= 3600*t;  // time in seconds
PHI= 1 - exp(-((k*ts)^2));   // failure fraction in mechanism #1
plot("nl", t, PHI, 'r--');
f= 1 - exp(-ts/3E5); // "pure" failure fraction in mechanism #2
F=conv(PHI,f,"same")/length(t); //"convoluted" failure fraction mechanism #2
plot(t,F,'b--');
legend('failure fraction #1','subsequent failure fraction #2',4);
xlabel('time (hours)');
ylabel('failure fraction');
title('Component failure evolution #1 is followed by #2');
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

[Scilab-users] convol

[Heinz Nabielek]

Friends:

I have an object consisting of many (~10,000) small components that can fail in 
a statistical way during long-term operation at extreme conditions.

My primary failure model is described by PHI(t) going monotonically from zero 
to one at times from t=0 to T. In the computer, this is realized in n timesteps.

Mechanism 2: There is a secondary failure mechanism that starts only when a 
component has failed by mechanism 1 and occurs after some delay. The delay 
function is f(t) and final expression for the evolution of the mechanism 2 
failure is given by:

F(t)  = Int_{0}^{T}  PHI(t) . f(T-t) . dt

a classical convolution Faltungsintegral.

In Scilab, I write

F= conv(PHI, f, "same")/ n;

and the resulting function F looks reasonable for most of the time. Under some 
conditions, however, I can have

F>PHI

at early stages, but this is physically impossible. Because it comes later, F 
must always be below PHI.
What do I do wrong?
Heinz


___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Re: [Scilab-users] (no subject)

[Heinz Nabielek]

On 10.01.2023, at 18:05, Federico Miyara  wrote:
> 
> 
> My knowledge of Matlab is freezed as of 2017, when I abandoned Matlab and 
> embraced Scilab, so might be some of my comments are outdated  
> 
> It is free, open source software
> 
> No licenses that eventually expire
> 
> Strings are better handled. In Matlab you can check for equality only if 
> their lengths are the same since they are just vectors of codes. In Scilab a 
> string is a true type for supporting many operations in a more intuitive way.
> 
> Variable names support extended characters
> 
> You can extract components from the result of any function: y = fft(x)(1:5)
> 
> The same applies to any field of structure
> 
> It is mostly compatible with Matlab and you can automatically convert Matlab 
> scripts to Scilab scripts, which allows to reuse Matlab code (it is not 
> perfect and manual correction is needed)
> 
> Insanely powerful control of graphics.
> 
> Impressive community to help with any doubts.



150% agreed and confirmed
Heinz
___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] x-axis of an Arrhenius diagram with °C scale

[Heinz Nabielek]

This time no question, but a solution: x-axis of an Arrhenius diagram with °C 
scale.
Heinz




T=[(500:100:1000)'; 1200;1400;1600;2000;2500]
TT=flipdim(1e4./(T+273.15),1)
a=gca();
a.x_ticks=tlist(["ticks","locations","labels"],.. 
 [TT],.. 
 ["${2500}$";.. 
 "${2000}$";.. 
 "${1600}$";.. 
 "${1400}$";.. 
 "${1200}$";.. 
 "${1000}$";.. 
 "${900}$";.. 
 "${800}$";.. 
 "${700}$";.. 
 "${600}$";.. 
 "${500}$";..
]);___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] LATEX: easy help?

[Heinz Nabielek]

Great. Thanks.

Is \: and \, equivalent for a blank?
Heinz



> On 27.12.2022, at 22:15, Stéphane Mottelet  wrote:
> 
> Hi,
> 
> The following should render the labels as you want:
> 
> xlabel('$\mbox{inverse temperature}\,[10^4/T(K)]$');
>  
> 
> ylabel('$\log_{10} D\,[m^2/s]$')
> ;
> 
> S.
> 
> 
> Le 24/12/2022 à 05:16, Heinz Nabielek a écrit :
>> 
>> xlabel('inverse temperature [1E4/T(K)]'); 
>> ylabel('log10 D [m2/s]'); 

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] LATEX: easy help?

[Heinz Nabielek]

I need subscripts and superscripts  log10D [m2/s]  and 104/T(K) labelling 
Scilab diagrams, but I am lost on LATEX: easy help?

Heinz

xlabel('inverse temperature [1E4/T(K)]'); 
ylabel('log10 D [m2/s]'); ___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Problems with cdfbet

[Heinz Nabielek]

N = number of particles investigated. Sorry for typo..

> On 05.11.2022, at 01:03, Heinz Nabielek  wrote:
> 
> 
> N = number of particles investiged
> 
> 

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Problems with cdfbet

[Heinz Nabielek]

On 04.11.2022, at 09:46, Dang Ngoc Chan, Christophe <christophe.d...@sidel.com> wrote:Hello Heinz,De : users <users-boun...@lists.scilab.org> De la part de Heinz NabielekEnvoyé : mardi 1 novembre 2022 22:02I need to compute the inverse of the cumulative beta function, but I just cannot handle cdfbet.From what I undestand inhttps://help.scilab.org/docs/6.1.1/en_US/cdfbet.htmlthe notation for the beta function in Scilab is :P = \int_{0}^{X} t^{A-1} * (1-t)^{B-1} dtWhat you want is X knowing [P, A, B] right ?The syntax seems to be[X,Y]=cdfbet("XY",A,B,P,Q)With Q = 1 - P.With your examples:--> [X,Y]=cdfbet("XY",1 ,10 ,0.95 ,0.05)X  = 0.2588656Y  = 0.7411344--> [X,Y]=cdfbet("XY",1 ,100 ,0.95 ,0.05)X  = 0.0295130Y  = 0.9704870RegardsChristophe Dang Ngoc ChanDear Christophe et Stéphane,great many thanks for the patient explanations.This way I was able to compute acceptance limits in the large-scale manufacture of small objects that are difficult to find in textbooks.GreetingsHeinz   ... I am certain that you can generate a much more efficient Scilab code.N=10^(3:9);
def=[300 100 30 10 3 1 0];
for i=1:7;
NN=N(i);
for j=1:7;
defects=def(j);
A(i,j)=cdfbet("XY", defects+1, NN+1-defects, .95,.05);
end;
plot2d(1e9,1,logflag='ll');xgrid();
plot(N', A,'--','LineWidth',3);
a=gca();a.font_size=3;
legend ('n=','300 defects','100 defects','30 defects','10 defects','3 defects','1 defect','0 defects');
xlabel('N = number of particles investiged','fontsize',3);
ylabel('one-sided upper 95% defect fraction','fontsize',3);
title ('Acceptance limit of defect fraction n/N at 95% confidence','fontsize',4);

Acceptance limit of defect fraction n over N at 95% confidence.pdf
Description: Adobe PDF document
__Dr Heinz NabielekSchüttelstrasse 77A/11A-1020 Wien, ÖsterreichTel +43 1 276 56 13cell +43 677 616 349 22heinznabie...@me.com___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] Problems with cdfbet

[Heinz Nabielek]

I need to compute the inverse of the cumulative beta function, but I just 
cannot handle cdfbet.

Simple calculations in EXCEL would be
BETA.INV(0.95,1,10) = 0.259 
BETA.INV(0.95,1,100) = 0.0295...

How do I express this in cdfbet?
Heinz

___
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab 6.1.1 release

[Heinz Nabielek]

(base) heinz@Heinzs-iMac ~ %  
/Applications/scilab-6.1.1.app/Contents/MacOS/scilab
JavaVM: Failed to load JVM: 
/Library/Java/JavaVirtualMachines/adoptopenjdk-8.jre/Contents/Home/libserver.dylib
JavaVM FATAL: Failed to load the jvm library.
JavaVM: Failed to load JVM: 
/Library/Java/JavaVirtualMachines/adoptopenjdk-8.jre/Contents/Home/libserver.dylib
JavaVM FATAL: Failed to load the jvm library.
Error in the creation of the Java VM: Unknown JNI error
Options:
0: -Djava.library.path=
1: 
-Djava.class.path=/Applications/scilab-6.1.1.app/Contents/share/scilab/modules/jvm/jar/org.scilab.modules.jvm.jar
2: 
-Djava.util.logging.config.file=/Applications/scilab-6.1.1.app/Contents/share/scilab/etc/logging.properties
3: -Djava.compiler=JIT
4: -Xmx256m
5: -Xrs
6: -Dapple.laf.useScreenMenuBar=true
7: -Djogamp.gluegen.UseTempJarCache=false

Scilab cannot open JVM library.




> On 20.07.2021, at 17:13, Stéphane Mottelet  wrote:
> 
>  /Applications/scilab-6.1.1.app/Contents/MacOS/scilab

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab 6.1.1 release

[Heinz Nabielek]

On 20.07.2021, at 16:55, Stéphane Mottelet  wrote:
> 
> Please remove jdk-13.0.1.jdk
> 
> sudo rm -rf /Library/Java/JavaVirtualMachines/jdk-13.0.1.jdk

DONE
> 
> Then retry to launch Scilab.

NO SCILAB START

> If Scilab does not start, remove the other one:
> 
> sudo rm -rf /Library/Java/JavaVirtualMachines/adoptopenjdk-8.jdk

DONE

> 
> Scilab should start and propose to download and install the adoptopenjdk JRE. 
> Say Yes.


I DID.

STILL NO SCILAB START
Just a short flash.
Heinz

(base) heinz@Heinzs-iMac ~ % ls /Library/Java/JavaVirtualMachines
adoptopenjdk-8.jdk  jdk-13.0.1.jdk
(base) heinz@Heinzs-iMac ~ % sudo rm -rf 
/Library/Java/JavaVirtualMachines/jdk-13.0.1.jdk

(base) heinz@Heinzs-iMac ~ % sudo rm -rf 
/Library/Java/JavaVirtualMachines/adoptopenjdk-8.jdk
(base) heinz@Heinzs-iMac ~ % ls /Library/Java/JavaVirtualMachines   
 

(base) heinz@Heinzs-iMac ~ % ls /Library/Java/JavaVirtualMachines
adoptopenjdk-8.jre
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab 6.1.1 release

[Heinz Nabielek]

(base) heinz@Heinzs-iMac ~ % ls /Library/Java/JavaVirtualMachines
adoptopenjdk-8.jdk  jdk-13.0.1.jdk

> On 20.07.2021, at 16:03, Stéphane Mottelet  wrote:
> 
> ls /Library/Java/JavaVirtualMachines

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab 6.1.1 release

[Heinz Nabielek]

yes

….  heinz.  ….

> On 20.07.2021, at 15:59, Stéphane Mottelet  wrote:
> 
> Did you have scilab-branch-6.1 working before ?
> 
>> Le 20/07/2021 à 15:56, Heinz Nabielek a écrit :
>> I downloaded scilab-6.1.1-x86_64.dmg onto my iMac with macOS 10.15.7 
>> (19H1217), but the beast does not start. I had this before, but forgot what 
>> needs to be done...
>> Heinz
>> 
>> Previously running scilab-branch-6.1.app
>> 
>>>> Dear Scilab users,
>>>> 
>>>> We are pleased to announce the release of Scilab 6.1.1 as a joint effort 
>>>> between Scilab contributors and the Scilab team at ESI Group.
>>>> 
>>>> Scilab 6.1.1 is the first revision of the 6.1 development branch and is 
>>>> compatible with the 6.1.0 release. It fixes 156 bugs and implements 
>>>> missing features from the 6.1.0 version especially:
>>>> 
>>>> * Improve readability on plot's title and displayed values
>>>> * Fix crashes and slow behaviors on functions
>>>> * Add missing arguments and behaviors to existing functions
>>>> * Reword and add examples on many help pages
>>>> 
>>>> For the complete list of changes and bugs fixed, please take a look at the 
>>>> CHANGES file https://help.scilab.org/CHANGES
>>>> 
>>>> Download this brand new version at https://www.scilab.org/download/6.1.1
>>>> 
>>>> --
>>>> Clément DAVID on behalf of the Scilab team.
>> ___
>> users mailing list
>> users@lists.scilab.org
>> http://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
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab 6.1.1 release

[Heinz Nabielek]

I downloaded scilab-6.1.1-x86_64.dmg onto my iMac with macOS 10.15.7 (19H1217), 
but the beast does not start. I had this before, but forgot what needs to be 
done...
Heinz

Previously running scilab-branch-6.1.app

>> Dear Scilab users,
>> 
>> We are pleased to announce the release of Scilab 6.1.1 as a joint effort 
>> between Scilab contributors and the Scilab team at ESI Group.
>> 
>> Scilab 6.1.1 is the first revision of the 6.1 development branch and is 
>> compatible with the 6.1.0 release. It fixes 156 bugs and implements missing 
>> features from the 6.1.0 version especially:
>> 
>> * Improve readability on plot's title and displayed values
>> * Fix crashes and slow behaviors on functions
>> * Add missing arguments and behaviors to existing functions
>> * Reword and add examples on many help pages
>> 
>> For the complete list of changes and bugs fixed, please take a look at the 
>> CHANGES file https://help.scilab.org/CHANGES
>> 
>> Download this brand new version at https://www.scilab.org/download/6.1.1
>> 
>> --
>> Clément DAVID on behalf of the Scilab team.

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] How can I make a movie that shows how the sphere rotates around the z-axis?

[Heinz Nabielek]

Thanks for all the help from the group and saving the animation works well.
However [had somebody noticed?], there was an error in the statistical 
formulation to generate a uniform random point distribution within the sphere. 
The more correct (but perhaps clumsy) version is below.
Heinz


R=23;
N=14500;
nmax=1.1*6*N/%pi;
// Random points in cube 2R*2R*2R
x=grand(nmax,1,'unf',-R,R);
y=grand(nmax,1,'unf',-R,R);
z=grand(nmax,1,'unf',-R,R);
R2=R*R;
// Exclude points outside of sphere with radius R
T=(x^2+y^2+z^2)>R2;
M=[x y z];
M(T,:)=[];
[r c]=size(M);
M(((N+1):r), :)=[];
x=M(:,1); y=M(:,2); z=M(:,3); 
f = figure();
f.background = 8;
s = linspace(5, 1, N);
scatter3d(x,y,z,msizes=1,s);
a = gca();
outgif = "Sphere01.gif";
idGif = animaGIF(gcf(), outgif, ,2);
for i = 20:60;
a.rotation_angles = [45, i];
idGif = animaGIF(gcf(), idGif); // Adds the current figure to the GIF stream
end;
animaGIF(idGif);// Closes the GIF stream




> On 02.07.2021, at 00:15, Heinz Nabielek  wrote:
> 
> Colleagues:
> 
> the code below generates 14,500 random points in a spherical volume.
> How can I make a movie that shows how the sphere rotates around the z-axis? 
> And how do I store the animation? animaGIF does not exist on my Scilab 6.1 in 
> my macOS 10.15.7
> 
> Heinz
> 
> 
> R=23;
> N=14500;
> r=grand(N,1,'unf',0,R);
> phi=grand(N,1,'unf',0,2*%pi);
> theta=grand(N,1,'unf',0,%pi);
> x=r.*cos(phi).*sin(theta);
> y=r.*sin(phi).*sin(theta);
> z=r.*cos(theta);
> scatter3d(x,y,z,msizes=4);
> gca().rotation_angles = [60, 60];

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] How can I make a movie that shows how the sphere rotates around the z-axis?

[Heinz Nabielek]

No. How to do it?

I have downloaded but then what?


/Users/heinz/Downloads/animaGIF_1.0_61-bin
-rw-rw-r--@1 heinz  staff1386 Feb 28  2020 DESCRIPTION
-rw-rw-r--@1 heinz  staff  407469 Jul 22  2019 animaGIF_1.0.pdf
-rw-rw-r--@1 heinz  staff 396 Jul 21  2019 changelog.txt
drwxrwxr-x@3 heinz  staff  96 Feb 28  2020 etc
drwxrwxr-x@3 heinz  staff  96 Feb 28  2020 jar
-rw-rw-r--@1 heinz  staff   22677 Jul 21  2019 license.txt
-rw-rw-r--@1 heinz  staff 594 Feb 28  2020 loader.sce
drwxrwxr-x@6 heinz  staff 192 Feb 28  2020 locales
drwxrwxr-x@   10 heinz  staff 320 Feb 28  2020 macros
-rw-rw-r--@1 heinz  staff 539 Jul 21  2019 readme.txt
drwxrwxr-x@5 heinz  staff 160 Feb 28  2020 tests
-rw-rw-r--@1 heinz  staff 990 Feb 28  2020 unloader.sce

> On 02.07.2021, at 23:25, Stéphane Mottelet  wrote:
> 
> Heinz did you manage the offline install of animaGIF ?
> 
> Le 02/07/2021 à 23:03, Heinz Nabielek a écrit :
>> Thanks. The rotating graphics works well, but on my Mac with Scilab 6.1,  
>> already the first line "scicv_Init()" bumps...
>> What to do?
>> Heinz
>> 
>> 
>> 
>>> On 02.07.2021, at 10:50, P M  wrote:
>>> 
>>>  the rotation can be realised by changing   "rotation_angles"  in a for 
>>> loop
>>> 
>>> .. creating the movie...maybe use scicv?
>>> 
>>> in a for loop
>>> - set rotation angle
>>> - display the graph
>>> - save the actual graph as an image
>>> - reload the image
>>> - add image to frame
>>> 
>>> 
>>> see code below:
>>> Only draw back here:   the avi-file could not be played afterwards..and had 
>>> only 6kB in size.
>>> 
>>> I guess, this is because ffmpeg is missing on my PC, but I am pretty sure 
>>> that some scicv expert could solve this issue.
>>> 
>>> BR
>>> Philipp
>>> 
>>> scicv_Init();
>>> 
>>> 
>>> 
>>> R=23;
>>> N=14500;
>>>  
>>> r=grand(N,1,'unf',0,R);
>>> phi=grand(N,1,'unf',0,2*%pi);
>>> theta=grand(N,1,'unf',0,%pi);
>>> x=r.*cos(phi).*sin(theta);
>>> y=r.*sin(phi).*sin(theta);
>>> z=r.*cos(theta);
>>> 
>>> 
>>> 
>>> f = figure();
>>> f.background = 8;
>>> scatter3d(x,y,z,msizes=4);
>>> a = gca();
>>> 
>>> 
>>> 
>>> // create the avi file
>>> videoWriter = new_VideoWriter("F:\testAvi.avi", CV_FOURCC('M', 'P', '4', 
>>> '2'), 25, [400, 400]);
>>> 
>>> 
>>> 
>>> for i = 0:360
>>> 
>>> a.rotation_angles = [60, i];
>>> 
>>>  // save the graph as image
>>> 
>>> xs2png(f.figure_id, 'F:\testImg.bmp');
>>> 
>>>  // read image
>>> 
>>> img = imread('F:\testImg.png', CV_LOAD_IMAGE_GRAYSCALE);
>>> 
>>> 
>>> // add image as frame to avi
>>> 
>>> VideoWriter_write(videoWriter, img);
>>> 
>>> 
>>> 
>>> end
>>> 
>>> 
>>> 
>>> // free avi from memory
>>> delete_VideoWriter(videoWriter);
>>> 
>>> 
>>> 
>>> 
>>> 
>>> Am Fr., 2. Juli 2021 um 00:22 Uhr schrieb Heinz Nabielek 
>>> :
>>> Colleagues:
>>> 
>>> the code below generates 14,500 random points in a spherical volume.
>>> How can I make a movie that shows how the sphere rotates around the z-axis? 
>>> And how do I store the animation? animaGIF does not exist on my Scilab 6.1 
>>> in my macOS 10.15.7
>>> 
>>> Heinz
>>> 
>>> 
>>> R=23;
>>> N=14500;
>>> r=grand(N,1,'unf',0,R);
>>> phi=grand(N,1,'unf',0,2*%pi);
>>> theta=grand(N,1,'unf',0,%pi);
>>> x=r.*cos(phi).*sin(theta);
>>> y=r.*sin(phi).*sin(theta);
>>> z=r.*cos(theta);
>>> scatter3d(x,y,z,msizes=4);
>>> gca().rotation_angles = [60, 60];
>>> 
>>> 
>>> __
>>> Dr Heinz Nabielek
>>> Schüttelstrasse 77A/11
>>> A-1020 Wien, Österreich
>>> Tel +43 1 276 56 13
>>> cell +43 677 616 349 22
>>> heinznabie...@me.com
>>> ___
>>> users mailing list
>>> users@lists.scilab.org
>>> http://lists.scilab.org/mailman/listinfo/users
>>> ___
>>> users mailing list
>>> users@lists.scilab.org
>>> http://lists.scilab.org/mailman/listinfo/users
>> ___
>> users mailing list
>> users@lists.scilab.org
>> http://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
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] How can I make a movie that shows how the sphere rotates around the z-axis?

[Heinz Nabielek]

Thanks. The rotating graphics works well, but on my Mac with Scilab 6.1,  
already the first line "scicv_Init()" bumps...
What to do?
Heinz



> On 02.07.2021, at 10:50, P M  wrote:
> 
>  the rotation can be realised by changing   "rotation_angles"  in a for 
> loop
> 
> .. creating the movie...maybe use scicv?
> 
> in a for loop
> - set rotation angle
> - display the graph
> - save the actual graph as an image
> - reload the image
> - add image to frame
> 
> 
> see code below:
> Only draw back here:   the avi-file could not be played afterwards..and had 
> only 6kB in size.
> 
> I guess, this is because ffmpeg is missing on my PC, but I am pretty sure 
> that some scicv expert could solve this issue.
> 
> BR
> Philipp
> 
> scicv_Init();
> 
> 
> 
> R=23;
> N=14500;
>  
> 
> r=grand(N,1,'unf',0,R);
> phi=grand(N,1,'unf',0,2*%pi);
> theta=grand(N,1,'unf',0,%pi);
> x=r.*cos(phi).*sin(theta);
> y=r.*sin(phi).*sin(theta);
> z=r.*cos(theta);
> 
> 
> 
> f = figure();
> f.background = 8;
> scatter3d(x,y,z,msizes=4);
> a = gca();
> 
> 
> 
> // create the avi file
> videoWriter = new_VideoWriter("F:\testAvi.avi", CV_FOURCC('M', 'P', '4', 
> '2'), 25, [400, 400]);
> 
> 
> 
> for i = 0:360
> 
> 
> a.rotation_angles = [60, i];
> 
> 
> 
> // save the graph as image
> 
> 
> xs2png(f.figure_id, 'F:\testImg.bmp');
> 
> 
> 
> // read image
> 
> 
> img = imread('F:\testImg.png', CV_LOAD_IMAGE_GRAYSCALE);
> 
> 
> 
> // add image as frame to avi
> 
> 
> VideoWriter_write(videoWriter, img);
> 
> 
> 
> end
> 
> 
> 
> // free avi from memory
> delete_VideoWriter(videoWriter);
> 
> 
> 
> 
> 
> Am Fr., 2. Juli 2021 um 00:22 Uhr schrieb Heinz Nabielek 
> :
> Colleagues:
> 
> the code below generates 14,500 random points in a spherical volume.
> How can I make a movie that shows how the sphere rotates around the z-axis? 
> And how do I store the animation? animaGIF does not exist on my Scilab 6.1 in 
> my macOS 10.15.7
> 
> Heinz
> 
> 
> R=23;
> N=14500;
> r=grand(N,1,'unf',0,R);
> phi=grand(N,1,'unf',0,2*%pi);
> theta=grand(N,1,'unf',0,%pi);
> x=r.*cos(phi).*sin(theta);
> y=r.*sin(phi).*sin(theta);
> z=r.*cos(theta);
> scatter3d(x,y,z,msizes=4);
> gca().rotation_angles = [60, 60];
> 
> 
> __
> Dr Heinz Nabielek
> Schüttelstrasse 77A/11
> A-1020 Wien, Österreich
> Tel +43 1 276 56 13
> cell +43 677 616 349 22
> heinznabie...@me.com
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] How can I make a movie that shows how the sphere rotates around the z-axis?

[Heinz Nabielek]

Colleagues:

the code below generates 14,500 random points in a spherical volume.
How can I make a movie that shows how the sphere rotates around the z-axis? And 
how do I store the animation? animaGIF does not exist on my Scilab 6.1 in my 
macOS 10.15.7

Heinz


R=23;
N=14500;
r=grand(N,1,'unf',0,R);
phi=grand(N,1,'unf',0,2*%pi);
theta=grand(N,1,'unf',0,%pi);
x=r.*cos(phi).*sin(theta);
y=r.*sin(phi).*sin(theta);
z=r.*cos(theta);
scatter3d(x,y,z,msizes=4);
gca().rotation_angles = [60, 60];


__
Dr Heinz Nabielek
Schüttelstrasse 77A/11
A-1020 Wien, Österreich
Tel +43 1 276 56 13
cell +43 677 616 349 22
heinznabie...@me.com
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] xfpoly: would wish that the filling colour could be made transparent.

[Heinz Nabielek]

On 10.04.2021, at 08:10, Antoine Monmayrant mailto:antoine.monmayr...@laas.fr>> wrote:
> 
> 
> On 09/04/2021 23:55, Heinz Nabielek wrote:
>> xfpoly does generally a good job for me, sometimes I would wish that the 
>> filling colour could be made transparent.
> This is a much needed improvement of the scilab graphical stack that 
> currently does not provide any transparency.
> I think this is not an easy improvement.
> At the moment, my workaround is to plot everything I need, export as svg and 
> than add the transparency I need in the svg using inkscape or directly 
> editing the svg file with a text editor...
>> Heinz
>> 
>> PS 1: Is there a new version of the 2011 BetweenCurves around?
> Er, no, it was just a dirty hack I needed for my own publications and I think 
> 2011 is the most recent one.
> I can try to see how to improve it if this can improve scilab...




I had initiated by log vs lin plot with 
>> plot2d([0 80],[1 100], style=0,logflag = "nl");

but BetweenCurves starts with its own plot and here I would not know, what to 
do...
Heinz

BTW, how do the French infection rates look like?

>> 
>> PS 2: Any suggestion to make my clumsy coding more elegent, is highly welcome
>> 
>> PS 2: BTW, since the recent lockdown, infection rates are coming down in 
>> Austria.
>> 
>> 
>> A=[12.62813.94217.07717.05415.59414.97614.796
>> 11.87513.44816.50417.71719.44716.09913.30213.762 
>>16.03219.49222.09820.42521.08720.64914.268
>> 19.40222.52526.86223.51427.60323.85115.83021.570 
>>28.68226.10929.97428.72724.70521.45827.087
>> 28.40133.67035.11933.96228.12021.30127.24437.467 
>>37.71539.49037.56930.48027.09838.36636.951
>> 35.09743.75939.299]';
>> d=(1:length(A))';
>> up=10^(d/53);
>> M=[ones(up) up];
>> aa=M\A;
>> B=inv(M'*M);
>> DD=(1:110)';
>> U=10^(DD/53);
>> MM=[ones(U) U];
>> yh = M*aa; //Fitted values yh to approximate measured y's
>> e=A-yh; //Errors or residuals
>> SSE=e'*e; //Sum of squared errors
>> ybar=mean(A); R2=1-SSE/sum((A-ybar)^2);
>> [m n]=size(M);
>> MSE = SSE/(m-n-1); //Mean square error
>> C=MSE*B  // covariance matrix
>> sig=sqrt(MSE);
>> seb=sqrt(diag(C));
>> [aa seb]
>> [n pp]=size(M);
>> CONF=.95; alpha=1-CONF;
>> ta2 = cdft('T',n-pp,1-alpha/2,alpha/2); //t-value for alpha/2
>> yhh= MM*aa;
>> p=sig*sqrt(diag(1+MM*B*MM'));
>> N=[yhh+ta2*p  yhh-ta2*p];
>> polyX = [DD;flipdim(DD,1)];
>> polyY = [N(:,1);flipdim(N(:,2),1)];
>> plot2d([0 80],[1 100], style=0,logflag = "nl");
>> xgrid;
>> xfpoly(polyX, polyY,6);
>> plot(DD,MM*aa,'g.-');
>> plot(d,A,'b.') ;
>> title('AUSTRIA daily infection rates per 100,000','fontsize',5);
>> xlabel('days since 1 Feb 2021','fontsize',3);
>> ylabel('number of infections per day per 100,000','fontsize',3);
>> legend('data from Johns Hopkins GitHub','95% confidence range','model 
>> prediction','AUSTRIA recorded',4);
>> 
>> 
>> 
> ___
> users mailing list
> users@lists.scilab.org <mailto:users@lists.scilab.org>
> http://lists.scilab.org/mailman/listinfo/users 
> <http://lists.scilab.org/mailman/listinfo/users>

___
users mailing list
users@lists.scilab.org <mailto:users@lists.scilab.org>
http://lists.scilab.org/mailman/listinfo/users 
<http://lists.scilab.org/mailman/listinfo/users>___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] xfpoly: would wish that the filling colour could be made transparent.

[Heinz Nabielek]

On 10.04.2021, at 08:10, Antoine Monmayrant  wrote:
> 
> 
> On 09/04/2021 23:55, Heinz Nabielek wrote:
>> xfpoly does generally a good job for me, sometimes I would wish that the 
>> filling colour could be made transparent.
> This is a much needed improvement of the scilab graphical stack that 
> currently does not provide any transparency.
> I think this is not an easy improvement.
> At the moment, my workaround is to plot everything I need, export as svg and 
> than add the transparency I need in the svg using inkscape or directly 
> editing the svg file with a text editor...
>> Heinz
>> 
>> PS 1: Is there a new version of the 2011 BetweenCurves around?
> Er, no, it was just a dirty hack I needed for my own publications and I think 
> 2011 is the most recent one.
> I can try to see how to improve it if this can improve scilab...




I had initiated by log vs lin plot with 
>> plot2d([0 80],[1 100], style=0,logflag = "nl");

but BetweenCurves starts with its own plot and here I would not know, what to 
do...
Heinz

BTW, how do the French infection rates look like?

>> 
>> PS 2: Any suggestion to make my clumsy coding more elegent, is highly welcome
>> 
>> PS 2: BTW, since the recent lockdown, infection rates are coming down in 
>> Austria.
>> 
>> 
>> A=[12.62813.94217.07717.05415.59414.97614.796
>> 11.87513.44816.50417.71719.44716.09913.30213.762 
>>16.03219.49222.09820.42521.08720.64914.268
>> 19.40222.52526.86223.51427.60323.85115.83021.570 
>>28.68226.10929.97428.72724.70521.45827.087
>> 28.40133.67035.11933.96228.12021.30127.24437.467 
>>37.71539.49037.56930.48027.09838.36636.951
>> 35.09743.75939.299]';
>> d=(1:length(A))';
>> up=10^(d/53);
>> M=[ones(up) up];
>> aa=M\A;
>> B=inv(M'*M);
>> DD=(1:110)';
>> U=10^(DD/53);
>> MM=[ones(U) U];
>> yh = M*aa; //Fitted values yh to approximate measured y's
>> e=A-yh; //Errors or residuals
>> SSE=e'*e; //Sum of squared errors
>> ybar=mean(A); R2=1-SSE/sum((A-ybar)^2);
>> [m n]=size(M);
>> MSE = SSE/(m-n-1); //Mean square error
>> C=MSE*B  // covariance matrix
>> sig=sqrt(MSE);
>> seb=sqrt(diag(C));
>> [aa seb]
>> [n pp]=size(M);
>> CONF=.95; alpha=1-CONF;
>> ta2 = cdft('T',n-pp,1-alpha/2,alpha/2); //t-value for alpha/2
>> yhh= MM*aa;
>> p=sig*sqrt(diag(1+MM*B*MM'));
>> N=[yhh+ta2*p  yhh-ta2*p];
>> polyX = [DD;flipdim(DD,1)];
>> polyY = [N(:,1);flipdim(N(:,2),1)];
>> plot2d([0 80],[1 100], style=0,logflag = "nl");
>> xgrid;
>> xfpoly(polyX, polyY,6);
>> plot(DD,MM*aa,'g.-');
>> plot(d,A,'b.') ;
>> title('AUSTRIA daily infection rates per 100,000','fontsize',5);
>> xlabel('days since 1 Feb 2021','fontsize',3);
>> ylabel('number of infections per day per 100,000','fontsize',3);
>> legend('data from Johns Hopkins GitHub','95% confidence range','model 
>> prediction','AUSTRIA recorded',4);
>> 
>> 
>> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] xfpoly: would wish that the filling colour could be made transparent.

[Heinz Nabielek]

xfpoly does generally a good job for me, sometimes I would wish that the filling colour could be made transparent.HeinzPS 1: Is there a new version of the 2011 BetweenCurves around?PS 2: Any suggestion to make my clumsy coding more elegent, is highly welcomePS 2: BTW, since the recent lockdown, infection rates are coming down in Austria.A=[12.62813.94217.07717.05415.59414.97614.79611.87513.44816.50417.71719.44716.09913.30213.76216.03219.49222.09820.42521.08720.64914.26819.40222.52526.86223.51427.60323.85115.83021.57028.68226.10929.97428.72724.70521.45827.08728.40133.67035.11933.96228.12021.30127.24437.46737.71539.49037.56930.48027.09838.36636.95135.09743.75939.299]';
d=(1:length(A))';
up=10^(d/53);
M=[ones(up) up];
aa=M\A;
B=inv(M'*M);
DD=(1:110)';
U=10^(DD/53);
MM=[ones(U) U];
yh = M*aa; //Fitted values yh to approximate measured y's
e=A-yh; //Errors or residuals
SSE=e'*e; //Sum of squared errors
ybar=mean(A); R2=1-SSE/sum((A-ybar)^2);
[m n]=size(M);
MSE = SSE/(m-n-1); //Mean square error
C=MSE*B  // covariance matrix
sig=sqrt(MSE);
seb=sqrt(diag(C));
[aa seb]
[n pp]=size(M);
CONF=.95; alpha=1-CONF;
ta2 = cdft('T',n-pp,1-alpha/2,alpha/2); //t-value for alpha/2
yhh= MM*aa;
p=sig*sqrt(diag(1+MM*B*MM'));
N=[yhh+ta2*p  yhh-ta2*p];
polyX = [DD;flipdim(DD,1)];
polyY = [N(:,1);flipdim(N(:,2),1)];
plot2d([0 80],[1 100], style=0,logflag = "nl");
xgrid;
xfpoly(polyX, polyY,6);
plot(DD,MM*aa,'g.-');
plot(d,A,'b.') ;
title('AUSTRIA daily infection rates per 100,000','fontsize',5);
xlabel('days since 1 Feb 2021','fontsize',3);
ylabel('number of infections per day per 100,000','fontsize',3);
legend('data from Johns Hopkins GitHub','95% confidence range','model prediction','AUSTRIA recorded',4);

AUSTRIA daily infection rates per 100,000 + confidence.pdf
Description: Adobe PDF document
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Can't get current directory

[Heinz Nabielek]

Perfect for me...



--> pwd
 ans  =  "/Users/heinz"

--> cd Desktop
 ans  =   "/Users/heinz/Desktop"

> On 24.01.2021, at 22:05, Gerard Awanou  wrote:
> 
> I am using the GUI. It starts with the message
> 
> Startup execution:
>   loading initial environment
> 
> Then, I try to go to a specific directory
> 
> cd Desktop
> 
> I then got the error message
> 
> "/Users/gerardawanou/Desktop"
> Can't get current directory.
> 
> I did give Scilab full disc acsess
> 
> On Sun, Jan 24, 2021 at 2:52 PM Heinz Nabielek  wrote:
> Gérard,
> 
> I added SciLab into Security/Privacy/Full Disk Access.
> 
> Works fine. But furious about the way the newer Mac OSses make life ever more 
> difficult..
> 
> Does anybody have a better suggestion?
> Heinz
> 
> 
> 
> > On 24.01.2021, at 21:03, Gerard Awanou  wrote:
> > 
> > Hi
> > 
> > I am running Scilab 6.1 on mac OS Catalina and I am getting the annoying 
> > message
> > "Can't get current directory".
> > 
> > How do I fix this?
> > 
> > ___
> > users mailing list
> > users@lists.scilab.org
> > http://lists.scilab.org/mailman/listinfo/users
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Can't get current directory

[Heinz Nabielek]

Gérard,

I added SciLab into Security/Privacy/Full Disk Access.

Works fine. But furious about the way the newer Mac OSses make life ever more 
difficult..

Does anybody have a better suggestion?
Heinz



> On 24.01.2021, at 21:03, Gerard Awanou  wrote:
> 
> Hi
> 
> I am running Scilab 6.1 on mac OS Catalina and I am getting the annoying 
> message
> "Can't get current directory".
> 
> How do I fix this?
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] plot(x,y,':') & dotted line : is this a bug?

[Heinz Nabielek]

Ok: perhaps not perfect dots

scf();plot(1:10,':'); scf();plot(1:100,':');


> On 14.09.2020, at 16:28, Antoine Monmayrant  
> wrote:
> 
>  scf();plot(1:10,':')

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] plot(x,y,':') & dotted line : is this a bug?

[Heinz Nabielek]

> On 14.09.2020, at 16:14, Antoine Monmayrant  
> wrote:
> 
>  plot(rand(1:10),':')


I get a perfect dotted line.
Heinz

  "Scilab Version: ""6.1.1.0"   
 
  "Operating System: "  "Mac OS X 10.15.6"  
 
  "Java version: "  "1.8.0_51"  
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] calendar does not work

[Heinz Nabielek]

Ah: it is that terrible Austrian dialect that they are speaking here in 
Vienna

--> getenv("LANG")
 ans  =

  "en_AT.UTF-8"

--> calendar
 
 AugÂ2020
 MonÂÂTueÂÂWedÂÂThuÂÂFriÂÂSatÂÂSun
 ÂÂ12
 Â3456789
 10ÂÂÂ11ÂÂÂ12ÂÂÂ13ÂÂÂ14ÂÂÂ15ÂÂÂ16
 17ÂÂÂ18ÂÂÂ19ÂÂÂ20ÂÂÂ21ÂÂÂ22ÂÂÂ23
 24ÂÂÂ25ÂÂÂ26ÂÂÂ27ÂÂÂ28ÂÂÂ29ÂÂÂ30
 31ÂÂ

Who do I change to en_US.UTF-8 permanently under zsh? All my attempts change it 
only momentarily
Heinz



> On 25.08.2020, at 13:30, Stéphane Mottelet  wrote:
> 
> Hello,
> 
> With english defaults under OSX I cannot reproduce your output :
> 
> --> getenv("LANG")
>  ans  =
> 
>   "en_US.UTF-8"
> 
> --> calendar
> 
>  Aug 2020
>  Mon  Tue  Wed  Thu  Fri  Sat  Sun
>12
>   3456789
>  10   11   12   13   14   15   16
>  17   18   19   20   21   22   23
>  24   25   26   27   28   29   30
>  31
> 
> S.
> 
> Le 19/08/2020 à 00:10, Heinz Nabielek a écrit :
>> --> calendar
>>AugÂ2020
>>  MonÂÂTueÂÂWedÂÂThuÂÂFriÂÂSatÂÂSun
>>  ÂÂ12
>>  Â3456789
>>  10ÂÂÂ11ÂÂÂ12ÂÂÂ13ÂÂÂ14ÂÂÂ15ÂÂÂ16
>>  17ÂÂÂ18ÂÂÂ19ÂÂÂ20ÂÂÂ21ÂÂÂ22ÂÂÂ23
>>  24ÂÂÂ25ÂÂÂ26ÂÂÂ27ÂÂÂ28ÂÂÂ29ÂÂÂ30
>>  31ÂÂ
>> 
>> What had happened? Scilab 6.1 under macOS Catalina 10.15.6
>> 
>> Heinz
>> ___
>> users mailing list
>> users@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
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] errors (uncertainties) in non-linear least-squares fitting parameters

[Heinz Nabielek]

On 24.08.2020, at 23:08, Rafael Guerra  wrote:
> 
> Hi Heinz,
>  
> For the regression errors, I am not an expert but from wikipedia or from 
> reference below, I would risk the following code (at your peril):
> https://pages.mtu.edu/~fmorriso/cm3215/UncertaintySlopeInterceptOfLeastSquaresFit.pdf
>  
> // Note: for degrees of freedom>=6, t-distribution ~2
> N = length(MW);
> mx = mean(MW);
> SSxx = sum((MW -mx).^2);
> Ea = diag(2*sig/sqrt(SSxx))  // take Ea diagonals; slope 95% confidence 
> Eb = diag(2*sig*sqrt(1/N+mx^2/SSxx)) // take Eb diagonals; intercept 95% 
> confidence
>  
> Concerning the least squares regression part, it seems the code may be 
> written more compactly using reglin:
>  
> [a,b,sig]=reglin(MW',Y') // simple least squares linear regression
> GG= a.*.xx' + repmat(b,size(xx'))
> plot(xx,GG','LineWidth',1);


Here is a little misunderstanding (my fault: I had not explained it).
I want all three straight lines to go simultaneously through one point at the 
negative x-axis. This is why I had to use a non-linear least-squares fit.

Heinz




___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] errors (uncertainties) in non-linear least-squares fitting parameters

[Heinz Nabielek]

I have successfully fitted straight lines to my avian mortality Monte-Carlo 
simulation model (function of turbine size and wind speed distributions):

f(x) = pi (x+p4) with p1 for 7, p2 for 6 and p3 for 5 m/s mean wind speed
   p1 =   0.3930457
   p2 =   0.3492537
   p3 =   0.2987269
   p4 =   4.058154 

But, sorry, I do not manage to derive the errors (uncertainties) in the 
parameters. I know I need the Jacobian, but have been unable to code it in 
Scilab.

Can you help?
Best greetings
Heinz

BTW, these bird fatalities are far too large when compared to observations: as 
a next modelling step, I need to introduce avoidance factors.



Avian Mortality.sce
Description: Binary data



Avian Mortality.pdf
Description: Adobe PDF document
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] calendar does not work

[Heinz Nabielek]

--> calendar
 
 AugÂ2020
 MonÂÂTueÂÂWedÂÂThuÂÂFriÂÂSatÂÂSun
 ÂÂ12
 Â3456789
 10ÂÂÂ11ÂÂÂ12ÂÂÂ13ÂÂÂ14ÂÂÂ15ÂÂÂ16
 17ÂÂÂ18ÂÂÂ19ÂÂÂ20ÂÂÂ21ÂÂÂ22ÂÂÂ23
 24ÂÂÂ25ÂÂÂ26ÂÂÂ27ÂÂÂ28ÂÂÂ29ÂÂÂ30
 31ÂÂ

What had happened? Scilab 6.1 under macOS Catalina 10.15.6

Heinz
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] the linear least-squares fitting problem

[Heinz Nabielek]

Dear Scilabers,for mathematicians, the linear least-squares fitting is a trivial problem, but as an experimental physicist I had problems finding the right correlations for predicting the Confidence Limits of the fitted functional relationship. For a trivial set of experimental data, I wrote a few lines of Scilab code and have 2 questions:- is the computation of confidence limits correct?- for large data sets and higher degree polynomials, can we make the computations more efficient and reliable?Thanks in advanceHeinz_x=[0.1269868 0.135477 0.221034 0.45 0.8350086 0.9057919 0.9133759 0.9688678]';
y=[0.0695843 0.0154644 0.0893982 0.5619441 1.3879476 1.5639274  1.6570076 1.9061488]';
plot(x,y,'m.','markersize',14);xgrid(1,1,1);
X=[ones(x) x x^2]; 		// matrix of independent vectors
b=X\y; 	// parameters of linear LSFIT with a 2nd degree polynomial
yh = X*b; 	//Fitted values yh to approximate measured y's
e=y-yh; 	//Errors or residuals
SSE=e'*e; 	//Sum of squared errors
ybar=mean(y); R2=1-SSE/sum((y-ybar)^2);
[m n]=size(X);
MSE = SSE/(m-n-1); 	//Mean square error
se=sqrt(MSE);
C=MSE*inv(X'*X);
seb=sqrt(diag(C));
CONF=.95; alpha=1-CONF;
ta2 = cdft('T',m-n,1-alpha/2,alpha/2); 	//t-value for alpha/2
xx=(0.1:.025:1)';  	// cover the whole abscissa range in detail
XX=[ones(xx) xx xx^2];
yhh=XX*b; 	// predicted extended LSFIT curve
plot(xx,yhh,'kd-','MarkerSize',8);
[mm n]=size(XX);
sY = []; sYp = []; 	//Terms involved in Confidence Interval for Y, Ypred
for i=1:mm; 
sY = [sY; sqrt(XX(i,:)*C*XX(i,:)')]; 		//  standard error in fit
sYp = [sYp; se*sqrt(1+XX(i,:)*(C/se)*XX(i,:)')]; //  standard error in spread in population distribution
end;
plot(xx, yhh+ta2*sYp,'rd-');
plot(xx, yhh+ta2*sY,'bd-');
plot(xx, yhh-ta2*sY,'gd-');
plot(xx, yhh-ta2*sYp,'md-');
title('LINEAR LEAST-SQUARES FIT WITH 95% CONFIDENCE BOUNDS ON FIT AND ON THE POPULATION SPREAD','fontsize',2);
xlabel('x-values','fontsize',5);
ylabel('y-values','fontsize',5);
legend('measurement data','LSFIT of 2nd order polynomial','upper 95% CL population','upper 95% CL fit','lower 95% CL fit','lower 95% CL population',2);

LINEAR LEAST-SQUARES FIT WITH 95% CONFIDENCE BOUNDS ON FIT AND ON THE POPULATION SPREAD.pdf
Description: Adobe PDF document
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] [Scilab-Dev] Scilab 6.1 for OSX is out !

[Heinz Nabielek]

Thanks a lot. Works perfect.
Heinz

> On 18.06.2020, at 16:52, Stéphane Mottelet  wrote:
> 
> Hello,
> 
> The branch-6.1 OSX build has been updated yesterday evening. To have an idea 
> of bug fixes and other improvements since the last build see the following 
> list on CR:
> 
> https://codereview.scilab.org/#/q/status:merged+since:2020-04-20
> 
> S.
> 
> Le 21/04/2020 à 19:21, Stéphane Mottelet a écrit :
>> Hi all,
>> 
>> Finally I managed to package a branch-6.1 version. "branch-6.1" means that 
>> you will have the actual features of 6.1 version plus the bug fixes and 
>> improvements provided by the team and contributors since its release. The 
>> archive is available to dowload at the following url:
>> 
>> https://www.utc.fr/~mottelet/scilab_for_macOS.html
>> 
>> Tests have been made under HighSierra, Mojave, Catalina and everything seems 
>> OK. Please use the ML to report eventual problems.
>> 
>> Enjoy !
>> 
> -- 
> Stéphane Mottelet
> Ingénieur de recherche
> EA 4297 Transformations Intégrées de la Matière Renouvelable
> Département Génie des Procédés Industriels
> Sorbonne Universités - Université de Technologie de Compiègne
> CS 60319, 60203 Compiègne cedex
> Tel : +33(0)344234688
> http://www.utc.fr/~mottelet
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Pong contest

[Heinz Nabielek]

From an iMac macOS 10.15.5 (19F101):

...SCALAR RESULT:  Mflops= 2.915 ;  Pi ~ 3.141592 
...ARRAY RESULT:  Mflops= 1.047 ;  Pi ~ 3.141592 

  Model Identifier: iMac14,3
  Processor Name:   Quad-Core Intel Core i7
  Processor Speed:  3.1 GHz
  Number of Processors: 1
  Total Number of Cores:4

Heinz

> On 07.06.2020, at 12:01, Claus Futtrup  wrote:
> 
> Hi Rafael
> 
> ...SCALAR RESULT:  Mflops= 2.737 ;  Pi ~ 3.141592 
> ...ARRAY RESULT:  Mflops= 1.284 ;  Pi ~ 3.141592
> 
> The machine is a two-year old Lenovo T480, Core i7 8th Gen. Built-in Intel 
> Graphics (no external card)
> 
> /Claus
> 
> On 07-06-2020 09:54, Rafael Guerra wrote:
>> Hello Stéphane et al.,
>>  
>> Would it make sense in this “Pong contest”, in addition to the OS used, to 
>> also provide some sort of benchmarking of machine CPU?
>>  
>> Fyi, find below simple test by Gershenfeld (1999), it computes Pi=4*atan(1) 
>> by two methods:
>> Scalar test:  Pi(N)~Sum{i=1:N; 0.5/((i-0.75)*(i-0.25))} 
>> Array test:   Pi(N)= Pi(N-1) + 0.5/((N-0.75)*(N-0.25))
>>  
>> // Based on benchmarking test in Gershenfeld Mathematical modeling book 
>> (1999)
>> // There are 5 floating point operations per step, total 5 Mflop with N=1e6
>> N= 1e6;
>> tic();
>> ps= 0;
>> for i=1:N
>>   ps= ps + 0.5/((i-0.75)*(i-0.25));
>> end
>> dt=toc();
>> printf('...SCALAR RESULT:  Mflops= %4.3f ;  Pi ~ %7.6f \n',5.0/dt,ps);
>>  
>> tic();
>> pv(1)= 0.5/((1-0.75)*(1-0.25));
>> for i=2:N
>>   pv(i)= pv(i-1) + 0.5/((i-0.75)*(i-0.25));
>> end
>> dt=toc();
>> printf('...ARRAY RESULT:  Mflops= %4.3f ;  Pi ~ %7.6f \n',5.0/dt,pv(N));
>>  
>>  
>> PS:
>> ...SCALAR RESULT:  Mflops= 3.489
>> ...ARRAY RESULT:  Mflops= 1.367
>>  
>> Regards,
>> Rafael
>>  
>>  
>> On 05-06-2020 14:48, Stéphane Mottelet wrote:
>> Hello all, 
>> 
>> As you may have noticed, there is now a little game in the Gaphics/Animation 
>> section of the demonstrations. As the speed and responsiveness seems to be 
>> similar under all platforms, it would be funny to launch a little contest. 
>> In the next Scilab version it would be interesting to implement a high score 
>> online record, but until then, you can answer this message with the same 
>> kind of screenshot I joined, showing you score and just taken after the game 
>> end with and with the last dialog on the figure. 
>> 
>> Cheers, 
>> 
>> 
>> ___
>> users mailing list
>> 
>> users@lists.scilab.org
>> http://lists.scilab.org/mailman/listinfo/users
> 
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Pong contest

[Heinz Nabielek]

Heinz

> On 05.06.2020, at 14:48, Stéphane Mottelet  wrote:
> 
> Hello all,
> 
> As you may have noticed, there is now a little game in the Gaphics/Animation 
> section of the demonstrations. As the speed and responsiveness seems to be 
> similar under all platforms, it would be funny to launch a little contest. In 
> the next Scilab version it would be interesting to implement a high score 
> online record, but until then, you can answer this message with the same kind 
> of screenshot I joined, showing you score and just taken after the game end 
> with and with the last dialog on the figure.
> 
> Cheers,
> 
> -- 
> Stéphane Mottelet
> Ingénieur de recherche
> EA 4297 Transformations Intégrées de la Matière Renouvelable
> Département Génie des Procédés Industriels
> Sorbonne Universités - Université de Technologie de Compiègne
> CS 60319, 60203 Compiègne cedex
> Tel : +33(0)344234688
> http://www.utc.fr/~mottelet
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Can Scilab compute the inverse of the regularized Incomplete Beta Function?

[Heinz Nabielek]

CONF=.95, N=1432, n=1
cdfbet("XY", n+1, N+1-n, CONF, 1-CONF)

is doing it just fine. Correct and with high precision...
h

> On 18.05.2020, at 19:08, Rafael Guerra  wrote:
> 
> Hi Heinz,
>  
> Fyi, the following site provides Matlab code that may be translated to Scilab:
> https://people.sc.fsu.edu/~jburkardt/m_src/asa109/asa109.html
> betain.m : incomplete Beta function ratio
> xinbta.m : inverse of the incomplete Beta function
>  
> The Scilab Atoms 'Distfun' package contains:
> distfun_betainc : Regularized Incomplete Beta function
>  
> Regards,
> Rafael
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Can Scilab compute the inverse of the regularized Incomplete Beta Function?

[Heinz Nabielek]

CONF=.95, N=1432, n=1
cdfbet("XY", n+1, N+1-n, CONF, 1-CONF)

is doing it just fine. Correct and with high precision...
h


On 18.05.2020, at 18:09, Tim Wescott  wrote:
> 
> So you have \beta(x, n+1, N+1-n) = 0.95, and you want to solve for x?
> 
> fsolve will do this for a single value of the confidence.  Is that
> sufficient?
> 
> On Sun, 2020-05-17 at 23:49 +0200, Heinz Nabielek wrote:
>> Dear SciLabers:
>> 
>> can Scilab compute the inverse of the regularized Incomplete Beta
>> Function?
>> 
>> Example: in unbiased sampling in Austria with sample size N=1432,
>> they detected n=1 infections.
>> Therefore, expected infected fraction = 0.000698324.
>> 
>> But this does not say much, because the sample size was small and the
>> "success" was extremely small (fortunately).
>> 
>> The standard procedure therefore is to derive the one-sided 95% upper
>> confidence limit:
>> CONF=0.95; N=1432; n=1:
>> One-sided 95% upper confidence limit fraction = BETA.INV(CONF, n+1,
>> N+1-n) = 0.003306121
>> 
>> How would I do that in Scilab?
>> Heinz

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Can Scilab compute the inverse of the regularized Incomplete Beta Function?

[Heinz Nabielek]

> On 18.05.2020, at 14:18, Dang Ngoc Chan, Christophe 
>  wrote:
> 
> Hello,
> 
>> De  Heinz Nabielek
>> Envoyé : dimanche 17 mai 2020 23:50
>> 
>> CONF=0.95; N=1432; n=1:
>> One-sided 95% upper confidence limit fraction = BETA.INV(CONF, n+1, N+1-n)
>> = 0.003306121
>> 
>> How would I do that in Scilab?
> 
> Would it be :
> 
> [X,Y]=cdfbet("XY", n+1, N+1-n, CONF, 1-CONF)
> X  =
> 
>   0.0033061


Written with more digits

0.0033061215, and EXCEL finds

0.003306121493. Very good agreement. I am glad I had asked.

And the HELP provides real support with the statement

Description
Calculates any one parameter of the beta distribution given values for the 
others (The beta density is proportional to t^(A-1) * (1-t)^(B-1).

Since my originating binomial distribution is obviously

t^n  *  (1-t)^(N-n),

the parameters for the Beta Function have to be

A=n+1
B=N+1-n

to obtain an upper conficende limit fraction = Beta.Inv(CONF, n+1, N+1-n)  !!!

Textbooks write Beta.Inv(CONF, n+1, N-n) since the last 50 years and it is 
wrong. And "R" is also using this wrong correlation in their binom.test

Great many thanks for all the help.
Heinz

__
Dr Heinz Nabielek
Schüttelstrasse 77A/11
A-1020 Wien, Österreich
Tel +43 1 276 56 13
cell +43 677 616 349 22
heinznabie...@me.com


___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] Can Scilab compute the inverse of the regularized Incomplete Beta Function?

[Heinz Nabielek]

Dear SciLabers:

can Scilab compute the inverse of the regularized Incomplete Beta Function?

Example: in unbiased sampling in Austria with sample size N=1432, they detected 
n=1 infections.
Therefore, expected infected fraction = 0.000698324.

But this does not say much, because the sample size was small and the "success" 
was extremely small (fortunately).

The standard procedure therefore is to derive the one-sided 95% upper 
confidence limit:
CONF=0.95; N=1432; n=1:
One-sided 95% upper confidence limit fraction = BETA.INV(CONF, n+1, N+1-n) = 
0.003306121

How would I do that in Scilab?
Heinz

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] The application "scilab-branch-6.1.app" is not open any more.

[Heinz Nabielek]

On 08.05.2020, at 08:46, Stéphane Mottelet  wrote:
> 
>  [x,y]=getversion() 

x,y]=getversion() 
 x  =   "scilab-branch-6.1"
 y  =   "GCC"  "x64"  "modelicac"  "release"  "Apr 21 2020"  "15:18:18"

Thanks for all the help
Heinz

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] The application "scilab-branch-6.1.app" is not open any more.

[Heinz Nabielek]

--> getversion()
 ans  =   "scilab-branch-6.1"

but it still says Created: Tuesday, 28. May 2019 at 19:28 after complete reload 
and re-installation.

In any case, everything works perfectly fine, I just a completely upset iMac 
for a few hours...

Heinz


> On 07.05.2020, at 21:08, Stéphane Mottelet  wrote:
> 
> Do you have the latest build (20/04/2020) ?
> 
> S.
> 
>> Le 7 mai 2020 à 20:13, Heinz Nabielek  a écrit :
>> 
>> Suddenly I cannot start. It says
>> The application "scilab-branch-6.1.app" is not open any more.
>> 
>> What to do?
>> 
>> Had worked nicely for a week
>> Heinz
>> 
>> System Version:macOS 10.15.4 (19E287)
>> Kernel Version:Darwin 19.4.0
>> ___
>> users mailing list
>> users@lists.scilab.org
>> https://antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

[Scilab-users] The application "scilab-branch-6.1.app" is not open any more.

[Heinz Nabielek]

Suddenly I cannot start. It says
The application "scilab-branch-6.1.app" is not open any more.

What to do?

Had worked nicely for a week
Heinz

  System Version:   macOS 10.15.4 (19E287)
  Kernel Version:   Darwin 19.4.0
___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Corona modelling

[Heinz Nabielek]

And how does it compare with Johns Hopkins?
Heinz



> On 07.05.2020, at 13:49, hdf  wrote:
> 
> and here are the results, with parameters given:  
> Sent from the Scilab users - Mailing Lists Archives mailing list archive at 
> Nabble.com.
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users

Re: [Scilab-users] Scilab 6.1 for OSX is out !

[Heinz Nabielek]

Thanks very much indeed. Works to perfection under macOS 10.15.4 (19E287).
Heinz



> On 21.04.2020, at 19:21, Stéphane Mottelet  wrote:
> 
> Hi all,
> 
> Finally I managed to package a branch-6.1 version. "branch-6.1" means that 
> you will have the actual features of 6.1 version plus the bug fixes and 
> improvements provided by the team and contributors since its release. The 
> archive is available to dowload at the following url:
> 
> https://www.utc.fr/~mottelet/scilab_for_macOS.html
> 
> Tests have been made under HighSierra, Mojave, Catalina and everything seems 
> OK. Please use the ML to report eventual problems.
> 
> Enjoy !
> 
> -- 
> Stéphane Mottelet
> Ingénieur de recherche
> EA 4297 Transformations Intégrées de la Matière Renouvelable
> Département Génie des Procédés Industriels
> Sorbonne Universités - Université de Technologie de Compiègne
> CS 60319, 60203 Compiègne cedex
> Tel : +33(0)344234688
> http://www.utc.fr/~mottelet
> 
> ___
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users

___
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users