ecterrab

14540 Reputation

24 Badges

20 years, 21 days

MaplePrimes Activity


These are answers submitted by ecterrab

You know, tetrads are defined up to an arbitrary Lorentz rotation of the tetrad system, so naturally, there are infinitely many ways of writing a tetrad correctly. Equation (27.22) in Stephani's book shows this:

So in order to help you with this, or fix something in the code if necessary, would you mind please to write, in the worksheet you attached, the null tetrad matrix you expect, as a matrix, and post the worksheet again.

Please note: after you write the matrix - say M - that is the tetrad you expect, before reposting your worksheet input this:

Tetrads:-IsTetrad(M);
                                       Type of tetrad: xxx
                                                    true

So I expect true, indicating that the matrix M is indeed a tetrad, and xxx indicating what type of tetrad is it, orthonormal or null. If you don't get true but still think your M is a tetrad, please post anyway and I will give it a second look.

Independent of all the above, note that: a) the database of solutions is accessible regarding the solutions, not the tetrads, presented or not in the book, and b) you do have, in the Tetrads package, several commands to depart from one tetrad and obtain another one, notably TransformTetrad.

Not only TransformTetrad can perform all the different transformations for tetrads mentioned in Stephani's book, but more: it can, automatically, compute a canonicalform for the tetrad, so that the Weyl scalars appear with a certain form, all this explained in ?TransformTetrad. This help page, by the way, is not a standard one but something close to a textbook section on the matter. The topic is advanced. I intended to present a self-complete and clear page as opposed to a very technical and undecipherable one.

If you were already aware of all that, please disregard these comments. Otherwise, there is a source of computational manipulations there that you may be missing and that allows you to work with any tetrad of your preference.

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

Hi

Good catch. The variable keeping track of intermediate simplifications was reinitialized by mistake at line 19, as shown by showstat. This problem is fixed and the fix available to everybody using Maple 2020 within the Maplesoft Physics Updates v.828 or newer.

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

Hi nm,

I realize the natural interest in having that piece of information. In the past, there was a Maple worksheet, updated with each new version, with that contents.

But then the Updates grew in coverage and scope. The task of producing fixes that are compatible with "the previous version" (that for you is the "current version") plus having new developments also consistent, plus testing a large number of things before uploading a new version of the Updates, all together takes an amount of time that leaves 0 for that description you ask now.

It is also the case that I am producing these Physics Updates during spare time, only because I like the idea; the fixes and developments are happening - so why not present them around the clock? Most fixes in new versions, or new features, are directly related to feedback by people, so someone at least knows what is new in the version being uploaded. Unfortunately, this project also makes me divert valuable time on the infrastructure that makes it possible (MapleCloud, etc.).

So yes, it would be nice, but it won't happen.

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

It is probably related to the large number of changes happening at this point. I am adding a test for this example.

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

Nice to see someone playing with the symgen routines; there are, for me, several beautiful ideas behind this code. This is my first work in differential equations, done in 1995 with Luiz Duarte and Luis da Mota in Brasil. On to the problem, the core idea in the symgen routines, that surprised several for their efficacy, is that we map differential equations into an algebraic problem for which, frequently, solutions exist that are computable without solving any differential equation. You realize the gain this represents. On the other hand, it is not the case that such solutions always exist.

Take your ODE

You see, the solution is easy to compute BUT requires solving differential equations and so are out of the scope of symgen. That is why DEtools:-symgen(ode, HINT = [f(x), g(x)*y]) doesn't succeed in this example. An important consideration is that, at the time the symgen routines were written, the capabilities of Maple's dsolve were rudimentary. Mapping an ODE into a problem that required solving differential equations was just of no use. It was actually the presentation of symgen that changed the game, triggering a full rewriting of Maple's dsolve in 1997, now using the symgen routines at its core.

Years later I wrote the symmetry routines of PDEtools, a package also originally written in 1995 in Brasil, with Katherina von Bulow . Basically all the symmetry commands in (O)DEtools have an equivalent in PDEtools. For symgen, the equivalent is Infinitesimals. At the time of writing it, however, Maple's dsolve was already a powerful command, so I coded Infinitesimals such that it tackles the problem for the infinitesimals as a differential equations problem. That is why PDEtools:-Infinitesimals(ode, HINT = [f(x), g(x)*y]) succeeds in this case.

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

