Forums before death by AOL, social media and spammers... "We can't have nice things"
|    sci.math.symbolic    |    Symbolic algebra discussion    |    10,432 messages    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
|    Message 8,833 of 10,432    |
|    clicliclic@freenet.de to Waldek Hebisch    |
|    Re: fyi, rebuild CAS integration tests,     |
|    25 Jul 15 18:46:14    |
      Waldek Hebisch schrieb:       >       > clicliclic@freenet.de wrote:       > >       > > I don't recall hearing any advocates here on sci.math.symbolic, and       > > am unaware of problems caused by the use of numerical       > > differentiation to check symbolic derivatives - when such a check       > > fails, the symbolic result must be wrong provided the numerics is up       > > to the task. In practice, I am using the fast numerical integration       > > available for plotting on Derive to check dubious portions of       > > antiderivatives.       >       > You wrote "provided the numerics is up to the task". In general       > you need to correctly handle branches is numerics. Simply       > "plugging in" expressions into numeric evaluator will not       > produce correct braches for the problem: brach cuts in numeric       > definitions tend to produce spurious discontinuities.              ... and accuracy control allows one to make any numerical evaluation       as stable as desired.              > >       > > I don't get the point of this. Complex functions - be they atomic or       > > composite expressions - must often be rendered single-valued to be       > > meaningful, and everybody concerned (Ahlfors and a reader; a user       > > and FriCAS) must agree about the branch assigned to them. When       > > selecting a 'principal' branch for a complex function, the       > > corresponding real function should preferably be subsumed and       > > unnecessary discontinuities be avoided - otherwise the choice is       > > mostly by convention, including principles like 'CCC'       > > (Counter-Clockwise Continuity). So there are preferences and there       > > is freedom. Whenever one's mathematics depends on the choice, it has       > > to be made clear what branch one is using.       >       > You ignore the fact that symbolic computations usualy do not       > depend on choice of branches: from symbolic point of view to       > branches are indistingushable. Once correct symbolic result       > is obtained, to solve analytic problem one may need to       > determine branches. Choice of branches is specific to given       > problem and frequently requires information not present in       > symbolic version. Of course, once you have concrete analytic       > problem you need clear choice of branches. But usually such       > choice can be separated from other considerations. In particular       > you can produce useful symbolic results without caring about       > branches.              Would z*dif(asec(z), z) = 1/sqrt(z^2-1)) or log(z) + log(w) = log(z*w)       be such a "useful symbolic result"?              While many symbolic manipulations involving complex functions do not       depend on the choice of branches, they may mix those branches so that       the results can be expected to apply only if the functions are viewed as       multi-valued - as involving all their branches simultaneously, of which       there are often infinitely many. This view is obviously incompatible       with a context-independent unambiguous evaluation of the functions, and       with the possible exception (?) of Axiom/FriCAS, I know of no computer       algebra system that takes it. Rather, systems prefer to make complex       functions single-valued symbolically as well as numerically, restricting       the atomic functions to immutable standard branches and defining the       branches of arbitrary expressions through functional composition.              This restriction entails both a loss and a gain. On the downside, many       functional relations of multi-valued complex analysis cease to hold       without qualification, and composite functions may easily (depending on       user attention on input, system attention on output) be specified such       that they carry convoluted branch structures; on the upside, the       restriction allows symbolic results to be checked numerically by the       user - in particular visually after plotting - in the manner familiar       from real analysis, and systems may also use it internally to check and       supplant their symbolic computations - for instance, they may accept       numerical in place of undecidable symbolic zero-equivalence (cf.       Mathematica).              You seem to believe the gain not to be worth the loss, but I don't see       how the present FriCAS can be taken to support this belief, nor do I see       interesting unique capabilities that could result from pursuing it.              > >       > > Each of the complex functions of a CAS, on the other hand, comes       > > with a specific branch built into the numerics, and its symbolics is       > > hopefully designed to comply - no provision is made for the user to       > > specify a new branch choice for any atomic or composite function, or       > > for the numerics and symbolics to accept and act on such a       > > specification. This leaves the preinstalled branches and their       > > functional composition as the only unambiguous branch assignment       > > practical for arbitrary expressions.       >       > Once you are dealing with arbitrary expressions numeric branches       > frequently are bogus. You simply can not use "arbitrary expression"       > as useful definition of branch behaviour. Simply, branch information       > must be carried "out of band". If you want symbolic representation       > one possibility is to specify domain and choice at a single point       > (frequently just specifying value is enough to identify branch).       > If expression is analytic and domain is simply connected, then       > this uniquely determines the function. Clearly, significant       > effert will be needed to make approach like this practical.              Pie in the sky.              Branches "can" be assigned through functional composition because one is       free to do so. Though 'nicer' choices will often exist, this is "useful"       because an unambiguous evaluation is imparted to any expression without       reference to separate information. Your "must" does therefore only       indicate a personal viewpoint, not a mathematical necessity. And your       separate step of analytic continuation from a specified germ could be       performed equally well when symbolic results systematically reflect       input branch cuts; the ambiguity characteristic of alternative branch       assignments would then enter through this optional step alone, while the       gain outlined would not be sacrificed without need.              All this will be above the head of the great majority of users with       nothing but real functions of a real variable on their mind. Yet such a       user may well happen to cobble together a differential equation which       leads to a simple integral whose solution by FriCAS 1.2.6 doesn't hold       everywhere (because the standard branches cause the complex plane to       break into incompatible pieces of which FriCAS adopts one for the entire       plane), and whether his solution is valid (cf. Timofeev Chapter 1,       Example 64 for a>0, x>0) or invalid (ditto for a<0, x<0) is a matter of       luck. The nuisance caused by having to find out that the solution       returned doesn't satisfy the differential equation symbolically, and       then perhaps trying to determine numerically if and where in a              [continued in next message]              --- SoupGate-Win32 v1.05        * Origin: you cannot sedate... all the things you hate (1:229/2)    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
(c) 1994, bbs@darkrealms.ca