1 // Copyright (c) 2017-2024, Lawrence Livermore National Security, LLC and other CEED contributors. 2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details. 3 // 4 // SPDX-License-Identifier: BSD-2-Clause 5 // 6 // This file is part of CEED: http://github.com/ceed 7 #pragma once 8 9 #include <ceed.h> 10 #include <ceed/backend.h> 11 #include <ceed/jit-source/hip/hip-types.h> 12 #include <hip/hip_runtime.h> 13 #if (HIP_VERSION >= 50200000) 14 #include <hipblas/hipblas.h> // IWYU pragma: export 15 #else 16 #include <hipblas.h> // IWYU pragma: export 17 #endif 18 19 typedef struct { 20 CeedScalar *h_array; 21 CeedScalar *h_array_borrowed; 22 CeedScalar *h_array_owned; 23 CeedScalar *d_array; 24 CeedScalar *d_array_borrowed; 25 CeedScalar *d_array_owned; 26 } CeedVector_Hip; 27 28 typedef struct { 29 hipModule_t module; 30 hipFunction_t ApplyNoTranspose, ApplyTranspose; 31 hipFunction_t ApplyUnsignedNoTranspose, ApplyUnsignedTranspose; 32 hipFunction_t ApplyUnorientedNoTranspose, ApplyUnorientedTranspose; 33 CeedInt num_nodes; 34 const CeedInt *h_offsets; 35 const CeedInt *h_offsets_borrowed; 36 const CeedInt *h_offsets_owned; 37 const CeedInt *d_offsets; 38 const CeedInt *d_offsets_borrowed; 39 const CeedInt *d_offsets_owned; 40 const CeedInt *d_t_offsets; 41 const CeedInt *d_t_indices; 42 const CeedInt *d_l_vec_indices; 43 const bool *h_orients; 44 const bool *h_orients_borrowed; 45 const bool *h_orients_owned; 46 const bool *d_orients; 47 const bool *d_orients_borrowed; 48 const bool *d_orients_owned; 49 const CeedInt8 *h_curl_orients; 50 const CeedInt8 *h_curl_orients_borrowed; 51 const CeedInt8 *h_curl_orients_owned; 52 const CeedInt8 *d_curl_orients; 53 const CeedInt8 *d_curl_orients_borrowed; 54 const CeedInt8 *d_curl_orients_owned; 55 } CeedElemRestriction_Hip; 56 57 typedef struct { 58 hipModule_t module; 59 hipFunction_t Interp; 60 hipFunction_t Grad; 61 hipFunction_t Weight; 62 CeedScalar *d_interp_1d; 63 CeedScalar *d_grad_1d; 64 CeedScalar *d_q_weight_1d; 65 } CeedBasis_Hip; 66 67 typedef struct { 68 hipModule_t module; 69 hipFunction_t Interp; 70 hipFunction_t InterpTranspose; 71 hipFunction_t Deriv; 72 hipFunction_t DerivTranspose; 73 hipFunction_t Weight; 74 CeedScalar *d_interp; 75 CeedScalar *d_grad; 76 CeedScalar *d_div; 77 CeedScalar *d_curl; 78 CeedScalar *d_q_weight; 79 } CeedBasisNonTensor_Hip; 80 81 typedef struct { 82 hipModule_t module; 83 const char *qfunction_name; 84 const char *qfunction_source; 85 hipFunction_t QFunction; 86 Fields_Hip fields; 87 void *d_c; 88 } CeedQFunction_Hip; 89 90 typedef struct { 91 void *h_data; 92 void *h_data_borrowed; 93 void *h_data_owned; 94 void *d_data; 95 void *d_data_borrowed; 96 void *d_data_owned; 97 } CeedQFunctionContext_Hip; 98 99 typedef struct { 100 hipModule_t module, module_point_block; 101 hipFunction_t LinearDiagonal; 102 hipFunction_t LinearPointBlock; 103 CeedElemRestriction diag_rstr, point_block_diag_rstr; 104 CeedVector elem_diag, point_block_elem_diag; 105 CeedEvalMode *d_eval_modes_in, *d_eval_modes_out; 106 CeedScalar *d_identity, *d_interp_in, *d_grad_in, *d_div_in, *d_curl_in; 107 CeedScalar *d_interp_out, *d_grad_out, *d_div_out, *d_curl_out; 108 } CeedOperatorDiag_Hip; 109 110 typedef struct { 111 hipModule_t module; 112 hipFunction_t LinearAssemble; 113 CeedInt block_size_x, block_size_y, elems_per_block; 114 CeedScalar *d_B_in, *d_B_out; 115 } CeedOperatorAssemble_Hip; 116 117 typedef struct { 118 CeedVector *e_vecs; // E-vectors, inputs followed by outputs 119 CeedVector *q_vecs_in; // Input Q-vectors needed to apply operator 120 CeedVector *q_vecs_out; // Output Q-vectors needed to apply operator 121 CeedInt num_inputs, num_outputs; 122 CeedInt num_active_in, num_active_out; 123 CeedVector *qf_active_in; 124 CeedOperatorDiag_Hip *diag; 125 CeedOperatorAssemble_Hip *asmb; 126 } CeedOperator_Hip; 127 128 CEED_INTERN int CeedGetHipblasHandle_Hip(Ceed ceed, hipblasHandle_t *handle); 129 130 CEED_INTERN int CeedVectorCreate_Hip(CeedSize n, CeedVector vec); 131 132 CEED_INTERN int CeedElemRestrictionCreate_Hip(CeedMemType mem_type, CeedCopyMode copy_mode, const CeedInt *offsets, const bool *orients, 133 const CeedInt8 *curl_orients, CeedElemRestriction rstr); 134 135 CEED_INTERN int CeedBasisCreateTensorH1_Hip(CeedInt dim, CeedInt P_1d, CeedInt Q_1d, const CeedScalar *interp_1d, const CeedScalar *grad_1d, 136 const CeedScalar *q_ref_1d, const CeedScalar *q_weight_1d, CeedBasis basis); 137 CEED_INTERN int CeedBasisCreateH1_Hip(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts, const CeedScalar *interp, 138 const CeedScalar *grad, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis); 139 CEED_INTERN int CeedBasisCreateHdiv_Hip(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts, const CeedScalar *interp, 140 const CeedScalar *div, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis); 141 CEED_INTERN int CeedBasisCreateHcurl_Hip(CeedElemTopology topo, CeedInt dim, CeedInt num_nodes, CeedInt num_qpts, const CeedScalar *interp, 142 const CeedScalar *curl, const CeedScalar *q_ref, const CeedScalar *q_weight, CeedBasis basis); 143 144 CEED_INTERN int CeedQFunctionCreate_Hip(CeedQFunction qf); 145 146 CEED_INTERN int CeedQFunctionContextCreate_Hip(CeedQFunctionContext ctx); 147 148 CEED_INTERN int CeedOperatorCreate_Hip(CeedOperator op); 149