Welcome Guest! To enable all features please Login or Register.

Notification

Icon
Error

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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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 13 time(s).

7219_CDF Gamma.png

7250_CDF Gamma.png
Never Give up. Great things take time.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
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,213
Man
Italy
Location: Italy

Was thanked: 1120 time(s) in 731 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 please consider a donation to SMath Studio; for personal contributions to me: 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: 216
Man
Sri Lanka
Location: Colombo

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

7219_Settings.png
7250_Integral Accuracy.png
Never Give up. Great things take time.
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,382
Man
Russian Federation

Was thanked: 1619 time(s) in 542 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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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.
Never Give up. Great things take time.
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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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 8 time(s).

CDF MathCAD.png

7251_CDFpng.png

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

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

Never Give up. Great things take time.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#7 Posted : 10 November 2019 17:01:56(UTC)
Jean Giraud


Rank: Advanced Member

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 4,060
Canada

Was thanked: 714 time(s) in 558 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 6 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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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 7 time(s).

Test_Gamma_7251.pdf (193kb) downloaded 5 time(s).
Never Give up. Great things take time.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#9 Posted : 11 November 2019 16:59:00(UTC)
Jean Giraud


Rank: Advanced Member

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 4,060
Canada

Was thanked: 714 time(s) in 558 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 2 time(s).
Offline Jean Giraud  
#10 Posted : 11 November 2019 19:14:35(UTC)
Jean Giraud


Rank: Advanced Member

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 4,060
Canada

Was thanked: 714 time(s) in 558 post(s)
... max PDF located

NDTMA Test Gamma.sm (23kb) downloaded 4 time(s).
Golden Ratio jmG.sm (36kb) downloaded 4 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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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 4 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

Never Give up. Great things take time.
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: 405
Uruguay

Was thanked: 218 time(s) in 130 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 8 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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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

Never Give up. Great things take time.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#14 Posted : 12 November 2019 18:07:05(UTC)
Jean Giraud


Rank: Advanced Member

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 4,060
Canada

Was thanked: 714 time(s) in 558 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 5 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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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 2 time(s).

Test Gamma [2] 7251 [6179].png
Never Give up. Great things take time.
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: 405
Uruguay

Was thanked: 218 time(s) in 130 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: 216
Man
Sri Lanka
Location: Colombo

Was thanked: 67 time(s) in 45 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.
Never Give up. Great things take time.
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: 216
Man
Sri Lanka
Location: Colombo

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

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

7251 Vectorize Error.png
Never Give up. Great things take time.
Best Regards
Eng. NDTM Amarasekera - Sri Lanka
Offline Jean Giraud  
#19 Posted : 12 November 2019 21:17:55(UTC)
Jean Giraud


Rank: Advanced Member

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 4,060
Canada

Was thanked: 714 time(s) in 558 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: Advanced Member

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 4,060
Canada

Was thanked: 714 time(s) in 558 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.