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
|
//Test inputs for ex15f:
X=1:3;Y=4:6;
//...and 3 examples of Scilab argument functions:
deff('z=f1(x,y)','z=x+y');
deff('z=f2(x,y)','z=x*x+y*y');
deff('z=f3(x,y)','z=cos(x+y)');
//reference values for Z calculated by feval
Zref1=feval(X,Y,f1);
Zref2=feval(X,Y,f2);
Zref3=feval(X,Y,f3);
// Calling ex15f with a argument function which is a Scilab function:
Z1=ex15f(X,Y,f1);
if norm(Z1-Zref1) > %eps then pause,end
Z2=ex15f(X,Y,f2);
if norm(Z2-Zref2) > %eps then pause,end
Z3=ex15f(X,Y,f3);
if norm(Z3-Zref3) > %eps then pause,end
// Calling ex15f with a argument function which is a Fortran function:
Z1=ex15f(X,Y,'f1f'); // f1f defined in ex15f.f
if norm(Z1 - Zref1) > %eps then pause,end
Z2=ex15f(X,Y,'f2f'); // f2f defined in ex15f.f
if norm(Z2 - Zref2) > %eps then pause,end
// making f3f.o and linking f3f dynamically with Scilab
if ~c_link('f3f') then
// see link-examples-so for link examples
ilib_for_link('f3f','f3f.o',[],'f','f3fmake','f3floader.sce');
exec f3floader.sce;
end
// Calling ex15f with the Fortran argument function fp3
Z3=ex15f(X,Y,'f3f');
if norm(Z3 - Zref3) > %eps then pause,end
|
