Line intersections

#21 Pubblicato: 16/02/2018 09:01:36

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Wrote
Just an exercise on Smath to test functionalities and capabilities.
Try the attached file and play with it.

Thanks
Marco
InterX.sm (85,93 KiB) scaricato 1441 volta(e).

Hi Marco. Here some comments about your great file. As you suggest, I play with it.
Best regards.
Alvaro.

InterX.r1.sm (77,75 KiB) scaricato 1281 volta(e).

1 utenti hanno apprezzato questo post

Andrey Ivashov 16/02/2018 12:30:00

#22 Pubblicato: 16/02/2018 11:55:14

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Wrote
Try the attached file and play with it.

1. Marco => still computing after 15 min ? !
2. Alvaro => red from top to bottom

#23 Pubblicato: 16/02/2018 12:24:05

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Wrote
Try the attached file and play with it.

Attached is about timing in work sheet.
Observe the sequential timing [black].
t0:=time(1) is in collapsed area.

Contour_8 Hypersurface.sm (1,17 MiB) scaricato 1272 volta(e).

#24 Pubblicato: 16/02/2018 12:31:14

Andrey Ivashov

2.343 apprezzamenti in 4.056 post.

Gruppo: Super Administrator

Wrote
Hi Marco. Here some comments about your great file. As you suggest, I play with it.
Best regards.
Alvaro.

1 min 53 sec on my windows tablet pc.

#25 Pubblicato: 16/02/2018 13:09:46

Fridel Selitsky

520 apprezzamenti in 451 post.

Gruppo: User

Wrote

1 min 53 sec on my windows tablet pc.

SMath Studio 0.99.6611 time(0)-to =83.735 s on my windows10

#26 Pubblicato: 16/02/2018 13:31:34

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: User

Wrote
Originally Posted by: Andrey Ivashov [url=/forum/yaf_postsm47126_Line-intersections.aspx#post47126][/url]

1 min 53 sec on my windows tablet pc.
SMath Studio 0.99.6611time(0)-to =83.735 son my windows10

Time taken in my win 10 laptop with SS 0.99.6619 shown.

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

#27 Pubblicato: 16/02/2018 13:54:42

Davide Carpi

1.443 apprezzamenti in 2.936 post.

Gruppo: Moderator

5 min 25 s (win 10 x64, intel i5 quad core 2.40GHz)

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

#28 Pubblicato: 16/02/2018 15:07:54

Davide Carpi

1.443 apprezzamenti in 2.936 post.

Gruppo: Moderator

Wrote
5 min 25 s (win 10 x64, intel i5 quad core 2.40GHz)

reduced to 30-35 seconds using ranged indices to replace stack/augment/submatrix (most heavy operations were in stack/augment under loops)

InterX.r1 - ranged indices.sm (87,61 KiB) scaricato 1248 volta(e).

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

1 utenti hanno apprezzato questo post

Andrey Ivashov 16/02/2018 15:13:00

#29 Pubblicato: 16/02/2018 15:11:46

Andrey Ivashov

2.343 apprezzamenti in 4.056 post.

Gruppo: Super Administrator

18 seconds on my env.

Great! Thank you Davide!

#30 Pubblicato: 16/02/2018 15:49:26

Andrey Ivashov

2.343 apprezzamenti in 4.056 post.

Gruppo: Super Administrator

BTW, I would suggest to use another version of plotG which is much faster then previous one (and simplier):

plotG(x,y,char,size,color):line(k:range(1,rows(x)),augment(x,y,el(r3,k):char,el(r4,k):size,el(r5,k):color),2,1)

plotG(xy,char,size,color):line(k:range(1,rows(xy)),augment(xy,el(r3,k):char,el(r4,k):size,el(r5,k):color),2,1)

Best regards.

1 utenti hanno apprezzato questo post

Davide Carpi 16/02/2018 20:54:00

#31 Pubblicato: 16/02/2018 16:18:55

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Hi. Just one more thing, I guess: the only one remaining matlab function, mwfind replaced. (Orignal work have many) If I'm not wrong, it is now 100% SMath native. But double the amazing time obtained by Davide. But it's just a "scalar" implementation, sure there are better ones.

Best regards.

Alvaro.

PD: I don't replace plotG because don't understand very well that syntax.

InterX.r1 - ranged indices.sm (90,75 KiB) scaricato 1195 volta(e).

#32 Pubblicato: 16/02/2018 16:45:48

Andrey Ivashov

2.343 apprezzamenti in 4.056 post.

Gruppo: Super Administrator

Wrote
I don't replace plotG because don't understand very well that syntax.

It is simple - ranged indexes used.

Just take a look at the speed:

Best regards.

#33 Pubblicato: 16/02/2018 16:59:17

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Wrote

Wrote
I don't replace plotG because don't understand very well that syntax.

It is simple - ranged indexes used.

Thanks for the explanation, now I see the speed increment, and I don't think that it could be that big difference before you show it.

But actually the plotG(4) routine was which confuse me. I prefer

plotG(xy,char,size,color):if(cols(xy)>1,plotG(col(xy,1),col(xy,2),char,size,color),plotG(Re(xy),Im(xy),char,size,color))

because I think that any good plot implementation must to use complexes but not couple of reals. So, if xy argument have only one column, one can assume that data are as Z = X + iY.

Best regards.

Alvaro.

#34 Pubblicato: 16/02/2018 17:08:15

Davide Carpi

1.443 apprezzamenti in 2.936 post.

Gruppo: Moderator

Citazione
because I think that any good plot implementation must to use complexes but not couple of reals. So, if xy argument have only one column, one can assume that data are as Z = X + iY.

Andrey misses this detail of your function about the complex values, but you can still do it; keep the if statement and replace the loops (the case with more than a column as shown by Andrey doesn't require the col() functions, this is redundant and time consuming... If you think that more than 2 columns might be provided, add the case and a single sumbatrix instead

)

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

1 utenti hanno apprezzato questo post

Alvaro 16/02/2018 17:45:00

#35 Pubblicato: 16/02/2018 18:24:21

Arie

94 apprezzamenti in 285 post.

Gruppo: User

Wrote
But double the amazing time obtained by Davide. But it's just a "scalar" implementation, sure there are better ones.

Here's my pure SMATH find function. with similar time to MWFIND.

InterX.r1 - ranged indices_alyles.sm (94,58 KiB) scaricato 1238 volta(e).

Feel free to join the SMath Studio Users Discord Channel: https://discord.gg/PayZpJW

4 utenti hanno apprezzato questo post

Andrey Ivashov 16/02/2018 18:41:00, sergio 17/02/2018 06:31:00, Davide Carpi 16/02/2018 20:57:00, Alvaro 16/02/2018 22:03:00

#36 Pubblicato: 17/02/2018 09:35:24

Martin Kraska

1.361 apprezzamenti in 2.275 post.

Gruppo: Moderator

Wrote
BTW, I would suggest to use another version of plotG which is much faster then previous one (and simplier):
plotG(x,y,char,size,color):line(k:range(1,rows(x)),augment(x,y,el(r3,k):char,el(r4,k):size,el(r5,k):color),2,1)
plotG(xy,char,size,color):line(k:range(1,rows(xy)),augment(xy,el(r3,k):char,el(r4,k):size,el(r5,k):color),2,1)
Best regards.

For such uses it would be good if the augment and stack functions could expand scalar arguments (or more generally any arguments not matching the maximum number of rows or cols) to uniformly populated vectors of the required size).

