generic/packages/bisect.asc

(* This file implements a bisection search *)

IMPORT Bisection FROM libbisect;
(* Bisection is the name of the init function *)

MODEL bisect__base;

   ninputs IS_A integer;
   x[1..ninputs] IS_A generic_real;
   x_par[1..ninputs] IS_A generic_real;
    
END bisect__base;   

MODEL bisect;

   b IS_A bisect__base;
   ninputs IS_A integer;
   x[1..ninputs] IS_A generic_real;
   x_par[1..ninputs] IS_A generic_real;

   (* wire up the terms *)
   
   ninputs, b.ninputs ARE_THE_SAME;
   FOR i IN [1..ninputs] CREATE
    x[i], b.x[i]         ARE_THE_SAME;
    x_par[i], b.x_par[i] ARE_THE_SAME;
   END;
   
   y[1..ninputs] IS_A generic_real; (* this is the result vector *)
   
END bisect;


MODEL test_bisect REFINES bisect;

   ninputs := 4;
   myset IS_A set OF symbol;
   myset := ['mary','jane','sue'];
   multipliers[myset]  IS_A generic_real;
   
   multipliers['mary'] := 1.0;
   multipliers['jane'] := 2.5;
   multipliers['sue']  := 3.9;      
   
   INITIALIZATION
   
   PROCEDURE set_values;
    x[1] := 2.3;
    x[2] := 0.5657;
    x[3] := 0.0;
    x[4] := 67.893;
   END set_values;
   
   PROCEDURE set_values_external;
     FOR i IN myset DO
        EXTERNAL do_set_values(x[1..ninputs],x_par[1..ninputs],
                           multipliers[i]);
     END;
   END set_values_external;
   
   PROCEDURE bisection;
     EXTERNAL do_bisection(x[1..ninputs],x_par[1..ninputs],y[1..ninputs]);   
   END bisection;
   
   PROCEDURE all;
     RUN set_values_external;
     RUN bisection;
   END all;
   
END test_bisect;
   
1	(* This file implements a bisection search *)
2
3	IMPORT Bisection FROM libbisect;
4	(* Bisection is the name of the init function *)
5
6	MODEL bisect__base;
7
8	ninputs IS_A integer;
9	x[1..ninputs] IS_A generic_real;
10	x_par[1..ninputs] IS_A generic_real;
11
12	END bisect__base;
13
14	MODEL bisect;
15
16	b IS_A bisect__base;
17	ninputs IS_A integer;
18	x[1..ninputs] IS_A generic_real;
19	x_par[1..ninputs] IS_A generic_real;
20
21	(* wire up the terms *)
22
23	ninputs, b.ninputs ARE_THE_SAME;
24	FOR i IN [1..ninputs] CREATE
25	x[i], b.x[i] ARE_THE_SAME;
26	x_par[i], b.x_par[i] ARE_THE_SAME;
27	END;
28
29	y[1..ninputs] IS_A generic_real; (* this is the result vector *)
30
31	END bisect;
32
33
34	MODEL test_bisect REFINES bisect;
35
36	ninputs := 4;
37	myset IS_A set OF symbol;
38	myset := ['mary','jane','sue'];
39	multipliers[myset] IS_A generic_real;
40
41	multipliers['mary'] := 1.0;
42	multipliers['jane'] := 2.5;
43	multipliers['sue'] := 3.9;
44
45	INITIALIZATION
46
47	PROCEDURE set_values;
48	x[1] := 2.3;
49	x[2] := 0.5657;
50	x[3] := 0.0;
51	x[4] := 67.893;
52	END set_values;
53
54	PROCEDURE set_values_external;
55	FOR i IN myset DO
56	EXTERNAL do_set_values(x[1..ninputs],x_par[1..ninputs],
57	multipliers[i]);
58	END;
59	END set_values_external;
60
61	PROCEDURE bisection;
62	EXTERNAL do_bisection(x[1..ninputs],x_par[1..ninputs],y[1..ninputs]);
63	END bisection;
64
65	PROCEDURE all;
66	RUN set_values_external;
67	RUN bisection;
68	END all;
69
70	END test_bisect;
71