Revision 760 branches/threaded2/src/nh99/mod_hess.cpp
mod_hess.cpp (revision 760)  

66  66 
} 
67  67 
#endif 
68  68 
} 
69  
70  69 
void function_minimizer::hess_routine_noparallel(void) 
71  70 
{ 
72 
int nvar=initial_params::nvarcalc();


73 
independent_variables x(1,nvar); 

71 
int nvar = initial_params::nvarcalc();


72 
independent_variables x(1, nvar);


74  73 
initial_params::xinit(x); 
75  74 
double delta=1.e5; 
76 
dvector g1(1,nvar); 

77 
dvector g2(1,nvar); 

78 
dvector hess(1,nvar); 

79 
dvector hess1(1,nvar); 

80 
dvector hess2(1,nvar); 

81 
double eps=.1;


75 
dvector g1(1, nvar);


76 
dvector g2(1, nvar);


77 
dvector hess(1, nvar);


78 
dvector hess1(1, nvar);


79 
dvector hess2(1, nvar);


80 
double eps = 0.1;


82  81 
gradient_structure::set_YES_DERIVATIVES(); 
83  82  
84 
dvector gbest(1,nvar); 

85 
gbest.fill_seqadd(1.e+50,0.); 

86 


87  83 
std::ostringstream oss; 
88  84 
std::thread::id this_thread_id = std::this_thread::get_id(); 
89  85 
oss << *ad_comm::adprogram_name << this_thread_id << ".hes"; 
...  ...  
161  157 
} 
162  158 
} 
163  159 
ofs << gradient_structure::Hybrid_bounded_flag; 
164 
dvector tscale(1,nvar); // need to get scale from somewhere


160 
dvector tscale(1,nvar); 

165  161 
initial_params::stddev_scale(tscale,x); 
166  162 
ofs << tscale; 
167  163 
Also available in: Unified diff