Lines Matching refs:param
84 PetscParam param; member
122 appctx.param.N = 10; /* order of the spectral element */ in main()
123 appctx.param.E = 8; /* number of elements */ in main()
124 appctx.param.L = 1.0; /* length of the domain */ in main()
125 appctx.param.mu = 0.00001; /* diffusion coefficient */ in main()
126 appctx.param.a = 0.0; /* advection speed */ in main()
128 appctx.param.steps = PETSC_INT_MAX; in main()
129 appctx.param.Tend = 0.01; in main()
132 PetscCall(PetscOptionsGetInt(NULL, NULL, "-N", &appctx.param.N, NULL)); in main()
133 PetscCall(PetscOptionsGetInt(NULL, NULL, "-E", &appctx.param.E, NULL)); in main()
135 PetscCall(PetscOptionsGetReal(NULL, NULL, "-Tend", &appctx.param.Tend, NULL)); in main()
136 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &appctx.param.mu, NULL)); in main()
137 PetscCall(PetscOptionsGetReal(NULL, NULL, "-a", &appctx.param.a, NULL)); in main()
138 appctx.param.Le = appctx.param.L / appctx.param.E; in main()
143 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.… in main()
144 …PetscCall(PetscDTGaussLobattoLegendreQuadrature(appctx.param.N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINE… in main()
145 appctx.SEMop.gll.n = appctx.param.N; in main()
146 lenglob = appctx.param.E * (appctx.param.N - 1); in main()
179 xs = xs / (appctx.param.N - 1); in main()
180 xm = xm / (appctx.param.N - 1); in main()
187 for (j = 0; j < appctx.param.N - 1; j++) { in main()
188 … x = (appctx.param.Le / 2.0) * (appctx.SEMop.gll.nodes[j] + 1.0) + appctx.param.Le * i; in main()
189 ind = i * (appctx.param.N - 1) + j; in main()
191 wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
192 if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
229 PetscCall(TSSetMaxSteps(appctx.ts, appctx.param.steps)); in main()
230 PetscCall(TSSetMaxTime(appctx.ts, appctx.param.Tend)); in main()
244 PetscCall(ComputeReference(appctx.ts, appctx.param.Tend, appctx.dat.reference, &appctx)); in main()
331 lenglob = appctx->param.E * (appctx->param.N - 1); in InitialConditions()
370 lenglob = appctx->param.E * (appctx->param.N - 1); in TrueSolution()
401 lenglob = appctx->param.E * (appctx->param.N - 1); in ComputeReference()
405 tc = -appctx->param.mu * (j + 1) * (j + 1) * 4.0 * PETSC_PI * PETSC_PI * t; in ComputeReference()
406 …oefficients[j] * PetscSinScalar(2 * (j + 1) * PETSC_PI * (xg[i] + appctx->param.a * t)) * PetscExp… in ComputeReference()
465 for (i = 0; i < appctx->param.N; i++) { in RHSLaplacian()
466 vv = -appctx->param.mu * 2.0 / appctx->param.Le; in RHSLaplacian()
467 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in RHSLaplacian()
473 …PetscCheck(appctx->param.N - 1 >= 1, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Polynomial order… in RHSLaplacian()
474 xs = xs / (appctx->param.N - 1); in RHSLaplacian()
475 xn = xn / (appctx->param.N - 1); in RHSLaplacian()
477 PetscCall(PetscMalloc1(appctx->param.N, &rowsDM)); in RHSLaplacian()
482 for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l; in RHSLaplacian()
483 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_… in RHSLaplacian()
516 for (i = 0; i < appctx->param.N; i++) { in RHSAdvection()
517 vv = -appctx->param.a; in RHSAdvection()
518 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in RHSAdvection()
524 …PetscCheck(appctx->param.N - 1 >= 1, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Polynomial order… in RHSAdvection()
525 xs = xs / (appctx->param.N - 1); in RHSAdvection()
526 xn = xn / (appctx->param.N - 1); in RHSAdvection()
528 PetscCall(PetscMalloc1(appctx->param.N, &rowsDM)); in RHSAdvection()
533 for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l; in RHSAdvection()
534 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_… in RHSAdvection()