Function initialize_solution

src/solution.cxx:105–148

}

void initialize_solution(Sol& solution, Prob& problem, Alg& algorithm, Workspace* workspace)
{
   int nphases = problem.nphases;
   int nparam;
   int i;

   solution.states      = new MatrixXd[nphases];           
   solution.controls    = new MatrixXd[nphases];           
   solution.nodes       = new MatrixXd[nphases];           
   solution.integrand_cost= new MatrixXd[nphases];         
   solution.parameters  = new MatrixXd[nphases];           
   solution.relative_errors     = new MatrixXd[nphases];   

   solution.dual.costates = new MatrixXd[nphases];         
   solution.dual.path     = new MatrixXd[nphases];         
   solution.dual.events   = new MatrixXd[nphases];         
   solution.dual.Hamiltonian = new MatrixXd[nphases];      
   solution.dual.linkages    = new MatrixXd;               
   solution.endpoint_cost    = new double[nphases];
   solution.integrated_cost  = new double[nphases];
   solution.problem     = &problem;


   for (i=0;i<nphases; i++)
   {
      nparam = problem.phase[i].nparameters;
      solution.parameters[i].resize(nparam,1);
   }

   solution.error_flag = 0;
   solution.error_msg = "";

   solution.mesh_stats = new MeshStats[ get_number_of_mesh_refinement_iterations(problem,algorithm)];
   for (i=0;i<get_number_of_mesh_refinement_iterations(problem,algorithm) ; i++)
   {
      solution.mesh_stats[i].n_con_evals = 0;
      solution.mesh_stats[i].n_obj_evals = 0;
      solution.mesh_stats[i].n_ode_rhs_evals = 0;
      solution.mesh_stats[i].n_jacobian_evals = 0;
      solution.mesh_stats[i].n_hessian_evals = 0;
   }

   return;
}

void resize_solution(Sol& solution, Prob& problem, Alg& algorithm)
{