Finite gap Jacobi matrices and the Schottky-Klein prime function

A covering map formalism for studying the spectral curves assocaited with finite gap Jacobi matrices is presented. We advocate a constructive function theoretic framework based on use of the Schottky-Klein prime function. The single gap, or genus-one, case is studied in explicit detail.


Introduction
This paper is motivated by recent work of Christiansen, Simon and Zinchenko [9] who describe a covering map formalism to define and study the isospectral torus associated with finite gap Jacobi matrices. In turn, that paper drew on earlier work of Sodin & Yuditskii [37] and Peherstorfer & Yuditskii [34] which concentrated on infinite gap sets. The focus in the present paper is the finite-gap situation, and our goals are limited to offering a complementary perspective to the covering map formalism of prior work. It is offered both for mathematical interest, and because it is believed that the approach here provides a framework for constructive computational methods.
In other work, the author has been interested in physical science problems involving arbitrary multiply connected domains, that is, finitely connected planar domains generally unpossessing of any geometrical symmetries. The mathematical approach has been to develop a constructive framework based on the Schottky uniformization of the so-called Schottky double -a compact Riemann surface naturally associated to any given multiply connected planar domain -and to make use of the Schottky-Klein prime function that lives on that surface as the basic functional building block for more complicated functions (we adopt the designation "Schottky-Klein prime function" following Baker [3]; see also [7] [8]). This idea has proven to be surprisingly fruitful. A survey of many other results, including applications to potential theory, the construction of multiply connected quadrature domains, the solution to a variety of problems in fluid dynamics, and the derivation of a multiply connected extension of the classical Schwarz-Christoffel conformal mapping formula is given in the review articles [10] [11]. Most recently, a connection between a problem in the field of vortex dynamics and the planar limit of random normal matrices has been pointed out [22]; both those application areas benefit from the function theoretic approach just outlined, either for the description of vortex patches ("V-states") or, in the random matrix case, the limiting shapes of the multi-support eigenvalue sets.
The central point of contact with finite-gap Jacobi matrices is the observation that the hyperelliptic curve associated with a given multi-support spectrum can be viewed as the Schottky double of the multiply connected domain exterior to a finite set of bands on the real axis. For general multiply connected domains, not necessarily enjoying any reflectional symmetry about the real axis, the Schottky double is made up of two "sides", often referred to as a front-side and a back-side, with an antiholomorphic involution relating them [28]. With an additional reflectional symmetry about the real axis the natural holomorphic coordinate on the back-side of the double becomes the same as that on the front-side and the Schottky double is a two-sheeted hyperelliptic Riemann surface which we call S following [9]. The two sheets of the surface will be denoted by S + (the "physical" sheet) and S − . Accordingly, the points at infinity on the two sheets will be denoted by ∞ + and ∞ − .
To fix ideas, a (one-sided) Jacobi matrix J as considered in [9] is one of the tridiagonal form The key mathematical results of this paper can be summarized concisely as follows. Let the essential spectrum of such a Jacobi matrix take the form e = e 0 ∪ e 2 ∪ · · · ∪ e l , which is a union of l + 1 disjoint closed intervals, or bands, with l gaps. The case l = 2 is shown schematically in Figure 1. It is known that the isospectral torus of a Jacobi matrix having such a spectrum can be parametrized in terms Fig. 1 A typical finite-gap spectrum on the real axis with l = 2.
of so-called minimal Herglotz functions [9]; the isospectral torus essentially encodes the set of Jacobi matrices having the same essential spectrum. We show here that a parametrization of a typical minimal Herglotz function m(z) associated to such a finite gap set is furnished, in terms of a parameter ζ, by the explicit formulas where ω(., .) is the Schottky-Klein prime function associated to a multiply connected circular domain D ζ consisting of the unit ζ-disc with l smaller circular discs excised, all with centres located on the real diameter (−1, 1). A and D are constants, the set {r j |j = 1, ..., l} give the l poles of the minimal Herglotz function, while {s j |j = 1, ..., l} is a set of its zeros. Another pole at the real point α corresponds to ∞ + . The exponential terms appearing in M (ζ) will be explained later in the paper. Based on these results we find the following explicit expressions that assist in the inverse problem of finding the Jacobi parameters associated to e, as explained later. Namely, where is meromorphic on S (it is automorphic under a Schottky group action associated to D ζ ). The set of zeros {s j |j = 1, ..., l} is given by the solution of a system of equations written down in §6. The parameter X ∞ is related to the logarithmic capacity of e while Y ∞ is geometrical parameter (explained later) associated with the spectral bands. Formulas for both are also available in terms of the S-K prime function and a related functionω(ζ, α) defined later in (28): The Schottky double setting has been employed previously by Aptekarev [2] in a study of orthogonal polynomials on a general set of contours, but he builds the associated function theory based on a Riemann theta function formalism. Here we eschew the latter and use instead the Schottky-Klein prime function as the basic building block in the belief that it is a natural way to build the requisite functions. Lukashov and Peherstorfer [25] (see also [7] [8]) have used a function theory similar to ours in their studies of orthogonal polynomials (although they use the designation Schottky-Burnside function to refer to what we call the Schottky-Klein prime function). Christiansen, Simon and Zinchenko [9] also avoid use of the Riemann theta function and develop their own construction based on representations of the requisite functions in terms of infinite Blaschke products. Additional discussion of the Schottky model and its applications in approximation theory is given in the monograph by Bogatyrëv [6].
It is worth clarifying how our approach differs from previous work since it allows us to underline the reason for writing this article. Like [9] we adopt a covering map formalism but our viewpoint differs in that it is the bands, not the gaps, that correspond to Möbius-identified circles in our uniformization domain. But the most important distinction is that, within this alternative model of the underlying curves, the function theory we employ to solve the inverse problem is different to that used in previous work: here we use the Schottky-Klein prime function which affords many constructive advantages. This is because the efficient computation of S-K prime functions has been a focus of much recent work by the author and co-workers who have prepared freely downloadable software, on the MATLAB platform, for its fast and effective computation [17] [23]. Further discussion of this is given in §11.

