| /libCEED/include/ceed/jit-source/magma/ |
| H A D | magma-common-tensor.h | 18 template <typename T, int LENGTH, int NUM_COMP> 19 static __device__ __inline__ void read_1d(const T *devptr, const int compstride, T *sBuffer[NUM_COM… in read_1d() 30 template <typename T, int LENGTH, int NUM_COMP> 31 static __device__ __inline__ void write_1d(T *sBuffer[NUM_COMP], T *devptr, const int compstride, c… in write_1d() 42 template <typename T, int LENGTH, int NUM_COMP> 43 static __device__ __inline__ void sum_1d(T *sBuffer[NUM_COMP], T *devptr, const int compstride, con… in sum_1d() 58 template <typename T, int P, int DIM_U, int NUM_COMP, int rU_SIZE, int i_DIM> 59 …c __device__ __inline__ void read_U_2d(const T *dU, const int compstride, T rU[DIM_U][NUM_COMP][rU… in read_U_2d() 94 template <typename T, int Q, int DIM_V, int NUM_COMP, int rV_SIZE, int i_DIM> 95 static __device__ __inline__ void read_V_2d(const T *dV, const int compstride, T rV[DIM_V][NUM_COMP… in read_V_2d() [all …]
|
| H A D | magma-common-nontensor.h | 19 template <typename T, int P, int Q, int BY> 20 …line__ void read_A_notrans_g2r_1D_nosync(const int tx, const int ty, const T *dA, T *sA, T rA[Q]) { in read_A_notrans_g2r_1D_nosync() 44 template <typename T, int P, int Q, int BY> 45 …inline__ void read_A_trans_g2r_1D_nosync(const int tx, const int ty, const T *dA, T *sA, T rA[Q]) { in read_A_trans_g2r_1D_nosync() 69 template <typename T, int P, int Q, int NB> 70 static __device__ __inline__ void read_B_g2s_1D_nosync(const int tx, const int n, const T *dB, T *s… in read_B_g2s_1D_nosync() 93 template <typename T, int P, int Q, int NB> 94 static __device__ __inline__ void write_C_r2g_1D_nosync(const int tx, const int n, T rC[NB], T *dC)… in write_C_r2g_1D_nosync() 112 template <typename T, int P, int Q, int NB> 113 static __device__ __inline__ void sum_C_r2g_1D_nosync(const int tx, const int n, T rC[NB], T *dC) { in sum_C_r2g_1D_nosync() [all …]
|
| H A D | magma-basis-grad-3d.h | 19 template <typename T, bool Add> 22 template <typename T> 23 struct magma_grad_3d_device_accumulate<T, true> { 24 static __device__ __inline__ void op(T &rV, const T &rTmp) { rV += rTmp; } 27 template <typename T> 28 struct magma_grad_3d_device_accumulate<T, false> { 29 static __device__ __inline__ void op(T &rV, const T &rTmp) { rV = rTmp; } 40 template <typename T, int DIM_U, int DIM_V, int NUM_COMP, int P, int Q, int rU_SIZE, int rV_SIZE, i… 41 static __device__ __inline__ void magma_grad_3d_device(const T *sTinterp, const T *sTgrad, T rU[DIM… 42 … T rV[DIM_V][NUM_COMP][rV_SIZE], const int tx, T rTmp, T *swork) { [all …]
|
| H A D | magma-basis-grad-2d.h | 18 template <typename T, bool Add> 21 template <typename T> 22 struct magma_grad_2d_device_accumulate<T, true> { 23 static __device__ __inline__ void op(T &rV, const T &rTmp) { rV += rTmp; } 26 template <typename T> 27 struct magma_grad_2d_device_accumulate<T, false> { 28 static __device__ __inline__ void op(T &rV, const T &rTmp) { rV = rTmp; } 39 template <typename T, int DIM_U, int DIM_V, int NUM_COMP, int P, int Q, int rU_SIZE, int rV_SIZE, i… 40 static __device__ __inline__ void magma_grad_2d_device(const T *sTinterp, const T *sTgrad, T rU[DIM… 41 … T rV[DIM_V][NUM_COMP][rV_SIZE], const int tx, T rTmp, T *swork) { [all …]
|
| H A D | magma-basis-interp-3d.h | 18 template <typename T, int DIM_U, int DIM_V, int NUM_COMP, int P, int Q, int rU_SIZE, int rV_SIZE> 19 …__device__ __inline__ void magma_interp_3d_device(const T *sT, T rU[DIM_U][NUM_COMP][rU_SIZE], T r… in magma_interp_3d_device() 20 T rTmp[Q], T *swork) { in magma_interp_3d_device() 37 T *sTmp = swork + batchid * (1 * Q); in magma_interp_3d_device() 53 T *sTmp = swork + batchid * (Q * P); // sTmp is input in magma_interp_3d_device() 68 T *sTmp = swork + batchid * (Q * Q); in magma_interp_3d_device() 80 T *sTmp = swork; in magma_interp_3d_device()
|
| H A D | magma-basis-interp-2d.h | 18 template <typename T, int DIM_U, int DIM_V, int NUM_COMP, int P, int Q, int rU_SIZE, int rV_SIZE> 19 …__device__ __inline__ void magma_interp_2d_device(const T *sT, T rU[DIM_U][NUM_COMP][rU_SIZE], T r… in magma_interp_2d_device() 20 T rTmp, T *swork) { in magma_interp_2d_device() 37 T *sTmp = swork + batchid * (1 * Q); in magma_interp_2d_device() 52 T *sTmp = swork + batchid * (Q * P); in magma_interp_2d_device()
|
| H A D | magma-basis-weight-1d.h | 14 template <typename T, int Q> 15 static __device__ __inline__ void magma_weight_1d_device(const T *sTweight, T *sV, const int tx) { in magma_weight_1d_device()
|
| H A D | magma-basis-weight-2d.h | 14 template <typename T, int DIM, int NUM_COMP, int Q, int i_DIM, int i_COMP> 15 static __device__ __inline__ void magma_weight_2d_device(const T *sTweight, T rV[DIM][NUM_COMP][Q],… in magma_weight_2d_device()
|
| H A D | magma-basis-weight-3d.h | 14 template <typename T, int DIM, int NUM_COMP, int Q, int i_DIM, int i_COMP> 15 static __device__ __inline__ void magma_weight_3d_device(const T *sTweight, T rV[DIM][NUM_COMP][Q],… in magma_weight_3d_device()
|
| H A D | magma-basis-grad-1d.h | 17 template <typename T, int DIM, int NUM_COMP, int P, int Q> 18 static __device__ __inline__ void magma_grad_1d_device(const T *sT, T *sU[NUM_COMP], T *sV[NUM_COMP… in magma_grad_1d_device() 30 T rv = 0.0; in magma_grad_1d_device()
|
| H A D | magma-basis-interp-1d.h | 17 template <typename T, int DIM, int NUM_COMP, int P, int Q> 18 static __device__ __inline__ void magma_interp_1d_device(const T *sT, T *sU[NUM_COMP], T *sV[NUM_CO… in magma_interp_1d_device() 30 T rv = 0.0; in magma_interp_1d_device()
|
| H A D | magma-basis-interp-deriv-nontensor.h | 13 template <typename T, int Q_COMP, int P, int Q, int NB> 56 template <typename T, int Q_COMP, int P, int Q, int NB> 102 template <typename T, int Q_COMP, int P, int Q, int NB> 148 template <typename T, int P, int Q, int NB> 184 template <typename T, int P, int Q, int NB> 220 template <typename T, int P, int Q, int NB>
|
| /libCEED/examples/nek/boxes/ |
| H A D | b1e.rea | 129 T IFFLOW 131 T IFTRAN 132 T T F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields) 133 F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh) 185 T COORDINATES 186 T VELOCITY 187 T PRESSURE 188 T TEMPERATURE
|
| /libCEED/julia/LibCEED.jl/test/ |
| H A D | buildmats.jl | 1 function build_mats_hdiv(qref, qweight, ::Type{T}=Float64) where {T} argument 3 interp = Array{T}(undef, dim, Q, P) 4 div = Array{T}(undef, Q, P) 42 function build_mats_hcurl(qref, qweight, ::Type{T}=Float64) where {T} argument 44 interp = Array{T}(undef, dim, Q, P) 45 curl = Array{T}(undef, 1, Q, P)
|
| /libCEED/examples/petsc/ |
| H A D | index.md | 90 …\bm x \rVert} \bm I_{(3\times 3)} - \frac{1}{\lVert \bm x \rVert^3} (\bm x \bm x^T)_{(3\times 3)} . 109 …T \frac{\partial\overset{\circ}{\bm{x}}}{\partial \bm{X}}_{(3\times2)} \right)^{-1} \frac{\partial… 121 …rtial \overset\circ{\bm x}} \left( \frac{\partial u}{\partial \overset\circ{\bm x}} \right)^T \, dS 122 …\partial \overset\circ{\bm x}} \right)^T}_{\bm g_{(2\times 2)}} \left(\frac{\partial u}{\partial … 129 …T \frac{\partial\overset{\circ}{\bm{x}}}{\partial \bm{X}} \right)^{-1}_{(2\times 2)} \frac{\partia… 136 …T \, dS = \int_{\Omega} \frac{\partial v}{\partial \bm X} \underbrace{\left(\frac{\partial \ov… 153 …on operators given by interpolation basis operations, $\bm{B}$, and $\bm{B}^T$, respectively, act …
|
| /libCEED/examples/fluids/problems/ |
| H A D | bc_freestream.c | 261 const CeedScalar T = 200; in TestRiemannHLL_fwd() local 263 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas->cv * T; in TestRiemannHLL_fwd() 266 const CeedScalar Y0_left[5] = {p, u[0], u[1], u[2], T}; in TestRiemannHLL_fwd() 267 const CeedScalar Y0_right[5] = {1.2 * p, 1.2 * u[0], 1.2 * u[1], 1.2 * u[2], 1.2 * T}; in TestRiemannHLL_fwd() 323 const CeedScalar T = 200; in TestRiemannHLLC_fwd() local 325 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas->cv * T; in TestRiemannHLLC_fwd() 328 const CeedScalar Y0_left[5] = {p, u[0], u[1], u[2], T}; in TestRiemannHLLC_fwd() 329 const CeedScalar Y0_right[5] = {1.2 * p, 1.2 * u[0], 1.2 * u[1], 1.2 * u[2], 1.2 * T}; in TestRiemannHLLC_fwd() 385 const CeedScalar T = 200; in TestComputeHLLSpeeds_Roe_fwd() local 387 const CeedScalar p = (HeatCapacityRatio(gas) - 1) * rho * gas->cv * T; in TestComputeHLLSpeeds_Roe_fwd() [all …]
|
| /libCEED/examples/fluids/qfunctions/ |
| H A D | eulervortex.h | 87 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local 92 T = 1 + delta_T; in Exact_Euler() 96 rho = pow(T / S_vortex, 1 / (gamma - 1.)); in Exact_Euler() 97 P = rho * T; in Exact_Euler() 136 T = 1. - S_bubble * exp(1. - r * r); in Exact_Euler() 137 rho = P / (R * T); in Exact_Euler() 144 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler() 149 T = 1. - S_bubble * exp(1. - r * r); in Exact_Euler() 150 rho = P / (R * T); in Exact_Euler() 159 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler() [all …]
|
| /libCEED/julia/LibCEED.jl/examples/ |
| H A D | ex3-volume.jl | 91 # Store diffusion components (J^T * J) 136 # The qdata contains the symmetric diffusion tensor (J^T*J) 137 # dv_i = sum_j (J^T*J)_{i,j} * du_j 170 # Compute the mesh volume using the mass+diffusion operator: vol = 1^T \cdot (M + K) \cdot 1 171 print("Computing the mesh volume using the formula: vol = 1^T * (M + K) * 1...")
|
| H A D | ex2-surface.jl | 138 # sa = 1^T \cdot abs( K \cdot x). 139 print("Computing the mesh surface area using the formula: sa = 1^T.|K.x| ...")
|
| H A D | ex1-volume.jl | 131 # Compute the mesh volume using the mass operator: vol = 1^T \cdot M \cdot 1 132 print("Computing the mesh volume using the formula: vol = 1^T.M.1 ...")
|
| /libCEED/examples/solids/ |
| H A D | index.md | 70 \bm{\epsilon} = \dfrac{1}{2}\left(\nabla \bm{u} + \nabla \bm{u}^T \right). 78 Q \bm \sigma(\bm \epsilon) Q^T = \bm \sigma(Q \bm \epsilon Q^T), 93 \gamma(\bm \epsilon) = \gamma(Q \bm \epsilon Q^T) 170 \diff \bm{\epsilon} = \dfrac{1}{2}\left( \nabla \diff \bm{u} + \nabla \diff \bm{u}^T \right) 256 \bm C = \bm F^T \bm F 262 …C - \bm I_3) = \frac 1 2 \Big( \nabla_X \bm u + (\nabla_X \bm u)^T + (\nabla_X \bm u)^T \nabla_X \… 432 … E}{\partial \bm F} \!:\! \diff \bm F = \frac 1 2 \Big( \diff \bm F^T \bm F + \bm F^T \diff \bm F … 476 …+ \lambda (\bm C^{-1} : \diff \bm E) \bm F^{-T} + 2(\mu - \lambda \log J) \bm F^{-T} \diff\bm E \,… 478 …+ \lambda (\bm F^{-T} : \diff \bm F) \bm F^{-T} + (\mu - \lambda \log J) \bm F^{-T} (\bm F^T \diff… 480 …lambda (\bm F^{-T} : \diff \bm F) \bm F^{-T} + (\mu - \lambda \log J) \Big( \diff \bm F\, \bm C^{-… [all …]
|
| /libCEED/examples/ |
| H A D | bps.md | 39 \langle v,u \rangle = \bm v^T M \bm u , \qquad \langle v,f\rangle = \bm v^T \bm b \,. 69 a(v,u) = \bm v^T K \bm u ,
|
| /libCEED/doc/sphinx/source/ |
| H A D | libCEEDapi.md | 51 …rder ($Q_3$) scalar continuous ($H^1$) elements, where we use the notions **T-vector**, **L-vector… 75 - True degrees of freedom/unknowns, **T-vector**: 81 > ```{image} ../../img/T-vector.svg 122 > - here X can be any of the T, L, E, or Q categories 123 > - for example, the mass matrix operator maps a **T-vector** to a **T'-vector** 124 > - the solutions vector is a **T-vector**, and the RHS vector is a **T'-vector** 125 …> - using the parallel prolongation operator, one can map the solution **T-vector** to a solution … 132 …> - the T in TA indicates that the data format represents an operator from a **T-vector** to a **T… 138 > - note that $\bm{P}$ is the operator mapping T-vectors to L-vectors. 159 For example, one can compute and store a global matrix on **T-vector** level. [all …]
|
| /libCEED/julia/LibCEED.jl/src/ |
| H A D | Quadrature.jl | 16 struct QuadratureMode{T} end argument
|
| H A D | UserQFunction.jl | 8 @inline function extract_context(ptr, ::Type{T}) where {T} argument 9 unsafe_load(Ptr{T}(ptr))
|