@Michael_Watson 

In the reviewed worksheet I posted with the previous reply, I do not see difficulties in eqn3 - what is the problem you are experiencing? Perhaps if you share a worksheet showing the issue?

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

@acer 

The edited version you posted looks excellent Acer!

Best

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

@VicenteAtal 

I know some bits on the internals of the Maple worksheet representation but this is beyond my knowledge - I forwarded your worksheet to the Maple GUI gurus.

Best

Edgardo

Hi

Perhaps you could post a worksheet showing, say, where you are with your problem and where is that you don't know how to proceed further?

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

@J F Ogilvie 

Yes, any of these two you mentioned would serve the purpose you have, but you need to reverse the order in the equations; i.e use 

BTW, the "powerhouses" for developments like this one are the packages, that require more and more functionality.

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

@J F Ogilvie 

Replacing assume by this new Physics:-Library:-Assume would not be appropriate at this point, after so many years of existence of assume. If we were to do that, worksheets or programs already developed by people and that rely on, or make use of the redefinition of variables performed by the old assume, might then require adjustements, because removing this inconvenient redefinition is actually the main feature of this new Physics:-Library:-Assume.

On the other hand, after a week using this new Physics:-Library:-Assume, the advantages are so obvious that Physics:-Library:-Assume could better be a user level command, instead of one somehow hidden two package-levels within Physics.

A similar situation happens with this new functionality, automaticsimplification, as well as the redefinesum and combinepowersofsamebase, all of them current computing-mode options of the Physics package, that introduce functionality very useful also outside the framework of Physics. 

Likewise, PDEtools:-declare, PDEtools:-Solve and PDEtools:-casesplit have highly regarded functionality useful beyond differential equations, as well as the commands for performing Union and Intersection of regions of the complex plane expressed algebraically, both currently within the MathematicalFunctions package.

But things also change with time. The originally internal solver of ODEtools, called `ODEtools/Solve` ended up as a useful addition to solve; the code under `ODEtools/first_integral` got renamed as int/homotopic and is nowadays an important part of int, the same way as internal DE code formerly developed to compute certain kind of elliptic integrals, and the current `int/parallel`, also former DE code for integration. In the same line, IntegrationTools:-Change is actually a wrapper around PDEtools:-dchange. Or more important, the internal DE code to manipulate assumptions related to DE parameters ended up as the skeleton of the current assuming command, the internal DE conversion routines for mathematical functions related to DE special function solutions are today the skeleton of the FunctionAdvisor, the routines for compacting DE solutions are now known as simplify(expresion, size), also now at the core of this new automaticsimplification of Physics, and etc. the list is large.

All this is to say that the Physics, PDEtools and MathematicalFunction packages are also, somehow, small powerhouses for these developments. That is good, as opposed to "bad because of happening within packages". Let's see the filled part of the cup. And when the development is very good it often ends up available as a user-level command in one way or another.

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

@Markiyan Hirnyk 

{1+I*y : y>=0} union {-1+I*y : y<=0} is not a branch cut of sqrt(1-z^2). Just clarifying. Having said that, yes, a command of the form is(..., branch_cut of ...) would be useful. All the tools for developing it are already in place; these are the hidden exports of MathematicalFunctions, `&Intersect`, `&Minus`, `&Union`, that I still need to document.

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

The reference you are mentioning actually corresponds to the work of the U.Bath group previous to Maple 17. The actual reference describing the Maple 17 implementation is:

"England, M., Cheb-Terrab, E., Bradford, R., Davenport, J. and Wilson, D. (2013) Branch Cuts in Maple 17. ACM Communications in Computer Algebra, 48 (1). pp. 24-27. ISSN 1932-2232"

The main reference for computer algebra implementations of the computation of branch cuts is anyway 

"A. Dingle and R.J. Fateman. Branch cuts in computer algebra. In Proceedings of ISSAC ’94, pages 250–257, 1994."

The method used by the U.Bath group, that is the departure point for the Maple 17 implementation, is in fact mentioned in this paper from 1994.

In both papers, the one from 1994 and the one from 2014, it is pointed out that the difficult part is the detection of spurious cuts (in order to not return them as if they were actual cuts). Despite the terminology, spurious cuts are easy to understand: they basically happen when different sub-expressions of one expression have cuts that intersect over some segments but with opposite signs, cancelling each other over those intersecting segments. Ie over those segments there are no branch cut. All the computer algebra implementations of computation of branch cuts previous to Maple 18 have trouble in detecting these intersecting segments.

