void verlet(const char* filename){
dictionary* ini;
ini = iniparser_load(filename);
iniparser_dump(ini, stdout);
FILE *file;
file = fopen(LOGFILE, "a");
iniparser_dump(ini, file);
fclose(file);
ReportMessage("\n");
t_key key = key_init(ini);
t_pSys pSys = pSys_init (ini);
t_opts opts = opts_init (ini);
t_pair p = t_pair_init(ini);
t_dump dump = dump_init(ini);
t_fix fix = fix_init(ini);
t_compute compute = compute_init(ini, &key);
init (&key, &pSys, &opts, &p, &compute);
if(key.dump) dump_run(&key, &pSys, &opts, &dump); // Make initial snapshot.
while(opts.thisIter < opts.targIters){ // Repeat until iterational limit.
verlet_iter (&key, &pSys, &opts, &p, &fix, &compute);
if(key.dump) dump_run(&key, &pSys, &opts, &dump);
resetStep(&pSys);
}
}
void ImplicitStepper::prepareForExternalSolve()
{
if( !m_usedNonlinearSolver )
{ // calls to solve everything if hasnt been done already
resetStep();
solveUnconstrained( m_params.m_alwaysUseNonLinear );
}
}
void DeltaProcess::returnToDebugger() {
resetStepping(); // reset dispatch table
resetStep(); // disable stepping
suspend(stopped);// stop!
}
