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
|
/* Copyright INRIA */
/************************************
* (very) simple example 1
* with a C wrapper
* -->link('ext1c.o','ext1cI','C');
* -->a=[1,2,3];b=[4,5,6];n=3;
* -->c=fort('ext1cI',n,a,b)
* c=a+b
************************************/
#include "../../routines/stack-c.h"
int ext14cI(fname)
char* fname;
{
int m1,n1,l1,m2,n2,l2,m3,n3,l3,n,l4;
int minlhs=1, minrhs=3, maxlhs=1, maxrhs=3;
Nbvars = 0;
CheckRhs(minrhs,maxrhs) ;
CheckLhs(minlhs,maxlhs) ;
GetRhsVar( 1, "c", &m1, &n1, &l1);
GetRhsVar( 2, "d", &m2, &n2, &l2);
GetRhsVar( 3, "d", &m3, &n3, &l3);
if ( m3*n3 != m2*n2)
{
sciprint("%s :Incompatible dimensions\r\n",fname);
Error(999); return(0);
}
CreateVar( 4, "d", &m2, &n2, &l4);
n=m3*n3;
ext14c(&n,stk(l2),stk(l3),stk(l4));
LhsVar(1) = 4;
PutLhsVar();
return(0);
}
int ext14c(n, a, b, c)
int *n;
double *a, *b, *c;
{
int k;
for (k = 0; k < *n; ++k)
{
c[k] = a[k] + b[k];
}
return(0);
}
|
