Carl Love

Carl Love

28070 Reputation

25 Badges

13 years, 31 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

@nm Yes, I was aware of this when I posted the code. I was waiting for someone, mostly likely you, to notice ;-).

Actually, what stopped me was that in order to fix this properly, I need to know what the degree of the following is supposed to be:

diff(y(x), x) = y(x)^(1/2)  # (1)

Is it 1 or 2? Clearly, the degree of

diff(y(x),x,x) = diff(y(x),x)^(1/2) + y(x) 

is 2. But with #(1), I don't know if the y(x) is to be considered a "zeroeth-order derivative". If yes, then the degree of #(1) is 2; if no, then it's 1.

@ThU Perhaps I'm being naive, but I interpretted the Question to mean that the OP wanted the 2s to be placed in the algebraic exponents, while the minus signs were to remain in the arguments to exp. Your code doesn't separate the minus signs from the 2s. 

If I am indeed being naive about the OP's intent, then your method does seem like a reasonable way to get the desired result if you can get to the expression before it's entered.

@Guy1 To "vote up", you need to click on the thumb icon in the header of the Answer (upper right, under the timestamp). Please do that.

I added some significant details to the Answer, so you should reread it; specifically, I made a distinction between "merger of animations" and "animation of a merger".

I have no idea how or why this Answer got turned into a separate Question. I had nothing to do with it.  Unfortunately, I have no way to put it back where it belongs: an Answer to "Expand and negative exponent".

@Rouben Rostamian  Your choice of file extension, .maple, is a bit unfortunate because it's used for the "meta-worksheets" that were introduced a few versions ago. (These are collections of worksheets, variable assignments, and associated files.) The most commonly used extension for plaintext Maple code is .mpl.

@Oliveira See the help page ?worksheet,reference,initialization

@radaar The command map doesn't do multiprocessing. Why not try gradually increasing the value of numcpus as I suggested? Start with kernelopts(numcpus= 2).

@Thomas Richard Yes, I only used expand for pedagogical reasons. For practical coding, I'd use orthopoly:-T(5,x) to get the polynomial in place of expand(cos(5*arccos(x))).

@radaar Yes, using Grid:-Map will use a lot of memory. From your example worksheet, it doesn't seem like that will be a problem. Just increase the number of processors until you either run out of processors or memory. The memory usage should be roughly a linear function of the number of processors. If it isn't, there's likely some bug that'll require expert attention.

@nm Your Answer is okay for the n=2 case, but has no relevance to the other cases. That's fine, but I want to alert the OP not to try to apply your Answer to the other cases, because doing so could cause confusion.

I don't see how this Answer deserves a vote up, although someone gave it one.

@radaar Use

kernelopts(numcpus= 8);

etc. It may be that you need to issue this command immediately after a restart. Or, maybe that bug has been fixed and you can issue it anytime.

@minhthien2016 Geez, you can't figure that out on your own? For 3-4-5, the solutions are pure rational, so there's no RootOfs, so SolA is a single solution set rather than a sequence of multiple solution sets, so SolA[1] needs to be changed to SolA.

@WmEllis That's great. But I wonder how it originally got set incorrectly.

If you measure "total time" / "number of solutions", your procedure is still much faster.

It should be possible to easily modify your procedure to return the first solution found, if that's wanted.

 

@WmEllis You should set this permanently, because it affects all output, not just matrices. You can do this from the menus: Tools => Options => Display => Output display; select either "Typeset Notation" or "2-D Math Notation"; close the menus with "Apply Globally".

If at any time you need to see output in plaintext ASCII characters (which is very useful for debugging), you can do

lprint(%);

to see just the last output in plaintext; or you can do 

interface(prettyprint= 0); 

to see it in plaintext until you reset it with

interface(prettyprint= 2); #or 3

A closely related setting is 

interface(typesetting= extended);

(My Goodness, why didn't they just make this prettyprint= 4?) In my personal opinion, this extended level is not fully implemented, and it's riddled with bugs. But it is user-adjustable for every type of expression to be output. The opposite of extended is standard.

First 261 262 263 264 265 266 267 Last Page 263 of 709