cahnHilliard.dat.3d 3.45 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
dimension of world: 3

% =================== MESH ================================
mesh_name: mesh
mesh->macro file name:      ./macro/fliese_20x20x2.per.3d
mesh->periodic file:        ./macro/fliese_20x20x2.per
mesh->global refinements:   3
mesh->check:                0

mesh->scale mesh: 1
x: 10
y: ${x}
z: ${x} / 10.0

mesh->dimension: [${x}, ${y}, ${z}]

#include "init/reinit.inc.2d"

% ============== USER-PARAMETER ==========================
ch->use reinit:            0
ch->epsilon:               0.5
ch->gamma:                 1
ch->double-well type:      1

% =========== OUTPUT ==============================================
ch->space->output->filename: ./output/ch3d_

% ==================== TIMESTEPS ===============================
adapt->timestep:      0.1
adapt->max timestep:  1e+10
adapt->min timestep:  1e-10
adapt->start time:    0.0
adapt->end time:      10000.0

% ============= PROBLEM-SPACES ==================================
ch->space->components:           2
ch->space->polynomial degree[0]: 2
ch->space->polynomial degree[1]: 2
ch->space->dim:                  3
ch->space->mesh:                 mesh

% ================== SOLVER ======================================
ch->space->solver:             fgmres
ch->space->solver->use old initial guess: 1
ch->space->solver->backend:            mtl
%ch->space->solver->petsc prefix:       ch
ch->space->solver->right precon:    ch2
ch->space->solver->orthogonalization:  1
ch->space->solver->restart:            30
ch->space->solver->ell:            3
ch->space->solver->max iteration:      2000
ch->space->solver->print cycle:        10
ch->space->solver->tolerance:          1.e-8
ch->space->solver->info:               10


ch->space->solver->right precon->subsolver M->solver: cg
ch->space->solver->right precon->subsolver M->solver->backend: mtl
ch->space->solver->right precon->subsolver M->solver->tolerance: 1.e-3
ch->space->solver->right precon->subsolver M->solver->max iteration: 5
ch->space->solver->right precon->subsolver M->solver->left precon: diag

ch->space->solver->right precon->subsolver MpL->solver: cg
ch->space->solver->right precon->subsolver MpL->solver->backend: mtl
ch->space->solver->right precon->subsolver MpL->solver->tolerance: 1.e-3
ch->space->solver->right precon->subsolver MpL->solver->max iteration: 5
ch->space->solver->right precon->subsolver MpL->solver->left precon: diag

ch->space->solver->right precon->subsolver S->solver->max iteration: 10
% 
% ch->space->solver->right precon->subsolver M->solver: cg
% ch->space->solver->right precon->subsolver M->solver->backend: mtl
% ch->space->solver->right precon->subsolver M->solver->tolerance: 1.e-3
% ch->space->solver->right precon->subsolver M->solver->max iteration: 5
% ch->space->solver->right precon->subsolver M->solver->left precon: diag
% 
% ch->space->solver->right precon->subsolver MpL->solver: cg
% ch->space->solver->right precon->subsolver MpL->solver->backend: mtl
% ch->space->solver->right precon->subsolver MpL->solver->tolerance: 1.e-3
% ch->space->solver->right precon->subsolver MpL->solver->max iteration: 20
% ch->space->solver->right precon->subsolver MpL->solver->left precon: diag

% =================== OUTPUT =========================================
ch->space->output->ParaView animation:        1
ch->space->output->ParaView format:           1
ch->space->output->write every i-th timestep: 10
ch->space->output->append index:              1
ch->space->output->index length:              9
ch->space->output->index decimals:            7