by **B Lautrup** from CERN Courier Volume 12, Number 3, March 1972 (pdf)

Theoretical physicists mostly use the computer facilities for numerical evaluation of theories in order to compare theory with experiment (or theory with theory). Usually this only puts moderate demands on the computer and its software. Most programs require only a library of special functions and virtually never use tapes or other ‘ complicated ‘ peripheral equipment. They are reasonably satisfied with a card reader, for input and a line printer for output. About 80% of the computing jobs’ from the CERN Theory Division are of this scale. They would in fact be even better off in an interactive environment with an ‘ advanced desk calculator’. The GAMMA system described above is, among other things, a sophisticated attempt to cater for these needs. Some theorists use computers on a grander scale. Numerical integration in many dimensions, phase shift analysis and multiparameter fits to phenomenological models lead to large central processor time consumption and require fast computers. This kind of computation is analogous to the last stage in the analysis of experimental results. A more exotic use of the computer by theorists is in the field of algebraic manipulations. A theoretical expression may not immediately submit to numerical evaluation but require a certain amount of algebraic simplification before numbers can be extracted. This is for instance the case in the evaluation of Feynman graphs in quantum electrodynamics and related fields. A Feynman graph leads to multi-dimensional integrals, over infinite regions, of complicated products of matrices, vectors and spi-nors. Not only is the integration region infinite but the integrand has singularities within this region. This prevents a direct numerical attack on the integral and other methods are needed. First of all the problem is to trans-B. Lautrup form the integration region into a finite domain. This is done in a series of well-defined steps reducing the integral to a standard form. In this process matrices and tensors are multiplied out, new integration variables introduced, vectors multiplied with each other to yield physical scalar quantities, the old integration variables integrated out, etc. Although this is quite simple to do by hand for uncomplicated graphs, it becomes extremely cumbersome for higher order graphs. The number of terms is so large that it is beyond the ability of the human mind to handle them. As the operations performed on each term are almost identical, the problem is mostly of a bookkeeping nature. This is an ideal task for the computer since it is able to do the same operation over and over again without tiring, or making errors. At CERN a large-scale algebraic manipulation program (SCHOON-SCHIP) capable of carrying out the reduction described above is available. It was written by a single person, M. Veltman from Utrecht, in the middle of the 1960s and has been continually expanded and modified since then. The program accepts input in the form of an expression containing gamma-matrices, vectors, scalars and even spinors. It allows for expansion, trace calculation, collection of identical terms, a variety of substitutions, limited factorisation, asymptotic expansion, etc. It has even been possible to perform several-dimensional analytic integration using the substitution facilities. Another program capable of carrying out the same kind of calculations will eventually become available at CERN. This is the language oriented program REDUCE written by A. Hearn from Salt Lake City. It is more flexible than SCHOONSCHIP but requires much more software. It has been implemented interactively on a PDP-10.

Image source: twitter