@Joe Riel Click on the text which appears not to be an active link. That'll take you to the post. There's another example of how the controls are not highlighted.

The new Notifications flag is disfunctional in several ways:

1. The notifications only apply to posts to which I am email subscribed. I don't email subscribe to anything because I don't want all that email! What I want a notification about is when someone replies to a thread that I've participated in. And I want that flag to turn orange automatically when I have a response, not requiring a page refresh!

2. The flag ceases to be orange after it has been opened, regardless of whether I have checked or dismissed the threads.

Please look at and study any StackExchange forum. Notice the vast amount of information that's on one screen? Notice the vast number of flags that update automatically? Anything that you can do to make the look and feel similar to StackExchange would be a great improvement.

@vv The randomize test that I concocted above proves beyond any reasonable doubt that the code intentionally makes a call to Maple's random number generator. So, whether some unintentional randomization is also occurring is irrelevant.

@Joe Riel Your getparams code is missing the closing )] on the seq command.

@annarita It can be done by doing a parametric plot3d for each of the six faces and then assembling the results with plots:-display:

Pz0:= plot3d([x,y,0], x= 0..1, y= 0..1, color= x^2+y^2-0^2):
Pz1:= plot3d([x,y,1], x= 0..1, y= 0..1, color= x^2+y^2-1^2):
Px0:= plot3d([0,y,z], y= 0..1, z= 0..1, color= 0^2+y^2-z^2):
Px1:= plot3d([1,y,z], y= 0..1, z= 0..1, color= 1^2+y^2-z^2):
Py0:= plot3d([x,0,z], x= 0..1, z= 0..1, color= x^2+0^2-z^2):
Py1:= plot3d([x,1,z], x= 0..1, z= 0..1, color= x^2+1^2-z^2):
plots:-display(P||(x,y,z)||(0,1));

This is the basic plot on the (one-dimensional) HUE color scale (the visual spectrum). Many variations are possible using the three-dimensional RGB or HSV color spaces. Additional options to control the grid, transparency, etc., can be added to the end of the display command.

@tomleslie Please reread what I wrote very carefully. It's only three lines. You are completely misinterpretting it. I said that in most modern languages a line break is treated the same as a space. I did not say that in most modern languages a space is treated as multiplication. Indeed, the only serious (*) general-purpose languages that I can recall where line break and space are treated differently are Fortran, Basic, Cobol, and Python and its descendents. (I'm sure that I'm missing several. I used to use Foxpro / dBase a lot, but I can't recall how the space-vs-line-break issue is handled there.) And of course I mean syntactic uses, not uses where the character---line break or space---is in a string.

(*) There's a non-serious language called Whitespace. You can look it up on Wikipedia.

@tomleslie In Maple (and I think most moderm languages), a line break is parsed just liked a space. When code using your indentation style is copy-and-pasted as 2-D Input, there is a line break, hence a space, hence an implicit multiplication, between the procedure names and their following left parenthesis.

@annarita What version of Maple are you using? You seem to "know" about the colorscheme option, yet your version of Maple doesn't seem to support it. When I run your worksheet, my version of the command produces the plot without error.

@Christian Wolinski Your use of @ contributes to the simplicity of the code. However, the terminal punctuation definitely needs to be a comma, not a semicolon. Your version of Maple doesn't use initializers in variable declarations. The following two procedures are different:

proc() local x:= 1, y:= 2; end proc:

proc() local x:= 1; y:= 2 end proc:

In the first, y is explicitly declared local. With the second, y is only implicitly declared local, with a warning message.

@vv Oh, your Answer is in Reply to YasH's Reply about Eigenvectors rather than being an answer to the OP's Question. It is confusing that you (or someone) classified it as an Answer.

@vv Has the posted code been changed? I can find no use of indices in either of the posted worksheets. I searched with the GUI's Find command (from the Edit menu). I searched both inside and outside of the Code Edit Regions.

@YasH Your test with randomize is very strong evidence that the variability of the order of the returned eigenvectors is not due to the code's use of some type of randomization. So, the variability must be due to some address-based storage issue, such as what's discussed about the indices command above. I use indices in nearly all of my programs where the order of the returned results doesn't matter because it's an extremely efficient command, so I'd guess that it's widely used in the library code also.

@vv Yes, indices are stored by some address scheme, and hence the order is session dependent. If you want the results of indices or entries to be returned in a consistent order, then use option indexorder:

indices(MyTable, indexorder);

and likewise for entries. The order will then be the same as if the indices (not the entries) were an ordinary Maple set.

@YasH That's a good question, and, yes, it's related. I don't know the answer, but it's not too difficult to figure out with a few experiments. For example, after a restart, do the four always come in the same order? If yes, that's a sign that random numbers are used in the algorithm. If no, that's a sign that the ordering is based on the addresses of certain objects, most likely vectors.

The idea that random numbers might be used in the algorithm can be further investigated by checking if the pre-determined sequence of random numbers changes after calling Eigenvectors.

@vv Vote up. My tests below show that a simple for-loop approach, taking the remainder at each step, is fastest, at least for word-sized integers.

In a post a few months ago, I showed that the divide-and-conquer approach was much faster than plain mul or a linear for-loop for plain multiplication. Apparently, plain mul uses a linear simple-loop approach. This probably prevents the use of the much faster FFT-based multiplication algorithm (which is only used when both factors are supra-word-length integers).

Taking the remaindering into account, the divide-and-conquer approach is still much faster than mul, but it doesn't beat the simple for loop. ModMul.mw

Download ModMul.mw