Rank: Advanced Member Groups: Registered, Advanced Member Joined: 06/04/2023(UTC) Posts: 345 Was thanked: 26 time(s) in 25 post(s)
|
The work with not quadrat matrix. How better solve? See please the attached sm-file. 1-Asteroid-5-Points.sm (26kb) downloaded 14 time(s).Edited by user 01 November 2023 19:02:45(UTC)
| Reason: Not specified
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
Originally Posted by: OchkovVF The work with not quadrat matrix. How better solve? As it looks, the blue planet and green voyager are tangent. Take care Valery ... Jean. 1-Asteroid-5-Points.sm (35kb) downloaded 2 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,361 Was thanked: 822 time(s) in 520 post(s)
|
Hi Valery. You can use dn_LinAlgLLS_SVD or dn_LinAlgLLS_COF or dn_LinAlgLLS_QRorLQ for non-square matrices, instead dn_LinAlgSolve. 1-Asteroid-5-Points.sm (28kb) downloaded 3 time(s).Best regards. Alvaro.
|
1 user thanked Razonar for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 06/04/2023(UTC) Posts: 345 Was thanked: 26 time(s) in 25 post(s)
|
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,361 Was thanked: 822 time(s) in 520 post(s)
|
A little more compact. Best regards. Alvaro.
|
1 user thanked Razonar for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 06/04/2023(UTC) Posts: 345 Was thanked: 26 time(s) in 25 post(s)
|
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,361 Was thanked: 822 time(s) in 520 post(s)
|
Originally Posted by: OchkovVF ... New problem with Ranks Hi Valery. SMath really needs to improve on linear algebra. As a workaround, you can redefine rank with a plugin Also, the problem with rank is not that it can't deal with complicated arrays, it can't handle simple ones either 10-2-Asteroid-7-points-Rank.sm (20kb) downloaded 5 time(s).Best regards- Alvaro.
|
2 users thanked Razonar for this useful post.
|
on 02/11/2023(UTC), on 03/11/2023(UTC)
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 06/04/2023(UTC) Posts: 345 Was thanked: 26 time(s) in 25 post(s)
|
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,361 Was thanked: 822 time(s) in 520 post(s)
|
Originally Posted by: OchkovVF Why? Hi. Due to numerical errors. These ranks are evaluated with numerical procedures. That means that for some cases they must necessarily fail. You can introduce some kind of tolerance in the function. Maple also fails in this example. This shows that the rank function that is executed is that of maple. However, I don't know how to call Maxima's rank function, it is not in the Takeover options. Here some notes about how to handle a little better the issue. 10-2-Asteroid-7-points-Rank.sm (36kb) downloaded 4 time(s).Best regards. Alvaro. Edited by user 03 November 2023 07:52:42(UTC)
| Reason: SMath with notes added.
|
1 user thanked Razonar for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Maxima gives the correct answer. Yet I think that maple is not to blame, because maple(rank(M)) sends 4 to maple. At least this happens with Maxima(rank(M)). In Maxima I have a workaround for such cases. First I convey the definition of M to Maxima and then call rank(M). The format "$...§" conveys everything inside the markup as verbatim expression to Maxima. without giving SMath the chance to interfere. Should not be too complicated to add rank() to the takeover candidates. 10-2-Asteroid-7-points-Rank _Kr.sm (41kb) downloaded 3 time(s). |
|
2 users thanked mkraska for this useful post.
|
on 04/11/2023(UTC), on 04/11/2023(UTC)
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,361 Was thanked: 822 time(s) in 520 post(s)
|
Hi Martin. Originally Posted by: mkraska Maxima gives the correct answer ... Nope. The issue is about a matrix C=augment( A,b ), where b is a column vector so, it must to have the same rank as A. Originally Posted by: mkraska ... Yet I think that maple is not to blame, because maple(rank(M)) sends 4 to maple. At least this happens with Maxima(rank(M)). Either. In some way maple holds the rank function and evaluates it on its own kernel. Actually, the example with the matrix M is for show that. Best regards. Alvaro.
|
1 user thanked Razonar for this useful post.
|
|
|
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.