1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
|
#include <string.h>
#include <malloc.h>
#include "../machine.h"
void C2F(maxcpl)(c,ex,la2,lp2,ma,mm,n,or,x,z)
int *ex,*la2,*lp2,*ma,*mm,*n,*or,*x;
double *c,*z;
{
int *alphi,*beta,*deg,*hp1,*hp2,*ihp1,*ihp2,*imin;
int *ind,*l1,*l2,*nlp2,*m1,*m2,*q,*s1,*s2,*sat;
double *p1w,*p2w,*pimin,*pivw,*r;
int i,ndim,n1dim;
int isize = sizeof(int);
int dsize = sizeof(double);
ndim = 3 * *n / 2 + 1;
n1dim = ndim + 1;
if ((alphi = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((beta = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((deg = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((hp1 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((hp2 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((ihp1 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((ihp2 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((imin = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((ind = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((l1 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((l2 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((nlp2 = (int *)malloc(n1dim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
for (i = 0; i < (*n + 1); i++) nlp2[i]=lp2[i];
if ((m1 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((m2 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((p1w = (double *)malloc(ndim * dsize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((p2w = (double *)malloc(ndim * dsize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((pimin = (double *)malloc(ndim * dsize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((pivw = (double *)malloc(ndim * dsize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((q = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((r = (double *)malloc(ndim * dsize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((s1 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((s2 = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
if ((sat = (int *)malloc(ndim * isize)) == NULL) {
cerro("Running out of memory");
return;
}
F2C(cplmax)(alphi,beta,c,deg,ex,hp1,hp2,ihp1,ihp2,imin,
ind,l1,l2,la2,nlp2,m1,m2,ma,mm,n,&n1dim,&ndim,or,
p1w,p2w,pimin,pivw,q,r,s1,s2,sat,x,z);
free(alphi); free(beta); free(deg); free(hp1); free(hp2);
free(ihp1); free(ihp2); free(imin); free(ind); free(l1);
free(l2); free(nlp2); free(m1); free(m2); free(p1w);
free(p2w); free(pimin); free(pivw); free(q); free(r);
free(s1); free(s2); free(sat);
}
|
