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.e-5;
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