| /libCEED/include/ceed/jit-source/sycl/ |
| H A D | sycl-gen-templates.h | 20 inline void loadMatrix(const CeedInt N, const CeedScalar *restrict d_B, CeedScalar *restrict B) { in loadMatrix() 21 const CeedInt item_id = get_local_linear_id(); in loadMatrix() 22 const CeedInt group_size = get_local_size(0) * get_local_size(1) * get_local_size(2); in loadMatrix() 23 for (CeedInt i = item_id; i < N; i += group_size) B[i] = d_B[i]; in loadMatrix() 33 inline void readDofsOffset1d(const CeedInt num_comp, const CeedInt strides_comp, const CeedInt P_1D… in readDofsOffset1d() 34 …const global CeedInt *restrict indices, const global CeedScalar *restrict d_u, private CeedScalar … in readDofsOffset1d() 35 const CeedInt item_id_x = get_local_id(0); in readDofsOffset1d() 36 const CeedInt elem = get_global_id(2); in readDofsOffset1d() 39 const CeedInt node = item_id_x; in readDofsOffset1d() 40 const CeedInt ind = indices[node + elem * P_1D]; in readDofsOffset1d() [all …]
|
| H A D | sycl-shared-basis-read-write-templates.h | 15 inline void loadMatrix(const CeedInt N, const CeedScalar *restrict d_B, CeedScalar *restrict B) { in loadMatrix() 16 const CeedInt item_id = get_local_linear_id(); in loadMatrix() 17 const CeedInt group_size = get_local_size(0) * get_local_size(1) * get_local_size(2); in loadMatrix() 18 for (CeedInt i = item_id; i < N; i += group_size) B[i] = d_B[i]; in loadMatrix() 28 inline void ReadElementStrided1d(const CeedInt NUM_COMP, const CeedInt P_1D, const CeedInt num_elem… in ReadElementStrided1d() 29 … const CeedInt strides_comp, const CeedInt strides_elem, global const CeedScalar *restrict d_u, in ReadElementStrided1d() 31 const CeedInt item_id_x = get_local_id(0); in ReadElementStrided1d() 32 const CeedInt elem = get_global_id(2); in ReadElementStrided1d() 35 const CeedInt node = item_id_x; in ReadElementStrided1d() 36 const CeedInt ind = node * strides_node + elem * strides_elem; in ReadElementStrided1d() [all …]
|
| H A D | sycl-shared-basis-tensor-templates.h | 19 inline void ContractX1d(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *restrict … in ContractX1d() 21 const CeedInt item_id_x = get_local_id(0); in ContractX1d() 28 for (CeedInt i = 0; i < P_1D; i++) { in ContractX1d() 38 inline void ContractTransposeX1d(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *… in ContractTransposeX1d() 40 const CeedInt item_id_x = get_local_id(0); in ContractTransposeX1d() 47 for (CeedInt i = 0; i < Q_1D; i++) { in ContractTransposeX1d() 57 inline void Interp1d(const CeedInt NUM_COMP, const CeedInt P_1D, const CeedInt Q_1D, private const … in Interp1d() 59 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in Interp1d() 67 inline void InterpTranspose1d(const CeedInt NUM_COMP, const CeedInt P_1D, const CeedInt Q_1D, priva… in InterpTranspose1d() 69 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in InterpTranspose1d() [all …]
|
| /libCEED/include/ceed/jit-source/cuda/ |
| H A D | cuda-shared-basis-read-write-templates.h | 17 for (CeedInt i = data.t_id; i < P * Q; i += blockDim.x * blockDim.y * blockDim.z) B[i] = d_B[i]; in LoadMatrix() 28 …eadElementStrided1d(SharedData_Cuda &data, const CeedInt elem, const CeedInt strides_node, const C… in ReadElementStrided1d() 29 … const CeedInt strides_elem, const CeedScalar *__restrict__ d_u, CeedScalar *r_u) { in ReadElementStrided1d() 31 const CeedInt node = data.t_id_x; in ReadElementStrided1d() 32 const CeedInt ind = node * strides_node + elem * strides_elem; in ReadElementStrided1d() 34 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in ReadElementStrided1d() 44 …iteElementStrided1d(SharedData_Cuda &data, const CeedInt elem, const CeedInt strides_node, const C… in WriteElementStrided1d() 45 … const CeedInt strides_elem, const CeedScalar *r_v, CeedScalar *d_v) { in WriteElementStrided1d() 47 const CeedInt node = data.t_id_x; in WriteElementStrided1d() 48 const CeedInt ind = node * strides_node + elem * strides_elem; in WriteElementStrided1d() [all …]
|
| H A D | cuda-gen-templates.h | 17 for (CeedInt i = data.t_id; i < P * Q; i += blockDim.x * blockDim.y * blockDim.z) B[i] = d_B[i]; in LoadMatrix() 28 …_device__ void ReadPoint(SharedData_Cuda &data, const CeedInt elem, const CeedInt p, const CeedInt… in ReadPoint() 29 … const CeedInt *__restrict__ indices, const CeedScalar *__restrict__ d_u, CeedScalar *r_u) { in ReadPoint() 30 const CeedInt ind = indices[p + elem * NUM_PTS]; in ReadPoint() 32 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in ReadPoint() 41 …device__ void WritePoint(SharedData_Cuda &data, const CeedInt elem, const CeedInt p, const CeedInt… in WritePoint() 42 … const CeedInt *__restrict__ indices, const CeedScalar *__restrict__ r_u, CeedScalar *d_u) { in WritePoint() 44 const CeedInt ind = indices[p + elem * NUM_PTS]; in WritePoint() 46 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in WritePoint() 60 inline __device__ void SetEVecStandard1d_Single(SharedData_Cuda &data, const CeedInt n, const CeedS… in SetEVecStandard1d_Single() [all …]
|
| H A D | cuda-ref-basis-tensor-at-points.h | 19 …for (CeedInt i = 2; i < Q_1D; i++) chebyshev_x[i] = 2 * x * chebyshev_x[i - 1] - chebyshev_x[i - 2… in ChebyshevPolynomialsAtPoint() 30 for (CeedInt i = 2; i < Q_1D; i++) { in ChebyshevDerivativeAtPoint() 43 extern "C" __global__ void InterpAtPoints(const CeedInt num_elem, const CeedScalar *__restrict__ ch… in InterpAtPoints() 44 … const CeedInt *__restrict__ points_per_elem, const CeedScalar *__restrict__ coords, in InterpAtPoints() 46 const CeedInt i = threadIdx.x; in InterpAtPoints() 54 for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) { in InterpAtPoints() 58 const CeedInt P = BASIS_P_1D; in InterpAtPoints() 59 const CeedInt Q = BASIS_Q_1D; in InterpAtPoints() 60 const CeedInt u_stride = BASIS_NUM_NODES; in InterpAtPoints() 61 const CeedInt v_stride = BASIS_NUM_PTS; in InterpAtPoints() [all …]
|
| H A D | cuda-ref-basis-tensor.h | 19 extern "C" __global__ void Interp(const CeedInt num_elem, const CeedInt is_transpose, const CeedSca… in Interp() 21 const CeedInt i = threadIdx.x; in Interp() 27 for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) { in Interp() 31 const CeedInt P = is_transpose ? BASIS_Q_1D : BASIS_P_1D; in Interp() 32 const CeedInt Q = is_transpose ? BASIS_P_1D : BASIS_Q_1D; in Interp() 33 const CeedInt stride_0 = is_transpose ? 1 : BASIS_P_1D; in Interp() 34 const CeedInt stride_1 = is_transpose ? BASIS_P_1D : 1; in Interp() 35 const CeedInt u_stride = is_transpose ? BASIS_NUM_QPTS : BASIS_NUM_NODES; in Interp() 36 const CeedInt v_stride = is_transpose ? BASIS_NUM_NODES : BASIS_NUM_QPTS; in Interp() 37 const CeedInt u_comp_stride = num_elem * (is_transpose ? BASIS_NUM_QPTS : BASIS_NUM_NODES); in Interp() [all …]
|
| H A D | cuda-ref-restriction-offset.h | 15 extern "C" __global__ void OffsetNoTranspose(const CeedInt *__restrict__ indices, const CeedScalar … in OffsetNoTranspose() 16 …for (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; … in OffsetNoTranspose() 17 const CeedInt ind = indices[node]; in OffsetNoTranspose() 18 const CeedInt loc_node = node % RSTR_ELEM_SIZE; in OffsetNoTranspose() 19 const CeedInt elem = node / RSTR_ELEM_SIZE; in OffsetNoTranspose() 21 for (CeedInt comp = 0; comp < RSTR_NUM_COMP; comp++) { in OffsetNoTranspose() 31 extern "C" __global__ void OffsetTranspose(const CeedInt *__restrict__ indices, const CeedScalar *_… in OffsetTranspose() 32 …for (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; … in OffsetTranspose() 33 const CeedInt ind = indices[node]; in OffsetTranspose() 34 const CeedInt loc_node = node % RSTR_ELEM_SIZE; in OffsetTranspose() [all …]
|
| H A D | cuda-ref-restriction-at-points.h | 16 extern "C" __global__ void AtPointsTranspose(const CeedInt *__restrict__ indices, const CeedInt *__… in AtPointsTranspose() 18 …for (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; … in AtPointsTranspose() 19 const CeedInt ind = indices[node]; in AtPointsTranspose() 20 const CeedInt loc_node = node % RSTR_ELEM_SIZE; in AtPointsTranspose() 21 const CeedInt elem = node / RSTR_ELEM_SIZE; in AtPointsTranspose() 24 for (CeedInt comp = 0; comp < RSTR_NUM_COMP; comp++) { in AtPointsTranspose() 30 extern "C" __global__ void AtPointsTranspose(const CeedInt *__restrict__ l_vec_indices, const CeedI… in AtPointsTranspose() 31 … const CeedInt *__restrict__ points_per_elem, const CeedInt *__restrict__ t_offsets, in AtPointsTranspose() 35 …for (CeedInt i = blockIdx.x * blockDim.x + threadIdx.x; i < RSTR_NUM_NODES; i += blockDim.x * grid… in AtPointsTranspose() 36 const CeedInt ind = l_vec_indices[i]; in AtPointsTranspose() [all …]
|
| /libCEED/include/ceed/jit-source/hip/ |
| H A D | hip-shared-basis-read-write-templates.h | 17 for (CeedInt i = data.t_id; i < P * Q; i += blockDim.x * blockDim.y * blockDim.z) B[i] = d_B[i]; in LoadMatrix() 28 …ReadElementStrided1d(SharedData_Hip &data, const CeedInt elem, const CeedInt strides_node, const C… in ReadElementStrided1d() 29 … const CeedInt strides_elem, const CeedScalar *__restrict__ d_u, CeedScalar *r_u) { in ReadElementStrided1d() 31 const CeedInt node = data.t_id_x; in ReadElementStrided1d() 32 const CeedInt ind = node * strides_node + elem * strides_elem; in ReadElementStrided1d() 34 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in ReadElementStrided1d() 44 …riteElementStrided1d(SharedData_Hip &data, const CeedInt elem, const CeedInt strides_node, const C… in WriteElementStrided1d() 45 … const CeedInt strides_elem, const CeedScalar *r_v, CeedScalar *d_v) { in WriteElementStrided1d() 47 const CeedInt node = data.t_id_x; in WriteElementStrided1d() 48 const CeedInt ind = node * strides_node + elem * strides_elem; in WriteElementStrided1d() [all …]
|
| H A D | hip-gen-templates.h | 17 for (CeedInt i = data.t_id; i < P * Q; i += blockDim.x * blockDim.y * blockDim.z) B[i] = d_B[i]; in LoadMatrix() 28 …__device__ void ReadPoint(SharedData_Hip &data, const CeedInt elem, const CeedInt p, const CeedInt… in ReadPoint() 29 … const CeedInt *__restrict__ indices, const CeedScalar *__restrict__ d_u, CeedScalar *r_u) { in ReadPoint() 30 const CeedInt ind = indices[p + elem * NUM_PTS]; in ReadPoint() 32 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in ReadPoint() 41 …_device__ void WritePoint(SharedData_Hip &data, const CeedInt elem, const CeedInt p, const CeedInt… in WritePoint() 42 … const CeedInt *__restrict__ indices, const CeedScalar *__restrict__ r_u, CeedScalar *d_u) { in WritePoint() 44 const CeedInt ind = indices[p + elem * NUM_PTS]; in WritePoint() 46 for (CeedInt comp = 0; comp < NUM_COMP; comp++) { in WritePoint() 60 inline __device__ void SetEVecStandard1d_Single(SharedData_Hip &data, const CeedInt n, const CeedSc… in SetEVecStandard1d_Single() [all …]
|
| H A D | hip-ref-basis-tensor-at-points.h | 19 …for (CeedInt i = 2; i < Q_1D; i++) chebyshev_x[i] = 2 * x * chebyshev_x[i - 1] - chebyshev_x[i - 2… in ChebyshevPolynomialsAtPoint() 30 for (CeedInt i = 2; i < Q_1D; i++) { in ChebyshevDerivativeAtPoint() 43 extern "C" __global__ void InterpAtPoints(const CeedInt num_elem, const CeedScalar *__restrict__ ch… in InterpAtPoints() 44 … const CeedInt *__restrict__ points_per_elem, const CeedScalar *__restrict__ coords, in InterpAtPoints() 46 const CeedInt i = threadIdx.x; in InterpAtPoints() 54 for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) { in InterpAtPoints() 58 const CeedInt P = BASIS_P_1D; in InterpAtPoints() 59 const CeedInt Q = BASIS_Q_1D; in InterpAtPoints() 60 const CeedInt u_stride = BASIS_NUM_NODES; in InterpAtPoints() 61 const CeedInt v_stride = BASIS_NUM_PTS; in InterpAtPoints() [all …]
|
| H A D | hip-ref-basis-tensor.h | 19 extern "C" __global__ void Interp(const CeedInt num_elem, const CeedInt is_transpose, const CeedSca… in Interp() 21 const CeedInt i = threadIdx.x; in Interp() 27 for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) { in Interp() 31 const CeedInt P = is_transpose ? BASIS_Q_1D : BASIS_P_1D; in Interp() 32 const CeedInt Q = is_transpose ? BASIS_P_1D : BASIS_Q_1D; in Interp() 33 const CeedInt stride_0 = is_transpose ? 1 : BASIS_P_1D; in Interp() 34 const CeedInt stride_1 = is_transpose ? BASIS_P_1D : 1; in Interp() 35 const CeedInt u_stride = is_transpose ? BASIS_NUM_QPTS : BASIS_NUM_NODES; in Interp() 36 const CeedInt v_stride = is_transpose ? BASIS_NUM_NODES : BASIS_NUM_QPTS; in Interp() 37 const CeedInt u_comp_stride = num_elem * (is_transpose ? BASIS_NUM_QPTS : BASIS_NUM_NODES); in Interp() [all …]
|
| H A D | hip-ref-restriction-offset.h | 15 extern "C" __global__ void OffsetNoTranspose(const CeedInt *__restrict__ indices, const CeedScalar … in OffsetNoTranspose() 16 …for (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; … in OffsetNoTranspose() 17 const CeedInt ind = indices[node]; in OffsetNoTranspose() 18 const CeedInt loc_node = node % RSTR_ELEM_SIZE; in OffsetNoTranspose() 19 const CeedInt elem = node / RSTR_ELEM_SIZE; in OffsetNoTranspose() 21 for (CeedInt comp = 0; comp < RSTR_NUM_COMP; comp++) { in OffsetNoTranspose() 31 extern "C" __global__ void OffsetTranspose(const CeedInt *__restrict__ indices, const CeedScalar *_… in OffsetTranspose() 32 …for (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; … in OffsetTranspose() 33 const CeedInt ind = indices[node]; in OffsetTranspose() 34 const CeedInt loc_node = node % RSTR_ELEM_SIZE; in OffsetTranspose() [all …]
|
| H A D | hip-ref-restriction-at-points.h | 16 extern "C" __global__ void AtPointsTranspose(const CeedInt *__restrict__ indices, const CeedInt *__… in AtPointsTranspose() 18 …for (CeedInt node = blockIdx.x * blockDim.x + threadIdx.x; node < RSTR_NUM_ELEM * RSTR_ELEM_SIZE; … in AtPointsTranspose() 19 const CeedInt ind = indices[node]; in AtPointsTranspose() 20 const CeedInt loc_node = node % RSTR_ELEM_SIZE; in AtPointsTranspose() 21 const CeedInt elem = node / RSTR_ELEM_SIZE; in AtPointsTranspose() 24 for (CeedInt comp = 0; comp < RSTR_NUM_COMP; comp++) { in AtPointsTranspose() 30 extern "C" __global__ void AtPointsTranspose(const CeedInt *__restrict__ l_vec_indices, const CeedI… in AtPointsTranspose() 31 … const CeedInt *__restrict__ points_per_elem, const CeedInt *__restrict__ t_offsets, in AtPointsTranspose() 35 …for (CeedInt i = blockIdx.x * blockDim.x + threadIdx.x; i < RSTR_NUM_NODES; i += blockDim.x * grid… in AtPointsTranspose() 36 const CeedInt ind = l_vec_indices[i]; in AtPointsTranspose() [all …]
|
| /libCEED/backends/sycl-ref/ |
| H A D | ceed-sycl-ref.hpp | 29 CeedInt num_nodes; 30 CeedInt num_elem; 31 CeedInt num_comp; 32 CeedInt elem_size; 33 CeedInt comp_stride; 34 CeedInt strides[3]; 35 CeedInt *h_offsets; 36 CeedInt *h_offsets_borrowed; 37 CeedInt *h_offsets_owned; 38 CeedInt *d_offsets; [all …]
|
| /libCEED/backends/cuda-ref/ |
| H A D | ceed-cuda-ref.h | 29 CeedInt num_nodes; 30 const CeedInt *h_offsets; 31 const CeedInt *h_offsets_borrowed; 32 const CeedInt *h_offsets_owned; 33 const CeedInt *d_offsets; 34 const CeedInt *d_offsets_borrowed; 35 const CeedInt *d_offsets_owned; 36 const CeedInt *d_t_offsets; 37 const CeedInt *d_t_indices; 38 const CeedInt *d_l_vec_indices; [all …]
|
| /libCEED/backends/hip-ref/ |
| H A D | ceed-hip-ref.h | 34 CeedInt num_nodes; 35 const CeedInt *h_offsets; 36 const CeedInt *h_offsets_borrowed; 37 const CeedInt *h_offsets_owned; 38 const CeedInt *d_offsets; 39 const CeedInt *d_offsets_borrowed; 40 const CeedInt *d_offsets_owned; 41 const CeedInt *d_t_offsets; 42 const CeedInt *d_t_indices; 43 const CeedInt *d_l_vec_indices; [all …]
|
| /libCEED/backends/avx/ |
| H A D | ceed-avx-tensor.c | 42 …TensorContract_Avx_Blocked(CeedTensorContract contract, CeedInt A, CeedInt B, CeedInt C, CeedInt J, in CeedTensorContract_Avx_Blocked() 43 … const CeedScalar *restrict t, CeedTransposeMode t_mode, const CeedInt add, in CeedTensorContract_Avx_Blocked() 44 … const CeedScalar *restrict u, CeedScalar *restrict v, const CeedInt JJ, const CeedInt CC) { in CeedTensorContract_Avx_Blocked() 45 CeedInt t_stride_0 = B, t_stride_1 = 1; in CeedTensorContract_Avx_Blocked() 52 for (CeedInt a = 0; a < A; a++) { in CeedTensorContract_Avx_Blocked() 54 for (CeedInt j = 0; j < (J / JJ) * JJ; j += JJ) { in CeedTensorContract_Avx_Blocked() 55 for (CeedInt c = 0; c < (C / CC) * CC; c += CC) { in CeedTensorContract_Avx_Blocked() 57 for (CeedInt jj = 0; jj < JJ; jj++) { in CeedTensorContract_Avx_Blocked() 58 … for (CeedInt cc = 0; cc < CC / 4; cc++) vv[jj][cc] = loadu(&v[(a * J + j + jj) * C + c + cc * 4]); in CeedTensorContract_Avx_Blocked() 60 for (CeedInt b = 0; b < B; b++) { in CeedTensorContract_Avx_Blocked() [all …]
|
| /libCEED/include/ |
| H A D | ceed-impl.h | 88 CeedInt num_vecs, max_vecs; 98 CeedInt num_view_tabs; 111 CeedInt num_rust_source_roots, max_rust_source_roots, num_rust_source_roots_readers; 112 CeedInt num_jit_source_roots, max_jit_source_roots, num_jit_source_roots_readers; 115 CeedInt num_jit_defines, max_jit_defines, num_jit_defines_readers; 121 …int (*ElemRestrictionCreate)(CeedMemType, CeedCopyMode, const CeedInt *, const bool *, const CeedI… 122 …int (*ElemRestrictionCreateAtPoints)(CeedMemType, CeedCopyMode, const CeedInt *, const bool *, con… 123 …int (*ElemRestrictionCreateBlocked)(CeedMemType, CeedCopyMode, const CeedInt *, const bool *, cons… 124 …int (*BasisCreateTensorH1)(CeedInt, CeedInt, CeedInt, const CeedScalar *, const CeedScalar *, cons… 125 …int (*BasisCreateH1)(CeedElemTopology, CeedInt, CeedInt, CeedInt, const CeedScalar *, const CeedSc… [all …]
|
| /libCEED/examples/ceed/ |
| H A D | ex2-surface.c | 48 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[3]); 49 …CartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[3], CeedInt degree, CeedInt num_comp,… 51 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[3], CeedInt mesh_degree, CeedVector … 52 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords); 57 CeedInt dim = 3; // dimension of the mesh in main() 58 CeedInt num_comp_x = 3; // number of x components in main() 59 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main() 60 CeedInt sol_degree = 4; // polynomial degree for the solution in main() 61 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() 62 CeedInt prob_size = -1; // approximate problem size in main() [all …]
|
| H A D | ex3-volume.c | 47 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim… 48 …artesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp… 50 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVecto… 51 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords); 56 CeedInt dim = 3; // dimension of the mesh in main() 57 CeedInt num_comp_x = 3; // number of x components in main() 58 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main() 59 CeedInt sol_degree = 4; // polynomial degree for the solution in main() 60 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() 61 CeedInt prob_size = -1; // approximate problem size in main() [all …]
|
| H A D | ex1-volume.c | 48 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim… 49 …artesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp… 51 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVecto… 52 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords); 57 CeedInt dim = 3; // dimension of the mesh in main() 58 CeedInt num_comp_x = 3; // number of x components in main() 59 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main() 60 CeedInt sol_degree = 4; // polynomial degree for the solution in main() 61 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() 62 CeedInt prob_size = -1; // approximate problem size in main() [all …]
|
| /libCEED/examples/rust-qfunctions/ |
| H A D | ex1-volume.c | 49 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim… 50 …artesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp… 52 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVecto… 53 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords); 58 CeedInt dim = 3; // dimension of the mesh in main() 59 CeedInt num_comp_x = 3; // number of x components in main() 60 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main() 61 CeedInt sol_degree = 4; // polynomial degree for the solution in main() 62 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main() 63 CeedInt prob_size = -1; // approximate problem size in main() [all …]
|
| /libCEED/backends/ref/ |
| H A D | ceed-ref.h | 21 const CeedInt *offsets; 22 const CeedInt *offsets_borrowed; 23 const CeedInt *offsets_owned; 30 …int (*Apply)(CeedElemRestriction, CeedInt, CeedInt, CeedInt, CeedInt, CeedInt, CeedTransposeMode, … 53 CeedInt *e_data_out_indices; 60 CeedInt num_inputs, num_outputs; 61 CeedInt qf_size_in, qf_size_out; 67 …estrictionCreate_Ref(CeedMemType mem_type, CeedCopyMode copy_mode, const CeedInt *offsets, const b… 70 CEED_INTERN int CeedBasisCreateTensorH1_Ref(CeedInt dim, CeedInt P_1d, CeedInt Q_1d, const CeedScal… 72 …_INTERN int CeedBasisCreateH1_Ref(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt n… [all …]
|