blob: adaa26d9d4dd2ace0006195303d2c9bb8faa9b06 (
plain)
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
|
/*
* Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
* Copyright (C) INRIA -
*
* This file must be used under the terms of the CeCILL.
* This source file is licensed as described in the file COPYING, which
* you should have received as part of this distribution. The terms
* are also available at
* http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
*
*/
#include "gentril.h"
#define TRIL(Type) {\
Type *DX;\
DX=(Type *)dx;\
kk = *k;\
nn = *n;\
l = 0;\
if (kk>=0) {\
l = l + (*m)*(kk+1);\
nn = nn - kk - 1;\
kk = -1;\
}\
ls = l;\
ll = -kk;\
for (j = 1; j <= nn; ++j) {\
if (ll > *m) ll = *m;\
for (i = 0; i < ll; ++i) \
DX[ls+i] = 0;\
ll++;\
ls = ls + *m;\
}\
}
int C2F(gentril)(int *typ, int *m, int *n, int *k, int *dx)
{
static int kk, l, ls, ll, j, i, nn;
switch (*typ)
{
case 1:
TRIL(integer1);
break;
case 2:
TRIL(integer2);
break;
case 4:
TRIL(int) ;
break;
case 11:
TRIL(unsigned char);
break;
case 12:
TRIL(unsigned short);
break;
case 14:
TRIL(unsigned int);
break;
}
return 0;
}
|