New in Maple 18, however, spurious cuts that are over the real axis are actually detected (ie are not returned). The algorithm implemented in Maple 18 is new, so not described in any of the two papers above. This algorithm was developed entirely at Maplesoft. Through rotations, this algorithm can be extended to cover straight segments anywhere over the complex plane, and there is a way (more complicated, but doable) to also cover non-straight segments.

I still need to find the time to do the writeup of the algorithm implemented in Maple 18, then present the resulting description for publication elsewhere, as well as to finish extending the implementation as said in the previous paragraph. Not so soon ... For now I'm rather busy with the Physics package, mainly new developments for general relativity and covering some gaps in functionality for particle physics.

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

@Andriy 

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

@nm 
Other people already requested the automation of the installation of Physics updates by means of a mini-installer. Tha would save having to unzip, read installation instructions, move files around etc. These requests were forwarded to the people who take care of that at Maplesoft.

Note anyway that there are two different things: having a mini-installer for the update and having a button to check for updates. The latter is more complicated in that it requires a change in the GUI ... perhaps in a future release. Also, currently, updates hare happening mostly everyday, so a button to 'check for updates' would be redundant. The relevant step ahead today, I think, would be to have a mini-installer for these updates.

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

Updates for Physics are available on the Maplesoft R&D Physics webpage around the clock, posted every week, and include bug fixes (e.g. those reported here in Mapleprimes) and also new developments, functionality.

Maple 18.02 includes all the changes that happened in the package till July/30, that are listed below in a linked PDF, as well as in a linked worksheet including examples for all the changes. 

Likewise all the changes that happened in the differential equation programs and the MathematicalFunctions package till July/30, made available for download along the year on the Maplesoft R&D DEs and Mathematical Functions webpage, are also included in Maple 18.02.

The update being distributed for Physics in its Maplesoft R&D webpage, however, includes a substantial number of further changes and new functionality not present in Maple 18.02. To mention but one, there is a new Tetrads subpackage for General Relativity. These further changes include all the fixes and new developments from August/1 till today, and run natively in the new Maple 18.02 just released.

NOTE: If you have been updating Physics till recently, installing Maple 18.02 will automatically reset the changes in Physics to July/30, which is not a problem: to activate all the changes from August/1 till today please update Physics again, as said downloading it from the Maplesoft R&D Physics webpage.

The links mentioned: a PDF describing the changes in Physics included in Maple 18.02 and Worksheet to reproduce the changes in Physics 18.02 described in the PDF.

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

@NOh 

I cannot explain why the website suddenly loses contents. It happened other times with other people too. The contents is there for a while, after uploading, and then, sometimes, suddenly, it is not there anymore. It was reported. I thought it was fixed. Unfortunately it is still happening.

I uploaded the worksheet again.

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

@guyp 

Yes, it would be interesting to have these issues fixed directly in sum instead of having to redefine sum. Note anyway that the redefinition mentioned, is already a core Maple procedure, and the whole Physics package is too. This issue is pretty relevant in Physics, where we keep entering these sums and not understanding why these unexpected results.

But there are more things to take into account: some people would prefer to see multi-index summation directly in the add command, while others would prefer to not change the premature evaluation because they take advantage of that (sum is a command with more than 30 years around), etc.

The issue is then debatable, and I personally would also prefer to just fix the standard and old Maple sum command for everybody, and just be sorry if someone was taking advantage of this premature evaluation. But for now what we have is this redefinition provided, it concretely fixes the core sum command (it does exactly what we would do if we were to fix this problem in sum), is provided as an official approach to the problem, and tries to address all these different things mentioned at once.

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

Hi NOh

I'm pretty busy in this moment - will give a look at this and return some feedback hopefully by the end of tomorrow. Meantime, it would help concretely when revising/optimizing your code if you could explain, before that review, in few words (say a paragraph per routine), what is what you are trying to achieve with each of these routines you have programmed.

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

Hi Age

Sounds like a good idea to me. Regarding the problem you mentioned in this first post, it is more a GUI issue, where % and (2) are not exactly the same thing (they should). Independent of that I tweaked the Physics code to avoid this GUI subtlety. So the two results are the same, equal to 1, and in that sense the problem is fixed (not the underlying GUI issue though). The adjustment can be downloaded as an update to the Physics package from the usual place, the Maplesoft R&D Physics webpage.

Independent of that: note there is a distinction between `.` and `*`, and when you compute algebraically the operation is `*`. The code takes care of this duality in some places but using `.` to perform algebraic (non-matricial and non Bra, Ket stuff) you will not always get the operation performed.

There were other posts in Mapleprimes regarding how to implement the algebra of Pauli matrices using commutators and anticommutators. The plan is to have these algebras known by the system, so that nobody needs to enter them, but that is still not implemented. This week and perhaps the next one I will still be working on tetrads.

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