Revision 1053 branches/threaded2/examples/threaded/multisimple/msimple1.tpl
msimple1.tpl (revision 1053)  

1  1 
//$Id: msimple1.tpl 742M 20130602 00:06:06Z (local) $ 
2 
GLOBALS_SECTION 

3 
int NSLAVES = 7; 

4  2  
5  3 
DATA_SECTION 
6  4 
!! ad_comm::change_datafile_name("msimple.dat"); 
7  5 
init_int nobs 
8 
int nslaves 

9 
!! nslaves = NSLAVES; 

10 
matrix Y(1,nslaves,1,nobs) 

11 
matrix x(1,nslaves,1,nobs) 

6 
init_int nrow 

7 
matrix Y(1,nrow,1,nobs) 

8 
matrix x(1,nrow,1,nobs) 

9 
number A 

10 
number B 

12  11  
13  12 
LOCAL_CALCS 
14 
double A = 2.0;


15 
double B = 4.0;


13 
A = 2.0; 

14 
B = 4.0; 

16  15 
random_number_generator rng(101); 
17  16 
dvector err(1,nobs); 
18 
for (int i = 1; i <= nslaves; i++)


17 
for (int i = 1; i <= nrow; i++)


19  18 
{ 
20  19 
x(i).fill_randu(rng); 
21  20 
x(i) *= 100.0; 
...  ...  
27  26 
PARAMETER_SECTION 
28  27 
init_number a 
29  28 
init_number b 
30 
vector ff(1,nslaves)


31 
matrix pred_Y(1,nslaves,1,nobs)


29 
vector ff(1,nrow)


30 
matrix pred_Y(1,nrow,1,nobs)


32  31 
objective_function_value f 
33  32  
34  33 
PRELIMINARY_CALCS_SECTION 
...  ...  
36  35 
b = 2.0; 
37  36  
38  37 
PROCEDURE_SECTION 
39 
for (int i=1;i<=nslaves;i++)


38 
for (int i=1;i<=nrow;i++)


40  39 
{ 
41  40 
pred_Y(i) = a * x(i) + b; 
42  41 
ff(i) = norm2(pred_Y(i)Y(i)); 
43 
ff(i) = nobs/2.*log(ff(i)); // make it a likelihood function so that 

44 
// covariance matrix is correct 

45  42 
} 
46  43 
f = sum(ff); 
44 
f = (nrow*nobs)/2.*log(f); // make it a likelihood function so that 

45 
// covariance matrix is correct 

46 
REPORT_SECTION 

47 
report << "A = " << A << "; B = " << B <<endl; 

48 
report << "a = " << a << "; b = " << b <<endl; 

47  49  
50 
Also available in: Unified diff