Line intersections

#21 Publicado: 16/2/2018 9:01:36

1.139 me gusta en 1.805 mensajes.

Grupo: 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) descargado 1437 vez(es).

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) descargado 1276 vez(es).

A 1 usuarios les gustó este mensaje

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

#22 Publicado: 16/2/2018 11:55:14

Jean Giraud

983 me gusta en 6.866 mensajes.

Grupo: 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 Publicado: 16/2/2018 12:24:05

Jean Giraud

983 me gusta en 6.866 mensajes.

Grupo: 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) descargado 1269 vez(es).

#24 Publicado: 16/2/2018 12:31:14

Andrey Ivashov

2.343 me gusta en 4.056 mensajes.

Grupo: 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 Publicado: 16/2/2018 13:09:46

Fridel Selitsky

520 me gusta en 451 mensajes.

Grupo: 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 Publicado: 16/2/2018 13:31:34

NDTM Amarasekera

130 me gusta en 352 mensajes.

Grupo: 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 Publicado: 16/2/2018 13:54:42

Davide Carpi

1.443 me gusta en 2.936 mensajes.

Grupo: 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 Publicado: 16/2/2018 15:07:54

Davide Carpi

1.443 me gusta en 2.936 mensajes.

Grupo: 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) descargado 1243 vez(es).

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

A 1 usuarios les gustó este mensaje

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

#29 Publicado: 16/2/2018 15:11:46

Andrey Ivashov

2.343 me gusta en 4.056 mensajes.

Grupo: Super Administrator

18 seconds on my env.

Great! Thank you Davide!

#30 Publicado: 16/2/2018 15:49:26

Andrey Ivashov

2.343 me gusta en 4.056 mensajes.

Grupo: 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.

A 1 usuarios les gustó este mensaje

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

#31 Publicado: 16/2/2018 16:18:55

Alvaro

1.139 me gusta en 1.805 mensajes.

Grupo: 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) descargado 1191 vez(es).

#32 Publicado: 16/2/2018 16:45:48

Andrey Ivashov

2.343 me gusta en 4.056 mensajes.

Grupo: 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 Publicado: 16/2/2018 16:59:17

Alvaro

1.139 me gusta en 1.805 mensajes.

Grupo: 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 Publicado: 16/2/2018 17:08:15

Davide Carpi

1.443 me gusta en 2.936 mensajes.

Grupo: Moderator

Citar
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

A 1 usuarios les gustó este mensaje

Alvaro 16/2/2018 17:45:00

#35 Publicado: 16/2/2018 18:24:21

Arie

94 me gusta en 285 mensajes.

Grupo: 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) descargado 1235 vez(es).

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

A 4 usuarios les gustó este mensaje

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

#36 Publicado: 17/2/2018 9:35:24

Martin Kraska

1.361 me gusta en 2.275 mensajes.

Grupo: 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

A 2 usuarios les gustó este mensaje

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

#37 Publicado: 17/2/2018 9:39:40

Andrey Ivashov

2.343 me gusta en 4.056 mensajes.

Grupo: Super Administrator

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

Best regards.

A 4 usuarios les gustó este mensaje

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

#38 Publicado: 17/2/2018 13:20:11

Radovan Omorjan

325 me gusta y 1 no me gusta en 2.052 mensajes.

Grupo: 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!"

A 1 usuarios les gustó este mensaje

frapuano 19/2/2018 7:14:00

#39 Publicado: 18/2/2018 22:58:48

Jean Giraud

983 me gusta en 6.866 mensajes.

Grupo: 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) descargado 1475 vez(es).

Maths InterseX [Marco].sm (51,7 KiB) descargado 1008 vez(es).

#40 Publicado: 19/2/2018 0:26:16

Alvaro

1.139 me gusta en 1.805 mensajes.

Grupo: 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) descargado 1192 vez(es).

Foro SMath

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