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 /* Common Block Declarations */
16
17 struct {
18 logical lmiss;
19 real rmiss;
20 integer npm;
21 } szbpm1_;
22
23 #define szbpm1_1 szbpm1_
24
25 struct {
26 char cmark[1];
27 } szbpm2_;
28
29 #define szbpm2_1 szbpm2_
30
31 struct {
32 logical lclip;
33 } szbtx3_;
34
35 #define szbtx3_1 szbtx3_
36
37 /* ----------------------------------------------------------------------- */
38 /* Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved. */
39 /* ----------------------------------------------------------------------- */
szpmzr_(integer * n,real * rpx,real * rpy)40 /* Subroutine */ int szpmzr_(integer *n, real *rpx, real *rpy)
41 {
42 /* System generated locals */
43 integer i__1, i__2;
44
45 /* Local variables */
46 static integer i__;
47 static logical lflag;
48 extern /* Subroutine */ int stepr2_(void), strpr2_(void);
49 static logical lclipz;
50 extern /* Subroutine */ int sztxzv_(real *, real *, char *, ftnlen);
51
52 /* Parameter adjustments */
53 --rpy;
54 --rpx;
55
56 /* Function Body */
57 lclipz = szbtx3_1.lclip;
58 szbtx3_1.lclip = FALSE_;
59 stepr2_();
60 i__1 = *n;
61 i__2 = szbpm1_1.npm;
62 for (i__ = 1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ += i__2) {
63 lflag = szbpm1_1.lmiss && (rpx[i__] == szbpm1_1.rmiss || rpy[i__] ==
64 szbpm1_1.rmiss);
65 if (! lflag) {
66 sztxzv_(&rpx[i__], &rpy[i__], szbpm2_1.cmark, (ftnlen)1);
67 }
68 /* L10: */
69 }
70 szbtx3_1.lclip = lclipz;
71 strpr2_();
72 return 0;
73 } /* szpmzr_ */
74
75