1 static char help[] = "Poisson Problem in 2d and 3d with simplicial finite elements.\n\
2 We solve the Poisson problem in a rectangular\n\
3 domain, using a parallel unstructured mesh (DMPLEX) to discretize it.\n\
4 This example supports discretized auxiliary fields (conductivity) as well as\n\
5 multilevel nonlinear solvers.\n\n\n";
6
7 /*
8 A visualization of the adaptation can be accomplished using:
9
10 -dm_adapt_view hdf5:$PWD/adapt.h5 -sol_adapt_view hdf5:$PWD/adapt.h5::append -dm_adapt_pre_view hdf5:$PWD/orig.h5 -sol_adapt_pre_view hdf5:$PWD/orig.h5::append
11
12 Information on refinement:
13
14 -info :~sys,vec,is,mat,ksp,snes,ts
15 */
16
17 #include <petscdmplex.h>
18 #include <petscdmadaptor.h>
19 #include <petscsnes.h>
20 #include <petscds.h>
21 #include <petscviewerhdf5.h>
22
23 typedef enum {
24 NEUMANN,
25 DIRICHLET,
26 NONE
27 } BCType;
28 typedef enum {
29 RUN_FULL,
30 RUN_EXACT,
31 RUN_TEST,
32 RUN_PERF
33 } RunType;
34 typedef enum {
35 COEFF_NONE,
36 COEFF_ANALYTIC,
37 COEFF_FIELD,
38 COEFF_NONLINEAR,
39 COEFF_BALL,
40 COEFF_CROSS,
41 COEFF_CHECKERBOARD_0,
42 COEFF_CHECKERBOARD_1
43 } CoeffType;
44
45 typedef struct {
46 RunType runType; /* Whether to run tests, or solve the full problem */
47 PetscBool jacobianMF; /* Whether to calculate the Jacobian action on the fly */
48 PetscBool showInitial, showSolution, restart, quiet, nonzInit;
49 /* Problem definition */
50 BCType bcType;
51 CoeffType variableCoefficient;
52 PetscErrorCode (**exactFuncs)(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx);
53 PetscBool fieldBC;
54 void (**exactFields)(PetscInt, PetscInt, PetscInt, const PetscInt[], const PetscInt[], const PetscScalar[], const PetscScalar[], const PetscScalar[], const PetscInt[], const PetscInt[], const PetscScalar[], const PetscScalar[], const PetscScalar[], PetscReal, const PetscReal[], PetscInt, const PetscScalar[], PetscScalar[]);
55 PetscBool bdIntegral; /* Compute the integral of the solution on the boundary */
56 /* Reproducing tests from SISC 40(3), pp. A1473-A1493, 2018 */
57 PetscInt div; /* Number of divisions */
58 PetscInt k; /* Parameter for checkerboard coefficient */
59 PetscInt *kgrid; /* Random parameter grid */
60 PetscBool rand; /* Make random assignments */
61 /* Solver */
62 PC pcmg; /* This is needed for error monitoring */
63 PetscBool checkksp; /* Whether to check the KSPSolve for runType == RUN_TEST */
64 } AppCtx;
65
zero(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)66 static PetscErrorCode zero(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
67 {
68 u[0] = 0.0;
69 return PETSC_SUCCESS;
70 }
71
ecks(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)72 static PetscErrorCode ecks(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
73 {
74 u[0] = x[0];
75 return PETSC_SUCCESS;
76 }
77
78 /*
79 In 2D for Dirichlet conditions, we use exact solution:
80
81 u = x^2 + y^2
82 f = 4
83
84 so that
85
86 -\Delta u + f = -4 + 4 = 0
87
88 For Neumann conditions, we have
89
90 -\nabla u \cdot -\hat y |_{y=0} = (2y)|_{y=0} = 0 (bottom)
91 -\nabla u \cdot \hat y |_{y=1} = -(2y)|_{y=1} = -2 (top)
92 -\nabla u \cdot -\hat x |_{x=0} = (2x)|_{x=0} = 0 (left)
93 -\nabla u \cdot \hat x |_{x=1} = -(2x)|_{x=1} = -2 (right)
94
95 Which we can express as
96
97 \nabla u \cdot \hat n|_\Gamma = {2 x, 2 y} \cdot \hat n = 2 (x + y)
98
99 The boundary integral of this solution is (assuming we are not orienting the edges)
100
101 \int^1_0 x^2 dx + \int^1_0 (1 + y^2) dy + \int^1_0 (x^2 + 1) dx + \int^1_0 y^2 dy = 1/3 + 4/3 + 4/3 + 1/3 = 3 1/3
102 */
quadratic_u_2d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)103 static PetscErrorCode quadratic_u_2d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
104 {
105 *u = x[0] * x[0] + x[1] * x[1];
106 return PETSC_SUCCESS;
107 }
108
quadratic_u_field_2d(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar uexact[])109 static void quadratic_u_field_2d(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar uexact[])
110 {
111 uexact[0] = a[0];
112 }
113
ball_u_2d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)114 static PetscErrorCode ball_u_2d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
115 {
116 const PetscReal alpha = 500.;
117 const PetscReal radius2 = PetscSqr(0.15);
118 const PetscReal r2 = PetscSqr(x[0] - 0.5) + PetscSqr(x[1] - 0.5);
119 const PetscReal xi = alpha * (radius2 - r2);
120
121 *u = PetscTanhScalar(xi) + 1.0;
122 return PETSC_SUCCESS;
123 }
124
cross_u_2d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)125 static PetscErrorCode cross_u_2d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
126 {
127 const PetscReal alpha = 50 * 4;
128 const PetscReal xy = (x[0] - 0.5) * (x[1] - 0.5);
129
130 *u = PetscSinReal(alpha * xy) * (alpha * PetscAbsReal(xy) < 2 * PETSC_PI ? 1 : 0.01);
131 return PETSC_SUCCESS;
132 }
133
f0_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])134 static void f0_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
135 {
136 f0[0] = 4.0;
137 }
138
f0_ball_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])139 static void f0_ball_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
140 {
141 PetscInt d;
142 const PetscReal alpha = 500., radius2 = PetscSqr(0.15);
143 PetscReal r2, xi;
144
145 for (d = 0, r2 = 0.0; d < dim; ++d) r2 += PetscSqr(x[d] - 0.5);
146 xi = alpha * (radius2 - r2);
147 f0[0] = (-2.0 * dim * alpha - 8.0 * PetscSqr(alpha) * r2 * PetscTanhReal(xi)) * PetscSqr(1.0 / PetscCoshReal(xi));
148 }
149
f0_cross_u_2d(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])150 static void f0_cross_u_2d(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
151 {
152 const PetscReal alpha = 50 * 4;
153 const PetscReal xy = (x[0] - 0.5) * (x[1] - 0.5);
154
155 f0[0] = PetscSinReal(alpha * xy) * (alpha * PetscAbsReal(xy) < 2 * PETSC_PI ? 1 : 0.01);
156 }
157
f0_checkerboard_0_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])158 static void f0_checkerboard_0_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
159 {
160 f0[0] = -20.0 * PetscExpReal(-(PetscSqr(x[0] - 0.5) + PetscSqr(x[1] - 0.5)));
161 }
162
f0_bd_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],const PetscReal n[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])163 static void f0_bd_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
164 {
165 PetscInt d;
166 for (d = 0, f0[0] = 0.0; d < dim; ++d) f0[0] += -n[d] * 2.0 * x[d];
167 }
168
169 /* gradU[comp*dim+d] = {u_x, u_y} or {u_x, u_y, u_z} */
f1_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f1[])170 static void f1_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f1[])
171 {
172 PetscInt d;
173 for (d = 0; d < dim; ++d) f1[d] = u_x[d];
174 }
175
176 /* < \nabla v, \nabla u + {\nabla u}^T >
177 This just gives \nabla u, give the perdiagonal for the transpose */
g3_uu(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,PetscReal u_tShift,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar g3[])178 static void g3_uu(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, PetscReal u_tShift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])
179 {
180 PetscInt d;
181 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0;
182 }
183
184 /*
185 In 2D for x periodicity and y Dirichlet conditions, we use exact solution:
186
187 u = sin(2 pi x)
188 f = -4 pi^2 sin(2 pi x)
189
190 so that
191
192 -\Delta u + f = 4 pi^2 sin(2 pi x) - 4 pi^2 sin(2 pi x) = 0
193 */
xtrig_u_2d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)194 static PetscErrorCode xtrig_u_2d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
195 {
196 *u = PetscSinReal(2.0 * PETSC_PI * x[0]);
197 return PETSC_SUCCESS;
198 }
199
f0_xtrig_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])200 static void f0_xtrig_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
201 {
202 f0[0] = -4.0 * PetscSqr(PETSC_PI) * PetscSinReal(2.0 * PETSC_PI * x[0]);
203 }
204
205 /*
206 In 2D for x-y periodicity, we use exact solution:
207
208 u = sin(2 pi x) sin(2 pi y)
209 f = -8 pi^2 sin(2 pi x)
210
211 so that
212
213 -\Delta u + f = 4 pi^2 sin(2 pi x) sin(2 pi y) + 4 pi^2 sin(2 pi x) sin(2 pi y) - 8 pi^2 sin(2 pi x) = 0
214 */
xytrig_u_2d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)215 static PetscErrorCode xytrig_u_2d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
216 {
217 *u = PetscSinReal(2.0 * PETSC_PI * x[0]) * PetscSinReal(2.0 * PETSC_PI * x[1]);
218 return PETSC_SUCCESS;
219 }
220
f0_xytrig_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])221 static void f0_xytrig_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
222 {
223 f0[0] = -8.0 * PetscSqr(PETSC_PI) * PetscSinReal(2.0 * PETSC_PI * x[0]);
224 }
225
226 /*
227 In 2D for Dirichlet conditions with a variable coefficient, we use exact solution:
228
229 u = x^2 + y^2
230 f = 6 (x + y)
231 nu = (x + y)
232
233 so that
234
235 -\div \nu \grad u + f = -6 (x + y) + 6 (x + y) = 0
236 */
nu_2d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)237 static PetscErrorCode nu_2d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
238 {
239 *u = x[0] + x[1];
240 return PETSC_SUCCESS;
241 }
242
checkerboardCoeff(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)243 static PetscErrorCode checkerboardCoeff(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
244 {
245 AppCtx *user = (AppCtx *)ctx;
246 PetscInt div = user->div;
247 PetscInt k = user->k;
248 PetscInt mask = 0, ind = 0, d;
249
250 PetscFunctionBeginUser;
251 for (d = 0; d < dim; ++d) mask = (mask + (PetscInt)(x[d] * div)) % 2;
252 if (user->kgrid) {
253 for (d = 0; d < dim; ++d) {
254 if (d > 0) ind *= dim;
255 ind += (PetscInt)(x[d] * div);
256 }
257 k = user->kgrid[ind];
258 }
259 u[0] = mask ? 1.0 : PetscPowRealInt(10.0, -k);
260 PetscFunctionReturn(PETSC_SUCCESS);
261 }
262
f0_analytic_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])263 void f0_analytic_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
264 {
265 f0[0] = 6.0 * (x[0] + x[1]);
266 }
267
268 /* gradU[comp*dim+d] = {u_x, u_y} or {u_x, u_y, u_z} */
f1_analytic_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f1[])269 void f1_analytic_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f1[])
270 {
271 PetscInt d;
272 for (d = 0; d < dim; ++d) f1[d] = (x[0] + x[1]) * u_x[d];
273 }
274
f1_field_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f1[])275 void f1_field_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f1[])
276 {
277 PetscInt d;
278 for (d = 0; d < dim; ++d) f1[d] = a[0] * u_x[d];
279 }
280
281 /* < \nabla v, \nabla u + {\nabla u}^T >
282 This just gives \nabla u, give the perdiagonal for the transpose */
g3_analytic_uu(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,PetscReal u_tShift,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar g3[])283 void g3_analytic_uu(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, PetscReal u_tShift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])
284 {
285 PetscInt d;
286 for (d = 0; d < dim; ++d) g3[d * dim + d] = x[0] + x[1];
287 }
288
g3_field_uu(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,PetscReal u_tShift,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar g3[])289 void g3_field_uu(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, PetscReal u_tShift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])
290 {
291 PetscInt d;
292 for (d = 0; d < dim; ++d) g3[d * dim + d] = a[0];
293 }
294
295 /*
296 In 2D for Dirichlet conditions with a nonlinear coefficient (p-Laplacian with p = 4), we use exact solution:
297
298 u = x^2 + y^2
299 f = 16 (x^2 + y^2)
300 nu = 1/2 |grad u|^2
301
302 so that
303
304 -\div \nu \grad u + f = -16 (x^2 + y^2) + 16 (x^2 + y^2) = 0
305 */
f0_analytic_nonlinear_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])306 void f0_analytic_nonlinear_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
307 {
308 f0[0] = 16.0 * (x[0] * x[0] + x[1] * x[1]);
309 }
310
311 /* gradU[comp*dim+d] = {u_x, u_y} or {u_x, u_y, u_z} */
f1_analytic_nonlinear_u(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f1[])312 void f1_analytic_nonlinear_u(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f1[])
313 {
314 PetscScalar nu = 0.0;
315 PetscInt d;
316 for (d = 0; d < dim; ++d) nu += u_x[d] * u_x[d];
317 for (d = 0; d < dim; ++d) f1[d] = 0.5 * nu * u_x[d];
318 }
319
320 /*
321 grad (u + eps w) - grad u = eps grad w
322
323 1/2 |grad (u + eps w)|^2 grad (u + eps w) - 1/2 |grad u|^2 grad u
324 = 1/2 (|grad u|^2 + 2 eps <grad u,grad w>) (grad u + eps grad w) - 1/2 |grad u|^2 grad u
325 = 1/2 (eps |grad u|^2 grad w + 2 eps <grad u,grad w> grad u)
326 = eps (1/2 |grad u|^2 grad w + grad u <grad u,grad w>)
327 */
g3_analytic_nonlinear_uu(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,PetscReal u_tShift,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar g3[])328 void g3_analytic_nonlinear_uu(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, PetscReal u_tShift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])
329 {
330 PetscScalar nu = 0.0;
331 PetscInt d, e;
332 for (d = 0; d < dim; ++d) nu += u_x[d] * u_x[d];
333 for (d = 0; d < dim; ++d) {
334 g3[d * dim + d] = 0.5 * nu;
335 for (e = 0; e < dim; ++e) g3[d * dim + e] += u_x[d] * u_x[e];
336 }
337 }
338
339 /*
340 In 3D for Dirichlet conditions we use exact solution:
341
342 u = 2/3 (x^2 + y^2 + z^2)
343 f = 4
344
345 so that
346
347 -\Delta u + f = -2/3 * 6 + 4 = 0
348
349 For Neumann conditions, we have
350
351 -\nabla u \cdot -\hat z |_{z=0} = (2z)|_{z=0} = 0 (bottom)
352 -\nabla u \cdot \hat z |_{z=1} = -(2z)|_{z=1} = -2 (top)
353 -\nabla u \cdot -\hat y |_{y=0} = (2y)|_{y=0} = 0 (front)
354 -\nabla u \cdot \hat y |_{y=1} = -(2y)|_{y=1} = -2 (back)
355 -\nabla u \cdot -\hat x |_{x=0} = (2x)|_{x=0} = 0 (left)
356 -\nabla u \cdot \hat x |_{x=1} = -(2x)|_{x=1} = -2 (right)
357
358 Which we can express as
359
360 \nabla u \cdot \hat n|_\Gamma = {2 x, 2 y, 2z} \cdot \hat n = 2 (x + y + z)
361 */
quadratic_u_3d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)362 static PetscErrorCode quadratic_u_3d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
363 {
364 *u = 2.0 * (x[0] * x[0] + x[1] * x[1] + x[2] * x[2]) / 3.0;
365 return PETSC_SUCCESS;
366 }
367
ball_u_3d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)368 static PetscErrorCode ball_u_3d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
369 {
370 const PetscReal alpha = 500.;
371 const PetscReal radius2 = PetscSqr(0.15);
372 const PetscReal r2 = PetscSqr(x[0] - 0.5) + PetscSqr(x[1] - 0.5) + PetscSqr(x[2] - 0.5);
373 const PetscReal xi = alpha * (radius2 - r2);
374
375 *u = PetscTanhScalar(xi) + 1.0;
376 return PETSC_SUCCESS;
377 }
378
quadratic_u_field_3d(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar uexact[])379 static void quadratic_u_field_3d(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar uexact[])
380 {
381 uexact[0] = a[0];
382 }
383
cross_u_3d(PetscInt dim,PetscReal time,const PetscReal x[],PetscInt Nc,PetscScalar * u,PetscCtx ctx)384 static PetscErrorCode cross_u_3d(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
385 {
386 const PetscReal alpha = 50 * 4;
387 const PetscReal xyz = (x[0] - 0.5) * (x[1] - 0.5) * (x[2] - 0.5);
388
389 *u = PetscSinReal(alpha * xyz) * (alpha * PetscAbsReal(xyz) < 2 * PETSC_PI ? (alpha * PetscAbsReal(xyz) > -2 * PETSC_PI ? 1.0 : 0.01) : 0.01);
390 return PETSC_SUCCESS;
391 }
392
f0_cross_u_3d(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],PetscInt numConstants,const PetscScalar constants[],PetscScalar f0[])393 static void f0_cross_u_3d(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
394 {
395 const PetscReal alpha = 50 * 4;
396 const PetscReal xyz = (x[0] - 0.5) * (x[1] - 0.5) * (x[2] - 0.5);
397
398 f0[0] = PetscSinReal(alpha * xyz) * (alpha * PetscAbsReal(xyz) < 2 * PETSC_PI ? (alpha * PetscAbsReal(xyz) > -2 * PETSC_PI ? 1.0 : 0.01) : 0.01);
399 }
400
bd_integral_2d(PetscInt dim,PetscInt Nf,PetscInt NfAux,const PetscInt uOff[],const PetscInt uOff_x[],const PetscScalar u[],const PetscScalar u_t[],const PetscScalar u_x[],const PetscInt aOff[],const PetscInt aOff_x[],const PetscScalar a[],const PetscScalar a_t[],const PetscScalar a_x[],PetscReal t,const PetscReal x[],const PetscReal n[],PetscInt numConstants,const PetscScalar constants[],PetscScalar * uint)401 static void bd_integral_2d(PetscInt dim, PetscInt Nf, PetscInt NfAux, const PetscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar *uint)
402 {
403 uint[0] = u[0];
404 }
405
ProcessOptions(MPI_Comm comm,AppCtx * options)406 static PetscErrorCode ProcessOptions(MPI_Comm comm, AppCtx *options)
407 {
408 const char *bcTypes[3] = {"neumann", "dirichlet", "none"};
409 const char *runTypes[4] = {"full", "exact", "test", "perf"};
410 const char *coeffTypes[8] = {"none", "analytic", "field", "nonlinear", "ball", "cross", "checkerboard_0", "checkerboard_1"};
411 PetscInt bc, run, coeff;
412
413 PetscFunctionBeginUser;
414 options->runType = RUN_FULL;
415 options->bcType = DIRICHLET;
416 options->variableCoefficient = COEFF_NONE;
417 options->fieldBC = PETSC_FALSE;
418 options->jacobianMF = PETSC_FALSE;
419 options->showInitial = PETSC_FALSE;
420 options->showSolution = PETSC_FALSE;
421 options->restart = PETSC_FALSE;
422 options->quiet = PETSC_FALSE;
423 options->nonzInit = PETSC_FALSE;
424 options->bdIntegral = PETSC_FALSE;
425 options->checkksp = PETSC_FALSE;
426 options->div = 4;
427 options->k = 1;
428 options->kgrid = NULL;
429 options->rand = PETSC_FALSE;
430
431 PetscOptionsBegin(comm, "", "Poisson Problem Options", "DMPLEX");
432 run = options->runType;
433 PetscCall(PetscOptionsEList("-run_type", "The run type", "ex12.c", runTypes, 4, runTypes[options->runType], &run, NULL));
434 options->runType = (RunType)run;
435 bc = options->bcType;
436 PetscCall(PetscOptionsEList("-bc_type", "Type of boundary condition", "ex12.c", bcTypes, 3, bcTypes[options->bcType], &bc, NULL));
437 options->bcType = (BCType)bc;
438 coeff = options->variableCoefficient;
439 PetscCall(PetscOptionsEList("-variable_coefficient", "Type of variable coefficient", "ex12.c", coeffTypes, 8, coeffTypes[options->variableCoefficient], &coeff, NULL));
440 options->variableCoefficient = (CoeffType)coeff;
441
442 PetscCall(PetscOptionsBool("-field_bc", "Use a field representation for the BC", "ex12.c", options->fieldBC, &options->fieldBC, NULL));
443 PetscCall(PetscOptionsBool("-jacobian_mf", "Calculate the action of the Jacobian on the fly", "ex12.c", options->jacobianMF, &options->jacobianMF, NULL));
444 PetscCall(PetscOptionsBool("-show_initial", "Output the initial guess for verification", "ex12.c", options->showInitial, &options->showInitial, NULL));
445 PetscCall(PetscOptionsBool("-show_solution", "Output the solution for verification", "ex12.c", options->showSolution, &options->showSolution, NULL));
446 PetscCall(PetscOptionsBool("-restart", "Read in the mesh and solution from a file", "ex12.c", options->restart, &options->restart, NULL));
447 PetscCall(PetscOptionsBool("-quiet", "Don't print any vecs", "ex12.c", options->quiet, &options->quiet, NULL));
448 PetscCall(PetscOptionsBool("-nonzero_initial_guess", "nonzero initial guess", "ex12.c", options->nonzInit, &options->nonzInit, NULL));
449 PetscCall(PetscOptionsBool("-bd_integral", "Compute the integral of the solution on the boundary", "ex12.c", options->bdIntegral, &options->bdIntegral, NULL));
450 if (options->runType == RUN_TEST) PetscCall(PetscOptionsBool("-run_test_check_ksp", "Check solution of KSP", "ex12.c", options->checkksp, &options->checkksp, NULL));
451 PetscCall(PetscOptionsInt("-div", "The number of division for the checkerboard coefficient", "ex12.c", options->div, &options->div, NULL));
452 PetscCall(PetscOptionsInt("-k", "The exponent for the checkerboard coefficient", "ex12.c", options->k, &options->k, NULL));
453 PetscCall(PetscOptionsBool("-k_random", "Assign random k values to checkerboard", "ex12.c", options->rand, &options->rand, NULL));
454 PetscOptionsEnd();
455 PetscFunctionReturn(PETSC_SUCCESS);
456 }
457
CreateBCLabel(DM dm,const char name[])458 static PetscErrorCode CreateBCLabel(DM dm, const char name[])
459 {
460 DM plex;
461 DMLabel label;
462
463 PetscFunctionBeginUser;
464 PetscCall(DMCreateLabel(dm, name));
465 PetscCall(DMGetLabel(dm, name, &label));
466 PetscCall(DMConvert(dm, DMPLEX, &plex));
467 PetscCall(DMPlexMarkBoundaryFaces(plex, 1, label));
468 PetscCall(DMDestroy(&plex));
469 PetscFunctionReturn(PETSC_SUCCESS);
470 }
471
CreateMesh(MPI_Comm comm,AppCtx * user,DM * dm)472 static PetscErrorCode CreateMesh(MPI_Comm comm, AppCtx *user, DM *dm)
473 {
474 PetscFunctionBeginUser;
475 PetscCall(DMCreate(comm, dm));
476 PetscCall(DMSetType(*dm, DMPLEX));
477 PetscCall(DMSetFromOptions(*dm));
478 {
479 char convType[256];
480 PetscBool flg;
481
482 PetscOptionsBegin(comm, "", "Mesh conversion options", "DMPLEX");
483 PetscCall(PetscOptionsFList("-dm_plex_convert_type", "Convert DMPlex to another format", "ex12", DMList, DMPLEX, convType, 256, &flg));
484 PetscOptionsEnd();
485 if (flg) {
486 DM dmConv;
487
488 PetscCall(DMConvert(*dm, convType, &dmConv));
489 if (dmConv) {
490 PetscCall(DMDestroy(dm));
491 *dm = dmConv;
492 }
493 PetscCall(DMSetFromOptions(*dm));
494 PetscCall(DMSetUp(*dm));
495 }
496 }
497 PetscCall(DMViewFromOptions(*dm, NULL, "-dm_view"));
498 if (user->rand) {
499 PetscRandom r;
500 PetscReal val;
501 PetscInt dim, N, i;
502
503 PetscCall(DMGetDimension(*dm, &dim));
504 N = PetscPowInt(user->div, dim);
505 PetscCall(PetscMalloc1(N, &user->kgrid));
506 PetscCall(PetscRandomCreate(PETSC_COMM_SELF, &r));
507 PetscCall(PetscRandomSetFromOptions(r));
508 PetscCall(PetscRandomSetInterval(r, 0.0, user->k));
509 PetscCall(PetscRandomSetSeed(r, 1973));
510 PetscCall(PetscRandomSeed(r));
511 for (i = 0; i < N; ++i) {
512 PetscCall(PetscRandomGetValueReal(r, &val));
513 user->kgrid[i] = 1 + (PetscInt)val;
514 }
515 PetscCall(PetscRandomDestroy(&r));
516 }
517 PetscFunctionReturn(PETSC_SUCCESS);
518 }
519
SetupProblem(DM dm,AppCtx * user)520 static PetscErrorCode SetupProblem(DM dm, AppCtx *user)
521 {
522 PetscDS ds;
523 DMLabel label;
524 PetscWeakForm wf;
525 const PetscReal *L;
526 const PetscInt id = 1;
527 PetscInt bd, dim;
528
529 PetscFunctionBeginUser;
530 PetscCall(DMGetDS(dm, &ds));
531 PetscCall(DMGetDimension(dm, &dim));
532 PetscCall(DMGetPeriodicity(dm, NULL, NULL, &L));
533 switch (user->variableCoefficient) {
534 case COEFF_NONE:
535 if (L && L[0]) {
536 if (L && L[1]) {
537 PetscCall(PetscDSSetResidual(ds, 0, f0_xytrig_u, f1_u));
538 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_uu));
539 } else {
540 PetscCall(PetscDSSetResidual(ds, 0, f0_xtrig_u, f1_u));
541 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_uu));
542 }
543 } else {
544 PetscCall(PetscDSSetResidual(ds, 0, f0_u, f1_u));
545 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_uu));
546 }
547 break;
548 case COEFF_ANALYTIC:
549 PetscCall(PetscDSSetResidual(ds, 0, f0_analytic_u, f1_analytic_u));
550 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_analytic_uu));
551 break;
552 case COEFF_FIELD:
553 PetscCall(PetscDSSetResidual(ds, 0, f0_analytic_u, f1_field_u));
554 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_field_uu));
555 break;
556 case COEFF_NONLINEAR:
557 PetscCall(PetscDSSetResidual(ds, 0, f0_analytic_nonlinear_u, f1_analytic_nonlinear_u));
558 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_analytic_nonlinear_uu));
559 break;
560 case COEFF_BALL:
561 PetscCall(PetscDSSetResidual(ds, 0, f0_ball_u, f1_u));
562 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_uu));
563 break;
564 case COEFF_CROSS:
565 switch (dim) {
566 case 2:
567 PetscCall(PetscDSSetResidual(ds, 0, f0_cross_u_2d, f1_u));
568 break;
569 case 3:
570 PetscCall(PetscDSSetResidual(ds, 0, f0_cross_u_3d, f1_u));
571 break;
572 default:
573 SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Invalid dimension %" PetscInt_FMT, dim);
574 }
575 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_uu));
576 break;
577 case COEFF_CHECKERBOARD_0:
578 PetscCall(PetscDSSetResidual(ds, 0, f0_checkerboard_0_u, f1_field_u));
579 PetscCall(PetscDSSetJacobian(ds, 0, 0, NULL, NULL, NULL, g3_field_uu));
580 break;
581 default:
582 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Invalid variable coefficient type %d", user->variableCoefficient);
583 }
584 switch (dim) {
585 case 2:
586 switch (user->variableCoefficient) {
587 case COEFF_BALL:
588 user->exactFuncs[0] = ball_u_2d;
589 break;
590 case COEFF_CROSS:
591 user->exactFuncs[0] = cross_u_2d;
592 break;
593 case COEFF_CHECKERBOARD_0:
594 user->exactFuncs[0] = zero;
595 break;
596 default:
597 if (L && L[0]) {
598 if (L && L[1]) {
599 user->exactFuncs[0] = xytrig_u_2d;
600 } else {
601 user->exactFuncs[0] = xtrig_u_2d;
602 }
603 } else {
604 user->exactFuncs[0] = quadratic_u_2d;
605 user->exactFields[0] = quadratic_u_field_2d;
606 }
607 }
608 if (user->bcType == NEUMANN) {
609 PetscCall(DMGetLabel(dm, "boundary", &label));
610 PetscCall(DMAddBoundary(dm, DM_BC_NATURAL, "wall", label, 1, &id, 0, 0, NULL, NULL, NULL, user, &bd));
611 PetscCall(PetscDSGetBoundary(ds, bd, &wf, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL));
612 PetscCall(PetscWeakFormSetIndexBdResidual(wf, label, id, 0, 0, 0, f0_bd_u, 0, NULL));
613 }
614 break;
615 case 3:
616 switch (user->variableCoefficient) {
617 case COEFF_BALL:
618 user->exactFuncs[0] = ball_u_3d;
619 break;
620 case COEFF_CROSS:
621 user->exactFuncs[0] = cross_u_3d;
622 break;
623 default:
624 user->exactFuncs[0] = quadratic_u_3d;
625 user->exactFields[0] = quadratic_u_field_3d;
626 }
627 if (user->bcType == NEUMANN) {
628 PetscCall(DMGetLabel(dm, "boundary", &label));
629 PetscCall(DMAddBoundary(dm, DM_BC_NATURAL, "wall", label, 1, &id, 0, 0, NULL, NULL, NULL, user, &bd));
630 PetscCall(PetscDSGetBoundary(ds, bd, &wf, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL));
631 PetscCall(PetscWeakFormSetIndexBdResidual(wf, label, id, 0, 0, 0, f0_bd_u, 0, NULL));
632 }
633 break;
634 default:
635 SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Invalid dimension %" PetscInt_FMT, dim);
636 }
637 /* Setup constants */
638 switch (user->variableCoefficient) {
639 case COEFF_CHECKERBOARD_0: {
640 PetscScalar constants[2];
641
642 constants[0] = user->div;
643 constants[1] = user->k;
644 PetscCall(PetscDSSetConstants(ds, 2, constants));
645 } break;
646 default:
647 break;
648 }
649 PetscCall(PetscDSSetExactSolution(ds, 0, user->exactFuncs[0], user));
650 /* Setup Boundary Conditions */
651 if (user->bcType == DIRICHLET) {
652 PetscCall(DMGetLabel(dm, "marker", &label));
653 if (!label) {
654 /* Right now, p4est cannot create labels immediately */
655 PetscCall(PetscDSAddBoundaryByName(ds, user->fieldBC ? DM_BC_ESSENTIAL_FIELD : DM_BC_ESSENTIAL, "wall", "marker", 1, &id, 0, 0, NULL, user->fieldBC ? (PetscVoidFn *)user->exactFields[0] : (PetscVoidFn *)user->exactFuncs[0], NULL, user, NULL));
656 } else {
657 PetscCall(DMAddBoundary(dm, user->fieldBC ? DM_BC_ESSENTIAL_FIELD : DM_BC_ESSENTIAL, "wall", label, 1, &id, 0, 0, NULL, user->fieldBC ? (PetscVoidFn *)user->exactFields[0] : (PetscVoidFn *)user->exactFuncs[0], NULL, user, NULL));
658 }
659 }
660 PetscFunctionReturn(PETSC_SUCCESS);
661 }
662
SetupMaterial(DM dm,DM dmAux,AppCtx * user)663 static PetscErrorCode SetupMaterial(DM dm, DM dmAux, AppCtx *user)
664 {
665 PetscErrorCode (*matFuncs[1])(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar u[], PetscCtx ctx) = {nu_2d};
666 PetscCtx ctx[1];
667 Vec nu;
668
669 PetscFunctionBegin;
670 ctx[0] = user;
671 if (user->variableCoefficient == COEFF_CHECKERBOARD_0) matFuncs[0] = checkerboardCoeff;
672 PetscCall(DMCreateLocalVector(dmAux, &nu));
673 PetscCall(PetscObjectSetName((PetscObject)nu, "Coefficient"));
674 PetscCall(DMProjectFunctionLocal(dmAux, 0.0, matFuncs, ctx, INSERT_ALL_VALUES, nu));
675 PetscCall(DMSetAuxiliaryVec(dm, NULL, 0, 0, nu));
676 PetscCall(VecDestroy(&nu));
677 PetscFunctionReturn(PETSC_SUCCESS);
678 }
679
SetupBC(DM dm,DM dmAux,AppCtx * user)680 static PetscErrorCode SetupBC(DM dm, DM dmAux, AppCtx *user)
681 {
682 PetscErrorCode (*bcFuncs[1])(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar u[], PetscCtx ctx);
683 Vec uexact;
684 PetscInt dim;
685
686 PetscFunctionBegin;
687 PetscCall(DMGetDimension(dm, &dim));
688 if (dim == 2) bcFuncs[0] = quadratic_u_2d;
689 else bcFuncs[0] = quadratic_u_3d;
690 PetscCall(DMCreateLocalVector(dmAux, &uexact));
691 PetscCall(DMProjectFunctionLocal(dmAux, 0.0, bcFuncs, NULL, INSERT_ALL_VALUES, uexact));
692 PetscCall(DMSetAuxiliaryVec(dm, NULL, 0, 0, uexact));
693 PetscCall(VecDestroy(&uexact));
694 PetscFunctionReturn(PETSC_SUCCESS);
695 }
696
SetupAuxDM(DM dm,PetscFE feAux,AppCtx * user)697 static PetscErrorCode SetupAuxDM(DM dm, PetscFE feAux, AppCtx *user)
698 {
699 DM dmAux, coordDM;
700
701 PetscFunctionBegin;
702 /* MUST call DMGetCoordinateDM() in order to get p4est setup if present */
703 PetscCall(DMGetCoordinateDM(dm, &coordDM));
704 if (!feAux) PetscFunctionReturn(PETSC_SUCCESS);
705 PetscCall(DMClone(dm, &dmAux));
706 PetscCall(DMSetCoordinateDM(dmAux, coordDM));
707 PetscCall(DMSetField(dmAux, 0, NULL, (PetscObject)feAux));
708 PetscCall(DMCreateDS(dmAux));
709 if (user->fieldBC) PetscCall(SetupBC(dm, dmAux, user));
710 else PetscCall(SetupMaterial(dm, dmAux, user));
711 PetscCall(DMDestroy(&dmAux));
712 PetscFunctionReturn(PETSC_SUCCESS);
713 }
714
SetupDiscretization(DM dm,AppCtx * user)715 static PetscErrorCode SetupDiscretization(DM dm, AppCtx *user)
716 {
717 DM plex, cdm = dm;
718 PetscFE fe, feAux = NULL;
719 PetscBool simplex;
720 PetscInt dim;
721
722 PetscFunctionBeginUser;
723 PetscCall(DMGetDimension(dm, &dim));
724 PetscCall(DMConvert(dm, DMPLEX, &plex));
725 PetscCall(DMPlexIsSimplex(plex, &simplex));
726 PetscCall(DMDestroy(&plex));
727 PetscCall(PetscFECreateDefault(PETSC_COMM_SELF, dim, 1, simplex, NULL, -1, &fe));
728 PetscCall(PetscObjectSetName((PetscObject)fe, "potential"));
729 if (user->variableCoefficient == COEFF_FIELD || user->variableCoefficient == COEFF_CHECKERBOARD_0) {
730 PetscCall(PetscFECreateDefault(PETSC_COMM_SELF, dim, 1, simplex, "mat_", -1, &feAux));
731 PetscCall(PetscObjectSetName((PetscObject)feAux, "coefficient"));
732 PetscCall(PetscFECopyQuadrature(fe, feAux));
733 } else if (user->fieldBC) {
734 PetscCall(PetscFECreateDefault(PETSC_COMM_SELF, dim, 1, simplex, "bc_", -1, &feAux));
735 PetscCall(PetscFECopyQuadrature(fe, feAux));
736 }
737 /* Set discretization and boundary conditions for each mesh */
738 PetscCall(DMSetField(dm, 0, NULL, (PetscObject)fe));
739 PetscCall(DMCreateDS(dm));
740 PetscCall(SetupProblem(dm, user));
741 while (cdm) {
742 PetscCall(SetupAuxDM(cdm, feAux, user));
743 if (user->bcType == DIRICHLET) {
744 PetscBool hasLabel;
745
746 PetscCall(DMHasLabel(cdm, "marker", &hasLabel));
747 if (!hasLabel) PetscCall(CreateBCLabel(cdm, "marker"));
748 }
749 PetscCall(DMCopyDisc(dm, cdm));
750 PetscCall(DMGetCoarseDM(cdm, &cdm));
751 }
752 PetscCall(PetscFEDestroy(&fe));
753 PetscCall(PetscFEDestroy(&feAux));
754 PetscFunctionReturn(PETSC_SUCCESS);
755 }
756
main(int argc,char ** argv)757 int main(int argc, char **argv)
758 {
759 DM dm; /* Problem specification */
760 SNES snes; /* nonlinear solver */
761 Vec u; /* solution vector */
762 Mat A, J; /* Jacobian matrix */
763 MatNullSpace nullSpace; /* May be necessary for Neumann conditions */
764 AppCtx user; /* user-defined work context */
765 JacActionCtx userJ; /* context for Jacobian MF action */
766 PetscReal error = 0.0; /* L_2 error in the solution */
767
768 PetscFunctionBeginUser;
769 PetscCall(PetscInitialize(&argc, &argv, NULL, help));
770 PetscCall(ProcessOptions(PETSC_COMM_WORLD, &user));
771 PetscCall(SNESCreate(PETSC_COMM_WORLD, &snes));
772 PetscCall(CreateMesh(PETSC_COMM_WORLD, &user, &dm));
773 PetscCall(SNESSetDM(snes, dm));
774 PetscCall(DMSetApplicationContext(dm, &user));
775
776 PetscCall(PetscMalloc2(1, &user.exactFuncs, 1, &user.exactFields));
777 PetscCall(SetupDiscretization(dm, &user));
778
779 PetscCall(DMCreateGlobalVector(dm, &u));
780 PetscCall(PetscObjectSetName((PetscObject)u, "potential"));
781
782 PetscCall(DMCreateMatrix(dm, &J));
783 if (user.jacobianMF) {
784 PetscInt M, m, N, n;
785
786 PetscCall(MatGetSize(J, &M, &N));
787 PetscCall(MatGetLocalSize(J, &m, &n));
788 PetscCall(MatCreate(PETSC_COMM_WORLD, &A));
789 PetscCall(MatSetSizes(A, m, n, M, N));
790 PetscCall(MatSetType(A, MATSHELL));
791 PetscCall(MatSetUp(A));
792 #if 0
793 PetscCall(MatShellSetOperation(A, MATOP_MULT, (PetscErrorCodeFn *)FormJacobianAction));
794 #endif
795
796 userJ.dm = dm;
797 userJ.J = J;
798 userJ.user = &user;
799
800 PetscCall(DMCreateLocalVector(dm, &userJ.u));
801 if (user.fieldBC) PetscCall(DMProjectFieldLocal(dm, 0.0, userJ.u, user.exactFields, INSERT_BC_VALUES, userJ.u));
802 else PetscCall(DMProjectFunctionLocal(dm, 0.0, user.exactFuncs, NULL, INSERT_BC_VALUES, userJ.u));
803 PetscCall(MatShellSetContext(A, &userJ));
804 } else {
805 A = J;
806 }
807
808 nullSpace = NULL;
809 if (user.bcType != DIRICHLET) {
810 PetscCall(MatNullSpaceCreate(PetscObjectComm((PetscObject)dm), PETSC_TRUE, 0, NULL, &nullSpace));
811 PetscCall(MatSetNullSpace(A, nullSpace));
812 }
813
814 PetscCall(DMPlexSetSNESLocalFEM(dm, PETSC_FALSE, &user));
815 PetscCall(SNESSetJacobian(snes, A, J, NULL, NULL));
816
817 PetscCall(SNESSetFromOptions(snes));
818
819 if (user.fieldBC) PetscCall(DMProjectField(dm, 0.0, u, user.exactFields, INSERT_ALL_VALUES, u));
820 else PetscCall(DMProjectFunction(dm, 0.0, user.exactFuncs, NULL, INSERT_ALL_VALUES, u));
821 if (user.restart) {
822 #if defined(PETSC_HAVE_HDF5)
823 PetscViewer viewer;
824 char filename[PETSC_MAX_PATH_LEN];
825
826 PetscCall(PetscOptionsGetString(NULL, NULL, "-dm_plex_filename", filename, sizeof(filename), NULL));
827 PetscCall(PetscViewerCreate(PETSC_COMM_WORLD, &viewer));
828 PetscCall(PetscViewerSetType(viewer, PETSCVIEWERHDF5));
829 PetscCall(PetscViewerFileSetMode(viewer, FILE_MODE_READ));
830 PetscCall(PetscViewerFileSetName(viewer, filename));
831 PetscCall(PetscViewerHDF5PushGroup(viewer, "/fields"));
832 PetscCall(VecLoad(u, viewer));
833 PetscCall(PetscViewerHDF5PopGroup(viewer));
834 PetscCall(PetscViewerDestroy(&viewer));
835 #endif
836 }
837 if (user.showInitial) {
838 Vec lv;
839 PetscCall(DMGetLocalVector(dm, &lv));
840 PetscCall(DMGlobalToLocalBegin(dm, u, INSERT_VALUES, lv));
841 PetscCall(DMGlobalToLocalEnd(dm, u, INSERT_VALUES, lv));
842 PetscCall(DMPrintLocalVec(dm, "Local function", 1.0e-10, lv));
843 PetscCall(DMRestoreLocalVector(dm, &lv));
844 }
845 if (user.runType == RUN_FULL || user.runType == RUN_EXACT) {
846 PetscErrorCode (*initialGuess[1])(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar u[], PetscCtx ctx) = {zero};
847
848 if (user.nonzInit) initialGuess[0] = ecks;
849 if (user.runType == RUN_FULL) PetscCall(DMProjectFunction(dm, 0.0, initialGuess, NULL, INSERT_VALUES, u));
850 PetscCall(VecViewFromOptions(u, NULL, "-guess_vec_view"));
851 PetscCall(SNESSolve(snes, NULL, u));
852 PetscCall(SNESGetSolution(snes, &u));
853 PetscCall(SNESGetDM(snes, &dm));
854
855 if (user.showSolution) {
856 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Solution\n"));
857 PetscCall(VecFilter(u, 3.0e-9));
858 PetscCall(VecView(u, PETSC_VIEWER_STDOUT_WORLD));
859 }
860 } else if (user.runType == RUN_PERF) {
861 Vec r;
862 PetscReal res = 0.0;
863
864 PetscCall(SNESGetFunction(snes, &r, NULL, NULL));
865 PetscCall(SNESComputeFunction(snes, u, r));
866 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Initial Residual\n"));
867 PetscCall(VecFilter(r, 1.0e-10));
868 PetscCall(VecNorm(r, NORM_2, &res));
869 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "L_2 Residual: %g\n", (double)res));
870 } else {
871 Vec r;
872 PetscReal res = 0.0, tol = 1.0e-11;
873
874 /* Check discretization error */
875 PetscCall(SNESGetFunction(snes, &r, NULL, NULL));
876 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Initial guess\n"));
877 if (!user.quiet) PetscCall(VecView(u, PETSC_VIEWER_STDOUT_WORLD));
878 PetscCall(DMComputeL2Diff(dm, 0.0, user.exactFuncs, NULL, u, &error));
879 if (error < tol) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "L_2 Error: < %2.1e\n", (double)tol));
880 else PetscCall(PetscPrintf(PETSC_COMM_WORLD, "L_2 Error: %g\n", (double)error));
881 /* Check residual */
882 PetscCall(SNESComputeFunction(snes, u, r));
883 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Initial Residual\n"));
884 PetscCall(VecFilter(r, 1.0e-10));
885 if (!user.quiet) PetscCall(VecView(r, PETSC_VIEWER_STDOUT_WORLD));
886 PetscCall(VecNorm(r, NORM_2, &res));
887 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "L_2 Residual: %g\n", (double)res));
888 /* Check Jacobian */
889 {
890 Vec b;
891
892 PetscCall(SNESComputeJacobian(snes, u, A, A));
893 PetscCall(VecDuplicate(u, &b));
894 PetscCall(VecSet(r, 0.0));
895 PetscCall(SNESComputeFunction(snes, r, b));
896 PetscCall(MatMult(A, u, r));
897 PetscCall(VecAXPY(r, 1.0, b));
898 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Au - b = Au + F(0)\n"));
899 PetscCall(VecFilter(r, 1.0e-10));
900 if (!user.quiet) PetscCall(VecView(r, PETSC_VIEWER_STDOUT_WORLD));
901 PetscCall(VecNorm(r, NORM_2, &res));
902 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Linear L_2 Residual: %g\n", (double)res));
903 /* check solver */
904 if (user.checkksp) {
905 KSP ksp;
906
907 if (nullSpace) PetscCall(MatNullSpaceRemove(nullSpace, u));
908 PetscCall(SNESComputeJacobian(snes, u, A, J));
909 PetscCall(MatMult(A, u, b));
910 PetscCall(SNESGetKSP(snes, &ksp));
911 PetscCall(KSPSetOperators(ksp, A, J));
912 PetscCall(KSPSolve(ksp, b, r));
913 PetscCall(VecAXPY(r, -1.0, u));
914 PetscCall(VecNorm(r, NORM_2, &res));
915 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "KSP Error: %g\n", (double)res));
916 }
917 PetscCall(VecDestroy(&b));
918 }
919 }
920 PetscCall(VecViewFromOptions(u, NULL, "-vec_view"));
921 {
922 Vec nu;
923
924 PetscCall(DMGetAuxiliaryVec(dm, NULL, 0, 0, &nu));
925 if (nu) PetscCall(VecViewFromOptions(nu, NULL, "-coeff_view"));
926 }
927
928 if (user.bdIntegral) {
929 DMLabel label;
930 PetscBdPointFn *func[1] = {bd_integral_2d};
931 PetscInt id = 1;
932 PetscScalar bdInt = 0.0;
933 PetscReal exact = 3.3333333333;
934
935 PetscCall(DMGetLabel(dm, "marker", &label));
936 PetscCall(DMPlexComputeBdIntegral(dm, u, label, 1, &id, func, &bdInt, NULL));
937 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Solution boundary integral: %.4g\n", (double)PetscAbsScalar(bdInt)));
938 PetscCheck(PetscAbsReal(PetscAbsScalar(bdInt) - exact) <= PETSC_SQRT_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Invalid boundary integral %g != %g", (double)PetscAbsScalar(bdInt), (double)exact);
939 }
940
941 PetscCall(MatNullSpaceDestroy(&nullSpace));
942 if (user.jacobianMF) PetscCall(VecDestroy(&userJ.u));
943 if (A != J) PetscCall(MatDestroy(&A));
944 PetscCall(MatDestroy(&J));
945 PetscCall(VecDestroy(&u));
946 PetscCall(SNESDestroy(&snes));
947 PetscCall(DMDestroy(&dm));
948 PetscCall(PetscFree2(user.exactFuncs, user.exactFields));
949 PetscCall(PetscFree(user.kgrid));
950 PetscCall(PetscFinalize());
951 return 0;
952 }
953
954 /*TEST
955 # 2D serial P1 test 0-4
956 test:
957 suffix: 2d_p1_0
958 requires: triangle
959 args: -run_type test -bc_type dirichlet -dm_plex_interpolate 0 -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
960
961 test:
962 suffix: 2d_p1_1
963 requires: triangle
964 args: -run_type test -bc_type dirichlet -petscspace_degree 1 -show_initial -dm_plex_print_fem 1 -cdm_dm_plex_coordinate_dim {{2 3}}
965
966 test:
967 suffix: 2d_p1_1b
968 requires: triangle
969 args: -run_type test -bc_type dirichlet -petscspace_degree 1 -show_initial -dm_plex_print_fem 1 -dm_refine 3 -dm_coord_space 0 \
970 -dm_plex_option_phases proj_ -cdm_proj_dm_plex_coordinate_dim 3 -proj_dm_coord_space \
971 -proj_dm_coord_remap -proj_dm_coord_map sinusoid -proj_dm_coord_map_params 0.1,1.,1.
972
973 test:
974 suffix: 2d_p1_2
975 requires: triangle
976 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -bc_type dirichlet -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
977
978 test:
979 suffix: 2d_p1_neumann_0
980 requires: triangle
981 args: -dm_coord_space 0 -run_type test -bc_type neumann -dm_plex_boundary_label boundary -petscspace_degree 1 -show_initial -dm_plex_print_fem 1 -dm_view ascii::ascii_info_detail
982
983 test:
984 suffix: 2d_p1_neumann_1
985 requires: triangle
986 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -bc_type neumann -dm_plex_boundary_label boundary -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
987
988 # 2D serial P2 test 5-8
989 test:
990 suffix: 2d_p2_0
991 requires: triangle
992 args: -run_type test -bc_type dirichlet -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
993
994 test:
995 suffix: 2d_p2_1
996 requires: triangle
997 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -bc_type dirichlet -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
998
999 test:
1000 suffix: 2d_p2_neumann_0
1001 requires: triangle
1002 args: -dm_coord_space 0 -run_type test -bc_type neumann -dm_plex_boundary_label boundary -petscspace_degree 2 -show_initial -dm_plex_print_fem 1 -dm_view ascii::ascii_info_detail
1003
1004 test:
1005 suffix: 2d_p2_neumann_1
1006 requires: triangle
1007 args: -dm_coord_space 0 -run_type test -dm_refine_volume_limit_pre 0.0625 -bc_type neumann -dm_plex_boundary_label boundary -petscspace_degree 2 -show_initial -dm_plex_print_fem 1 -dm_view ascii::ascii_info_detail
1008
1009 test:
1010 suffix: bd_int_0
1011 requires: triangle
1012 args: -run_type test -bc_type dirichlet -petscspace_degree 2 -bd_integral -dm_view -quiet
1013
1014 test:
1015 suffix: bd_int_1
1016 requires: triangle
1017 args: -run_type test -dm_refine 2 -bc_type dirichlet -petscspace_degree 2 -bd_integral -dm_view -quiet
1018
1019 # 3D serial P1 test 9-12
1020 test:
1021 suffix: 3d_p1_0
1022 requires: ctetgen
1023 args: -run_type test -dm_plex_dim 3 -bc_type dirichlet -dm_plex_interpolate 0 -petscspace_degree 1 -show_initial -dm_plex_print_fem 1 -dm_view
1024
1025 test:
1026 suffix: 3d_p1_1
1027 requires: ctetgen
1028 args: -run_type test -dm_plex_dim 3 -bc_type dirichlet -petscspace_degree 1 -show_initial -dm_plex_print_fem 1 -dm_view
1029
1030 test:
1031 suffix: 3d_p1_2
1032 requires: ctetgen
1033 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -bc_type dirichlet -petscspace_degree 1 -show_initial -dm_plex_print_fem 1 -dm_view
1034
1035 test:
1036 suffix: 3d_p1_neumann_0
1037 requires: ctetgen
1038 args: -run_type test -dm_plex_dim 3 -bc_type neumann -dm_plex_boundary_label boundary -petscspace_degree 1 -snes_fd -show_initial -dm_plex_print_fem 1 -dm_view
1039
1040 # Analytic variable coefficient 13-20
1041 test:
1042 suffix: 13
1043 requires: triangle
1044 args: -run_type test -variable_coefficient analytic -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1045 test:
1046 suffix: 14
1047 requires: triangle
1048 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -variable_coefficient analytic -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1049 test:
1050 suffix: 15
1051 requires: triangle
1052 args: -run_type test -variable_coefficient analytic -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1053 test:
1054 suffix: 16
1055 requires: triangle
1056 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -variable_coefficient analytic -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1057 test:
1058 suffix: 17
1059 requires: ctetgen
1060 args: -run_type test -dm_plex_dim 3 -variable_coefficient analytic -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1061
1062 test:
1063 suffix: 18
1064 requires: ctetgen
1065 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -variable_coefficient analytic -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1066
1067 test:
1068 suffix: 19
1069 requires: ctetgen
1070 args: -run_type test -dm_plex_dim 3 -variable_coefficient analytic -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1071
1072 test:
1073 suffix: 20
1074 requires: ctetgen
1075 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -variable_coefficient analytic -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1076
1077 # P1 variable coefficient 21-28
1078 test:
1079 suffix: 21
1080 requires: triangle
1081 args: -run_type test -variable_coefficient field -petscspace_degree 1 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1082
1083 test:
1084 suffix: 22
1085 requires: triangle
1086 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -variable_coefficient field -petscspace_degree 1 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1087
1088 test:
1089 suffix: 23
1090 requires: triangle
1091 args: -run_type test -variable_coefficient field -petscspace_degree 2 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1092
1093 test:
1094 suffix: 24
1095 requires: triangle
1096 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -variable_coefficient field -petscspace_degree 2 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1097
1098 test:
1099 suffix: 25
1100 requires: ctetgen
1101 args: -run_type test -dm_plex_dim 3 -variable_coefficient field -petscspace_degree 1 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1102
1103 test:
1104 suffix: 26
1105 requires: ctetgen
1106 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -variable_coefficient field -petscspace_degree 1 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1107
1108 test:
1109 suffix: 27
1110 requires: ctetgen
1111 args: -run_type test -dm_plex_dim 3 -variable_coefficient field -petscspace_degree 2 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1112
1113 test:
1114 suffix: 28
1115 requires: ctetgen
1116 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -variable_coefficient field -petscspace_degree 2 -mat_petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1117
1118 # P0 variable coefficient 29-36
1119 test:
1120 suffix: 29
1121 requires: triangle
1122 args: -run_type test -variable_coefficient field -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1123
1124 test:
1125 suffix: 30
1126 requires: triangle
1127 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -variable_coefficient field -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1128
1129 test:
1130 suffix: 31
1131 requires: triangle
1132 args: -run_type test -variable_coefficient field -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1133
1134 test:
1135 requires: triangle
1136 suffix: 32
1137 args: -run_type test -dm_refine_volume_limit_pre 0.0625 -variable_coefficient field -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1138
1139 test:
1140 requires: ctetgen
1141 suffix: 33
1142 args: -run_type test -dm_plex_dim 3 -variable_coefficient field -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1143
1144 test:
1145 suffix: 34
1146 requires: ctetgen
1147 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -variable_coefficient field -petscspace_degree 1 -show_initial -dm_plex_print_fem 1
1148
1149 test:
1150 suffix: 35
1151 requires: ctetgen
1152 args: -run_type test -dm_plex_dim 3 -variable_coefficient field -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1153
1154 test:
1155 suffix: 36
1156 requires: ctetgen
1157 args: -run_type test -dm_plex_dim 3 -dm_refine_volume_limit_pre 0.0125 -variable_coefficient field -petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1158
1159 # Full solve 39-44
1160 test:
1161 suffix: 39
1162 requires: triangle !single
1163 args: -run_type full -dm_refine_volume_limit_pre 0.015625 -petscspace_degree 2 -pc_type gamg -pc_gamg_esteig_ksp_type cg -pc_gamg_esteig_ksp_max_it 10 -snes_rtol 1.0e-6 -ksp_rtol 1.0e-7 -ksp_monitor -ksp_converged_reason -snes_monitor_short -snes_converged_reason ::ascii_info_detail
1164 test:
1165 suffix: 40
1166 requires: triangle !single
1167 args: -run_type full -dm_refine_volume_limit_pre 0.015625 -variable_coefficient nonlinear -petscspace_degree 2 -pc_type svd -ksp_rtol 1.0e-10 -snes_monitor_short -snes_converged_reason ::ascii_info_detail
1168 test:
1169 suffix: 41
1170 requires: triangle !single
1171 args: -run_type full -dm_refine_volume_limit_pre 0.03125 -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 2 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -dm_refine_hierarchy 1 -snes_view -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short
1172 test:
1173 suffix: 42
1174 requires: triangle !single
1175 args: -run_type full -dm_refine_volume_limit_pre 0.0625 -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -dm_refine_hierarchy 2 -snes_view -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short -fas_levels_2_snes_type newtonls -fas_levels_2_pc_type svd -fas_levels_2_ksp_rtol 1.0e-10 -fas_levels_2_snes_atol 1.0e-11 -fas_levels_2_snes_monitor_short
1176 test:
1177 suffix: 43
1178 requires: triangle !single
1179 nsize: 2
1180 args: -run_type full -dm_refine_volume_limit_pre 0.03125 -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 2 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -dm_refine_hierarchy 1 -snes_view -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short
1181
1182 test:
1183 suffix: 44
1184 requires: triangle !single
1185 nsize: 2
1186 args: -run_type full -dm_refine_volume_limit_pre 0.0625 -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -dm_refine_hierarchy 2 -dm_plex_print_fem 0 -snes_view -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short -fas_levels_2_snes_type newtonls -fas_levels_2_pc_type svd -fas_levels_2_ksp_rtol 1.0e-10 -fas_levels_2_snes_atol 1.0e-11 -fas_levels_2_snes_monitor_short
1187
1188 # These tests use a loose tolerance just to exercise the PtAP operations for MATIS and multiple PCBDDC setup calls inside PCMG
1189 testset:
1190 requires: triangle !single
1191 nsize: 3
1192 args: -run_type full -petscspace_degree 1 -dm_mat_type is -pc_type mg -mg_coarse_pc_type bddc -pc_mg_galerkin pmat -ksp_rtol 1.0e-2 -snes_converged_reason -dm_refine_hierarchy 2 -snes_max_it 4
1193 test:
1194 suffix: gmg_bddc
1195 filter: sed -e "s/CONVERGED_FNORM_RELATIVE iterations 3/CONVERGED_FNORM_RELATIVE iterations 4/g"
1196 args: -mg_levels_pc_type jacobi
1197 test:
1198 filter: sed -e "s/iterations [0-4]/iterations 4/g"
1199 suffix: gmg_bddc_lev
1200 args: -mg_levels_pc_type bddc
1201
1202 # VTU viewer with empty processes
1203 test:
1204 requires: !complex
1205 suffix: vtu_empty
1206 args: -quiet -run_type test -dm_plex_simplex 0 -dm_plex_box_faces 2,2 -vec_view vtk:test.vtu:vtk_vtu -petscspace_degree 1 -petscpartitioner_type simple
1207
1208 # Restarting
1209 testset:
1210 suffix: restart
1211 requires: hdf5 triangle !complex
1212 args: -run_type test -bc_type dirichlet -petscspace_degree 1
1213 test:
1214 args: -dm_view hdf5:sol.h5 -vec_view hdf5:sol.h5::append
1215 test:
1216 args: -dm_plex_filename sol.h5 -dm_plex_name box -restart
1217
1218 # Periodicity
1219 test:
1220 suffix: periodic_0
1221 requires: triangle
1222 args: -run_type full -bc_type dirichlet -petscspace_degree 1 -snes_converged_reason ::ascii_info_detail
1223
1224 test:
1225 requires: !complex
1226 suffix: periodic_1
1227 args: -quiet -run_type test -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -dm_plex_box_bd periodic,periodic -vec_view vtk:test.vtu:vtk_vtu -petscspace_degree 1 -dm_refine 1
1228
1229 # 2D serial P1 test with field bc
1230 test:
1231 suffix: field_bc_2d_p1_0
1232 requires: triangle
1233 args: -run_type test -bc_type dirichlet -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1234
1235 test:
1236 suffix: field_bc_2d_p1_1
1237 requires: triangle
1238 args: -run_type test -dm_refine 1 -bc_type dirichlet -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1239
1240 test:
1241 suffix: field_bc_2d_p1_neumann_0
1242 requires: triangle
1243 args: -run_type test -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1244
1245 test:
1246 suffix: field_bc_2d_p1_neumann_1
1247 requires: triangle
1248 args: -run_type test -dm_refine 1 -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1249
1250 # 3D serial P1 test with field bc
1251 test:
1252 suffix: field_bc_3d_p1_0
1253 requires: ctetgen
1254 args: -run_type test -dm_plex_dim 3 -bc_type dirichlet -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1255
1256 test:
1257 suffix: field_bc_3d_p1_1
1258 requires: ctetgen
1259 args: -run_type test -dm_plex_dim 3 -dm_refine 1 -bc_type dirichlet -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1260
1261 test:
1262 suffix: field_bc_3d_p1_neumann_0
1263 requires: ctetgen
1264 args: -run_type test -dm_plex_dim 3 -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1265
1266 test:
1267 suffix: field_bc_3d_p1_neumann_1
1268 requires: ctetgen
1269 args: -run_type test -dm_plex_dim 3 -dm_refine 1 -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 1 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1270
1271 # 2D serial P2 test with field bc
1272 test:
1273 suffix: field_bc_2d_p2_0
1274 requires: triangle
1275 args: -run_type test -bc_type dirichlet -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1276
1277 test:
1278 suffix: field_bc_2d_p2_1
1279 requires: triangle
1280 args: -run_type test -dm_refine 1 -bc_type dirichlet -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1281
1282 test:
1283 suffix: field_bc_2d_p2_neumann_0
1284 requires: triangle
1285 args: -run_type test -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1286
1287 test:
1288 suffix: field_bc_2d_p2_neumann_1
1289 requires: triangle
1290 args: -run_type test -dm_refine 1 -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1291
1292 # 3D serial P2 test with field bc
1293 test:
1294 suffix: field_bc_3d_p2_0
1295 requires: ctetgen
1296 args: -run_type test -dm_plex_dim 3 -bc_type dirichlet -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1297
1298 test:
1299 suffix: field_bc_3d_p2_1
1300 requires: ctetgen
1301 args: -run_type test -dm_plex_dim 3 -dm_refine 1 -bc_type dirichlet -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1302
1303 test:
1304 suffix: field_bc_3d_p2_neumann_0
1305 requires: ctetgen
1306 args: -run_type test -dm_plex_dim 3 -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1307
1308 test:
1309 suffix: field_bc_3d_p2_neumann_1
1310 requires: ctetgen
1311 args: -run_type test -dm_plex_dim 3 -dm_refine 1 -bc_type neumann -dm_plex_boundary_label boundary -field_bc -petscspace_degree 2 -bc_petscspace_degree 2 -show_initial -dm_plex_print_fem 1
1312
1313 # Full solve simplex: Convergence
1314 test:
1315 suffix: 3d_p1_conv
1316 requires: ctetgen
1317 args: -run_type full -dm_plex_dim 3 -dm_refine 1 -bc_type dirichlet -petscspace_degree 1 \
1318 -snes_convergence_estimate -convest_num_refine 1 -pc_type lu
1319
1320 # Full solve simplex: PCBDDC
1321 test:
1322 suffix: tri_bddc
1323 requires: triangle !single
1324 nsize: 5
1325 args: -run_type full -petscpartitioner_type simple -dm_refine 2 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -dm_mat_type is -pc_type bddc -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0
1326
1327 # Full solve simplex: PCBDDC
1328 test:
1329 suffix: tri_parmetis_bddc
1330 requires: triangle !single parmetis
1331 nsize: 4
1332 args: -run_type full -petscpartitioner_type parmetis -dm_refine 2 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -dm_mat_type is -pc_type bddc -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0
1333
1334 testset:
1335 args: -run_type full -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -petscpartitioner_type simple -dm_refine 2 -bc_type dirichlet -petscspace_degree 2 -dm_mat_type is -pc_type bddc -ksp_type gmres -snes_monitor_short -ksp_monitor_short -snes_view -petscspace_poly_tensor -pc_bddc_corner_selection -ksp_rtol 1.e-9 -pc_bddc_use_edges 0
1336 nsize: 5
1337 output_file: output/ex12_quad_bddc.out
1338 filter: sed -e "s/aijcusparse/aij/g" -e "s/aijviennacl/aij/g" -e "s/factorization: cusparse/factorization: petsc/g"
1339 test:
1340 requires: !single
1341 suffix: quad_bddc
1342 test:
1343 requires: !single cuda
1344 suffix: quad_bddc_cuda
1345 args: -mat_is_localmat_type aijcusparse -pc_bddc_dirichlet_pc_factor_mat_solver_type cusparse -pc_bddc_neumann_pc_factor_mat_solver_type cusparse
1346 test:
1347 requires: !single viennacl
1348 suffix: quad_bddc_viennacl
1349 args: -mat_is_localmat_type aijviennacl
1350
1351 # Full solve simplex: ASM
1352 test:
1353 suffix: tri_q2q1_asm_lu
1354 requires: triangle !single
1355 args: -run_type full -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -pc_type asm -pc_asm_type restrict -pc_asm_blocks 4 -sub_pc_type lu -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0
1356
1357 test:
1358 suffix: tri_q2q1_msm_lu
1359 requires: triangle !single
1360 args: -run_type full -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -pc_type asm -pc_asm_type restrict -pc_asm_local_type multiplicative -pc_asm_blocks 4 -sub_pc_type lu -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0
1361
1362 test:
1363 suffix: tri_q2q1_asm_sor
1364 requires: triangle !single
1365 args: -run_type full -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -pc_type asm -pc_asm_type restrict -pc_asm_blocks 4 -sub_pc_type sor -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0
1366
1367 test:
1368 suffix: tri_q2q1_msm_sor
1369 requires: triangle !single
1370 args: -run_type full -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -ksp_rtol 1.0e-9 -pc_type asm -pc_asm_type restrict -pc_asm_local_type multiplicative -pc_asm_blocks 4 -sub_pc_type sor -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0
1371
1372 # Full solve simplex: FAS
1373 test:
1374 suffix: fas_newton_0
1375 requires: triangle !single
1376 args: -run_type full -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 2 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -dm_refine_hierarchy 1 -snes_view -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short
1377
1378 test:
1379 suffix: fas_newton_1
1380 requires: triangle !single
1381 args: -run_type full -dm_refine_hierarchy 3 -petscspace_degree 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type lu -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_snes_linesearch_type basic -fas_levels_ksp_rtol 1.0e-10 -fas_levels_snes_monitor_short
1382 filter: sed -e "s/total number of linear solver iterations=14/total number of linear solver iterations=15/g"
1383
1384 test:
1385 suffix: fas_ngs_0
1386 requires: triangle !single
1387 args: -run_type full -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 2 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -dm_refine_hierarchy 1 -snes_view -fas_levels_1_snes_type ngs -fas_levels_1_snes_monitor_short
1388
1389 # These two tests are broken because DMPlexComputeInjectorFEM() only works for regularly refined meshes
1390 test:
1391 suffix: fas_newton_coarse_0
1392 requires: pragmatic triangle
1393 TODO: broken
1394 args: -run_type full -variable_coefficient nonlinear -petscspace_degree 1 \
1395 -dm_refine 2 -dm_coarsen_hierarchy 1 -dm_plex_hash_location -dm_adaptor pragmatic \
1396 -snes_type fas -snes_fas_levels 2 -snes_converged_reason ::ascii_info_detail -snes_monitor_short -snes_view \
1397 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -fas_coarse_snes_linesearch_type basic \
1398 -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short
1399
1400 test:
1401 suffix: mg_newton_coarse_0
1402 requires: triangle pragmatic
1403 TODO: broken
1404 args: -run_type full -petscspace_degree 1 \
1405 -dm_refine 3 -dm_coarsen_hierarchy 3 -dm_plex_hash_location -dm_adaptor pragmatic \
1406 -snes_atol 1.0e-8 -snes_rtol 0.0 -snes_monitor_short -snes_converged_reason ::ascii_info_detail -snes_view \
1407 -ksp_type richardson -ksp_atol 1.0e-8 -ksp_rtol 0.0 -ksp_norm_type unpreconditioned -ksp_monitor_true_residual \
1408 -pc_type mg -pc_mg_levels 4 \
1409 -mg_levels_ksp_type gmres -mg_levels_pc_type ilu -mg_levels_ksp_max_it 10
1410
1411 # Test cgns writer for ranks with no elements
1412 test:
1413 suffix: cgns
1414 nsize: 5
1415 requires: cgns
1416 args: -quiet -run_type test -dm_plex_simplex 0 -petscspace_degree 1 -dm_plex_box_faces 2,2 -vec_view cgns:test.cgns -dm_refine 0 -petscpartitioner_type simple
1417
1418 # Full solve tensor
1419 test:
1420 suffix: tensor_plex_2d
1421 args: -run_type test -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -dm_refine_hierarchy 2
1422
1423 test:
1424 suffix: tensor_p4est_2d
1425 requires: p4est
1426 args: -run_type test -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -dm_forest_initial_refinement 2 -dm_forest_minimum_refinement 0 -dm_plex_convert_type p4est
1427
1428 test:
1429 suffix: tensor_plex_3d
1430 args: -run_type test -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -dm_plex_dim 3 -dm_refine_hierarchy 1 -dm_plex_box_faces 2,2,2
1431
1432 test:
1433 suffix: tensor_p4est_3d
1434 requires: p4est
1435 args: -run_type test -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -dm_forest_initial_refinement 1 -dm_forest_minimum_refinement 0 -dm_plex_dim 3 -dm_plex_convert_type p8est -dm_plex_box_faces 2,2,2
1436
1437 test:
1438 suffix: p4est_test_q2_conformal_serial
1439 requires: p4est
1440 args: -run_type test -petscspace_degree 2 -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2
1441
1442 test:
1443 suffix: p4est_test_q2_conformal_parallel
1444 requires: p4est
1445 nsize: 7
1446 args: -run_type test -petscspace_degree 2 -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -petscpartitioner_type simple
1447
1448 test:
1449 suffix: p4est_test_q2_conformal_parallel_parmetis
1450 requires: parmetis p4est
1451 nsize: 4
1452 args: -run_type test -petscspace_degree 2 -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -petscpartitioner_type parmetis
1453
1454 test:
1455 suffix: p4est_test_q2_nonconformal_serial
1456 requires: p4est
1457 filter: grep -v "CG or CGNE: variant"
1458 args: -run_type test -petscspace_degree 2 -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1459
1460 test:
1461 suffix: p4est_test_q2_nonconformal_parallel
1462 requires: p4est
1463 filter: grep -v "CG or CGNE: variant"
1464 nsize: 7
1465 args: -run_type test -petscspace_degree 2 -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type simple
1466
1467 test:
1468 suffix: p4est_test_q2_nonconformal_parallel_parmetis
1469 requires: parmetis p4est
1470 nsize: 4
1471 args: -run_type test -petscspace_degree 2 -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type parmetis
1472
1473 test:
1474 suffix: p4est_exact_q2_conformal_serial
1475 requires: p4est !single !complex !__float128
1476 args: -run_type exact -petscspace_degree 2 -fas_levels_snes_atol 1.e-10 -snes_max_it 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type none -fas_coarse_ksp_type preonly -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type none -fas_levels_ksp_type preonly -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2
1477
1478 test:
1479 suffix: p4est_exact_q2_conformal_parallel
1480 requires: p4est !single !complex !__float128
1481 nsize: 4
1482 args: -run_type exact -petscspace_degree 2 -fas_levels_snes_atol 1.e-10 -snes_max_it 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type none -fas_coarse_ksp_type preonly -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type none -fas_levels_ksp_type preonly -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2
1483
1484 test:
1485 suffix: p4est_exact_q2_conformal_parallel_parmetis
1486 requires: parmetis p4est !single
1487 nsize: 4
1488 args: -run_type exact -petscspace_degree 2 -fas_levels_snes_linesearch_type basic -fas_levels_snes_atol 1.e-10 -snes_max_it 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_snes_converged_reason -fas_coarse_pc_type none -fas_coarse_ksp_type preonly -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type none -fas_levels_ksp_type preonly -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -petscpartitioner_type parmetis
1489
1490 test:
1491 suffix: p4est_exact_q2_nonconformal_serial
1492 requires: p4est
1493 args: -run_type exact -petscspace_degree 2 -fas_levels_snes_atol 1.e-10 -snes_max_it 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type none -fas_coarse_ksp_type preonly -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type none -fas_levels_ksp_type preonly -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1494
1495 test:
1496 suffix: p4est_exact_q2_nonconformal_parallel
1497 requires: p4est
1498 nsize: 7
1499 args: -run_type exact -petscspace_degree 2 -fas_levels_snes_atol 1.e-10 -snes_max_it 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type none -fas_coarse_ksp_type preonly -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type none -fas_levels_ksp_type preonly -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type simple
1500
1501 test:
1502 suffix: p4est_exact_q2_nonconformal_parallel_parmetis
1503 requires: parmetis p4est
1504 nsize: 4
1505 args: -run_type exact -petscspace_degree 2 -fas_levels_snes_atol 1.e-10 -snes_max_it 1 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type none -fas_coarse_ksp_type preonly -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type none -fas_levels_ksp_type preonly -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type parmetis
1506
1507 test:
1508 suffix: p4est_full_q2_nonconformal_serial
1509 requires: p4est !single
1510 filter: grep -v "variant HERMITIAN"
1511 args: -run_type full -petscspace_degree 2 -snes_max_it 20 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type jacobi -fas_coarse_ksp_type cg -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type jacobi -fas_levels_ksp_type cg -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1512
1513 test:
1514 suffix: p4est_full_q2_nonconformal_parallel
1515 requires: p4est !single
1516 filter: grep -v "variant HERMITIAN"
1517 nsize: 7
1518 args: -run_type full -petscspace_degree 2 -snes_max_it 20 -snes_type fas -snes_fas_levels 3 -fas_coarse_pc_type jacobi -fas_coarse_ksp_type cg -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type jacobi -fas_levels_ksp_type cg -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type simple
1519
1520 test:
1521 suffix: p4est_full_q2_nonconformal_parallel_bddcfas
1522 requires: p4est !single
1523 filter: grep -v "variant HERMITIAN"
1524 nsize: 7
1525 args: -run_type full -petscspace_degree 2 -snes_max_it 20 -snes_type fas -snes_fas_levels 3 -dm_mat_type is -fas_coarse_pc_type bddc -fas_coarse_ksp_type cg -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type bddc -fas_levels_ksp_type cg -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type simple
1526
1527 test:
1528 suffix: p4est_full_q2_nonconformal_parallel_bddc
1529 requires: p4est !single
1530 filter: grep -v "variant HERMITIAN"
1531 nsize: 7
1532 args: -run_type full -petscspace_degree 2 -snes_max_it 20 -snes_type newtonls -dm_mat_type is -pc_type bddc -ksp_type cg -snes_monitor_short -snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type simple
1533
1534 test:
1535 TODO: broken
1536 suffix: p4est_fas_q2_conformal_serial
1537 requires: p4est !complex !__float128
1538 args: -run_type full -variable_coefficient nonlinear -petscspace_degree 2 -snes_max_it 20 -snes_type fas -snes_fas_levels 3 -pc_type jacobi -ksp_type gmres -fas_coarse_pc_type svd -fas_coarse_ksp_type gmres -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type svd -fas_levels_ksp_type gmres -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_refine_hierarchy 3
1539
1540 test:
1541 TODO: broken
1542 suffix: p4est_fas_q2_nonconformal_serial
1543 requires: p4est
1544 args: -run_type full -variable_coefficient nonlinear -petscspace_degree 2 -snes_max_it 20 -snes_type fas -snes_fas_levels 3 -pc_type jacobi -ksp_type gmres -fas_coarse_pc_type jacobi -fas_coarse_ksp_type gmres -fas_coarse_ksp_monitor_true_residual -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_snes_type newtonls -fas_levels_pc_type jacobi -fas_levels_ksp_type gmres -fas_levels_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1545
1546 test:
1547 suffix: fas_newton_0_p4est
1548 requires: p4est !single !__float128
1549 args: -run_type full -variable_coefficient nonlinear -petscspace_degree 1 -snes_type fas -snes_fas_levels 2 -fas_coarse_pc_type svd -fas_coarse_ksp_rtol 1.0e-10 -fas_coarse_snes_monitor_short -snes_monitor_short -fas_coarse_snes_linesearch_type basic -snes_converged_reason ::ascii_info_detail -snes_view -fas_levels_1_snes_type newtonls -fas_levels_1_pc_type svd -fas_levels_1_ksp_rtol 1.0e-10 -fas_levels_1_snes_monitor_short -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1550
1551 # Full solve simplicial AMR
1552 test:
1553 suffix: tri_p1_adapt_init_pragmatic
1554 requires: pragmatic
1555 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_initial 1 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor pragmatic
1556
1557 test:
1558 suffix: tri_p2_adapt_init_pragmatic
1559 requires: pragmatic
1560 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_initial 1 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor pragmatic
1561
1562 test:
1563 suffix: tri_p1_adapt_init_mmg
1564 requires: mmg
1565 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_initial 1 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor mmg
1566
1567 test:
1568 suffix: tri_p2_adapt_init_mmg
1569 requires: mmg
1570 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_initial 1 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor mmg
1571
1572 test:
1573 suffix: tri_p1_adapt_seq_pragmatic
1574 requires: pragmatic
1575 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_sequence 2 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor pragmatic
1576
1577 test:
1578 suffix: tri_p2_adapt_seq_pragmatic
1579 requires: pragmatic
1580 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_sequence 2 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor pragmatic
1581
1582 test:
1583 suffix: tri_p1_adapt_seq_mmg
1584 requires: mmg
1585 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_sequence 2 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor mmg
1586
1587 test:
1588 suffix: tri_p2_adapt_seq_mmg
1589 requires: mmg
1590 args: -run_type exact -dm_refine 5 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -snes_adapt_sequence 2 -adaptor_target_num 4000 -dm_plex_metric_h_max 0.5 -dm_adaptor mmg
1591
1592 test:
1593 suffix: tri_p1_adapt_analytic_pragmatic
1594 requires: pragmatic
1595 args: -run_type exact -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient cross -snes_adapt_initial 4 -adaptor_target_num 500 -dm_plex_metric_h_min 0.0001 -dm_plex_metric_h_max 0.05 -dm_adaptor pragmatic
1596 output_file: output/empty.out
1597
1598 test:
1599 suffix: tri_p2_adapt_analytic_pragmatic
1600 requires: pragmatic
1601 args: -run_type exact -dm_refine 3 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient cross -snes_adapt_initial 4 -adaptor_target_num 500 -dm_plex_metric_h_min 0.0001 -dm_plex_metric_h_max 0.05 -dm_adaptor pragmatic
1602 output_file: output/empty.out
1603
1604 test:
1605 suffix: tri_p1_adapt_analytic_mmg
1606 requires: mmg
1607 args: -run_type exact -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient cross -snes_adapt_initial 4 -adaptor_target_num 500 -dm_plex_metric_h_max 0.5 -dm_adaptor mmg
1608 output_file: output/empty.out
1609
1610 test:
1611 suffix: tri_p2_adapt_analytic_mmg
1612 requires: mmg
1613 args: -run_type exact -dm_refine 3 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient cross -snes_adapt_initial 4 -adaptor_target_num 500 -dm_plex_metric_h_max 0.5 -dm_adaptor mmg
1614 output_file: output/empty.out
1615
1616 test:
1617 suffix: tri_p1_adapt_uniform_pragmatic
1618 requires: pragmatic tetgen
1619 nsize: 2
1620 args: -run_type full -dm_plex_box_faces 8,8,8 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient none -snes_converged_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence 3 -adaptor_target_num 400 -dm_plex_metric_h_max 0.5 -dm_plex_dim 3 -dm_adaptor pragmatic
1621 timeoutfactor: 2
1622
1623 test:
1624 suffix: tri_p2_adapt_uniform_pragmatic
1625 requires: pragmatic tetgen
1626 nsize: 2
1627 args: -run_type full -dm_plex_box_faces 8,8,8 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient none -snes_converged_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence 1 -adaptor_target_num 400 -dm_plex_metric_h_max 0.5 -dm_plex_dim 3 -dm_adaptor pragmatic
1628 timeoutfactor: 1
1629
1630 test:
1631 suffix: tri_p1_adapt_uniform_mmg
1632 requires: mmg tetgen
1633 args: -run_type full -dm_plex_box_faces 4,4,4 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient none -snes_converged_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence 3 -adaptor_target_num 400 -dm_plex_metric_h_max 0.5 -dm_plex_dim 3 -dm_adaptor mmg
1634 timeoutfactor: 2
1635
1636 test:
1637 suffix: tri_p2_adapt_uniform_mmg
1638 requires: mmg tetgen
1639 TODO: broken
1640 args: -run_type full -dm_plex_box_faces 4,4,4 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient none -snes_converged_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence 1 -adaptor_target_num 400 -dm_plex_metric_h_max 0.5 -dm_plex_dim 3 -dm_adaptor mmg
1641 timeoutfactor: 1
1642
1643 test:
1644 suffix: tri_p1_adapt_uniform_parmmg
1645 requires: parmmg tetgen
1646 nsize: 2
1647 args: -run_type full -dm_plex_box_faces 8,8,8 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient none -snes_converged_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence 3 -adaptor_target_num 400 -dm_plex_metric_h_max 0.5 -dm_plex_dim 3 -dm_adaptor parmmg
1648 timeoutfactor: 2
1649
1650 test:
1651 suffix: tri_p2_adapt_uniform_parmmg
1652 requires: parmmg tetgen
1653 nsize: 2
1654 args: -run_type full -dm_plex_box_faces 8,8,8 -bc_type dirichlet -petscspace_degree 2 -variable_coefficient none -snes_converged_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence 1 -adaptor_target_num 400 -dm_plex_metric_h_max 0.5 -dm_plex_dim 3 -dm_adaptor parmmg
1655 timeoutfactor: 1
1656
1657 # Full solve tensor AMR
1658 test:
1659 suffix: quad_q1_adapt_0
1660 requires: p4est
1661 args: -run_type exact -dm_plex_simplex 0 -dm_plex_convert_type p4est -bc_type dirichlet -petscspace_degree 1 -variable_coefficient ball -snes_converged_reason ::ascii_info_detail -pc_type lu -dm_forest_initial_refinement 4 -snes_adapt_initial 1 -dm_view
1662 filter: grep -v DM_
1663
1664 test:
1665 suffix: amr_0
1666 nsize: 5
1667 args: -run_type test -petscpartitioner_type simple -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -dm_refine 1
1668
1669 test:
1670 suffix: amr_1
1671 requires: p4est !complex
1672 args: -run_type test -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -dm_plex_convert_type p4est -dm_p4est_refine_pattern center -dm_forest_maximum_refinement 5 -dm_view vtk:amr.vtu:vtk_vtu -vec_view vtk:amr.vtu:vtk_vtu:append
1673
1674 test:
1675 suffix: p4est_solve_bddc
1676 requires: p4est !complex
1677 args: -run_type full -variable_coefficient nonlinear -nonzero_initial_guess 1 -petscspace_degree 2 -snes_max_it 20 -snes_type newtonls -dm_mat_type is -pc_type bddc -ksp_type cg -snes_monitor_short -ksp_monitor -snes_linesearch_type bt -snes_converged_reason -snes_view -dm_plex_simplex 0 -petscspace_poly_tensor -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -petscpartitioner_type simple -pc_bddc_detect_disconnected
1678 nsize: 4
1679
1680 test:
1681 suffix: p4est_solve_fas
1682 requires: p4est
1683 args: -run_type full -variable_coefficient nonlinear -nonzero_initial_guess 1 -petscspace_degree 2 -snes_max_it 10 -snes_type fas -snes_linesearch_type bt -snes_fas_levels 3 -fas_coarse_snes_type newtonls -fas_coarse_snes_linesearch_type basic -fas_coarse_ksp_type cg -fas_coarse_pc_type jacobi -fas_coarse_snes_monitor_short -fas_levels_snes_max_it 4 -fas_levels_snes_type newtonls -fas_levels_snes_linesearch_type bt -fas_levels_ksp_type cg -fas_levels_pc_type jacobi -fas_levels_snes_monitor_short -fas_levels_cycle_snes_linesearch_type bt -snes_monitor_short -snes_converged_reason -snes_view -dm_plex_simplex 0 -petscspace_poly_tensor -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1684 nsize: 4
1685 TODO: identical machine two runs produce slightly different solver trackers
1686
1687 test:
1688 suffix: p4est_convergence_test_1
1689 requires: p4est
1690 args: -quiet -run_type test -petscspace_degree 1 -dm_plex_simplex 0 -petscspace_poly_tensor -dm_plex_convert_type p4est -dm_forest_minimum_refinement 2 -dm_forest_initial_refinement 2 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash
1691 nsize: 4
1692
1693 # Serial tests with GLVis visualization
1694 test:
1695 suffix: glvis_2d_tet_p1
1696 args: -quiet -run_type test -bc_type dirichlet -petscspace_degree 1 -vec_view glvis: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_plex_boundary_label marker -dm_plex_gmsh_periodic 0 -dm_coord_space 0
1697 test:
1698 suffix: glvis_2d_tet_p2
1699 args: -quiet -run_type test -bc_type dirichlet -petscspace_degree 2 -vec_view glvis: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_plex_boundary_label marker -dm_plex_gmsh_periodic 0 -dm_coord_space 0
1700 test:
1701 suffix: glvis_2d_hex_p1
1702 args: -quiet -run_type test -bc_type dirichlet -petscspace_degree 1 -vec_view glvis: -dm_plex_simplex 0 -dm_refine 1 -dm_coord_space 0
1703 test:
1704 suffix: glvis_2d_hex_p2
1705 args: -quiet -run_type test -bc_type dirichlet -petscspace_degree 2 -vec_view glvis: -dm_plex_simplex 0 -dm_refine 1 -dm_coord_space 0
1706 test:
1707 suffix: glvis_2d_hex_p2_p4est
1708 requires: p4est
1709 args: -quiet -run_type test -bc_type dirichlet -petscspace_degree 2 -vec_view glvis: -dm_plex_simplex 0 -dm_plex_convert_type p4est -dm_forest_minimum_refinement 0 -dm_forest_initial_refinement 1 -dm_forest_maximum_refinement 4 -dm_p4est_refine_pattern hash -viewer_glvis_dm_plex_enable_ncmesh
1710 test:
1711 suffix: glvis_2d_tet_p0
1712 args: -run_type exact -guess_vec_view glvis: -nonzero_initial_guess 1 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_plex_boundary_label marker -petscspace_degree 0 -dm_coord_space 0 -pc_type jacobi
1713 test:
1714 suffix: glvis_2d_hex_p0
1715 args: -run_type exact -guess_vec_view glvis: -nonzero_initial_guess 1 -dm_plex_box_faces 5,7 -dm_plex_simplex 0 -petscspace_degree 0 -dm_coord_space 0 -pc_type jacobi
1716
1717 # PCHPDDM tests
1718 testset:
1719 nsize: 4
1720 requires: hpddm slepc !single defined(PETSC_HAVE_DYNAMIC_LIBRARIES) defined(PETSC_USE_SHARED_LIBRARIES)
1721 args: -run_type test -run_test_check_ksp -quiet -petscspace_degree 1 -petscpartitioner_type simple -bc_type none -dm_plex_simplex 0 -pc_type hpddm -pc_hpddm_levels_1_sub_pc_type lu -pc_hpddm_levels_1_eps_nev 2 -pc_hpddm_coarse_p 1 -pc_hpddm_coarse_pc_type svd -ksp_rtol 1.e-10 -pc_hpddm_levels_1_st_pc_factor_shift_type INBLOCKS -ksp_converged_reason
1722 test:
1723 suffix: quad_singular_hpddm
1724 args: -dm_plex_box_faces 6,7
1725 test:
1726 requires: p4est
1727 suffix: p4est_singular_2d_hpddm
1728 args: -dm_plex_convert_type p4est -dm_forest_minimum_refinement 1 -dm_forest_initial_refinement 3 -dm_forest_maximum_refinement 3
1729 test:
1730 requires: p4est
1731 suffix: p4est_nc_singular_2d_hpddm
1732 args: -dm_plex_convert_type p4est -dm_forest_minimum_refinement 1 -dm_forest_initial_refinement 1 -dm_forest_maximum_refinement 3 -dm_p4est_refine_pattern hash
1733 testset:
1734 nsize: 4
1735 requires: hpddm slepc triangle !single defined(PETSC_HAVE_DYNAMIC_LIBRARIES) defined(PETSC_USE_SHARED_LIBRARIES)
1736 args: -run_type full -petscpartitioner_type simple -dm_refine 2 -bc_type dirichlet -petscspace_degree 2 -ksp_type gmres -ksp_gmres_restart 100 -pc_type hpddm -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0 -pc_type hpddm -pc_hpddm_levels_1_sub_pc_type lu -pc_hpddm_levels_1_eps_nev 4 -pc_hpddm_coarse_p 2 -pc_hpddm_coarse_pc_type redundant -ksp_rtol 1.e-1
1737 test:
1738 args: -pc_hpddm_coarse_mat_type baij -options_left no
1739 suffix: tri_hpddm_reuse_baij
1740 test:
1741 requires: !complex
1742 suffix: tri_hpddm_reuse
1743 testset:
1744 nsize: 4
1745 requires: hpddm slepc !single defined(PETSC_HAVE_DYNAMIC_LIBRARIES) defined(PETSC_USE_SHARED_LIBRARIES)
1746 args: -run_type full -petscpartitioner_type simple -dm_plex_box_faces 7,5 -dm_refine 2 -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 2 -ksp_type gmres -ksp_gmres_restart 100 -pc_type hpddm -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0 -pc_type hpddm -pc_hpddm_levels_1_sub_pc_type lu -pc_hpddm_levels_1_eps_nev 4 -pc_hpddm_coarse_p 2 -pc_hpddm_coarse_pc_type redundant -ksp_rtol 1.e-1
1747 test:
1748 args: -pc_hpddm_coarse_mat_type baij -options_left no
1749 suffix: quad_hpddm_reuse_baij
1750 test:
1751 requires: !complex
1752 suffix: quad_hpddm_reuse
1753 testset:
1754 nsize: 4
1755 requires: hpddm slepc !single defined(PETSC_HAVE_DYNAMIC_LIBRARIES) defined(PETSC_USE_SHARED_LIBRARIES)
1756 args: -run_type full -petscpartitioner_type simple -dm_plex_box_faces 7,5 -dm_refine 2 -dm_plex_simplex 0 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -pc_type hpddm -snes_monitor_short -ksp_monitor_short -snes_converged_reason ::ascii_info_detail -ksp_converged_reason -snes_view -show_solution 0 -pc_type hpddm -pc_hpddm_levels_1_sub_pc_type lu -pc_hpddm_levels_1_eps_threshold_absolute 0.1 -pc_hpddm_coarse_p 2 -pc_hpddm_coarse_pc_type redundant -ksp_rtol 1.e-1
1757 test:
1758 args: -pc_hpddm_coarse_mat_type baij -options_left no
1759 suffix: quad_hpddm_reuse_threshold_baij
1760 test:
1761 requires: !complex
1762 suffix: quad_hpddm_reuse_threshold
1763 testset:
1764 nsize: 4
1765 requires: hpddm slepc parmetis !single defined(PETSC_HAVE_DYNAMIC_LIBRARIES) defined(PETSC_USE_SHARED_LIBRARIES)
1766 filter: sed -e "s/linear solver iterations=17/linear solver iterations=16/g"
1767 args: -run_type full -petscpartitioner_type parmetis -dm_refine 3 -bc_type dirichlet -petscspace_degree 1 -ksp_type gmres -ksp_gmres_restart 100 -pc_type hpddm -snes_monitor_short -snes_converged_reason ::ascii_info_detail -snes_view -show_solution 0 -pc_type hpddm -pc_hpddm_levels_1_sub_pc_type icc -pc_hpddm_levels_1_eps_nev 20 -pc_hpddm_coarse_p 2 -pc_hpddm_coarse_pc_type redundant -ksp_rtol 1.e-10 -dm_plex_filename ${PETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_plex_boundary_label marker -pc_hpddm_levels_1_sub_pc_factor_levels 3 -variable_coefficient ball -dm_plex_gmsh_periodic 0 -fp_trap 0
1768 test:
1769 args: -pc_hpddm_coarse_mat_type baij -options_left no
1770 filter: grep -v " total: nonzeros=" | grep -v " rows=" | sed -e "s/total number of linear solver iterations=[1-2][4-7]/total number of linear solver iterations=16/g"
1771 suffix: tri_parmetis_hpddm_baij
1772 test:
1773 filter: grep -v " total: nonzeros=" | grep -v " rows=" | sed -e "s/total number of linear solver iterations=[1-2][4-7]/total number of linear solver iterations=16/g"
1774 requires: !complex
1775 suffix: tri_parmetis_hpddm
1776
1777 # 2D serial P1 tests for adaptive MG
1778 test:
1779 suffix: 2d_p1_adaptmg_0
1780 requires: triangle
1781 args: -petscpartitioner_type simple -dm_refine_hierarchy 3 -dm_plex_box_faces 4,4 -bc_type dirichlet -petscspace_degree 1 \
1782 -variable_coefficient checkerboard_0 -mat_petscspace_degree 0 -div 16 -k 3 \
1783 -snes_max_it 1 -ksp_converged_reason \
1784 -ksp_rtol 1e-8 -pc_type mg
1785 test:
1786 suffix: 2d_p1_adaptmg_1
1787 requires: triangle bamg
1788 args: -petscpartitioner_type simple -dm_refine_hierarchy 3 -dm_plex_box_faces 4,4 -bc_type dirichlet -petscspace_degree 1 \
1789 -variable_coefficient checkerboard_0 -mat_petscspace_degree 0 -div 16 -k 3 \
1790 -snes_max_it 1 -ksp_converged_reason \
1791 -ksp_rtol 1e-8 -pc_type mg -pc_mg_galerkin -pc_mg_adapt_interp_coarse_space eigenvector -pc_mg_adapt_interp_n 1 \
1792 -pc_mg_mesp_ksp_type richardson -pc_mg_mesp_ksp_richardson_self_scale -pc_mg_mesp_ksp_max_it 100 -pc_mg_mesp_pc_type none
1793 test:
1794 suffix: 2d_p1_adaptmg_gdsw
1795 requires: triangle
1796 nsize: 4
1797 args: -petscpartitioner_type simple -dm_refine 3 -dm_plex_box_faces 4,4 -bc_type dirichlet -petscspace_degree 1 \
1798 -variable_coefficient checkerboard_0 -mat_petscspace_degree 0 -div 16 -k 3 \
1799 -snes_max_it 1 -ksp_converged_reason \
1800 -ksp_rtol 1e-8 -pc_type mg -pc_mg_galerkin -pc_mg_adapt_interp_coarse_space gdsw -pc_mg_levels 2 -mg_levels_pc_type asm -dm_mat_type {{aij is}}
1801
1802 test:
1803 suffix: 2d_p1_adaptmg_agdsw
1804 requires: triangle mumps
1805 nsize: 4
1806 args: -petscpartitioner_type simple -dm_refine 3 -dm_plex_box_faces 4,4 -bc_type dirichlet -petscspace_degree 1 \
1807 -variable_coefficient checkerboard_0 -mat_petscspace_degree 0 -div 16 -k 3 \
1808 -snes_max_it 1 -ksp_converged_reason \
1809 -ksp_rtol 1e-8 -pc_type mg -pc_mg_galerkin -pc_mg_adapt_interp_coarse_space gdsw -pc_mg_levels 2 -mg_levels_pc_type asm -dm_mat_type is -mg_levels_gdsw_tolerance 0.1 -mg_levels_gdsw_pseudo_pc_type qr
1810
1811 test:
1812 suffix: p4est_2d_asm
1813 requires: p4est
1814 nsize: 4
1815 args: -run_type test -run_test_check_ksp -quiet -petscspace_degree 1 -petscpartitioner_type simple -bc_type none -dm_plex_simplex 0 \
1816 -pc_type asm -ksp_converged_reason -dm_plex_convert_type p4est -dm_forest_minimum_refinement 1 -dm_forest_initial_refinement 3 -dm_forest_maximum_refinement 5 \
1817 -pc_asm_dm_subdomains -dm_p4est_refine_pattern hash -dm_plex_dd_overlap 1 -sub_pc_type lu
1818
1819 TEST*/
1820