ecterrab

14727 Reputation

24 Badges

20 years, 330 days

MaplePrimes Activity


These are replies submitted by ecterrab

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

@itsme 

Thanks for your new post; the problem was inadvertently introduced in March during one more round of improvements in SortProducts, and is fixed in the Maplesoft Physics Updates v.1015 or newer (for Maple 2021), uploaded earlier today. Attached is what I receive, with all the terms correctly sorted.

physics_sorting_products_(reviewed).mw

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

@nm 

This kind of problem is tackled using differential algebra elimination. There are three computational implementations of that: DEtools:-rifsimp, DifferentialAlgebra and, more recently, DifferentialThomas. The three use different algorithms to accomplish, basically, the same. The algorithms are extended to handle mathematical functions using the approach explained in PDEtools:-dpolyform; i.e.: representing the (non-rational) mathematical functions by (rational) differential equations and auxiliary variables.

This approach is extremely powerful but, depending on the problem, it may hang. The example you are asking about, however, does not hang if a) you tackle it using DifferentialThomas (256 seconds) or more appropriately b) you use convert/exp (0.6 seconds):

Since the output of independentof could be expressed with any mathematical functions, the convert/exp step should be applied by default, directly within PDEtools:-Solve. I will test this idea and if it works as I expect put it in one of the next Maplesoft Physics Updates.

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

@ArashMhasani 

You are using a program you wrote, Ldiff, that I copy and paste here:

Ldiff := proc(f, h, x)
    local n, j, ld;
    if not type(f, 'vector') then
        ERROR("ldiff: 1st argument must be a vector")
    else
        n := linalg[vectdim](f)
    fi;
    if type(h, 'table') then
        if not (linalg[vectdim](h) = 1) then
            ERROR("2nd argument must be a scalar valued function")
        fi
    fi;
    if not (type(x, 'vector') and linalg[vectdim](x) = n) then
        lprint("ldiff: 3th argument is the variable vector and must be \
            defined as ");
        lprint("       an unassigned vector, wrong type or number of el\
            ements, ");
        lprint("       see your model description");
        lprint("  ");
        ERROR("see last comment")
    fi;
    ld := 0;
    for j to n do
        ld := ld + diff(h, x[j])*f[j]
    od;
    simplify(normal(ld))
end

I noticed that, removing the call to simplify, it takes 263 seconds to compute Ldiff(f,z[6],x), the input you say in your worksheet that hangs, resulting on an expression of length 167189212 (see ?length). If you also remove the call to normal,  Ldiff(f,z[6],x) is computed in less than 1 sec and results in an expression of length 846930. It is then unclear the advantage of calling normal. Also, simplify includes normalization, so I would suggest you to just remove the call to normal, and in this particular case of Ldiff(f,z[6],x) remove also the call to simplify.

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

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