| /libCEED/examples/solids/qfunctions/ |
| H A D | manufactured-force.h | 43 CeedScalar x = coords[i + 0 * Q], y = coords[i + 1 * Q], z = coords[i + 2 * Q]; in SetupMMSForce() local 48 …force[i + 0 * Q] = (-(E * (cos(x * 2.0) * cos(y * 3.0) * exp(z * 4.0) * 4.0 - cos(z * 4.0) * sin(y… in SetupMMSForce() 50 …(E * (cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * (4.5) + sin(x * 2.0) * sin(z * 4.0) * exp(y * 3… in SetupMMSForce() 52 … (E * nu * cos(x * 2.0) * cos(y * 3.0) * exp(z * 4.0) * 8.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce() 53 … (E * nu * sin(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 6.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce() 54 …(E * cos(z * 4.0) * sin(y * 3.0) * exp(x * 2.0) * (nu - 1.0) * 4.0) / ((nu * 2.0 - 1.0) * (nu + 1.… in SetupMMSForce() 58 …force[i + 1 * Q] = (-(E * (cos(y * 3.0) * cos(z * 4.0) * exp(x * 2.0) * 3.0 - cos(x * 2.0) * sin(z… in SetupMMSForce() 60 …(E * (cos(x * 2.0) * sin(z * 4.0) * exp(y * 3.0) * 8.0 + sin(x * 2.0) * sin(y * 3.0) * exp(z * 4.0… in SetupMMSForce() 62 … (E * nu * cos(y * 3.0) * cos(z * 4.0) * exp(x * 2.0) * 6.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce() 63 … (E * nu * sin(x * 2.0) * sin(y * 3.0) * exp(z * 4.0) * 12.0) / ((nu * 2.0 - 1.0) * (nu + 1.0)) - in SetupMMSForce() [all …]
|
| /libCEED/tests/ |
| H A D | t130-vector.c | 16 CeedVector x, y; in main() local 23 CeedGetWorkVector(ceed, 20, &y); in main() 24 if (y != x_copy) printf("failed to return same work vector"); in main() 25 CeedRestoreWorkVector(ceed, &y); in main() 30 CeedVector x, y; in main() local 37 CeedGetWorkVector(ceed, 30, &y); in main() 38 if (y == x_copy) printf("failed to return new work vector"); in main() 39 CeedRestoreWorkVector(ceed, &y); in main()
|
| H A D | t122-vector.c | 14 CeedVector x, y, w; in main() local 22 CeedVectorCreate(ceed, len, &y); in main() 29 CeedVectorSetArray(y, CEED_MEM_HOST, CEED_COPY_VALUES, array); in main() 33 CeedVectorPointwiseMult(w, x, y); in main() 45 CeedVectorPointwiseMult(w, w, y); in main() 73 CeedVectorSyncArray(y, type); in main() 75 CeedVectorPointwiseMult(y, y, y); in main() 76 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 84 CeedVectorRestoreArrayRead(y, &read_array); in main() 87 CeedVectorDestroy(&y); in main()
|
| H A D | t125-vector.c | 14 CeedVector x, y; in main() local 21 CeedVectorCreate(ceed, len, &y); in main() 27 CeedVectorSetArray(y, CEED_MEM_HOST, CEED_COPY_VALUES, array); in main() 33 CeedVectorSyncArray(y, type); in main() 35 CeedVectorAXPBY(y, -0.5, 1.0, x); in main() 40 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 48 CeedVectorRestoreArrayRead(y, &read_array); in main() 52 CeedVectorDestroy(&y); in main()
|
| H A D | t121-vector.c | 14 CeedVector x, y; in main() local 21 CeedVectorCreate(ceed, len, &y); in main() 27 CeedVectorSetArray(y, CEED_MEM_HOST, CEED_COPY_VALUES, array); in main() 33 CeedVectorSyncArray(y, type); in main() 35 CeedVectorAXPY(y, -0.5, x); in main() 40 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 48 CeedVectorRestoreArrayRead(y, &read_array); in main() 52 CeedVectorDestroy(&y); in main()
|
| H A D | t127-vector.c | 10 CeedVector x, y; in main() local 16 CeedVectorCreate(ceed, len, &y); in main() 38 CeedVectorSetValue(y, 0.0); in main() 39 CeedVectorCopyStrided(x, start, -1, step, y); in main() 43 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 53 CeedVectorRestoreArrayRead(y, &read_array); in main() 57 CeedVectorDestroy(&y); in main()
|
| H A D | t200-elemrestriction-f.f90 | 7 integer x,y variable 42 call ceedvectorcreate(ceed,2*ne,y,err); 43 call ceedvectorsetvalue(y,0.d0,err); 44 call ceedelemrestrictionapply(r,ceed_notranspose,x,y,& 47 call ceedvectorgetarrayread(y,ceed_mem_host,yy,yoffset,err) 56 call ceedvectorrestorearrayread(y,yy,yoffset,err) 59 call ceedvectordestroy(y,err)
|
| H A D | t201-elemrestriction-f.f90 | 7 integer x,y variable 38 call ceedvectorcreate(ceed,2*ne,y,err) 39 call ceedvectorsetvalue(y,0.d0,err) 40 call ceedelemrestrictionapply(r,ceed_notranspose,x,y,& 43 call ceedvectorgetarrayread(y,ceed_mem_host,yy,yoffset,err) 59 call ceedvectorrestorearrayread(y,yy,yoffset,err) 62 call ceedvectordestroy(y,err)
|
| H A D | t105-vector.c | 9 CeedVector x, y; in main() local 15 CeedVectorCreate(ceed, len, &y); in main() 26 CeedVectorSetArray(y, CEED_MEM_DEVICE, CEED_COPY_VALUES, (CeedScalar *)read_array); in main() 32 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 36 CeedVectorRestoreArrayRead(y, &read_array); in main() 40 CeedVectorDestroy(&y); in main()
|
| H A D | t103-vector.c | 9 CeedVector x, y; in main() local 15 CeedVectorCreate(ceed, len, &y); in main() 27 CeedVectorSetArray(y, CEED_MEM_HOST, CEED_COPY_VALUES, (CeedScalar *)read_array); in main() 33 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 37 CeedVectorRestoreArrayRead(y, &read_array); in main() 41 CeedVectorDestroy(&y); in main()
|
| H A D | t105-vector-f.f90 | 7 integer x,y,i,n variable 22 call ceedvectorcreate(ceed,n,y,err) 31 call ceedvectorsetarray(y,ceed_mem_device,ceed_copy_values,b,boffset,err) 33 call ceedvectorgetarrayread(y,ceed_mem_host,c,coffset,err) 42 call ceedvectorrestorearrayread(y,c,coffset,err) 45 call ceedvectordestroy(y,err)
|
| H A D | t103-vector-f.f90 | 7 integer x,y,i,n variable 21 call ceedvectorcreate(ceed,n,y,err) 31 call ceedvectorsetarray(y,ceed_mem_host,ceed_copy_values,xx,xoffset,err) 34 call ceedvectorgetarrayread(y,ceed_mem_host,yy,yoffset,err) 45 call ceedvectorrestorearrayread(y,yy,yoffset,err) 47 call ceedvectordestroy(y,err)
|
| H A D | t106-vector.c | 9 CeedVector x, y; in main() local 16 CeedVectorCreate(ceed, len, &y); in main() 24 CeedVectorSetArray(y, CEED_MEM_HOST, CEED_COPY_VALUES, initial_array); in main() 29 CeedVectorGetArrayRead(y, CEED_MEM_DEVICE, &read_array); in main() 31 CeedVectorRestoreArrayRead(y, &read_array); in main() 39 CeedVectorDestroy(&y); in main()
|
| H A D | t207-elemrestriction.c | 10 CeedVector x, y; in main() local 18 CeedVectorCreate(ceed, 2 * (num_elem + 1), &y); in main() 19 CeedVectorSetValue(y, 0.0); // Allocates array in main() 46 CeedElemRestrictionApply(elem_restriction, CEED_TRANSPOSE, x, y, CEED_REQUEST_IMMEDIATE); in main() 52 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &y_array); in main() 66 CeedVectorRestoreArrayRead(y, &y_array); in main() 70 CeedVectorDestroy(&y); in main()
|
| H A D | t200-elemrestriction.c | 9 CeedVector x, y; in main() local 23 CeedVectorCreate(ceed, num_elem * 2, &y); in main() 30 CeedElemRestrictionApply(elem_restriction, CEED_NOTRANSPOSE, x, y, CEED_REQUEST_IMMEDIATE); in main() 34 CeedVectorGetArrayRead(y, CEED_MEM_HOST, &read_array); in main() 38 CeedVectorRestoreArrayRead(y, &read_array); in main() 42 CeedVectorDestroy(&y); in main()
|
| H A D | t106-vector-f.f90 | 7 integer x,y,i,n variable 22 call ceedvectorcreate(ceed,n,y,err) 34 call ceedvectorsetarray(y,ceed_mem_host,ceed_use_pointer,b,boffset,err) 37 call ceedvectorsetarray(y,ceed_mem_device,ceed_copy_values,c,coffset,err) 40 call ceedvectorsyncarray(y,ceed_mem_host,err) 51 call ceedvectordestroy(y,err)
|
| H A D | t208-elemrestriction-f.f90 | 7 integer x,y variable 44 call ceedvectorcreate(ceed,blksize*elemsize,y,err); 45 call ceedvectorsetvalue(y,0.d0,err) 48 call ceedelemrestrictionapplyblock(r,1,ceed_notranspose,x,y,& 50 call ceedvectorgetarrayread(y,ceed_mem_host,yy,yoffset,err) 70 call ceedvectorrestorearrayread(y,yy,yoffset,err) 74 call ceedelemrestrictionapplyblock(r,1,ceed_transpose,y,x,& 94 call ceedvectordestroy(y,err)
|
| /libCEED/examples/solids/src/ |
| H A D | boundary.c | 25 PetscScalar y = coords[1]; in BCMMS() local 30 u[0] = exp(2 * x) * sin(3 * y) * cos(4 * z) / 1e8 * load_increment; in BCMMS() 31 u[1] = exp(3 * y) * sin(4 * z) * cos(2 * x) / 1e8 * load_increment; in BCMMS() 32 u[2] = exp(4 * z) * sin(2 * x) * cos(3 * y) / 1e8 * load_increment; in BCMMS() 44 PetscScalar y = coords[1]; in BCClamp() local 56 c_0 = clampMax[6] * M_PI, c_1 = clampMax[7] * M_PI, cx = kx * x + ky * y + kz * z, in BCClamp() 61 …u[0] = lx + s * (-kz * y + ky * z) + (1 - c) * (-(ky * ky + kz * kz) * x + kx * ky * y + kx * kz *… in BCClamp() 62 …u[1] = ly + s * (kz * x + -kx * z) + (1 - c) * (kx * ky * x + -(kx * kx + kz * kz) * y + ky * kz *… in BCClamp() 63 …u[2] = lz + s * (-ky * x + kx * y) + (1 - c) * (kx * kz * x + ky * kz * y + -(kx * kx + ky * ky) *… in BCClamp()
|
| /libCEED/python/tests/ |
| H A D | test-2-elemrestriction.py | 37 y = ceed.Vector(2 * num_elem) 38 y.set_value(0) 40 r.apply(x, y) 42 with y.array_read() as y_array: 63 y = ceed.Vector(2 * num_elem) 64 y.set_value(0) 66 r.apply(x, y) 68 with y.array_read() as y_array: 96 y = ceed.Vector(num_blk * blk_size * elem_size) 97 y.set_value(0) [all …]
|
| H A D | test-1-vector.py | 102 y = ceed.Vector(n) 108 y.set_array(x_array, cmode=libceed.USE_POINTER) 110 with y.array_read() as y_array: 149 y = ceed.Vector(n) 155 y.set_array(arr, memtype=libceed.MEM_DEVICE) 158 with y.array_read() as b: 240 y = ceed.Vector(n) 244 y.set_array(a, cmode=libceed.COPY_VALUES) 246 y.axpy(-0.5, x) 247 with y.array() as b: [all …]
|
| /libCEED/include/ceed/jit-source/cuda/ |
| H A D | cuda-shared-basis-tensor-at-points.h | 29 data.t_id_y = threadIdx.y; in InterpAtPoints() 31 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpAtPoints() 60 …edInt point_loop_bound = (blockDim.x * blockDim.y) * ceil(1.0 * BASIS_NUM_PTS / (blockDim.x * bloc… in InterpAtPoints() 62 …for (CeedInt i = threadIdx.x + threadIdx.y * blockDim.x; i < point_loop_bound; i += blockDim.x * b… in InterpAtPoints() 85 data.t_id_y = threadIdx.y; in InterpTransposeAtPoints() 87 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTransposeAtPoints() 117 …edInt point_loop_bound = (blockDim.x * blockDim.y) * ceil(1.0 * BASIS_NUM_PTS / (blockDim.x * bloc… in InterpTransposeAtPoints() 119 …for (CeedInt i = threadIdx.x + threadIdx.y * blockDim.x; i < point_loop_bound; i += blockDim.x * b… in InterpTransposeAtPoints() 155 data.t_id_y = threadIdx.y; in InterpTransposeAddAtPoints() 157 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTransposeAddAtPoints() [all …]
|
| H A D | cuda-shared-basis-nontensor.h | 23 data.t_id_y = threadIdx.y; in Interp() 25 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in Interp() 50 data.t_id_y = threadIdx.y; in InterpTranspose() 52 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTranspose() 77 data.t_id_y = threadIdx.y; in InterpTransposeAdd() 79 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in InterpTransposeAdd() 106 data.t_id_y = threadIdx.y; in Grad() 108 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in Grad() 133 data.t_id_y = threadIdx.y; in GradTranspose() 135 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in GradTranspose() [all …]
|
| /libCEED/include/ceed/jit-source/hip/ |
| H A D | hip-shared-basis-tensor-at-points.h | 30 data.t_id_y = threadIdx.y; in __launch_bounds__() 32 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() 61 …edInt point_loop_bound = (blockDim.x * blockDim.y) * ceil(1.0 * BASIS_NUM_PTS / (blockDim.x * bloc… in __launch_bounds__() 63 …for (CeedInt i = threadIdx.x + threadIdx.y * blockDim.x; i < point_loop_bound; i += blockDim.x * b… in __launch_bounds__() 86 data.t_id_y = threadIdx.y; in __launch_bounds__() 88 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() 118 …edInt point_loop_bound = (blockDim.x * blockDim.y) * ceil(1.0 * BASIS_NUM_PTS / (blockDim.x * bloc… in __launch_bounds__() 120 …for (CeedInt i = threadIdx.x + threadIdx.y * blockDim.x; i < point_loop_bound; i += blockDim.x * b… in __launch_bounds__() 156 data.t_id_y = threadIdx.y; in __launch_bounds__() 158 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() [all …]
|
| H A D | hip-shared-basis-nontensor.h | 24 data.t_id_y = threadIdx.y; in __launch_bounds__() 26 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() 51 data.t_id_y = threadIdx.y; in __launch_bounds__() 53 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() 78 data.t_id_y = threadIdx.y; in __launch_bounds__() 80 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() 108 data.t_id_y = threadIdx.y; in __launch_bounds__() 110 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() 135 data.t_id_y = threadIdx.y; in __launch_bounds__() 137 data.t_id = threadIdx.x + threadIdx.y * blockDim.x + threadIdx.z * blockDim.y * blockDim.x; in __launch_bounds__() [all …]
|
| /libCEED/interface/ |
| H A D | ceed-vector.c | 770 int CeedVectorAXPY(CeedVector y, CeedScalar alpha, CeedVector x) { in CeedVectorAXPY() argument 776 CeedCall(CeedVectorGetLength(y, &length_y)); in CeedVectorAXPY() 778 CeedCheck(length_x == length_y, CeedVectorReturnCeed(y), CEED_ERROR_UNSUPPORTED, in CeedVectorAXPY() 782 …CeedCheck(x != y, CeedVectorReturnCeed(y), CEED_ERROR_UNSUPPORTED, "Cannot use same vector for x a… in CeedVectorAXPY() 785 CeedCheck(has_valid_array_x, CeedVectorReturnCeed(y), CEED_ERROR_BACKEND, in CeedVectorAXPY() 787 CeedCall(CeedVectorHasValidArray(y, &has_valid_array_y)); in CeedVectorAXPY() 788 CeedCheck(has_valid_array_y, CeedVectorReturnCeed(y), CEED_ERROR_BACKEND, in CeedVectorAXPY() 794 CeedCall(CeedVectorGetCeed(y, &ceed_y)); in CeedVectorAXPY() 800 CeedCheck(ceed_parent_x == ceed_parent_y, CeedVectorReturnCeed(y), CEED_ERROR_INCOMPATIBLE, in CeedVectorAXPY() 810 if (y->AXPY) { in CeedVectorAXPY() [all …]
|