1 /* -- translated by f2c (version 20100827).
2 You must link the resulting object file with libf2c:
3 on Microsoft Windows system, link with libf2c.lib;
4 on Linux or Unix systems, link with .../path/to/libf2c.a -lm
5 or, if you install libf2c.a in a standard place, with -lf2c -lm
6 -- in that order, at the end of the command line, as in
7 cc *.o -lf2c -lm
8 Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
9
10 http://www.netlib.org/f2c/libf2c.zip
11 */
12
13 #include "libtinyf2c.h"
14
15 /* ********************************************************************** */
16 /* FORWARD LEGENDRE TRANSFORMATION (UPPER LEVEL) */
17 /* ********************************************************************** */
18 /* Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved. */
19 /* ----------------------------------------------------------------------- */
shlfwu_(integer * mm,integer * jm,integer * m,integer * isw,real * wm,real * sm,real * sd,real * xs,real * xc,real * pm,real * wy,real * wx,real * px,real * r__,real * z__)20 /* Subroutine */ int shlfwu_(integer *mm, integer *jm, integer *m, integer *
21 isw, real *wm, real *sm, real *sd, real *xs, real *xc, real *pm, real
22 *wy, real *wx, real *px, real *r__, real *z__)
23 {
24 /* System generated locals */
25 integer sm_offset, wm_offset, pm_dim1, pm_offset, px_dim2, px_offset,
26 wy_dim1, wy_offset, wx_dim1, wx_offset, z_dim1, z_dim2, z_offset;
27
28 /* Local variables */
29 extern /* Subroutine */ int shppma_(integer *, integer *, integer *, real
30 *, real *, real *, real *), shlfwm_(integer *, integer *, integer
31 *, integer *, real *, real *, real *, real *, real *, real *,
32 real *, real *, real *);
33
34 /* Parameter adjustments */
35 --r__;
36 z_dim1 = *jm;
37 z_dim2 = *jm - 0 + 1;
38 z_offset = 1 + z_dim1 * (0 + z_dim2);
39 z__ -= z_offset;
40 px_dim2 = *jm - 0 + 1;
41 px_offset = 1 + 2 * (0 + px_dim2 * 0);
42 px -= px_offset;
43 wx_dim1 = *jm - 0 + 1;
44 wx_offset = 0 + wx_dim1;
45 wx -= wx_offset;
46 wy_dim1 = *jm - 0 + 1;
47 wy_offset = 0 + wy_dim1;
48 wy -= wy_offset;
49 pm_dim1 = *mm + 1 - 0 + 1;
50 pm_offset = 0 + pm_dim1 * 0;
51 pm -= pm_offset;
52 wm_offset = -(*jm);
53 wm -= wm_offset;
54 sm_offset = *m;
55 sm -= sm_offset;
56
57 /* Function Body */
58 shppma_(mm, jm, m, &pm[pm_offset], xs, &px[px_offset], &r__[1]);
59 shlfwm_(mm, jm, m, isw, &wm[wm_offset], &sm[sm_offset], sd, xc, &pm[
60 pm_offset], &wy[wy_offset], &wx[wx_offset], &r__[1], &z__[
61 z_offset]);
62 return 0;
63 } /* shlfwu_ */
64
65