Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


2 Pages12>
Options
Go to last post Go to first unread
Offline ndtma  
#1 Posted : 08 November 2019 11:08:51(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Thank you Andrey for the Latest Beta 7250.
But, the attached SS File yields an error due to different result of
integration in 7250 compared to result in 7219.

SPI_Calc_Plain.sm (112kb) downloaded 24 time(s).

7219_CDF Gamma.png

7250_CDF Gamma.png
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline Davide Carpi  
#2 Posted : 08 November 2019 12:13:29(UTC)
Davide Carpi


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 13/01/2012(UTC)
Posts: 2,639
Man
Italy
Location: Italy

Was thanked: 1323 time(s) in 873 post(s)
Do you have the same settings about integrals' accuracy in both the versions?

2019-11-08 10_11_42-SMath Studio - [SPI_Calc_Plain.sm].png 2019-11-08 10_10_32-SMath Studio - [SPI_Calc_Plain.sm_].png

If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects
Offline ndtma  
#3 Posted : 08 November 2019 14:20:08(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Yes, same settings in both versions.
But the difference is in Eqn 2

7219_Settings.png
7250_Integral Accuracy.png
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Andrey Ivashov  
#4 Posted : 08 November 2019 18:15:57(UTC)
Andrey Ivashov


Rank: Administration

Groups: Developers, Registered, Knovel Developers, Administrators, Advanced Member
Joined: 11/07/2008(UTC)
Posts: 1,616
Man
Russian Federation

Was thanked: 1978 time(s) in 666 post(s)
Thank you very much! It was a real calculation error I've introduced in 7250.

Fixed in 7251 (will be available very soon).

UPDATE: Released here: https://en.smath.com/for...1--08-November-2019.aspx

Best regards, Andrey Ivashov.

Edited by user 08 November 2019 18:28:11(UTC)  | Reason: Not specified

thanks 1 user thanked Andrey Ivashov for this useful post.
on 08/11/2019(UTC)
Offline ndtma  
#5 Posted : 10 November 2019 06:15:32(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
I deleted my previous post for more corrections.
Still I get errors in SS 7251.
I will come back soon.
Sorry for the inconvenience.
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline ndtma  
#6 Posted : 10 November 2019 09:59:21(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Integration results of Mathcad attached (using first 10 elements of X
and Eqns 1, 2 and 3.)

SPI_Calc_Plain.sm (126kb) downloaded 19 time(s).

CDF MathCAD.png

7251_CDFpng.png

Rainfal_Data.xlsx (487kb) downloaded 15 time(s).

Edited by user 10 November 2019 18:54:20(UTC)  | Reason: Forgot to attach the rainfall data file (Excel)

Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#7 Posted : 10 November 2019 17:01:56(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: ndtma Go to Quoted Post
Integration results of Mathcad attached (using first 10 elements of X
and Eqns 1, 2 and 3.)

What I understand from your previous reports is that
the Simpson integrator is kaput again in latest SS version.
What more is kaput ? SS 6179 does not open your attachment.
In the mean time, can you doctor rkfixed your document ?
Cheers ... Jean

NDTMA.PNG

NDTMA.sm (55kb) downloaded 17 time(s).

Offline ndtma  
#8 Posted : 11 November 2019 08:31:56(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
In the mean time, can you doctor rkfixed your document ?
Cheers ... Jean


Thank you. I will try my best ASAP.
In the mean time following files are attached.

Test_Gamma_7251.sm (36kb) downloaded 18 time(s).

Test_Gamma_7251.pdf (193kb) downloaded 15 time(s).
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#9 Posted : 11 November 2019 16:59:00(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: ndtma Go to Quoted Post
Thank you. I will try my best ASAP.
In the mean time following files are attached.

Here is doctored sanity MCD.
Don't forget Smath Gamma(x) is only globally 4 decimals,
immaterial in your CDF integration. If you need to check
the peak location of PDF, we have the Golden ratio algo.
Enjoy ... Jean.

NDTMA test.PNG

NDTMA Test Gamma.sm (46kb) downloaded 12 time(s).
Offline Jean Giraud  
#10 Posted : 11 November 2019 19:14:35(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
... max PDF located

NDTMA Test Gamma.sm (23kb) downloaded 14 time(s).
Golden Ratio jmG.sm (36kb) downloaded 15 time(s).
Offline ndtma  
#11 Posted : 12 November 2019 13:07:16(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Sincerely, I was trying to show that results of Gamma Function CDF
in 7251 were not compatible with 7219.
I wish to add that similar working in 7219 reasonably tally with MCD and EXCEL,
while 7251 is not so. I believe that may be the reason why I got errors using
SS ICDF.normal function in SPI calculations (some CDF values>1 ??).

Anyway, I am attaching the following files and images for further examination please.
Many thanks to all for their comments.

Fimal Gamma Results.zip (1,170kb) downloaded 14 time(s).

Test Gamma EXCEL.png

Test_Gamma_7219.png

Test_Gamma_7251.png

Test_Gamma_MCD.png

Edited by user 12 November 2019 16:00:12(UTC)  | Reason: Not specified

Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Razonar  
#12 Posted : 12 November 2019 16:16:46(UTC)
Razonar


