ecterrab

14742 Reputation

24 Badges

20 years, 354 days

MaplePrimes Activity


These are replies submitted by ecterrab

@burgarth 

I couldn't download the PhysRevLett paper you mentioned. Would you mind please sending it via email or posting it here? Regarding the question, after fixing the mentioned weakness in the code, I wanted to point out that what you call the boundary terms (double sum with a KroneckerDelta inside) cannot be simplified because that would be wrong, resulting in an additional term that should be equal to 0.

In advance to taking a look at the PhysRevLett paper, I have an idea, mainly around splitting the sum from 1 to d into a sum from 1 to d-1 plus the d term, which in this case is equal to 0. The output would be as you see in my answer above after doing that replacement manually. That would be a good improvement, independent of your example, and would always be mathematically correct.

Now on your questions:

1) Yes, you can tell Maple about the dimension of the related Hilbert space. See the help page ?Physics,Setup, the section on quantumbasisdimension. But the simplifier of sums is not looking there at this moment. From this conversation, it would be good to make the simplifier look there, independent of the solution I mentioned in the previous paragraph.

2) To tell Maple that c[0] = 0, ... the simplest way is to assign it, as in c[0] := 0, c[d+1] = 0, then in the sum defining H use c[n], but again the Simplifier is not looking for 0 at the top of the sum. To indicate they are real, Setup realobjects is correct; also Assume(c[0], real) is also correct. You can always check whether the system understood your assumption with the getassumptions command, e.g. getassumptions(c[0]) after Assume(c[0],real); by the way in the context of Assume, c[0]::real means the same as C[0],real.

I will try the solution I have in mind (second paragraph) and write here again today or tomorrow.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Carl Love 
I see. It seems I read too quickly his post.

Hi @Carl Love 
I see you are answering me. What I said is not that "a moderator deleted" anything, only that "I don't know if this question was accidentally deleted, but I don't see it anymore"; then that "Maybe @nm could say something ?"

By the way, I don't recall @nm deleting one of his own undisputable bug reports. Much less within only 24 hours of posting it. I also have in high regard @nm's questions, and he is also well aware that bugs he reports are always fixed, quickly, as it happened with this one. I don't know what happened here.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

I'm curious about this too. Maybe @nm could say something? I see that his comments in the Unboxing Maple 2022 post are also absent at this moment.

One of the most important mathematical formulations in human history is that of the Standard Model in particle physics. It describes all the elementary particles (leptons like the electron, quarks, bosons as the Higgs or the photon), which in different arrangements, form all the observable particles in nature. The formulation is not just a tremendous theoretical achievement that rendered Nobel prizes but also a practical one. Basically, all the measurements performed in the particle accelerators at CERN and the Fermilab take this mathematical, abstract formulation as the starting point. However, for computer algebra systems, the complexity of the model is somewhat extreme: is not only the number of terms in the corresponding Lagrangian impressively large but also the mathematical properties of each of these objects that represented an insurmountable challenge for a long time. With hacks of different kinds, the representation of only portions of the Standard Model was possible with minimal computational capabilities.

Hidden among the novelties of Maple 2022, a breakthrough in computer algebra is the introduction, for the first time, of a representation for the whole Standard Model. This representation is fully computable, including the accessory commands to calculate related scattering amplitudes  (the essence of the computations behind particle collision experiments) and related Feynman integrals. This is a remarkable achievement in computational physics. And from the educational point of view, it brings one more brick of knowledge from "the dark side" of the moon into"the bright side." Making the Standard Model computations be at the tip of one's fingers completely transforms the possible experience we can have with the underlying knowledge.

This new development is illustrated in the Mapleprimes post The Standard Model of Particle Physics in Maple 2022.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

 

sys := [a*x+y, b*x+y+1, c*x+2*y]

PDEtools:-casesplit(sys, [{x, y}, {a, b, c}])

`casesplit/ans`([y = a/(-a+b), x = 1/(a-b), c = 2*a], [])

(1)

As explained in the help page for casesplit , the input above means: "express{x, y} in terms of {a, b, c}, plus equations that involve only the latter". That is what you request when passing to Groebner:-Basis the argument lexdeg([x, y], [a, b, c]).

 

There are however five important differences (additional features):

1. 

Within casesplit, if you enclose the variables as a set as I did above, casesplit uses its own algorithms to decide which ordering within each set is more convenient.

2. 

casesplit handles, in equal footing, algebraic and or differential equations in the system tackled, that can also contain inequations, and functions of different number of variables.

3. 

casesplit can handle the imaginary unit and mathematical functions, so non-polynomial objects. It does that using the approach explained in dpolyform .

4. 

You have options to include or discard singular solutions, or plot a tree of cases.

5. 

You can choose among three different Maple engines for performing elimination using from Groebner basis to reduced involutive forms and triangularization equivalent to it. These engines are rifsimp , DifferentialAlgebra , and DifferentialThomas .

 

Maple's Groebner basis and RegularChains packages are excellent packages by all means. But I still prefer the simplicity of casesplit's input and output plus the features mentioned above. By the way PDEtools:-casesplit is the main command used by both dsolve and pdsolve to uncouple systems of differential equations, or to compute singular solutions to a single ODE.

Download casesplit.mw

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@nm 
No, it is not expected that simplifying/size checks the syntax of the functions entering an expression; as explained in its help page, it only does a structural analysis of it to shorten its length. Thus, simplify(eq, size) will return the shortening it achieves without any error message. In contrast, simplify(eq) will explicitly call the routine to simplify integrals among others, and the syntax will get checked at some point, resulting in the error message you see.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Kevin Dragnet 

Indeed there is something wrong with the help. It is now tracked in the database of issues, to appear fixed in the upcoming release. Meantime you can access the page directly, either entering ?Mini-Course, or from the Physics help browser colum that appears on the left, as per this image:

 

I will give a look at the question you mentioned about Alias later today.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

Hi
Entering ?Mini-Course I get the right help page (the mini-course). This page is also linked in all the Physics commands' help pages. It might be that one of them has a wrong hyperlink ... Exactly how are you trying to enter this page?

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@lemelinm 

It actually appears exactly as you show. Use the Vectors:-Component for the lhs and the LeviCivita and d_ for the rhs.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Carl Love 
:):):) !!! I sometimes have these gaffes! Indeed the documentation of this package is good, I think but indeed this help page for ChangeCoordinates is just missing! 

The command was introduced in Maple 2015, and it is mentioned there in the page ?updates,Maple2015,Physics, in the section for the Vectors package, also with examples of its use - the basics for its help page. I will move that material into a new help page filling the gap.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

Hi, would you mind please uploading a Maple worksheet with your input/output/expectation? Also, explain please what you do mean by "we want Maple to reconstruct them" Reconstruct what? H? e? What is the meaning of zeta[2] [e 0 , e 1 ]? 

To upload your worksheet please use the Green arrow. Uploading your worksheet helps to help you, avoiding others having to retype or copy and paste things that may not be correct Maple syntax.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@vv 

Indeed, I agree with @Rouben Rostamian, yours is beautiful and inspiring tackle!

Best!

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@snpa 
Would you mind please uploading a worksheet with your question, presented, as in "After having entered this and that (shown) I enter "this" and expect XXX but receive YYY".  In that way, the question is sufficiently clear and helping you becomes simpler. To upload the worksheet please use the green arrow you see when replying or posting a question.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft.

@michalkvasnicka 

This issue in MapleCloud is now fixed, so the problem is resolved.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

First 12 13 14 15 16 17 18 Last Page 14 of 65