SMath Studio Forum
»
SMath Studio
»
Bugs & Problems
»
Different Results of Integration in 7219 & 7250 (Gamma Function)
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 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). |
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka |
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 13/01/2012(UTC) Posts: 2,639 Location: Italy Was thanked: 1323 time(s) in 873 post(s)
|
Do you have the same settings about integrals' accuracy in both the versions? |
If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 Location: Colombo Was thanked: 124 time(s) in 82 post(s)
|
Yes, same settings in both versions. But the difference is in Eqn 2 |
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka |
|
|
|
Rank: Administration Groups: Developers, Registered, Knovel Developers, Administrators, Advanced Member Joined: 11/07/2008(UTC) Posts: 1,616 Was thanked: 1978 time(s) in 666 post(s)
|
|
1 user thanked Andrey Ivashov for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 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 |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 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). 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 |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: ndtma 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.sm (55kb) downloaded 17 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 Location: Colombo Was thanked: 124 time(s) in 82 post(s)
|
Originally Posted by: Jean Giraud 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 |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: ndtma 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 Gamma.sm (46kb) downloaded 12 time(s).
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 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).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 |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,356 Was thanked: 815 time(s) in 516 post(s)
|
Hi. You can try to implement your own algo's for gamma and int functions. Test_Gamma_7251.sm (104kb) downloaded 20 time(s).Best regards. Alvaro.
|
2 users thanked Razonar for this useful post.
|
on 12/11/2019(UTC), on 12/11/2019(UTC)
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 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 |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: ndtma 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).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 Location: Colombo Was thanked: 124 time(s) in 82 post(s)
|
Originally Posted by: Jean Giraud 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). |
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,356 Was thanked: 815 time(s) in 516 post(s)
|
Originally Posted by: ndtma
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.
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 Location: Colombo Was thanked: 124 time(s) in 82 post(s)
|
Originally Posted by: Jean Giraud 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 |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 05/06/2014(UTC) Posts: 348 Location: Colombo Was thanked: 124 time(s) in 82 post(s)
|
|
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka |
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: ndtma 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
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: Razonar Hi Nugegodage. I absolutely agree Thanks Alvaro, Never Give up. Great things take time.
|
|
|
|
SMath Studio Forum
»
SMath Studio
»
Bugs & Problems
»
Different Results of Integration in 7219 & 7250 (Gamma Function)
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.