The Schottky double of a domain
It is known, by a multiply connected extension of the Riemann mapping theorem [27], that any (M + 1)-connected domain (for M ≥ 0) is conformally equivalent to the unit ζ-disc with M smaller circular discs excised. Let such a domain be denoted D ζ . To model the so-called Schottky double of D ζ one can reflect, in the unit circle, the circular boundaries of the interior discs to produce M additional circles in |ζ| > 1. By this same antiholomorphic reflection of points in D ζ , one produces a precise copy, or "back-side", of the original multiply connected domain outside the unit ζ-disc. We then identify each circle inside |ζ| < 1 with its reflection in |ζ| > 1 by means of a holomorphic Möbius transformation. For an (M + 1)-connected domain there will be M such Möbius transformations. The basic idea of this construction lies at the heart of the Schottky model of algebraic curves and, by forming the union of D ζ with its back-side (its reflection in the unit circle) and identifying the boundary circles as just described, we produce a model of the Schottky double. Gustafsson [28] and Varchenko & Etingof [39] have discussed the socalled Hele-Shaw free boundary problem in a multiply connected domain by considering the Schottky double of the domain. Aptekarev [2] used it in his study of orthogonal polynomials on a general set of contours.
In more detail, suppose the M smaller circular discs inside |ζ| < 1 have centres {δ j |j = 1, ..., M } and radii {q j |j = 1, ..., M }. The data {δ j , q j |j = 1, ..., M } will be called the conformal moduli of D ζ and they can be uniquely associated to a given set of spectral bands up to specification of an automorphism of the unit disc. Let the unit circle be denoted C 0 and let the M interior circular boundaries be denoted by Figure 2 shows a schematic in the quadruply con- It is straightforward to check that for points on the circle C k We define the reflection of a point ζ in the circle C k by φ k (ζ). Then for k = 1, ..., M introduce the Möbius transformation θ k (ζ) defined by It follows from (9) and (7) that It is straightforward to show that θ k (ζ) maps C k holomorphically onto C k as illustrated schematically in Figure 2.
The set Θ consisting of all functional compositions of the maps {θ k (ζ)|k = 1, ..., M } and their inverses, including the identity map, is an example of a classical Schottky group [3] [4]. We refer to the maps {θ k (ζ)|k = 1, ..., M } and their inverses as the generators of the group Θ. A fundamental region of Θ is a connected region whose images under all maps in Θ tessellate the complement of the limit set of the group. Let us define F as the region consisting of D ζ and its reflection in C 0 , i.e., the 2M -connected region bounded by {C k , C k |k = 1, ..., M }. Then F is a fundamental region for the group Θ. [ where [v k (ζ)] Cj and [v k (ζ)] C j denote respectively the changes in v k (ζ) on traversing C j and C j with the interior of F on the right. δ jk denotes the Kronecker delta function. Then, for j, k = 1, ..., M , and modulo the addition of integers dependent on the choice of integration paths for for some {τ jk |j, k = 1, ..., M } which are constants (i.e., they are independent of ζ, but they depend on the conformal moduli). It is established in [29] that there exists a unique function X(ζ, α) defined by the properties: (i) X(ζ, α) is single-valued and analytic in F (except possibly at infinity if this point is contained in F ). (ii) X(ζ, α) has a second-order zero at each of the points θ(α), θ ∈ Θ.
(iii) It is normalized such that The Schottky-Klein prime function (henceforth, S-K prime function) is then defined as the square root of X(ζ, α). There is a similar characterization for the prime function itself subject only to an ambiguity in sign, a matter which has recently been discussed by Bogatyrëv [5].
For the special class of Schottky groups considered here it is also possible to establish [16] [40] that the corresponding S-K prime function and its Schwarz conjugate function satisfy the additional functional relation where the Schwarz conjugate h(ζ) of an analytic function h(ζ) is defined by where θ i (ζ) is taken to be all elements in the set, denoted by Θ , comprising all elements of the group Θ excluding the identity and all inverses. While this formula has been used in applications [13] [16] to compute the prime function it can be inefficient to implement numerically especially as the connectivity of the domains grows. It is for this reason that the author, with collaborators, has developed novel numerical algorithms for the evaluation of the prime function that do not rely on this classical infinite product formula [17] [23]. Indeed an important fact, with both theoretical and numerical implications (cf. Crowdy & Marshall [17]) is that all the symmetries of Schottky doubles of planar domains imply that the relevant integrals of the first-kind {v k (ζ)|k = 1, ..., M } are linear combinations of the analytic extensions of the harmonic measures of the domain (defined below). Crowdy & Marshall [17] used this fact to devise a novel and effective numerical algorithm to compute both the {v k (ζ)|k = 1, ..., M } -and hence, in effect, the Abel map -as well as the S-K prime function itself, and without recourse to the (often slowly convergent) infinite product formula (16). The algorithms of [17] have now been improved even further [23]. In more detail, it is known [17] that Im[v j (ζ)] = 0 on |ζ| = 1 so that by analytic continuation off this circle Now introduce the modified functions where {Ω j |j = 1, ..., l} are what we refer to as the harmonic measures [17] of the domain: for any j = 1, ..., M , the function Ω j is harmonic everywhere in the fundamental region and satisfies the boundary conditions for some symmetric invertible matrix Q jk where, in fact [17], where τ jk appeared earlier in (12).

Potential theory
In their development Christiansen, Simon and Zinchenko [9] make connections between the inverse Jacobi problem and the potential theory associated with the spectral bands. To make contact with those observations we introduce a few additional ideas.
For an M -connected domain D ζ there are M + 1 so-called modified Green's functions {G j (ζ, α)|j = 0, 1, . . . M }. The function G j (ζ, α) has a logarithmic singularity in D ζ at α, vanishes on boundary C j , and has a constant value on all other boundaries of D ζ [35]; those constants are determined by the conditions that where s denotes boundary arclength and ∂/∂n is the normal derivative. These modified Green's functions are distinct from the usual potential theoretic firsttype Green's function which has zero imaginary part on all boundaries (and does not, in general, satisfy (22)). Within our Schottky model of the hyperelliptic curve the analytic extensioñ G 0 (ζ, α) of a modified Green's function G 0 (ζ, α) for our domain has relevance to the inverse Jacobi problem. It was demonstrated in [16] [21] that the analytic extensionG 0 (ζ, α) of a multiply connected domain D ζ has, in terms of the associated S-K prime function, the concise representatioñ Preimage circular domain D ζ in the ζ plane: the unit circle with two excised circles with centres δ 1 and δ 2 on the real axis. Their radii are q 1 and q 2 . The real point α maps to ∞ + ; its reflection 1/α, which sits outside the unit disc, corresponds to ∞ − .
It follows from (23) that and it is this fact that will be useful later.
5 Conformal mapping to l + 1 real slits All of the mathematical technology we have introduced so far pertains to any multiply connected domain of connectivity M + 1. Henceforth we set M = l, to accord with the notation of [9] where the spectrum of a Jacobi matrix is taken to comprise l + 1 bands on the real axis with l finite gaps. Only now do we restrict to the special class of domains D ζ for which all the circles {C j |j = 1, . . . , l} are centred on the real interval (−1, 1). We wish to construct a conformal map from D ζ to the unbounded region exterior to the l+1 bands e and where the circular boundary C j is transplanted to the band [α j , β j ] for j = 0, 1, . . . , l. Portions of the real ζ-axis inside the unit disc will be transplanted to the gaps with some point α chosen to map to z = ∞. By the freedoms of the Riemann mapping theorem the point α can be chosen as we like provided it lies on the real ζ-axis outside the circles {C k |k = 1, ..., l} and inside C 0 .
The required conformal mapping is provided by an instance of the so-called "radial slit mapping" formula written down in terms of the S-K prime function by Crowdy and Marshall [12]; see also [25] where a similar results is recorded in a different context. For the present application we require α, β ∈ R. It is easy to show from the symmetry property (15) of the S-K prime function that (25) in fact maps the other half of F (i.e. the reflection of D ζ in the unit circle) onto another copy of the plane exterior to the slits. The explicit form (25) means that once a domain D ζ is fixed, so that the conformal modulus set {q j , δ j |j = 1, ..., l} and the parameters A, α and β are known, formulas for the finite gap set e follows directly from the formulas Conversely, if the finite gap set e is given then it is a simple numerical matter to inverse the explicit nonlinear relations (26) to find the relevant conformal moduli. Furthermore, from a Laurent expansion of (25) near ζ = α we have that, locally, for some constants X ∞ and Y ∞ . Indeed, if we define the functionω(ζ, α) via then we find the formulas (6). The parameter X ∞ is related to the conformal radius or logarithmic capacity of the domain. Immediately we observe that, near z = ∞ or ζ = α, The function (25) is invariant under the action of the Schottky group Θ, i.e., it is automorphic with respect to this group action. This can be argued directly from its required properties as the relevant conformal mapping function where the image of each circle C k (for k = 0, 1, ..., l) is required to lie on the real zaxis. Alternatively, given the explicit form (25), it can also be verified directly on use of the properties (14) of the S-K prime function.
Although we do not make use of them here we record that there are at least two alternative ways to write the conformal map taking D ζ to S + in terms of the associated S-K prime function. We could also write it as the "parallel slit map" [12] whereG 0 (ζ; α) is the analytic extension of the modified Green's function given in (23) and B and C are real constants. (This provides an interesting link between the covering map and the potential theory associated with the spectrum). This alternative form is also discussed in [12]. The representation where E and F are real constants, is also viable. Formula (31) has been used in an application studying fluid motion through gaps in walls [19].

Minimal Herglotz functions
From its characterization given in Theorem 6.1 of [9] we know that a minimal Herglotz function m(z) (a) is a rational function on the hyperelliptic curve S; (b) it has a simple pole at ∞ − and one more in each of the l gaps; (c) it also has a simple zero at ∞ + with asymptotic behaviour where a 1 and b 1 are Jacobi parameters. It also has one more zero in each of the l gaps.
The procedure of "coefficient stripping" [9] corresponds to the following: with the function m 1 (z), say, yielding parameters a 1 and b 1 successfully found, the parameters a 2 and b 2 can be determined by constructing another minimal Herglotz function m 2 (z), say, with all the same properties of m 1 (z) but where the simple zeros of m 1 (z) in the gaps become the simple poles of next generation Herglotz function m 2 (z). By iterating this construction the full set of Jacobi parameters can be generated.
To effect this construction within our model we introduce If m(z) is rational on S then the function M (ζ) must be automorphic with respect to the Schottky group Θ. It follows from the general theory for automorphic functions of this kind [3] [13] that we can write for some set of integers {n k |k = 1, ..., l} and some constant D where ω(., .) is the relevant Schottky-Klein prime function with the poles {1/α, r j |j = 1, ..., l} and zeros {α, s j |j = 1, ..., l} chosen to satisfy a set of automorphicity conditions [13] -these are just the conditions of Abel's theorem [3]. It is clear from (34) that M (ζ) has the required simple pole at ζ = 1/α, corresponding to z = ∞ − , and that it vanishes at ζ = α, i.e., at z = ∞ + . Each point ζ = r j is on the real axis in F and projects down to the gap [β j−1 , α j ] for j = 1, . . . , l.
On use of (14) these l conditions can be expressed as n k τ kj , j = 1, ..., l. (35) Observe that the first factor in the representation (34) for the minimal Herglotz function can be written, on use of (24), in the form whereD is some constant andG 0 (ζ, α) is the modified Green's function of (23). A similar expression involving the potential theoretic Green's function is found using the alternative approach of [9]. Note that finding the zero set {s j } for a given set of poles from the expressions in (35) is an example of a Jacobi inversion problem. We proceed to present our construction under the assumption that its solution exists and is unique. Finally, we also make a comment here on the integrals of the first kind for the classes of domains to which we have now restricted. For such domains we can always pick the real part of v j (ζ), for each j = 1, ..., l, for ζ on the real axis to be zero, so that Together, (17) and (37) imply that v j (1/ζ) = −v j (ζ), j = 1, ..., l.
7 Inverse problem for the Jacobi parameters We can write where Since both M (ζ) and f (ζ) are automorphic with respect to the Schottky group action then so is N (ζ). A Taylor expansion of N (ζ) near ζ = α gives or, on use of (29), On substitution of this into (39) we find Enforcing the condition that On substitution of this result into (43), we arrive at The logarithmic derivatives arising naturally in these expressions imply that they can be rewritten more explicitly as equations (4) which are independent of the prefactor AD in the formula (40) for N (ζ) which is why we have set AD = 1 in (5).

Periodic Jacobi matrices
To examine the nature of any conditions necessary for the determination of a periodic Jacobi matrix we now study a specific example. Modifications for the general case will be obvious from what follows. Suppose that we seek a p = 4 periodic Jacobi matrix in the l = 2 gap case. Suppose that the initial poles are at {r 1 , x 1 } with zeros, following from Abel's theorem, at {s 1 , y 1 } (in addition to the pole at 1/α and zero at α that are present at every stage). Then, from (35), we have for some integers {n where, since the zeros of the last iteration (aside from the one at α) will now become the poles of the next minimal Herglotz function (aside from the one at 1/α), and since the matrix is periodic by assumption, the sequence must repeat itself. The poles and zeros (48) satisfy, respectively, the Abel conditions and and for some {n (m) jk ∈ Z|j, k = 1, 2; m = 1, 2, 3, 4}. Addition of (47) and (49)-(51) produces a cancellation of most of the terms leaving, for j = 1, 2, where is some integer. But (38) allows us to deduce that (52) is k τ kj , for j = 1, ..., l.
But the matrix Q jk is symmetric and invertible [17] [35] so a comparison of (55) and (56) leads tõ which, on equating imaginary parts, implies that we must have for some integer m k . Thus, the condition for a p-periodic Jacobi matrix is that the l harmonic measures, when evaluated at infinity (i.e., at ζ = α), must be rational numbers of the form (58). Such matters are explored in more detail in [9] within the covering map formalism considered there; other related discussion is to be found in Corollary 5 of [25]. The above treatment aims to give a flavour of how these considerations translate to the Schottky uniformization used here.

Two spectral bands
An instructive exercise is to study the single gap case (genus one) within this Schottky model because the construction can be performed analytically, and quite explicitly. It is then possible to see both the construction based on our formulation in action, as well as many features of the underlying general theory [9]. All the required function theory can be derived from elementary considerations without invoking any special results from classical elliptic or theta function theory. The genus-one case has also been considered, from different perspectives, by other authors [1] [33].
For the concentric annulus ρ < |ζ| < 1 the Schottky-Klein prime function is where P (ζ) can be defined by the infinite product This is precisely the form given by (16) in this doubly connected case where the Schottky group has just a single generator (and its inverse). As for the prime function itself, we suppress the dependence of P (ζ) on the conformal modulus ρ. The product (60) can be shown to converge by standard methods and it is elementary to verify, directly from the infinite product definition, that P (ζ) satisfies Indeed (60) is proportional to the infinite product given in (16) when the Schottky group Θ is generated by θ 1 (ζ) = ρ 2 ζ and its inverse. The second of the relations (61) is essentially just (14) with the function v 1 (ζ) defined by v 1 (ζ) = 1 2πi log ζ, It is clear from (62), and its definition given earlier, that the single harmonic measure in this case is In passing, we note that the Schottky-Klein prime function (59) can be related to the first Jacobi theta function (see [10] [11]). The conformal mapping from the concentric annulus ρ < |ζ| < 1 to two bands on the real axis is for some real α with ρ < α < 1 and real parameters R, ρ and β. The choice of these 4 real parameters corresponds to the freedom to pick the four real end-points of the two spectral bands.
Period-2 solution: Let us take the first minimal Herglotz function be where the simple pole r 1 satisfies −1/ρ < r 1 < −ρ so that it sits in the gap. The condition for automorphicity, i.e., the requirement that M 1 (ζ) satisfies is easily seen to be α 2 s 1 ρ 2 r 1 = 1, or For a period-2 solution the second minimal Herglotz function is then where, in order that this is automorphic, we require Together (67) and (69) imply the condition which is consistent with the requirement that ρ < α < 1.
On taking a logarithm condition (70) is which, by (63), is equivalent to Period-3 solution: Let the first minimal Herglotz function be with −1/ρ < r 1 < −ρ. Again we must enforce the condition (67) for automorphicity. Now suppose the second minimal Herglotz function is then, in order that it is automorphic, we require Finally, for a period-3 solution the third minimal Herglotz function is where, for its automorphicity, we require It follows from (67), (75) and (77) that we must have which is consistent with the requirement that ρ < α < 1. On taking logarithms this is equivalent, by (63), to But there is another option. Instead of choosing M 2 (ζ) as in (80) we could alternatively pick implying the automorphicity condition Then, with M 3 (ζ) given by (76), conditions (67), (81) and (77) imply which is consistent with the requirement that ρ < α < 1. On taking logarithms this is equivalent, by (63), to 10 Two spectral bands of equal length The case of two spectral bands of equal length is of interest and has recently been studied, in a more general context, by Eichinger, Puchhammer & Yuditskii [24]. Our formulation can be applied to this special case. The condition (70) for a period-two Jacobi matrix is consistent with a spectrum comprising two equal intervals. To see this, consider the covering map given by It is readily verified that this map satisfies implying the required symmetry between the images of |ζ| = ρ and |ζ| = 1. For two spectral bands on the real axis between [−(λ + 2) 1/2 , −(λ − 2) 1/2 ] and [(λ − 2) 1/2 , (λ + 2) 1/2 ] where λ > 2 we must pick ρ such that and, for given λ, this is easily solved by Newton's method. Once ρ is determined R follows explicitly: With the covering map determined by the geometry of the spectral bands in this way it can be shown that where we defineP Also, with where we have used the expressions for M 1 (ζ) and M 2 (ζ) from (65) and (68) then it follows directly that where which, from (61), can be shown to satisfy On use of the results (easily checked) that K(−1) = 1/2, K(±ρ) = 0 we find With a second differentiation we find where we define It can be confirmed numerically that, for any choice of r 1 , In summary we find the following structure of period-2 Jacobi matrices over two equal intervals: where a and b are real parameters given by b = X ∞ √ ρ K(−1) + K(−ρ) + K(1/( √ ρr 1 )) − 1 − 2K(ρ) − K( √ ρ/r 1 ) , For a given spectrum, the parameters ρ, X ∞ and Y ∞ will be set leaving the parameters a and b in (99) dependent in the explicit way given above on the single real parameter r 1 . Indeed we are also able to confirm numerically that, for any r 1 , the formulas for a and b given in (99) satisfy the relation established in an appendix to [24]. Just as P (ζ) can be related to the first Jacobi theta function the functions K(ζ) and L(ζ) relate to the Weierstrass ζ and ℘ functions respectively so our results (99) can be viewed as expressions of the Jacobi parameters in terms of those special functions.

Discussion
This paper has presented the basis of a constructive approach to the inverse problem for finite-gap Jacobi matrices based on a Schottky model of the underlying hyperelliptic spectral curve and use of the Schottky-Klein prime function. The latter function can now be calculated with great speed and computational efficiency by means of novel numerical algorithms (see [17] [23] [20])) which are not based on the calculation of infinite products over the Schottky group action, or of any infinite Blaschke products. While, for the groups relevant here, all those products are known to converge (for example by some of Beardon's results -see [9] for references), they are not numerically efficient as a basis for constructive methods. The single gap case can be performed quite explicitly within our approach on use of a convenient infinite product formula for the Schottky-Klein prime function associated with a concentric annulus. The construction for more spectral bands can be performed entirely analogously, and with great numerical efficiency given the new software routines to evaluate Schottky-Klein prime functions on the Schottky doubles of planar domains [23].
A closing remark is that our explicit form (4) for the Jacobi parameters solving the inverse problem for a given spectrum is interesting. The appearance of a second derivative of log N (ζ) is highly reminiscent of similar formulas in the isospectral theory of the KdV equation with the automorphic function N (ζ) ostensibly playing the role of a tau function.