kusano 7d535a
1c1
kusano 7d535a
< /* mc64ad.f -- translated by f2c (version 20100827).
kusano 7d535a
---
kusano 7d535a
> /* mc64ad.f -- translated by f2c (version 20061008).
kusano 7d535a
13c13
kusano 7d535a
< #include "slu_ddefs.h"
kusano 7d535a
---
kusano 7d535a
> #include "superlu_ddefs.h"
kusano 7d535a
15,16d14
kusano 7d535a
< #define abs(a) ((a) > 0) ? (a) : -(a)
kusano 7d535a
< #define min(a,b) ((a) < (b)) ? (a) : (b)
kusano 7d535a
22c20,22
kusano 7d535a
< 
kusano 7d535a
---
kusano 7d535a
> /* extern double dlamch_(char *, int); */
kusano 7d535a
> #define abs(a) ((a) > 0) ? (a) : -(a)
kusano 7d535a
> #define min(a,b) ((a) < (b)) ? (a) : (b)
kusano 7d535a
51c51
kusano 7d535a
< /* Subroutine */ int mc64id_(int_t *icntl)
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64id(int_t *icntl)
kusano 7d535a
53c53
kusano 7d535a
<     int_t i__;
kusano 7d535a
---
kusano 7d535a
>     static int_t i__;
kusano 7d535a
115c115
kusano 7d535a
< /* Subroutine */ int mc64ad_(int_t *job, int_t *n, int_t *ne, int_t *
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64ad(int_t *job, int_t *n, int_t *ne, int_t *
kusano 7d535a
128,131c128,131
kusano 7d535a
<     int_t i__, j, k;
kusano 7d535a
<     double fact, rinf;
kusano 7d535a
<     extern /* Subroutine */ int_t mc21ad_(int_t *, int_t *, int_t *, 
kusano 7d535a
< 	    int_t *, int_t *, int_t *, int_t *, int_t *), mc64bd_(
kusano 7d535a
---
kusano 7d535a
>     static int_t i__, j, k;
kusano 7d535a
>     static double fact, rinf;
kusano 7d535a
>     extern /* Subroutine */ int_t mc21ad(int_t *, int_t *, int_t *, 
kusano 7d535a
> 	    int_t *, int_t *, int_t *, int_t *, int_t *), mc64bd(
kusano 7d535a
133,140c133,140
kusano 7d535a
< 	    *, int_t *, int_t *, int_t *, int_t *, int_t *, double *),
kusano 7d535a
< 	    mc64rd_(int_t *, int_t *, int_t *, int_t *, double *),
kusano 7d535a
< 	    mc64sd_(int_t *, int_t *, int_t *, int_t *, double *, int_t *,
kusano 7d535a
< 	    int_t *, int_t *, int_t *, 
kusano 7d535a
< 	    int_t *, int_t *, int_t *, int_t *, int_t *),
kusano 7d535a
< 	    mc64wd_(int_t *, int_t *, int_t *, int_t *, double *, int_t *,
kusano 7d535a
< 	    int_t *, int_t *, int_t *, int_t *, int_t *, int_t *,
kusano 7d535a
< 	    double *, double *);
kusano 7d535a
---
kusano 7d535a
> 	    *, int_t *, int_t *, int_t *, int_t *, int_t *, 
kusano 7d535a
> 	    double *), mc64rd(int_t *, int_t *, int_t *, int_t *,
kusano 7d535a
> 	     double *), mc64sd(int_t *, int_t *, int_t *, int_t *
kusano 7d535a
> 	    , double *, int_t *, int_t *, int_t *, int_t *, 
kusano 7d535a
> 	    int_t *, int_t *, int_t *, int_t *, int_t *), mc64wd(
kusano 7d535a
> 	    int_t *, int_t *, int_t *, int_t *, double *, int_t 
kusano 7d535a
> 	    *, int_t *, int_t *, int_t *, int_t *, int_t *, int_t 
kusano 7d535a
> 	    *, double *, double *);
kusano 7d535a
341,344c341
kusano 7d535a
< 	if (icntl[1] >= 0) {
kusano 7d535a
< 	    printf(" ****** Error in MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 		   " because JOB = %d\n", *job);
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9001) INFO(1),'JOB',JOB */
kusano 7d535a
351,354c348
kusano 7d535a
< 	if (icntl[1] >= 0) {
kusano 7d535a
< 	    printf(" ****** Error in MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 		   " because N = %d\n", *job);
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9001) INFO(1),'N',N */
kusano 7d535a
361,364c355
kusano 7d535a
< 	if (icntl[1] >= 0) {
kusano 7d535a
< 	    printf(" ****** Error in MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 		   " because NE = %d\n", *job);
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9001) INFO(1),'NE',NE */
kusano 7d535a
386,389c377
kusano 7d535a
< 	if (icntl[1] >= 0) {
kusano 7d535a
< 	    printf(" ****** Error in MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 		   " LIW too small, must be at least %8d\n", k);
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9004) INFO(1),K */
kusano 7d535a
410,413c398
kusano 7d535a
< 	    if (icntl[1] >= 0) {
kusano 7d535a
< 		printf(" ****** Error in MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 		       " LDW too small, must be at least %8d\n", k);
kusano 7d535a
< 	    }
kusano 7d535a
---
kusano 7d535a
> /*          IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9005) INFO(1),K */
kusano 7d535a
433,437c418
kusano 7d535a
< 		    if (icntl[1] >= 0) {
kusano 7d535a
< 			printf(" ****** Error in MC64A/AD. INFO(1) = %2d",
kusano 7d535a
< 			       info[1], " Column %8d", j, 
kusano 7d535a
< 			       " contains an entry with invalid row index %8d\n", i__);
kusano 7d535a
< 		    }
kusano 7d535a
---
kusano 7d535a
> /*              IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9006) INFO(1),J,I */
kusano 7d535a
444,447c425
kusano 7d535a
< 		    if (icntl[1] >= 0) {
kusano 7d535a
< 			printf(" ****** Error in MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 			       "        Column %8d", j,
kusano 7d535a
< 			       " contains two or more entries with row index %8d\n", i__);
kusano 7d535a
---
kusano 7d535a
> /*              IF (ICNTL(1).GE.0) WRITE(ICNTL(1),9007) INFO(1),J,I */
kusano 7d535a
459,480c437,440
kusano 7d535a
< 	printf("  ****** Input parameters for MC64A/AD: JOB = %8d,"
kusano 7d535a
< 	       " N = %d, NE = %8d\n", *job, *n, *ne);
kusano 7d535a
< 	printf(" IP(1:N+1)   = ");
kusano 7d535a
< 	for (j=1; j<=(*n+1); ++j) {
kusano 7d535a
< 	    printf("%8d", ip[j]);
kusano 7d535a
< 	    if (j%8 == 0) printf("\n");
kusano 7d535a
< 	}
kusano 7d535a
< 	printf("\n IRN(1:NE) = ");
kusano 7d535a
< 	for (j=1; j<=(*ne); ++j) {
kusano 7d535a
< 	    printf("%8d", irn[j]);
kusano 7d535a
< 	    if (j%8 == 0) printf("\n");
kusano 7d535a
< 	}
kusano 7d535a
< 	printf("\n");
kusano 7d535a
< 
kusano 7d535a
< 	if (*job > 1) {
kusano 7d535a
< 	    printf(" A(1:NE)     = ");
kusano 7d535a
< 	    for (j=1; j<=(*ne); ++j) {
kusano 7d535a
< 		printf("%f14.4", a[j]);
kusano 7d535a
< 		if (j%4 == 0) printf("\n");
kusano 7d535a
< 	    }
kusano 7d535a
< 	    printf("\n");
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        WRITE(ICNTL(3),9020) JOB,N,NE */
kusano 7d535a
> /*        WRITE(ICNTL(3),9021) (IP(J),J=1,N+1) */
kusano 7d535a
> /*        WRITE(ICNTL(3),9022) (IRN(J),J=1,NE) */
kusano 7d535a
> /*        IF (JOB.GT.1) WRITE(ICNTL(3),9023) (A(J),J=1,NE) */
kusano 7d535a
496c456
kusano 7d535a
< 	mc21ad_(n, &irn[1], ne, &ip[1], &iw[1], &cperm[1], num, &iw[*n + 1]);
kusano 7d535a
---
kusano 7d535a
> 	mc21ad(n, &irn[1], ne, &ip[1], &iw[1], &cperm[1], num, &iw[*n + 1]);
kusano 7d535a
502c462
kusano 7d535a
< 	mc64bd_(n, ne, &ip[1], &irn[1], &a[1], &cperm[1], num, &iw[1], &iw[*n 
kusano 7d535a
---
kusano 7d535a
> 	mc64bd(n, ne, &ip[1], &irn[1], &a[1], &cperm[1], num, &iw[1], &iw[*n 
kusano 7d535a
516c476
kusano 7d535a
< 	mc64rd_(n, ne, &ip[1], &iw[1], &dw[1]);
kusano 7d535a
---
kusano 7d535a
> 	mc64rd(n, ne, &ip[1], &iw[1], &dw[1]);
kusano 7d535a
518c478
kusano 7d535a
< 	mc64sd_(n, ne, &ip[1], &iw[1], &dw[1], &cperm[1], num, &iw[*ne + 1], &
kusano 7d535a
---
kusano 7d535a
> 	mc64sd(n, ne, &ip[1], &iw[1], &dw[1], &cperm[1], num, &iw[*ne + 1], &
kusano 7d535a
543c503
kusano 7d535a
< 	mc64wd_(n, ne, &ip[1], &irn[1], &dw[(*n << 1) + 1], &cperm[1], num, &
kusano 7d535a
---
kusano 7d535a
> 	mc64wd(n, ne, &ip[1], &irn[1], &dw[(*n << 1) + 1], &cperm[1], num, &
kusano 7d535a
578c538
kusano 7d535a
< 	mc64wd_(n, ne, &ip[1], &irn[1], &dw[*n * 3 + 1], &cperm[1], num, &iw[
kusano 7d535a
---
kusano 7d535a
> 	mc64wd(n, ne, &ip[1], &irn[1], &dw[*n * 3 + 1], &cperm[1], num, &iw[
kusano 7d535a
610,613c570
kusano 7d535a
< 	if (icntl[2] >= 0) {
kusano 7d535a
< 	    printf(" ****** Warning from MC64A/AD. INFO(1) = %2d", info[1],
kusano 7d535a
< 		   " The matrix is structurally singular.\n");
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        IF (ICNTL(2).GE.0) WRITE(ICNTL(2),9011) INFO(1) */
kusano 7d535a
617,620c574
kusano 7d535a
< 	if (icntl[2] >= 0) {
kusano 7d535a
< 	    printf(" ****** Warning from MC64A/AD. INFO(1) = %2d\n", info[1],
kusano 7d535a
< 		   "        Some scaling factors may be too large.\n");
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        IF (ICNTL(2).GE.0) WRITE(ICNTL(2),9012) INFO(1) */
kusano 7d535a
624,644c578,584
kusano 7d535a
< 	printf(" ****** Output parameters for MC64A/AD: INFO(1:2)  = %8d%8d\n",
kusano 7d535a
< 	       info[1], info[2]);
kusano 7d535a
< 	printf(" NUM        = ", *num);
kusano 7d535a
< 	printf(" CPERM(1:N) = ");
kusano 7d535a
< 	for (j=1; j<=*n; ++j) {
kusano 7d535a
< 	    printf("%8d", cperm[j]);
kusano 7d535a
< 	    if (j%8 == 0) printf("\n");
kusano 7d535a
< 	}
kusano 7d535a
< 	if (*job == 5) {
kusano 7d535a
< 	    printf("\n DW(1:N)    = ");
kusano 7d535a
< 	    for (j=1; j<=*n; ++j) {
kusano 7d535a
< 		printf("%11.3f", dw[j]);
kusano 7d535a
< 		if (j%5 == 0) printf("\n");
kusano 7d535a
< 	    }
kusano 7d535a
< 	    printf("\n DW(N+1:2N) = ");
kusano 7d535a
< 	    for (j=1; j<=*n; ++j) {
kusano 7d535a
< 		printf("%11.3f", dw[*n+j]);
kusano 7d535a
< 		if (j%5 == 0) printf("\n");
kusano 7d535a
< 	    }
kusano 7d535a
< 	    printf("\n");
kusano 7d535a
< 	}
kusano 7d535a
---
kusano 7d535a
> /*        WRITE(ICNTL(3),9030) (INFO(J),J=1,2) */
kusano 7d535a
> /*        WRITE(ICNTL(3),9031) NUM */
kusano 7d535a
> /*        WRITE(ICNTL(3),9032) (CPERM(J),J=1,N) */
kusano 7d535a
> /*        IF (JOB.EQ.5) THEN */
kusano 7d535a
> /*          WRITE(ICNTL(3),9033) (DW(J),J=1,N) */
kusano 7d535a
> /*          WRITE(ICNTL(3),9034) (DW(N+J),J=1,N) */
kusano 7d535a
> /*        ENDIF */
kusano 7d535a
648a589,604
kusano 7d535a
> /* L9001: */
kusano 7d535a
> /* L9004: */
kusano 7d535a
> /* L9005: */
kusano 7d535a
> /* L9006: */
kusano 7d535a
> /* L9007: */
kusano 7d535a
> /* L9011: */
kusano 7d535a
> /* L9012: */
kusano 7d535a
> /* L9020: */
kusano 7d535a
> /* L9021: */
kusano 7d535a
> /* L9022: */
kusano 7d535a
> /* L9023: */
kusano 7d535a
> /* L9030: */
kusano 7d535a
> /* L9031: */
kusano 7d535a
> /* L9032: */
kusano 7d535a
> /* L9033: */
kusano 7d535a
> /* L9034: */
kusano 7d535a
652c608
kusano 7d535a
< /* Subroutine */ int_t mc64bd_(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64bd(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
661,678c617,635
kusano 7d535a
<     int_t i__, j, k;
kusano 7d535a
<     double a0;
kusano 7d535a
<     int_t i0, q0;
kusano 7d535a
<     double ai, di;
kusano 7d535a
<     int_t ii, jj, kk;
kusano 7d535a
<     double bv;
kusano 7d535a
<     int_t up;
kusano 7d535a
<     double dq0;
kusano 7d535a
<     int_t kk1, kk2;
kusano 7d535a
<     double csp;
kusano 7d535a
<     int_t isp, jsp, low;
kusano 7d535a
<     double dnew;
kusano 7d535a
<     int_t jord, qlen, idum, jdum;
kusano 7d535a
<     double rinf;
kusano 7d535a
<     extern /* Subroutine */ int mc64dd_(int_t *, int_t *, int_t *, 
kusano 7d535a
< 	    double *, int_t *, int_t *), mc64ed_(int_t *, int_t *,
kusano 7d535a
< 	     int_t *, double *, int_t *, int_t *), mc64fd_(int_t *
kusano 7d535a
< 	    , int_t *, int_t *, int_t *, double *, int_t *, int_t *);
kusano 7d535a
---
kusano 7d535a
>     static int_t i__, j, k;
kusano 7d535a
>     static double a0;
kusano 7d535a
>     static int_t i0, q0;
kusano 7d535a
>     static double ai, di;
kusano 7d535a
>     static int_t ii, jj, kk;
kusano 7d535a
>     static double bv;
kusano 7d535a
>     static int_t up;
kusano 7d535a
>     static double dq0;
kusano 7d535a
>     static int_t kk1, kk2;
kusano 7d535a
>     static double csp;
kusano 7d535a
>     static int_t isp, jsp, low;
kusano 7d535a
>     static double dnew;
kusano 7d535a
>     static int_t jord, qlen, idum, jdum;
kusano 7d535a
>     static double rinf;
kusano 7d535a
>     extern /* Subroutine */ int_t mc64dd(int_t *, int_t *, int_t *, 
kusano 7d535a
> 	    double *, int_t *, int_t *), mc64ed(int_t *, int_t *,
kusano 7d535a
> 	     int_t *, double *, int_t *, int_t *), mc64fd(int_t *
kusano 7d535a
> 	    , int_t *, int_t *, int_t *, double *, int_t *, 
kusano 7d535a
> 	    int_t *);
kusano 7d535a
892c849
kusano 7d535a
< 		    mc64dd_(&i__, n, &q[1], &d__[1], &l[1], &c__1);
kusano 7d535a
---
kusano 7d535a
> 		    mc64dd(&i__, n, &q[1], &d__[1], &l[1], &c__1);
kusano 7d535a
914c871
kusano 7d535a
< 		    mc64ed_(&qlen, n, &q[1], &d__[1], &l[1], &c__1);
kusano 7d535a
---
kusano 7d535a
> 		    mc64ed(&qlen, n, &q[1], &d__[1], &l[1], &c__1);
kusano 7d535a
967c924
kusano 7d535a
< 			    mc64fd_(&l[i__], &qlen, n, &q[1], &d__[1], &l[1], 
kusano 7d535a
---
kusano 7d535a
> 			    mc64fd(&l[i__], &qlen, n, &q[1], &d__[1], &l[1], 
kusano 7d535a
979c936
kusano 7d535a
< 			mc64dd_(&i__, n, &q[1], &d__[1], &l[1], &c__1);
kusano 7d535a
---
kusano 7d535a
> 			mc64dd(&i__, n, &q[1], &d__[1], &l[1], &c__1);
kusano 7d535a
1079c1036
kusano 7d535a
< /* Subroutine */ int_t mc64dd_(int_t *i__, int_t *n, int_t *q, double 
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64dd(int_t *i__, int_t *n, int_t *q, double 
kusano 7d535a
1086,1087c1043,1044
kusano 7d535a
<     double di;
kusano 7d535a
<     int_t qk, pos, idum, posk;
kusano 7d535a
---
kusano 7d535a
>     static double di;
kusano 7d535a
>     static int_t qk, pos, idum, posk;
kusano 7d535a
1157c1114
kusano 7d535a
< /* Subroutine */ int_t mc64ed_(int_t *qlen, int_t *n, int_t *q, 
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64ed(int_t *qlen, int_t *n, int_t *q, 
kusano 7d535a
1164,1166c1121,1123
kusano 7d535a
<     int_t i__;
kusano 7d535a
<     double di, dk, dr;
kusano 7d535a
<     int_t pos, idum, posk;
kusano 7d535a
---
kusano 7d535a
>     static int_t i__;
kusano 7d535a
>     static double di, dk, dr;
kusano 7d535a
>     static int_t pos, idum, posk;
kusano 7d535a
1252c1209
kusano 7d535a
< /* Subroutine */ int_t mc64fd_(int_t *pos0, int_t *qlen, int_t *n, 
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64fd(int_t *pos0, int_t *qlen, int_t *n, 
kusano 7d535a
1259,1261c1216,1218
kusano 7d535a
<     int_t i__;
kusano 7d535a
<     double di, dk, dr;
kusano 7d535a
<     int_t qk, pos, idum, posk;
kusano 7d535a
---
kusano 7d535a
>     static int_t i__;
kusano 7d535a
>     static double di, dk, dr;
kusano 7d535a
>     static int_t qk, pos, idum, posk;
kusano 7d535a
1390c1347
kusano 7d535a
< /* Subroutine */ int_t mc64rd_(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64rd(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
1397,1401c1354,1358
kusano 7d535a
<     int_t j, k, r__, s;
kusano 7d535a
<     double ha;
kusano 7d535a
<     int_t hi, td, mid, len, ipj;
kusano 7d535a
<     double key;
kusano 7d535a
<     int_t last, todo[50], first;
kusano 7d535a
---
kusano 7d535a
>     static int_t j, k, r__, s;
kusano 7d535a
>     static double ha;
kusano 7d535a
>     static int_t hi, td, mid, len, ipj;
kusano 7d535a
>     static double key;
kusano 7d535a
>     static int_t last, todo[50], first;
kusano 7d535a
1538c1495
kusano 7d535a
< /* Subroutine */ int_t mc64sd_(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64sd(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
1547,1550c1504,1507
kusano 7d535a
<     int_t i__, j, k, l, ii, mod, cnt, num;
kusano 7d535a
<     double bval, bmin, bmax, rinf;
kusano 7d535a
<     int_t nval, wlen, idum1, idum2, idum3;
kusano 7d535a
<     extern /* Subroutine */ int_t mc64qd_(int_t *, int_t *, int_t *, 
kusano 7d535a
---
kusano 7d535a
>     static int_t i__, j, k, l, ii, mod, cnt, num;
kusano 7d535a
>     static double bval, bmin, bmax, rinf;
kusano 7d535a
>     static int_t nval, wlen, idum1, idum2, idum3;
kusano 7d535a
>     extern /* Subroutine */ int_t mc64qd(int_t *, int_t *, int_t *, 
kusano 7d535a
1552c1509
kusano 7d535a
< 	    mc64ud_(int_t *, int_t *, int_t *, int_t *, int_t *, 
kusano 7d535a
---
kusano 7d535a
> 	    mc64ud(int_t *, int_t *, int_t *, int_t *, int_t *, 
kusano 7d535a
1555a1513
kusano 7d535a
> 
kusano 7d535a
1627c1585
kusano 7d535a
<     mc64ud_(&cnt, &mod, n, &irn[1], ne, &ip[1], &len[1], &fc[1], &iw[1], numx,
kusano 7d535a
---
kusano 7d535a
>     mc64ud(&cnt, &mod, n, &irn[1], ne, &ip[1], &len[1], &fc[1], &iw[1], numx,
kusano 7d535a
1706c1664
kusano 7d535a
< 		mc64qd_(&ip[1], &lenl[1], &len[1], &w[1], &wlen, &a[1], &nval,
kusano 7d535a
---
kusano 7d535a
> 		mc64qd(&ip[1], &lenl[1], &len[1], &w[1], &wlen, &a[1], &nval,
kusano 7d535a
1767c1725
kusano 7d535a
< 	    mc64qd_(&ip[1], &len[1], &lenh[1], &w[1], &wlen, &a[1], &nval, &
kusano 7d535a
---
kusano 7d535a
> 	    mc64qd(&ip[1], &len[1], &lenh[1], &w[1], &wlen, &a[1], &nval, &
kusano 7d535a
1806c1764
kusano 7d535a
< 	mc64ud_(&cnt, &mod, n, &irn[1], ne, &ip[1], &len[1], &fc[1], &iw[1], &
kusano 7d535a
---
kusano 7d535a
> 	mc64ud(&cnt, &mod, n, &irn[1], ne, &ip[1], &len[1], &fc[1], &iw[1], &
kusano 7d535a
1854c1812
kusano 7d535a
< /* Subroutine */ int_t mc64qd_(int_t *ip, int_t *lenl, int_t *lenh, 
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64qd(int_t *ip, int_t *lenl, int_t *lenh, 
kusano 7d535a
1862,1865c1820,1823
kusano 7d535a
<     int_t j, k, s;
kusano 7d535a
<     double ha;
kusano 7d535a
<     int_t ii, pos;
kusano 7d535a
<     double split[10];
kusano 7d535a
---
kusano 7d535a
>     static int_t j, k, s;
kusano 7d535a
>     static double ha;
kusano 7d535a
>     static int_t ii, pos;
kusano 7d535a
>     static double split[10];
kusano 7d535a
1949c1907
kusano 7d535a
< /* Subroutine */ int mc64ud_(int_t *id, int_t *mod, int_t *n, int_t *
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64ud(int_t *id, int_t *mod, int_t *n, int_t *
kusano 7d535a
1958,1959c1916,1917
kusano 7d535a
<     int_t i__, j, k, j1, ii, kk, id0, id1, in1, in2, nfc, num0, num1, num2, 
kusano 7d535a
< 	    jord, last;
kusano 7d535a
---
kusano 7d535a
>     static int_t i__, j, k, j1, ii, kk, id0, id1, in1, in2, nfc, num0, num1,
kusano 7d535a
> 	     num2, jord, last;
kusano 7d535a
2028c1986
kusano 7d535a
< /* Integers ID0+1 to ID0+N are unique numbers for call ID to MC64U/UD, */
kusano 7d535a
---
kusano 7d535a
> /* Ints ID0+1 to ID0+N are unique numbers for call ID to MC64U/UD, */
kusano 7d535a
2152c2110
kusano 7d535a
< /* Subroutine */ int_t mc64wd_(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
---
kusano 7d535a
> /* Subroutine */ int_t mc64wd(int_t *n, int_t *ne, int_t *ip, int_t *
kusano 7d535a
2161,2175c2119,2133
kusano 7d535a
<     int_t i__, j, k, i0, k0, k1, k2, q0;
kusano 7d535a
<     double di;
kusano 7d535a
<     int_t ii, jj, kk;
kusano 7d535a
<     double vj;
kusano 7d535a
<     int_t up;
kusano 7d535a
<     double dq0;
kusano 7d535a
<     int_t kk1, kk2;
kusano 7d535a
<     double csp;
kusano 7d535a
<     int_t isp, jsp, low;
kusano 7d535a
<     double dmin__, dnew;
kusano 7d535a
<     int_t jord, qlen, jdum;
kusano 7d535a
<     double rinf;
kusano 7d535a
<     extern /* Subroutine */ int_t mc64dd_(int_t *, int_t *, int_t *, 
kusano 7d535a
< 	    double *, int_t *, int_t *), mc64ed_(int_t *, int_t *,
kusano 7d535a
< 	     int_t *, double *, int_t *, int_t *), mc64fd_(int_t *
kusano 7d535a
---
kusano 7d535a
>     static int_t i__, j, k, i0, k0, k1, k2, q0;
kusano 7d535a
>     static double di;
kusano 7d535a
>     static int_t ii, jj, kk;
kusano 7d535a
>     static double vj;
kusano 7d535a
>     static int_t up;
kusano 7d535a
>     static double dq0;
kusano 7d535a
>     static int_t kk1, kk2;
kusano 7d535a
>     static double csp;
kusano 7d535a
>     static int_t isp, jsp, low;
kusano 7d535a
>     static double dmin__, dnew;
kusano 7d535a
>     static int_t jord, qlen, jdum;
kusano 7d535a
>     static double rinf;
kusano 7d535a
>     extern /* Subroutine */ int_t mc64dd(int_t *, int_t *, int_t *, 
kusano 7d535a
> 	    double *, int_t *, int_t *), mc64ed(int_t *, int_t *,
kusano 7d535a
> 	     int_t *, double *, int_t *, int_t *), mc64fd(int_t *
kusano 7d535a
2432c2390
kusano 7d535a
< 		mc64dd_(&i__, n, &q[1], &d__[1], &l[1], &c__2);
kusano 7d535a
---
kusano 7d535a
> 		mc64dd(&i__, n, &q[1], &d__[1], &l[1], &c__2);
kusano 7d535a
2454c2412
kusano 7d535a
< 		mc64ed_(&qlen, n, &q[1], &d__[1], &l[1], &c__2);
kusano 7d535a
---
kusano 7d535a
> 		mc64ed(&qlen, n, &q[1], &d__[1], &l[1], &c__2);
kusano 7d535a
2508c2466
kusano 7d535a
< 			    mc64fd_(&l[i__], &qlen, n, &q[1], &d__[1], &l[1], 
kusano 7d535a
---
kusano 7d535a
> 			    mc64fd(&l[i__], &qlen, n, &q[1], &d__[1], &l[1], 
kusano 7d535a
2519c2477
kusano 7d535a
< 			mc64dd_(&i__, n, &q[1], &d__[1], &l[1], &c__2);
kusano 7d535a
---
kusano 7d535a
> 			mc64dd(&i__, n, &q[1], &d__[1], &l[1], &c__2);