The AbelMap command computes the Abel map between two points P and P_0 on a Riemann surface R of genus g, that is a g-tuple of complex numbers. The jth element of the Abel map is the integral of the jth normalized holomorphic differential integrated along a path from P to P_0.

The Riemann surface is entered as F; an irreducible, square-free polynomial in x and y. Floating point numbers are not allowed as coefficients of F. Algebraic numbers are allowed. Curves of arbitrary finite genus with arbitrary singularities are allowed.

The points P and P_0 are entered as , where a and b are constants, and r is an integer. If r < 0, that is, if entering one of the points for , then a = 0.

The differentials are normalized such that the jth differential integrated around the kth cycle, as given by algcurves[homology], is Kronecker delta (j, k).

Note: The Abel map will almost always be computed along with other objects associated with some polynomial F, such as the Riemann matrix. It is imperative that the order of the differential be the same for each of the objects, and at each stage of the calculation. As no order is imposed by algcurves[differentials], make sure to compute AbelMap and, for instance algcurves[periodmatrix], without a restart (or quit) in between.

This command is based on code written by Bernard Deconinck, Michael A. Nivala, and Matthew S. Patterson.