File: genmdiag.c

package info (click to toggle)
scilab 4.0-12
  • links: PTS
  • area: non-free
  • in suites: etch, etch-m68k
  • size: 100,640 kB
  • ctags: 57,333
  • sloc: ansic: 377,889; fortran: 242,862; xml: 179,819; tcl: 42,062; sh: 10,593; ml: 9,441; makefile: 4,377; cpp: 1,354; java: 621; csh: 260; yacc: 247; perl: 130; lex: 126; asm: 72; lisp: 30
file content (60 lines) | stat: -rw-r--r-- 934 bytes parent folder | download | duplicates (2) 
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
 
#include "../machine.h"
typedef signed char integer1;
typedef short integer2;
/*     Copyright INRIA */
#define Min(x,y)	(((x)<(y))?(x):(y))
#define Max(x,y)	(((x)>(y))?(x):(y))

#define MDIAG(Type) {\
Type *DX;\
Type *DY;\
    DX=(Type *)dx;\
    --DX;\
    DY=(Type *)dy;\
    if (*k>=0)\
      ix=1+(*k)*(*m);\
    else\
      ix=1-(*k);\
    for (i = 0; i < mn; ++i) \
       DY[i] = DX[ix+i*(*m+1)];\
 }


int C2F(genmdiag)(typ,m,n,k,dx,dy) 
integer *n;
integer *m;
integer *k;
integer *typ;
int *dx;
int *dy;
{
  static integer i, ix, mn;

  if (*k >= 0) 
    mn=Max(0,Min(*m,*n-*k));
  else
    mn=Max(0,Min(*m+*k,*n));

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