In particular, dsolve(equation, useint) will request dsolve to try to compute the integrals, regardless of their complexity or the result potentially be a wallpaper. Analogously, the useInt option requests not to attempt computing any integral. The output is a solution in terms of integrals. The call dsolve(equation) - neither useint nor useInt around - uses an internal algorithm that will attempt computing the integral most of the times but not always, and also discard wallpaper integrations preferring the compact uncomputed integral. That algorithm is pretty efficient, but it may end up avoiding the computation of an integral in a case where the integration can be performed smoothly.

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

@Carl Love,  nm

I wrote this help page and program. I wrote "work as well", that is what I meant. Not at all "works well". The page actually says, as you quoted, that for 1st order ODEs by using this algorithm you map the original problem into an equivalent one, so no gain, i.e. in that case it does not work well. I am a non-native-English speaker but I don't see what would be wrong with this use of "work as well" in this sentence of this help page.

Regarding the hang, as the help page tells, in the first step the algorithm constructed another ODE "equivalently difficult to solve", whose solution involves integrals with wallpaper integrands involving unresolvable RootOfs. The system hangs when trying to perform those integrations in order to move to the second step. You see the flow entering debug(`symgen/formal`). I don't see a problem here. The use of formal for 1st order ODEs is just inappropriate, even if allowed mainly for pedagogical purposes.

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

The problem with mathit is that it removes blank spaces, while in Maple, it is valid to have "symbols" with blank spaces. The textit approach resolved that. There are, of course, other options: add an extra \; after textit, or use the alternatives to work around the removal of blank spaces.

I will think about it.

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

One of the bugs in the old latex command is, precisely, what you are showing: the imaginary unit is I but appears translated as i. The fix consists of: respect your choice.

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

Thanks for the report; this one is fixed in the Maplesoft Physics Updates v.786 or newer. That version also added one more Usexxx variable: Latex:-UseColor := true (default value). With that, for example, Latex(Int(f(x),x)) gets translate using gray, but if you don't want colours set Latex:-UseColor := false.

Regarding one other question you posted: a) now, Typesetting settings do not get changed unless you with(Physics). Besides, you can also assign Latex:-UseTypesettingCurrentSettings := true for that same purpose. And b) assuming does not use Physics:-Assume unless you with(Physics), or you explicitly enter Physics:-Setup(assumingusesAssume = true). For an example of the convenience of using Physics:-Assume instead of assume, see the explanation and examples in ?assuming,details.

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

Update to Maplesoft Physics Updates v.779 or newer, then put Physics:-Latex:-UseTypesettingCurrentSettings := true at the top of your scripts, or in a worksheet before using Physics exports (e.g. Latex), or in your .mapleinit file read when you open Maple. Then your "current" Typesetting settings (whatever they are) will not be "extended", when you load Physics, or when you use Physics exports like Latex and execute a product that requires initializing the package.

Details.

In Maple, for a long time, the default was interface(typesetting = standard). Since some years it is interface(typesetting = extended). I like that. For a release or two interface(typesetting = extended) implied on Typesetting:-EnableTypesetRule(Typesetting:-SpecialFunctionRules); Typesetting:-Settings(':-typesetdot' = true); and Typesetting:-Settings(':-typesetprime' = true). And physicists enjoyed that, up to the feedback I receive,

