void test_high_order_interpolation() {
bcon bc_child_doc[] = { "k10", "v10", "k11", "v11", BEND };
bcon bc_parent_doc[] = { "k0", "v0", "k1", BPD(bc_child_doc), "k2", "v2", BEND };
test_bson_from_bcon( bc_parent_doc, BCON_OK, BSON_VALID );
bcon bc_child_array[] = { "k10", "v10", "k11", "v11", BEND };
bcon bc_parent_doc_array[] = { "k0", "v0", "k1", BPA(bc_child_array), "k2", "v2", BEND };
test_bson_from_bcon( bc_parent_doc_array, BCON_OK, BSON_VALID );
}
int main () {
int N=10; // number of hidden layers
int N_dat = 102; // number of training data
int c=0;
float *w_ji;
float *w_kj;
float *y_pj;
//float In [34] ;
//float d_pk [34] ;
float y_pk;
float b=0.0;
y_pj= (float*)malloc(N * sizeof(float));
w_ji = (float*)malloc(N* sizeof(float));;//memory allocation for hidden w
w_kj = (float*)malloc(N* sizeof(float));;//memory allocation for output w
float In [192] = {11.800009,16.400026,17.400030,16.400026,12.600012,16.200026,1.800000,13.600016,14.200018,15.000021,30.400080,17.200029,13.400015,17.800032,15.200022,13.600016,
13.200014,13.600016,14.000017,13.600016,14.200018,17.000029,8.999998,11.400007,13.600016,13.600016,29.000074,13.600016,32.800072,17.800032,13.600016,11.400007,11.600008,8.999998,
11.800009,16.400026,13.600016,8.999998,13.600016,11.200006,21.200045,14.000017,20.000040,14.000017,2.000000,32.000084,32.000084,32.200081,27.600069,26.000063,1.800000,1.800000,17.000029,
1.800000,25.400061,12.200010,1.800000,14.200018,16.400026,11.800009,13.600016,13.600016,13.400015,1.800000,11.800009,44.799889,13.600016,16.600027,11.200006,1.800000,14.200018,24.200056,
32.200081,1.800000,32.200081,16.200026,13.800016,33.200066,13.600016,16.400026,35.000038,14.200018,13.600016,16.400026,1.800000,11.800009,12.200010,12.200010,15.000021,1.800000,17.000029,
16.800028,13.600016,1.800000,1.800000,17.200029,32.200081,13.800016,14.200018,13.600016,14.200018,16.400026,17.200029,32.400078,12.200010,1.800000,1.800000,11.600008,13.600016,13.600016,11.000006,
13.600016,32.600075,32.400078,1.800000,14.200018,2.000000,11.600008,1.800000,13.600016,32.000084,17.200029,13.600016,39.399971,14.000017,2.000000,33.200066,15.000021,14.000017,11.600008,32.200081,
13.400015,1.800000,33.600060,19.600039,1.800000,14.200018,17.400030,17.200029,13.600016,17.200029,17.200029,14.200018,13.600016,14.400019,16.400026,17.200029,16.600027,16.400026,13.800016,13.600016,
6.999996,17.200029,16.400026,32.200081,17.200029,16.400026,14.000017,2.000000,16.400026,17.600031,32.600075,13.800016,17.600031,20.600042,21.000044,2.000000,14.000017,15.000021,13.600016,17.400030,
14.000017,20.600042,15.000021,11.000006,9.199999,13.600016,12.200010,14.200018,14.200018,13.600016,14.200018,32.200081,23.800055,13.400015,32.200081,17.200029,14.200018,13.600016,17.600031,17.200029,20.800043};
float d_pk [192] ;
for (c;c<192;c++) {
if (In [c] < 11.6) d_pk[c]=-1;
else d_pk[c]=1;
}
//In = (float*)malloc(N_dat * sizeof(float));
//d_pk = (float*)malloc(N_dat * sizeof(float));//memory allocation for desired values
float eta = 0.000000002; // 0.0000272
//input and desired output generation
/*int c=1;
for (c;c<N_dat;c++) {
In[c]=(c);
d_pk[c]=1.8 + 2.0* sqrt(c)/(10.0) + log(c);
}
* */
//Training
//BPA (N, N_dat,eta, In,d_pk, float *w_ji, float *w_kj,b)
BPA (N,N_dat,eta, In,d_pk, w_ji, w_kj,b);
//testing
int c2=0;
for (c2;c2<192;c2++) {
printf("%f \t", In [c2]);
//In[c2]= 2.0+c2;
ANN (N,In[c2], w_ji, w_kj,b, y_pj,y_pk);
printf ("%f \t", y_pk);
//printf ("%f \n", 1.8 +2.0*sqrt(In[c2])/10.0 + log(In[c2]));
printf ("%f \n", d_pk[c2]);
}
c=0;
for (c;c<N;c++) {
printf ("%f \t", w_ji[c]);
printf ("%f \n", w_kj[c]);
}
printf("b= %f \n",b);
return 0;
}