Rank: Advanced Member

Groups: Registered
Joined: 28/08/2014(UTC)
Posts: 1,356
Uruguay

Was thanked: 815 time(s) in 516 post(s)
Hi. You can try to implement your own algo's for gamma and int functions.

Clipboard01.gif

Test_Gamma_7251.sm (104kb) downloaded 20 time(s).

Best regards.
Alvaro.
thanks 2 users thanked Razonar for this useful post.
on 12/11/2019(UTC),  on 12/11/2019(UTC)
Offline ndtma  
#13 Posted : 12 November 2019 17:53:09(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Hi,Thank you. Your suggestion is highly commendable and accepted. Bravo!

But, very often, as practicing Engineers, we are inclined to use built-in
functions as much as possible in our day to day work in the design office.

The SPI calculation was such an effort I came across.

Edited by user 12 November 2019 17:56:48(UTC)  | Reason: Not specified

Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#14 Posted : 12 November 2019 18:07:05(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: ndtma Go to Quoted Post
Sincerely, I was trying to show that results of Gamma Function PDF and CDF
in 7251 were not compatible with 7219.
I wish to add that similar working in 7219 reasonably tally with MCD and EXCEL,
while 7251 is not so. I believe that may be the reason why I got errors using
SS ICDF.normal function in SPI calculations (some CDF values

As you say ... hard to doctor this patient.
Your download opened by default on my 6179 carries the bug.
Could it be the vectorize operator is the cause ?
If so: that would be quite a bug !
Cheers ... Jean

Test_Gamma_7251 [6179].sm (23kb) downloaded 16 time(s).

NDTMA 7251_6179.PNG
Offline ndtma  
#15 Posted : 12 November 2019 19:39:00(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
As you say ... hard to doctor this patient.
Your download opened by default on my 6179 carries the bug.
Could it be the vectorize operator is the cause ?
If so: that would be quite a bug !
Cheers ... Jean


It worked! Thanks, but how?

Test Gamma [2] 7251 [6179].sm (26kb) downloaded 13 time(s).

Test Gamma [2] 7251 [6179].png
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Razonar  
#16 Posted : 12 November 2019 19:43:07(UTC)
Razonar


Rank: Advanced Member

Groups: Registered
Joined: 28/08/2014(UTC)
Posts: 1,356
Uruguay

Was thanked: 815 time(s) in 516 post(s)
Originally Posted by: ndtma Go to Quoted Post


But, very often, as practicing Engineers, we are inclined to use built-in
functions as much as possible in our day to day work in the design office.



Hi Nugegodage. I absolutely agree. Actually, I post for try to find where the bug is.

I include gamma too, because I think that, like happen in matlab and maple, plugins can be written in native smath code, not in c++. Maybe in this case, gamma and other special functions are not the best example for do that, but this is only meanwhile they are in plugins, not in smath native functions.

In my point of view, plugins are good for graphics, some objects (like formatted cells, writer regions), communication with external programs (maple, maxima) and some other things.

Best regards.
Alvaro.
Offline ndtma  
#17 Posted : 12 November 2019 20:32:37(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
Could it be the vectorize operator is the cause ?


I have edited part of my original "SPI_Calc_Plain.sm" worksheet, and it
is evident the bug due to the vectorize operator.
Will come back soon with the edited version.

Thank you very much.
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline ndtma  
#18 Posted : 12 November 2019 20:52:43(UTC)
ndtma


Rank: Advanced Member

Groups: Registered
Joined: 05/06/2014(UTC)
Posts: 348
Man
Sri Lanka
Location: Colombo

Was thanked: 124 time(s) in 82 post(s)
This confirms Jean's guess about vectorize operator.

TEST-SPI_Calc_Plain.sm (53kb) downloaded 12 time(s).

7251 Vectorize Error.png
Look within!... The secret is inside you.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#19 Posted : 12 November 2019 21:17:55(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: ndtma Go to Quoted Post
It worked! Thanks, but how?

Answer is from observation(sss..)
My work sheet is clean and thus opened/calculated your 7251.
Your 7251 is corrupted and my 6179 executed the vectorize bug.
Mathematica has the expression shadow in context.
Read this: 6179 does not come with ZedGraph, one day some collab
posted a ZedGraph, 6179 did open and now my 6179 has ZedGraph.
SS 6346 had two versions: the official and un-official.
The un-official did lot more than the published version.
Detail about Gamma(x)like erf(x): under the integration operator,
they both need be vectorized. Maybe, the fact that I renamed/recast
it may have the virtue of doctoring. Though your Gamma(x) does
not play apparent role in the integration, it is still there as
a reluctant item to evaluate.
Now, if you recast your Gamma(x) Γ(x) in your original 7251
what happens bug/no_bug ? Like the English say: square the corner.

Cheers ... Jean
Offline Jean Giraud  
#20 Posted : 12 November 2019 21:23:46(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: Razonar Go to Quoted Post
Hi Nugegodage. I absolutely agree

Thanks Alvaro,
Never Give up. Great things take time.
Users browsing this topic
2 Pages12>
Forum Jump  
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.