But then those three Typesetting settings were removed from "interface(typesetting = extended)", which remained the default. After asking for opinions, we end up maintaining these three settings active when you initialize Physics (using 'with' or using its exports) and interface(typesetting = extended) (or the equivalent Physics command Setup(mathematicalnotation = true). Since then It's been several years, and it is the first time someone using Physics complains about using the prime or the dot settings - naturally because they are of no use in your case.

In brief, a possible solution to this is as implemented now: set your Typesetting settings as you prefer then enter Physics:-Latex:-UseTypesettingCurrentSettings := true, before using any Physics command or before loading the package, and your Typesetting settings will not be extended or changed.

As for any other issue you noticed (e.g. itsme seems to have crossed with other things but didn't provide examples), please recall the package is not finished, your feedback is very welcome, either describe them in an email with attached worksheet to physics@maplesoft.com, or post that worksheet here (Green arrow please). Thanks.

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


Both the equation and sol are translated properly ... this is what I see:

 

 

restart; sol := dsolve(t*(t-2)^2*(diff(diff(y(t), t), t))+t*(diff(y(t), t))+y(t) = 0, y(t))

y(t) = _C1*HeunC(1/2, 1, -1, -1/2, 3/2, -2/(t-2))+_C2*HeunC(1/2, 1, -1, -1/2, 3/2, -2/(t-2))*(Int(exp(1/(t-2))/HeunC(1/2, 1, -1, -1/2, 3/2, -2/(t-2))^2, t))

(2)

ode := t*(t-2)^2*(diff(diff(y(t), t), t))+t*(diff(y(t), t))+y(t) = 0

t*(t-2)^2*(diff(diff(y(t), t), t))+t*(diff(y(t), t))+y(t) = 0

(3)

with(Physics)

Latex(ode, output = string)

"t \left(t -2\right)^{2} \Mapleoverset{\ldots}{y}\left(t \right)+t \Mapleoverset{.}{y}\left(t \right)+y \left(t \right) = 0"

(4)

Latex(sol, output = string)

"y \left(t \right) = \textit{\_C1} \textit{HC}\left({\frac{1}{2}, 1, -1, -\frac{1}{2}, \frac{3}{2}, -\frac{2}{t -2}}\right)+\textit{\_C2} \textit{HC}\left({\frac{1}{2}, 1, -1, -\frac{1}{2}, \frac{3}{2}, -\frac{2}{t -2}}\right) \left(\int \frac{{\rm e}^{\frac{1}{t -2}}}{\textit{HC}\left({\frac{1}{2}, 1, -1, -\frac{1}{2}, \frac{3}{2}, -\frac{2}{t -2}}\right)^{2}}d t \right)"

(5)

``


Note anyway that the package is not finished.

About \Mapleoverset and other stuff: among the core elements of this new Physics:-Latex there are: to take full advantage of everything already implemented in the Maple Typesetting package, don't be shy in using CTAN official LaTeX packages, and naturally use the Maple own .sty file, maplestd2e.sty. This is how the equation looks when latex processed:

Because all the translation to latex is done taking advantage of the Maple Typesetting package, you can also control some of this translation using Typesetting settings, for example to not use the dot to represent time derivative (see the Typesetting help pages). The same applies to translating mathematical functions.

Regarding using CTAN official LaTeX packages, it is since 1997 that I am Editor at Computer Physics Communications. I have seen enough papers taking advantage of existing advanced typesetting. Perhaps for that reason, in Physics:-Latex, the idea is to really get the most, not just cover gaps in the old latex.

Download xyz.mw

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

Hi
Indeed, the three numbers refer to chapter, equation and case. I am correcting that in the help page. Regarding your other question, I cannot reproduce your problem. This is what I get:

You see the Coordinates there. To understand what could be the problem at your end, It is better if you upload a worksheet showing the input and output; for that purpose you can use the green arrow you see when you post a question or Edit your previous post.

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

Hi

Input ?Physics,Updates and you will see links to most of the pages and posts produced every year, several are related to general relativity. Under Physics Maple 2017 updates, the first link (itemized as "x.") is General Relativity: classification of solutions to Einstein's equations and the Tetrads package. Click it and you will see a couple of paragrams followed by Examples, the second reads "Equivalence for Schwarzschild metric (spherical and Kruskal coordinates)". The contents of that example answer your question. See also the toolbar: you can click the icon that opents that help page in a worksheet, then reproduce the computation pressing enter or running the whole worksheet (that `!!!` icon), and change the metrics of the example by the metrics you have in mind, the adjust things here and there to make this sort-of-template serve your purpose.

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

A PDEAdvisor command, although part of the design of PDEtools, never took off: besides "linear PDEs" or those that admit travelling wave solutions, the classification is scarce. There is, however, the infolevel: set it as in infolevel[pdsolve] := 3 (or 2) and you will see significant and relevant information regarding the type of problem and the methods used to tackle - eventually solve - it.

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

First 19 20 21 22 23 24 25 Last Page 21 of 59