int
main(int argc,char **argv)
{
anthy_context_t ac;
FILE *fp;
struct input cur_input;
struct res_db *db;
struct condition cond;
cur_input.serial = 0;
cur_input.str = 0;
init_condition(&cond);
parse_args(&cond, argc, argv);
db = create_db();
read_db(db, expdata);
printf("./test_anthy --help to print usage.\n");
print_run_env();
fp = fopen(testdata, "r");
if (!fp) {
printf("failed to open %s.\n", testdata);
return 0;
}
ac = init_lib(cond.use_utf8);
/* ファイルを読んでいくループ */
while (!read_file(fp, &cur_input)) {
if (check_cond(&cond, &cur_input)) {
set_string(&cond, db, &cur_input, ac);
}
}
anthy_release_context(ac);
anthy_quit();
if (cond.ask) {
/* ユーザに聞く */
ask_results(db);
}
show_stat(db);
save_db(expdata, db);
return 0;
}
int main(int argc, char* argv[])
{
// Load the mesh.
MeshSharedPtr mesh(new Mesh);
MeshReaderH2D mloader;
mloader.load("square.mesh", mesh);
// Perform initial mesh refinements.
for(int i = 0; i < INIT_GLOB_REF_NUM; i++)
mesh->refine_all_elements();
mesh->refine_towards_boundary("Bdy", INIT_BDY_REF_NUM);
// Initialize boundary conditions.
CustomEssentialBCNonConst bc_essential("Bdy");
EssentialBCs<double> bcs(&bc_essential);
// Create an H1 space with default shapeset.
SpaceSharedPtr<double> space(new H1Space<double>(mesh, &bcs, P_INIT));
int ndof = space->get_num_dofs();
// Initialize previous iteration solution for the Picard's method.
MeshFunctionSharedPtr<double> init_condition(new ConstantSolution<double>(mesh, INIT_COND_CONST));
// Initialize the weak formulation.
CustomNonlinearity lambda(alpha);
Hermes2DFunction<double> src(-heat_src);
CustomWeakFormPicard wf(init_condition, &lambda, &src);
// Initialize the Picard solver.
PicardSolver<double> picard(&wf, space);
picard.set_max_allowed_iterations(PICARD_MAX_ITER);
logger.info("Default tolerance");
try
{
picard.solve(init_condition);
}
catch(std::exception& e)
{
std::cout << e.what();
}
int iter_1 = picard.get_num_iters();
logger.info("Adjusted tolerance without Anderson");
// Perform the Picard's iteration (Anderson acceleration on by default).
picard.clear_tolerances();
picard.set_tolerance(PICARD_TOL, Hermes::Solvers::SolutionChangeAbsolute);
try
{
picard.solve(init_condition);
}
catch(std::exception& e)
{
std::cout << e.what();
}
int iter_2 = picard.get_num_iters();
// Translate the coefficient vector into a Solution.
MeshFunctionSharedPtr<double> sln(new Solution<double>);
Solution<double>::vector_to_solution(picard.get_sln_vector(), space, sln);
logger.info("Default tolerance without Anderson and good initial guess");
PicardSolver<double> picard2(&wf, space);
picard2.set_tolerance(1e-3, Hermes::Solvers::SolutionChangeRelative);
picard2.set_max_allowed_iterations(PICARD_MAX_ITER);
try
{
picard2.solve(sln);
}
catch(std::exception& e)
{
std::cout << e.what();
}
int iter_3 = picard2.get_num_iters();
logger.info("Default tolerance with Anderson and good initial guess");
picard2.use_Anderson_acceleration(true);
picard2.set_num_last_vector_used(PICARD_NUM_LAST_ITER_USED);
picard2.set_anderson_beta(PICARD_ANDERSON_BETA);
try
{
picard2.solve(sln);
}
catch(std::exception& e)
{
std::cout << e.what();
}
int iter_4 = picard2.get_num_iters();
#ifdef SHOW_OUTPUT
// Visualise the solution and mesh.
ScalarView s_view("Solution", new WinGeom(0, 0, 440, 350));
s_view.show_mesh(false);
s_view.show(sln);
OrderView o_view("Mesh", new WinGeom(450, 0, 420, 350));
o_view.show(space);
// Wait for all views to be closed.
View::wait();
#endif
bool success = Testing::test_value(iter_1, 14, "# of iterations 1", 1); // Tested value as of September 2013.
success = Testing::test_value(iter_2, 12, "# of iterations 2", 1) & success; // Tested value as of September 2013.
success = Testing::test_value(iter_3, 3, "# of iterations 3", 1) & success; // Tested value as of September 2013.
success = Testing::test_value(iter_4, 3, "# of iterations 4", 1) & success; // Tested value as of September 2013.
if(success)
{
printf("Success!\n");
return 0;
}
else
{
printf("Failure!\n");
return -1;
}
}
