Lines Matching refs:chebyshev_x
16 inline __device__ void ChebyshevPolynomialsAtPoint(const CeedScalar x, CeedScalar *chebyshev_x) { in ChebyshevPolynomialsAtPoint() argument
17 chebyshev_x[0] = 1.0; in ChebyshevPolynomialsAtPoint()
18 chebyshev_x[1] = 2 * x; in ChebyshevPolynomialsAtPoint()
19 …for (CeedInt i = 2; i < Q_1D; i++) chebyshev_x[i] = 2 * x * chebyshev_x[i - 1] - chebyshev_x[i - 2… in ChebyshevPolynomialsAtPoint()
24 CeedScalar chebyshev_x[3]; in ChebyshevDerivativeAtPoint() local
26 chebyshev_x[1] = 1.0; in ChebyshevDerivativeAtPoint()
27 chebyshev_x[2] = 2 * x; in ChebyshevDerivativeAtPoint()
31 chebyshev_x[(i + 1) % 3] = 2 * x * chebyshev_x[(i + 0) % 3] - chebyshev_x[(i + 2) % 3]; in ChebyshevDerivativeAtPoint()
32 …chebyshev_dx[i] = 2 * x * chebyshev_dx[i - 1] + 2 * chebyshev_x[(i + 0) % 3] - chebyshev_… in ChebyshevDerivativeAtPoint()
46 CeedScalar chebyshev_x[Q_1D]; in InterpAtPoints1d() local
49 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in InterpAtPoints1d()
57 r_V[comp] += chebyshev_x[i] * data.slice[i]; in InterpAtPoints1d()
68 CeedScalar chebyshev_x[Q_1D]; in InterpTransposeAtPoints1d() local
70 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in InterpTransposeAtPoints1d()
78 …atomicAdd(&data.slice[comp * Q_1D + (i + data.t_id_x) % Q_1D], chebyshev_x[(i + data.t_id_x) % Q_1… in InterpTransposeAtPoints1d()
93 CeedScalar chebyshev_x[Q_1D]; in GradAtPoints1d() local
95 ChebyshevDerivativeAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradAtPoints1d()
104 r_V[comp] += chebyshev_x[i] * data.slice[i]; in GradAtPoints1d()
115 CeedScalar chebyshev_x[Q_1D]; in GradTransposeAtPoints1d() local
117 ChebyshevDerivativeAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradTransposeAtPoints1d()
125 …atomicAdd(&data.slice[comp * Q_1D + (i + data.t_id_x) % Q_1D], chebyshev_x[(i + data.t_id_x) % Q_1… in GradTransposeAtPoints1d()
147 CeedScalar chebyshev_x[Q_1D]; in InterpAtPoints2d() local
154 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in InterpAtPoints2d()
158 buffer[i] += chebyshev_x[j] * data.slice[j + i * Q_1D]; in InterpAtPoints2d()
162 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in InterpAtPoints2d()
164 r_V[comp] += chebyshev_x[i] * buffer[i]; in InterpAtPoints2d()
177 CeedScalar chebyshev_x[Q_1D]; in InterpTransposeAtPoints2d() local
183 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in InterpTransposeAtPoints2d()
187 buffer[i] = chebyshev_x[i] * r_u; in InterpTransposeAtPoints2d()
190 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in InterpTransposeAtPoints2d()
198 …x < Q_1D && data.t_id_y < Q_1D) atomicAdd(&data.slice[jj + ii * Q_1D], chebyshev_x[jj] * buffer[ii… in InterpTransposeAtPoints2d()
216 CeedScalar chebyshev_x[Q_1D]; in GradAtPoints2d() local
224 if (dim == 0) ChebyshevDerivativeAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradAtPoints2d()
225 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradAtPoints2d()
229 buffer[i] += chebyshev_x[j] * data.slice[j + i * Q_1D]; in GradAtPoints2d()
233 if (dim == 1) ChebyshevDerivativeAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradAtPoints2d()
234 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradAtPoints2d()
236 r_V[comp + dim * NUM_COMP] += chebyshev_x[i] * buffer[i]; in GradAtPoints2d()
250 CeedScalar chebyshev_x[Q_1D]; in GradTransposeAtPoints2d() local
257 if (dim == 1) ChebyshevDerivativeAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradTransposeAtPoints2d()
258 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradTransposeAtPoints2d()
262 buffer[i] = chebyshev_x[i] * r_u; in GradTransposeAtPoints2d()
265 if (dim == 0) ChebyshevDerivativeAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradTransposeAtPoints2d()
266 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradTransposeAtPoints2d()
274 …x < Q_1D && data.t_id_y < Q_1D) atomicAdd(&data.slice[jj + ii * Q_1D], chebyshev_x[jj] * buffer[ii… in GradTransposeAtPoints2d()
297 CeedScalar chebyshev_x[Q_1D]; in InterpAtPoints3d() local
300 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[2], chebyshev_x); in InterpAtPoints3d()
301 const CeedScalar z = chebyshev_x[k]; in InterpAtPoints3d()
309 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in InterpAtPoints3d()
313 buffer[i] += chebyshev_x[j] * data.slice[j + i * Q_1D]; in InterpAtPoints3d()
317 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in InterpAtPoints3d()
319 r_V[comp] += chebyshev_x[i] * buffer[i] * z; in InterpAtPoints3d()
333 CeedScalar chebyshev_x[Q_1D]; in InterpTransposeAtPoints3d() local
336 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[2], chebyshev_x); in InterpTransposeAtPoints3d()
337 const CeedScalar z = chebyshev_x[k]; in InterpTransposeAtPoints3d()
344 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in InterpTransposeAtPoints3d()
348 buffer[i] = chebyshev_x[i] * r_u * z; in InterpTransposeAtPoints3d()
351 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in InterpTransposeAtPoints3d()
359 …x < Q_1D && data.t_id_y < Q_1D) atomicAdd(&data.slice[jj + ii * Q_1D], chebyshev_x[jj] * buffer[ii… in InterpTransposeAtPoints3d()
378 CeedScalar chebyshev_x[Q_1D]; in GradAtPoints3d() local
381 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[2], chebyshev_x); in GradAtPoints3d()
382 const CeedScalar z = chebyshev_x[k]; in GradAtPoints3d()
384 ChebyshevDerivativeAtPoint<Q_1D>(r_X[2], chebyshev_x); in GradAtPoints3d()
385 const CeedScalar dz = chebyshev_x[k]; in GradAtPoints3d()
395 if (dim == 0) ChebyshevDerivativeAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradAtPoints3d()
396 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradAtPoints3d()
400 buffer[i] += chebyshev_x[j] * data.slice[j + i * Q_1D]; in GradAtPoints3d()
404 if (dim == 1) ChebyshevDerivativeAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradAtPoints3d()
405 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradAtPoints3d()
409 r_V[comp + dim * NUM_COMP] += chebyshev_x[i] * buffer[i] * zz; in GradAtPoints3d()
424 CeedScalar chebyshev_x[Q_1D]; in GradTransposeAtPoints3d() local
427 ChebyshevPolynomialsAtPoint<Q_1D>(r_X[2], chebyshev_x); in GradTransposeAtPoints3d()
428 const CeedScalar z = chebyshev_x[k]; in GradTransposeAtPoints3d()
430 ChebyshevDerivativeAtPoint<Q_1D>(r_X[2], chebyshev_x); in GradTransposeAtPoints3d()
431 const CeedScalar dz = chebyshev_x[k]; in GradTransposeAtPoints3d()
440 if (dim == 1) ChebyshevDerivativeAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradTransposeAtPoints3d()
441 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[1], chebyshev_x); in GradTransposeAtPoints3d()
446 buffer[i] = chebyshev_x[i] * r_u * zz; in GradTransposeAtPoints3d()
449 if (dim == 0) ChebyshevDerivativeAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradTransposeAtPoints3d()
450 else ChebyshevPolynomialsAtPoint<Q_1D>(r_X[0], chebyshev_x); in GradTransposeAtPoints3d()
458 …x < Q_1D && data.t_id_y < Q_1D) atomicAdd(&data.slice[jj + ii * Q_1D], chebyshev_x[jj] * buffer[ii… in GradTransposeAtPoints3d()