summaryrefslogtreecommitdiff
path: root/modules/integer/src/c/gentriu.c
blob: 6c9dae3363635f72f86174ddc4d2350c9ba5e11e (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
66
67
/*
* 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 "gentriu.h"


#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)(int *typ, int *m, int *n, int *k, int *dx)
{
    static int kk, l, ls, ll, j, i, nn;

    switch (*typ)
    {
        case 1:
            TRIU(integer1);
            break;
        case 2:
            TRIU(integer2);
            break;
        case 4:
            TRIU(int) ;
            break;
        case 11:
            TRIU(unsigned char);
            break;
        case 12:
            TRIU(unsigned short);
            break;
        case 14:
            TRIU(unsigned int);
            break;
    }
    return 0;
}