Different Results of Integration in 7219 & 7250 (Gamma Function)

#1 Опубликовано: 08.11.2019 05:08:51

130 сообщений из 352 понравились пользователям.

Группа: User

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 (111 КиБ) скачан 45 раз(а).

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

#2 Опубликовано: 08.11.2019 06:13:29

Davide Carpi

1415 сообщений из 2872 понравились пользователям.

Группа: Moderator

Do you have the same settings about integrals' accuracy in both the versions?

Файл не найден. Файл не найден.

If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects

#3 Опубликовано: 08.11.2019 08:20:08

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

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

#4 Опубликовано: 08.11.2019 12:15:57

Andrey Ivashov

2269 сообщений из 3729 понравились пользователям.

Группа: Super Administrator

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/forum/yaf_postst17461_Beta--SMath-Studio-0-99-7251--08-November-2019.aspx

Best regards, Andrey Ivashov.

1 пользователям понравился этот пост

Martin_B 08.11.2019 13:54:00

#5 Опубликовано: 10.11.2019 00:15:32

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

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

#6 Опубликовано: 10.11.2019 03:59:21

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

Integration results of Mathcad attached (using first 10 elements of X
and Eqns 1, 2 and 3.)

SPI_Calc_Plain.sm (125 КиБ) скачан 37 раз(а).

Rainfal_Data.xlsx (486 КиБ) скачан 29 раз(а).

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

#7 Опубликовано: 10.11.2019 11:01:56

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Wrote
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 (54 КиБ) скачан 31 раз(а).

#8 Опубликовано: 11.11.2019 02:31:56

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

Wrote
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 (35 КиБ) скачан 41 раз(а).

Test_Gamma_7251.pdf (192 КиБ) скачан 38 раз(а).

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

#9 Опубликовано: 11.11.2019 10:59:00

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Wrote
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 (45 КиБ) скачан 31 раз(а).

#10 Опубликовано: 11.11.2019 13:14:35

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

... max PDF located

NDTMA Test Gamma.sm (22 КиБ) скачан 37 раз(а).
Golden Ratio jmG.sm (35 КиБ) скачан 31 раз(а).

#11 Опубликовано: 12.11.2019 07:07:16

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

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 МиБ) скачан 33 раз(а).

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

#12 Опубликовано: 12.11.2019 10:16:46

Alvaro Diaz Falconi

992 сообщений из 1674 понравились пользователям.

Группа: User

Hi. You can try to implement your own algo's for gamma and int functions.

Test_Gamma_7251.sm (103 КиБ) скачан 38 раз(а).

Best regards.
Alvaro.

2 пользователям понравился этот пост

frapuano 12.11.2019 10:31:00, NDTM Amarasekera 12.11.2019 11:54:00

#13 Опубликовано: 12.11.2019 11:53:09

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

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.

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

#14 Опубликовано: 12.11.2019 12:07:05

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Wrote
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 (22 КиБ) скачан 41 раз(а).

#15 Опубликовано: 12.11.2019 13:39:00

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

Wrote
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 (25 КиБ) скачан 31 раз(а).

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

#16 Опубликовано: 12.11.2019 13:43:07

Alvaro Diaz Falconi

992 сообщений из 1674 понравились пользователям.

Группа: User

Wrote

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.

#17 Опубликовано: 12.11.2019 14:32:37

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

Wrote
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

#18 Опубликовано: 12.11.2019 14:52:43

NDTM Amarasekera

130 сообщений из 352 понравились пользователям.

Группа: User

This confirms Jean's guess about vectorize operator.

TEST-SPI_Calc_Plain.sm (52 КиБ) скачан 32 раз(а).

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

#19 Опубликовано: 12.11.2019 15:17:55

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Wrote
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

#20 Опубликовано: 12.11.2019 15:23:46

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Wrote
Hi Nugegodage. I absolutely agree

Thanks Alvaro,
Never Give up. Great things take time.

Форум SMath

Different Results of Integration in 7219 & 7250 (Gamma Function) - Сообщения