Revision 752 branches/threaded2/tests/threaded_hessian_routine/function_minimizerhess_routine_threaded.cpp
function_minimizerhess_routine_threaded.cpp (revision 752)  

19  19  
20  20 
void function_minimizer::hess_routine_threaded(void) 
21  21 
{ 
22  
22  23 
int nvar=initial_params::nvarcalc(); // get the number of active parameters 
23  24 
//if (adjm_ptr) set_labels_for_hess(nvar); 
24  25 
independent_variables x(1,nvar); 
...  ...  
35  36 
gradient_structure::set_YES_DERIVATIVES(); 
36  37 
gbest.fill_seqadd(1.e+50,0.); 
37  38  
39 
/* 

38  40 
adstring tmpstring="admodel.hesthreaded"; 
39  41 
if (ad_comm::wd_flag) 
40  42 
tmpstring = ad_comm::adprogram_name + ".hesthreaded"; 
41  43 
uostream ofs((char*)tmpstring); 
42  44  
43  45 
ofs << nvar; 
46 
*/ 

47 
dvariable vf; 

48 
if (true) return; 

44  49 
{ 
45  50 
{ 
46  51 
dvariable vf=0.0; 
52 
if (true) return; 

47  53 
vf=initial_params::reset(dvar_vector(x)); 
48  54 
*objective_function_value::pobjfun=0.0; 
49  55 
pre_userfunction(); 
...  ...  
118  124 
hess2=(g1g2)/(sdelta1sdelta2); 
119  125 
hess=(eps2*hess1hess2) /(eps21.); 
120  126  
127 
/* 

121  128 
ofs << hess; 
122  129 
//if (adjm_ptr) ad_update_hess_stats_report(nvar,i); 
130 
*/ 

123  131 
} 
124  132 
} 
133 
/* 

125  134 
ofs << gradient_structure::Hybrid_bounded_flag; 
126  135 
dvector tscale(1,nvar); // need to get scale from somewhere 
127 
/*int check=*/initial_params::stddev_scale(tscale,x);


136 
initial_params::stddev_scale(tscale,x); 

128  137 
ofs << tscale; 
138 
*/ 

129  139 
} 
Also available in: Unified diff