Home
last modified time | relevance | path

Searched refs:N (Results 1 – 22 of 22) sorted by relevance

/honee/qfunctions/
H A Dutils_eigensolver_jacobi.h30 CEED_QFUNCTION_HELPER CeedInt MaxEntryRow(const CeedScalar *A, CeedInt N, CeedInt i) { in MaxEntryRow() argument
32 for (CeedInt j = i + 2; j < N; j++) in MaxEntryRow()
33 if (fabs(A[i * N + j]) > fabs(A[i * N + j_max])) j_max = j; in MaxEntryRow()
44 CEED_QFUNCTION_HELPER void MaxEntry(const CeedScalar *A, CeedInt N, CeedInt *max_idx_row, CeedInt *… in MaxEntry() argument
47 CeedScalar max_entry = fabs(A[*i_max * N + *j_max]); in MaxEntry()
48 for (CeedInt i = 1; i < N - 1; i++) { in MaxEntry()
50 if (fabs(A[i * N + j]) > max_entry) { in MaxEntry()
51 max_entry = fabs(A[i * N + j]); in MaxEntry()
67 CEED_QFUNCTION_HELPER void CalcRot(const CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedScalar… in CalcRot() argument
69 CeedScalar A_jj_ii = (A[j * N + j] - A[i * N + i]); in CalcRot()
[all …]
H A Dutils.h27 CEED_QFUNCTION_HELPER void ScaleN(CeedScalar *u, const CeedScalar alpha, const CeedInt N) { in ScaleN() argument
28 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] *= alpha; in ScaleN()
32 CEED_QFUNCTION_HELPER void SetValueN(CeedScalar *u, const CeedScalar alpha, const CeedInt N) { in SetValueN() argument
33 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] = alpha; in SetValueN()
37 …pyN(const CeedScalar *x, CeedScalar *y, const CeedInt N) { CeedPragmaSIMD for (CeedInt i = 0; i < … in CopyN() argument
43 CEED_QFUNCTION_HELPER CeedScalar DotN(const CeedScalar *u, const CeedScalar *v, const CeedInt N) { in DotN() argument
45 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) output += u[i] * v[i]; in DotN()
50 CEED_QFUNCTION_HELPER void AXPY(CeedScalar alpha, const CeedScalar *x, CeedScalar *y, CeedInt N) { in AXPY() argument
51 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) y[i] += alpha * x[i]; in AXPY()
81 CEED_QFUNCTION_HELPER void MatVecNM(const CeedScalar *A, const CeedScalar *x, const CeedInt N, cons… in MatVecNM() argument
[all …]
H A Ddensitycurrent.h21 CeedScalar N; member
85 const CeedScalar N = context->N; in Exact_DC() local
107 const CeedScalar theta = theta0 * exp(Square(N) * z / g) + delta_theta; in Exact_DC()
110 const CeedScalar Pi = 1. + Square(g) * (exp(-Square(N) * z / g) - 1.) / (cp * theta0 * Square(N)); in Exact_DC()
H A Ddifferential_filter.h74 …LHS_N(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out, const CeedInt N) { in DifferentialFilter_LHS_N() argument
86 CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) { in DifferentialFilter_LHS_N()
129 …dq[k] = Grad_q[0 * N + j][i] * dXdx[0][k] + Grad_q[1 * N + j][i] * dXdx[1][k] + Grad_q[2 * N + j][… in DifferentialFilter_LHS_N()
134 Grad_v[k * N + j][i] = wdetJ * dq_dXdx_a[k]; in DifferentialFilter_LHS_N()
H A Dblasius.h33 CEED_QFUNCTION_HELPER void ChebyshevEval(int N, const double *Tf, double x, double eta_max, double … in ChebyshevEval() argument
46 for (int i = 3; i < N; i++) { in ChebyshevEval()
69 CeedInt N = blasius->n_cheb; in BlasiusSolution() local
79 ChebyshevEval(N, blasius->Tf_cheb, X, blasius->eta_max, f); in BlasiusSolution()
80 ChebyshevEval(N - 1, blasius->Th_cheb, X, blasius->eta_max, h); in BlasiusSolution()
H A Dmass.h19 …ass_N(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out, const CeedInt N) { in Mass_N() argument
24 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) v[j]… in Mass_N()
H A Dadvection.h210 CEED_QFUNCTION_HELPER void StatePhysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, N… in StatePhysicalGradientFromReference_ND() argument
213 switch (N) { in StatePhysicalGradientFromReference_ND()
218 …i[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[… in StatePhysicalGradientFromReference_ND()
230 …i[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[… in StatePhysicalGradientFromReference_ND()
231 grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k]; in StatePhysicalGradientFromReference_ND()
H A Dsgs_dd_model.h33 CEED_QFUNCTION_HELPER void LeakyReLU(CeedScalar *x, const CeedScalar alpha, const CeedInt N) { in LeakyReLU() argument
34 for (CeedInt i = 0; i < N; i++) x[i] *= (x[i] < 0 ? alpha : 1.); in LeakyReLU()
H A Dshocktube.h21 CeedScalar N; member
/honee/problems/
H A Dblasius.c21 PetscInt N = blasius->n_cheb; in CompressibleBlasiusResidual() local
30 Th = Tf + N; in CompressibleBlasiusResidual()
34 ChebyshevEval(N, Tf, -1., blasius->eta_max, f); in CompressibleBlasiusResidual()
39 ChebyshevEval(N, Tf, 1., blasius->eta_max, f); in CompressibleBlasiusResidual()
42 for (int i = 0; i < N - 3; i++) { in CompressibleBlasiusResidual()
43 ChebyshevEval(N, Tf, blasius->X[i], blasius->eta_max, f); in CompressibleBlasiusResidual()
44 ChebyshevEval(N - 1, Th, blasius->X[i], blasius->eta_max, h); in CompressibleBlasiusResidual()
56 …r[N + 2 + i] = (mu_rho_tilde[0] * h[2] + mu_rho_tilde[1] * h[1]) + Pr * f[0] * h[1] + Pr * (gamma … in CompressibleBlasiusResidual()
60 ChebyshevEval(N - 1, Th, -1., blasius->eta_max, h); in CompressibleBlasiusResidual()
61 r[N] = h[0] - blasius->T_wall / S_infty.Y.temperature; in CompressibleBlasiusResidual()
[all …]
H A Ddensitycurrent.c33 CeedScalar N = 0.01; // 1/s in NS_DENSITY_CURRENT() local
48 PetscCall(PetscOptionsScalar("-N", "Brunt-Vaisala frequency", NULL, N, &N, NULL)); in NS_DENSITY_CURRENT()
73 N *= (1. / units->second); in NS_DENSITY_CURRENT()
79 density_current_ctx->N = N; in NS_DENSITY_CURRENT()
H A Dsgs_dd_model.c583 static PetscErrorCode TransposeMatrix(const PetscScalar *A, PetscScalar *B, const PetscInt N, const… in TransposeMatrix() argument
585 for (PetscInt i = 0; i < N; i++) { in TransposeMatrix()
587 B[j * N + i] = A[i * M + j]; in TransposeMatrix()
/honee/examples/meshes/
H A DMakefile9 cylinder-q1-n08.msh: GMSH_FLAGS = -order 1 -setnumber N 08 -setnumber Rb 0.7
10 cylinder-q1-n12.msh: GMSH_FLAGS = -order 1 -setnumber N 12 -setnumber Rb 0.8
11 cylinder-q1-n20.msh: GMSH_FLAGS = -order 1 -setnumber N 20
12 cylinder-q1-n40.msh: GMSH_FLAGS = -order 1 -setnumber N 40
13 cylinder-q2-n20.msh: GMSH_FLAGS = -order 2 -setnumber N 20
14 cylinder-q2-n40.msh: GMSH_FLAGS = -order 2 -setnumber N 40
15 cylinder-q1-n12-recirc.msh: GMSH_FLAGS = -order 1 -setnumber N 12 -setnumber H 1 -setnumber xL 1 -s…
H A Dcylinder.geo27 N = {20, Min 2, Max 80, Step 2,
40 Nx1 = N + 1; Rx1 = 1.00;
41 Nx2 = Ceil(N * 9 / 11 * xR / (2*xL)) + 1; Rx2 = 1.00;
42 Ny = N + 1; Ry = 2.00;
43 Nb = Ceil(RadialRatio * N) + 1;
44 Nc = N + 1; Rc = 1.00;
47 h = 2 * H / N;
/honee/src/
H A Dhonee-mass.c19 PetscErrorCode HoneeMassQFunctionCreate(Ceed ceed, CeedInt N, CeedInt q_data_size, CeedQFunction *q… in HoneeMassQFunctionCreate() argument
21 switch (N) { in HoneeMassQFunctionCreate()
54 …FMT " does not exist.\nA new mass qfunction can be easily added; see source code for pattern.", N); in HoneeMassQFunctionCreate()
57 PetscCallCeed(ceed, CeedQFunctionAddInput(*qf, "u", N, CEED_EVAL_INTERP)); in HoneeMassQFunctionCreate()
59 PetscCallCeed(ceed, CeedQFunctionAddOutput(*qf, "v", N, CEED_EVAL_INTERP)); in HoneeMassQFunctionCreate()
60 PetscCallCeed(ceed, CeedQFunctionSetUserFlopsEstimate(*qf, N)); in HoneeMassQFunctionCreate()
H A Dhonee-meshtransform.c59 …neeMeshTransform_PlateMesh(MPI_Comm comm, DM dm, PetscReal growth, PetscInt N, PetscReal refine_he… in HoneeMeshTransform_PlateMesh() argument
89 PetscReal dy1 = refine_height * (growth - 1) / (pow(growth, N) - 1); in HoneeMeshTransform_PlateMesh()
92 PetscReal logdy = (log(domain_max[1]) - log(refine_height)) / (faces[1] - N); in HoneeMeshTransform_PlateMesh()
100 if (y_box_index <= N) { in HoneeMeshTransform_PlateMesh()
104 PetscInt j = y_box_index - N; in HoneeMeshTransform_PlateMesh()
H A Dhonee-file.c148 PetscInt length, N; in HoneeLoadBinaryVec() local
150 PetscCall(VecGetSize(Q, &N)); in HoneeLoadBinaryVec()
162 …PetscCheck(length == N, comm, PETSC_ERR_ARG_INCOMP, "File Vec '%s' has length %" PetscInt_FMT " bu… in HoneeLoadBinaryVec()
163 filename, length, N); in HoneeLoadBinaryVec()
/honee/doc/
H A Dtheory.md525 Here, we define the number of wavemodes $N$, set of random numbers $ \{\bm{\sigma}^n, \bm{d}^n, \ph…
537 \kappa^n = \kappa_{\min} (1 + \alpha)^{n-1} \ , \quad \forall n=1, 2, ... , N
585 N;
606 rand --> N --Calc--> kn;
622 …and.dat` file is the table of the random number set, $\{\bm{\sigma}^n, \bm{d}^n, \phi^n\}_{n=1}^N$.
633 | $ \{\bm{\sigma}^n, \bm{d}^n, \phi^n\}_{n=1}^N$ | RN Set | No | Yes |
641 | $\{\kappa^n\}_{n=1}^N$ | k^n | No | Yes |
H A Dreferences.bib54 author = {Papanastasiou, TC and Malamataris, N and Ellwood K},
H A Druntime_options.md445 - N/A
H A Dexamples.md168 * - `-N`
/honee/include/
H A Dnavierstokes.h318 PetscErrorCode HoneeMassQFunctionCreate(Ceed ceed, CeedInt N, CeedInt q_data_size, CeedQFunction *q…