acer

30767 Reputation

29 Badges

19 years, 6 days
Ontario, Canada

Social Networks and Content at Maplesoft.com

MaplePrimes Activity


These are replies submitted by acer

Hi, sorry for adding this here; my link to contact a Member via message is not working.

Another Question by you (about Oct 2nd) was just deleted. I had an Answer in it relating to Explore. (I found your Explore usage followup goal -- about programmatic access to current parameter values -- so interesting that I wonder whether it might even make for a new Explore option, or at least a Help-page Example.)

Anyway, that whole Question thread has disappeared. Did you delete it, deliberately? If not, I might ask admin to reinstate it.

@Carl Love There are some chunks of right-click popup (or right-panel) menus which are purely GUI-generated.

Much of the remaining right-click popup menu items are pure GUI items, with all the items generated by the ContextMenu package appearing in the right-panel in modern versions.

The ContextMenu package is used to generate most of the menus of actions that appear in the right-panel in modern versions.

But even in the right-panel there are some items (appearing above the CM items, but below SubexpressionMenu items) which are GUI generated/related.

@salim-barzani Please don't spawn yet another Question thread on variation of the Explore use I just showed above.

Instead, you could just put your followup in a Reply here.

@Paras31 Regarding your followup query 1) immediately above, I would first observe that even for,
   plots:-odeplot(sols, [t, U[67](t)], t = 0 .. 10, color = "green")
I get a more convincing and plausible plot if I called dsolve with stepsize=0.5e-5 rather than 0.5e-4.

So I am suspicious that the results even out as far as t=100 might be inaccurate. (I find it difficult to investigate, because it's time consuming. I worry about numeric error effects, limitations of the forced method, and who knows possibly even stroboscopic mirages.)

The smaller forced stepsize has another effect; it pushes out further the time at which that "cannot compute further..." error message might be a problem. But it also makes the computation slower.

On my Maple 2024.1 on Linux I can recover some timing performance here by also passing the compile=true option to dsolve.

So, for fun, you might adjust along these lines: v1_1_1_ac1b.mw

An animation from that, which looks plausible: (This is 200 frames. at 400 frames it's even more convincing. So this seems to be reasonably accurate out at least to t=10.)

ps. There's an old joke about computing. "Speed, robustness, flexibility. Pick two." Here I worry that your situation might be, "Speed, accuracy, ability to compute. Pick one-and-a-half."

Having said that, you haven't actually stated outright 1) whether you think your results were already all accurate (I doubt it), 2) whether you're actually striving for accuracy out for large t, or not, and 3) whether you're forcing rkf45 so that it might be more sensibly compared to solvers in other software, and so on. You haven't provided much in the way of explicit requirements. (For all I know, you might be studying the "apparently chaotic" effects of accrued error here.)

@Paras31 For your query 2) it looks like your two calls to plot U[67] are not the same. The second attempt is not a repeat of the first. It uses a different range.

The first has a range that in 2D Input looks like,
   0...58
which as 2D Input parses to 0 .. 0.58
But your second has,
   0..58
which is the range from 0 to 58.

Please don't spawn yet another wholly separate Question thread on the code shown here.

If you have a followup query (such as on the performance of code given here), then you can ask about that here in a Reply.

Duplicate Question threads get flagged as such and may be deleted.

AFAIK there is no limit (outside of memory and cpu resources).

I have spawned many concurrent commandline-interface (CLI) instances, without problem, on each of Windows, Linux, and OSX.

But, might there not be an optimal number, depending say on the number of cores/(hyper)threads?

@salim-barzani It's not clear what you mean by, "parameter search automatically".

If you want to play around with the effects of the parameter values on the plot, then you could try the Explore command.

graph_in_simple_shape_acc.mw

But I don't see how "parameter search" is meaningful without an explicit goal for which you're searching. What is your goal?

It is unclear what question you might be asking.

@Paras31 It's confusing if your worksheet assign any numeric value to that name.

Your equations appear to contain the name,

   `#mover(mi("δ",fontstyle = "normal"),mo("ˆ"))`

(which pretty-prints like delta with a circumflex accent above it).

But that doesn't seem to be assigned any numeric value. Have I missed that assignment? If not, then what is its value?

@AHSAN You can add any column you want.

What's your explicit definition of the "residual" here? Use that.

@jrive I am going to submit a (likely duplicate) bug report on the following weakness, which characterizes an underlying weakness in simplify.

restart;

Len := expr -> [MmaTranslator:-Mma:-LeafCount(expr),
                length(expr),
                `simplify/size/size`(expr)]:

 

W := (K*(R + Q) + S)/K;

(K*(R+Q)+S)/K

simplify(W);
Len(%);

(K*(R+Q)+S)/K

[11, 33, 27]

simplify(W,size);
Len(%);

(K*(R+Q)+S)/K

[11, 33, 27]

expand(W);
Len(%);

Q+R+S/K

[8, 21, 21]

Download simp_weakness.mw

@jrive If F (the denominator, say) is obtained as the common multiplicative term that one wants to distribute through the sum (numerator), then I've showed that either of the following attain your target,

   map(`/`, Zin*F, F)
or,
   expand(Zin, F, op(Zin*F))


I picked off F as the whole term in the denominator. And of course that involved no up front check that it were indeed a "common" (multiplicative) factor that could be usefully cancelled via distribution through the numerator sum.

ps. The reason that I originally went with collect is that it sometimes happens that a judicious collect call gets a more compact result than does simplify(...,size). I have a hard-working compactifier oracle that looks for such, by trying different variable choices and followup coefficient-action (3nd arg).

The form of your target expression (and my interest in finding a shorter/shortest form) made me think trying collect first. In contrast, the form of your input expression leads more to using a distributed multiplication or a controlled expansion -- but one has to "see" that it'd work.

1 2 3 4 5 6 7 Last Page 1 of 564