Then the above example wouldn't even need a helper function as you could write

augment(x,y,"*",10,"red"

Technische Mechanik mit SMath Studio: https://link.springer.com/book/10.1007/978-3-658-50592-9

2 utenti hanno apprezzato questo post

Davide Carpi 17/02/2018 15:23:00, Radovan Omorjan 17/02/2018 12:43:00

#37 Pubblicato: 17/02/2018 09:39:40

Andrey Ivashov

2.343 apprezzamenti in 4.056 post.

Gruppo: Super Administrator

Very good suggestion, Martin! I think I could do this.

Best regards.

4 utenti hanno apprezzato questo post

Davide Carpi 17/02/2018 15:23:00, frapuano 17/02/2018 10:20:00, Radovan Omorjan 17/02/2018 12:43:00, Martin Kraska 17/02/2018 09:41:00

#38 Pubblicato: 17/02/2018 13:20:11

Radovan Omorjan

325 apprezzamenti e 1 disapprezzamenti in 2.052 post.

Gruppo: Moderator

Wrote
For such uses it would be good if the augment and stack functions could expand scalar arguments (or more generally any arguments not matching the maximum number of rows or cols) to uniformly populated vectors of the required size).

Then the above example wouldn't even need a helper function as you could write
augment(x,y,"*",10,"red"

I also vote for that augment() and stack () functionality.

In that case, we could even have "plot" functions (with the same functionality as augment())


plot(x,y)
plot(x,y,"*"
plot(x,y,"*",10)
plot(x,y,"*",10,"red"

This efforts are, of course, not justified if the micro Mathematics plot style would be introduced.

Regards,
Radovan

When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"

1 utenti hanno apprezzato questo post

frapuano 19/02/2018 07:14:00

#39 Pubblicato: 18/02/2018 22:58:48

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Hello Marco,

I'm trying for the intersection of the two Trifolium, .
blue => native ... red => rotated.
They are created in "Page1 Trifolium.sm
Then recalled for the InterX(T1,T2)
Nothing happens >>> 1 hr ? [Win 7 Starter Edition ... 1.66 GHz]
Is it supposed to work ? Thanks in advance.

Page1 Trifolium.sm (32,5 KiB) scaricato 1482 volta(e).

Maths InterseX [Marco].sm (51,7 KiB) scaricato 1011 volta(e).

#40 Pubblicato: 19/02/2018 00:26:16

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Hi Jean.

You are using an old version of the program. Try this one.

Notice that the point (0,0) is found several times by the algorithm. It could be interesting to check if it get the correct number of crosses at the origin.

Best regards.

Alvaro.

InterX.r1 - ranged indices_alyles.sm (88,24 KiB) scaricato 1195 volta(e).

Forum SMath

Line intersections - An example of translation of a function from Matlab to Smath - Messaggi