Lines Matching refs:param
68 PetscParam param; member
98 appctx.param.N = 10; /* order of the spectral element */ in main()
99 appctx.param.E = 10; /* number of elements */ in main()
100 appctx.param.L = 4.0; /* length of the domain */ in main()
101 appctx.param.mu = 0.01; /* diffusion coefficient */ in main()
103 appctx.param.steps = PETSC_INT_MAX; in main()
104 appctx.param.Tend = 4; in main()
106 PetscCall(PetscOptionsGetInt(NULL, NULL, "-N", &appctx.param.N, NULL)); in main()
107 PetscCall(PetscOptionsGetInt(NULL, NULL, "-E", &appctx.param.E, NULL)); in main()
108 PetscCall(PetscOptionsGetReal(NULL, NULL, "-Tend", &appctx.param.Tend, NULL)); in main()
109 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &appctx.param.mu, NULL)); in main()
110 appctx.param.Le = appctx.param.L / appctx.param.E; in main()
113 …PetscCheck((appctx.param.E % size) == 0, PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONG, "Number of element… in main()
118 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.… in main()
119 …PetscCall(PetscDTGaussLobattoLegendreQuadrature(appctx.param.N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINE… in main()
120 appctx.SEMop.gll.n = appctx.param.N; in main()
121 lenglob = appctx.param.E * (appctx.param.N - 1); in main()
149 xs = xs / (appctx.param.N - 1); in main()
150 xm = xm / (appctx.param.N - 1); in main()
157 for (j = 0; j < appctx.param.N - 1; j++) { in main()
158 … x = (appctx.param.Le / 2.0) * (appctx.SEMop.gll.nodes[j] + 1.0) + appctx.param.Le * i; in main()
159 ind = i * (appctx.param.N - 1) + j; in main()
161 wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
162 if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j]; in main()
208 PetscCall(TSSetMaxSteps(appctx.ts, appctx.param.steps)); in main()
209 PetscCall(TSSetMaxTime(appctx.ts, appctx.param.Tend)); in main()
267 …param.mu * PETSC_PI * PetscSinScalar(PETSC_PI * xg[i]) * PetscExpReal(-appctx->param.mu * PETSC_PI… in TrueSolution()
340 for (i = 0; i < appctx->param.N; i++) { in MatMult_Laplacian()
341 vv = -appctx->param.mu * 2.0 / appctx->param.Le; in MatMult_Laplacian()
342 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in MatMult_Laplacian()
347 PetscCall(PetscBLASIntCast(appctx->param.N, &n)); in MatMult_Laplacian()
348 …for (j = xs; j < xs + xn; j += appctx->param.N - 1) PetscCallBLAS("BLASgemv", BLASgemv_("N", &n, &… in MatMult_Laplacian()
383 PetscCall(PetscBLASIntCast(appctx->param.N, &n)); in MatMult_Advection()
384 …for (j = xs; j < xs + xn; j += appctx->param.N - 1) PetscCallBLAS("BLASgemv", BLASgemv_("N", &n, &… in MatMult_Advection()
431 …PetscCheck(appctx->param.N > 1, PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONG, "Spectral element order sho… in RHSMatrixLaplaciangllDM()
434 for (i = 0; i < appctx->param.N; i++) { in RHSMatrixLaplaciangllDM()
435 vv = -appctx->param.mu * 2.0 / appctx->param.Le; in RHSMatrixLaplaciangllDM()
436 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in RHSMatrixLaplaciangllDM()
442 xs = xs / (appctx->param.N - 1); in RHSMatrixLaplaciangllDM()
443 xn = xn / (appctx->param.N - 1); in RHSMatrixLaplaciangllDM()
445 PetscCall(PetscMalloc1(appctx->param.N, &rowsDM)); in RHSMatrixLaplaciangllDM()
450 for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l; in RHSMatrixLaplaciangllDM()
451 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_… in RHSMatrixLaplaciangllDM()
503 xs = xs / (appctx->param.N - 1); in RHSMatrixAdvectiongllDM()
504 xn = xn / (appctx->param.N - 1); in RHSMatrixAdvectiongllDM()
506 PetscCall(PetscMalloc1(appctx->param.N, &rowsDM)); in RHSMatrixAdvectiongllDM()
508 for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l; in RHSMatrixAdvectiongllDM()
509 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_… in RHSMatrixAdvectiongllDM()