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
|
#include "../machine.h"
typedef char integer1;
typedef short integer2;
/* Copyright INRIA */
#define TRIU(Type) {\
Type *DX;\
DX=(Type *)dx;\
kk = *k;\
nn = *n;\
l = 0;\
if (kk>0) {\
for (i = 0; i < (*m)*kk; ++i) \
DX[i] = 0;\
l = l+ (*m)*kk;\
nn = nn - kk;\
kk = 0;}\
ls = l +1 - kk;\
ll = *m -1 + kk;\
for (j = 1; j <= nn; ++j) {\
if (ll > 0) {\
for (i = 0; i < ll; ++i) \
DX[ls+i] = 0;\
ll--;\
ls = ls + *m +1;\
}\
}\
}
int C2F(gentriu)(typ,m,n,k,dx)
integer *n;
integer *m;
integer *k;
integer *typ;
int *dx;
{
static integer kk, l, ls, ll, j, i, nn, c0=0, c1=1;
switch (*typ) {
case 1:
TRIU(integer1);
break;
case 2:
TRIU(integer2);
break;
case 4:
TRIU(integer);
break;
case 11:
TRIU(unsigned char);
break;
case 12:
TRIU(unsigned short);
break;
case 14:
TRIU(unsigned int);
break;
}
}
|
