1 // Copyright (c) 2017-2022, 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 8 #include <ceed.h> 9 #include <ceed/backend.h> 10 #include <ceed/jit-source/hip/hip-types.h> 11 #include <hip/hip_runtime.h> 12 13 #include "../hip/ceed-hip-common.h" 14 #include "../hip/ceed-hip-compile.h" 15 #include "ceed-hip-ref-qfunction-load.h" 16 #include "ceed-hip-ref.h" 17 18 //------------------------------------------------------------------------------ 19 // Apply QFunction 20 //------------------------------------------------------------------------------ 21 static int CeedQFunctionApply_Hip(CeedQFunction qf, CeedInt Q, CeedVector *U, CeedVector *V) { 22 Ceed ceed; 23 CeedCallBackend(CeedQFunctionGetCeed(qf, &ceed)); 24 25 // Build and compile kernel, if not done 26 CeedCallBackend(CeedQFunctionBuildKernel_Hip_ref(qf)); 27 28 CeedQFunction_Hip *data; 29 CeedCallBackend(CeedQFunctionGetData(qf, &data)); 30 Ceed_Hip *ceed_Hip; 31 CeedCallBackend(CeedGetData(ceed, &ceed_Hip)); 32 CeedInt num_input_fields, num_output_fields; 33 CeedCallBackend(CeedQFunctionGetNumArgs(qf, &num_input_fields, &num_output_fields)); 34 const int blocksize = ceed_Hip->opt_block_size; 35 36 // Read vectors 37 for (CeedInt i = 0; i < num_input_fields; i++) { 38 CeedCallBackend(CeedVectorGetArrayRead(U[i], CEED_MEM_DEVICE, &data->fields.inputs[i])); 39 } 40 for (CeedInt i = 0; i < num_output_fields; i++) { 41 CeedCallBackend(CeedVectorGetArrayWrite(V[i], CEED_MEM_DEVICE, &data->fields.outputs[i])); 42 } 43 44 // Get context data 45 CeedCallBackend(CeedQFunctionGetInnerContextData(qf, CEED_MEM_DEVICE, &data->d_c)); 46 47 // Run kernel 48 void *args[] = {&data->d_c, (void *)&Q, &data->fields}; 49 CeedCallBackend(CeedRunKernel_Hip(ceed, data->QFunction, CeedDivUpInt(Q, blocksize), blocksize, args)); 50 51 // Restore vectors 52 for (CeedInt i = 0; i < num_input_fields; i++) { 53 CeedCallBackend(CeedVectorRestoreArrayRead(U[i], &data->fields.inputs[i])); 54 } 55 for (CeedInt i = 0; i < num_output_fields; i++) { 56 CeedCallBackend(CeedVectorRestoreArray(V[i], &data->fields.outputs[i])); 57 } 58 59 // Restore context 60 CeedCallBackend(CeedQFunctionRestoreInnerContextData(qf, &data->d_c)); 61 62 return CEED_ERROR_SUCCESS; 63 } 64 65 //------------------------------------------------------------------------------ 66 // Destroy QFunction 67 //------------------------------------------------------------------------------ 68 static int CeedQFunctionDestroy_Hip(CeedQFunction qf) { 69 CeedQFunction_Hip *data; 70 CeedCallBackend(CeedQFunctionGetData(qf, &data)); 71 Ceed ceed; 72 CeedCallBackend(CeedQFunctionGetCeed(qf, &ceed)); 73 if (data->module) CeedCallHip(ceed, hipModuleUnload(data->module)); 74 CeedCallBackend(CeedFree(&data)); 75 76 return CEED_ERROR_SUCCESS; 77 } 78 79 //------------------------------------------------------------------------------ 80 // Create QFunction 81 //------------------------------------------------------------------------------ 82 int CeedQFunctionCreate_Hip(CeedQFunction qf) { 83 Ceed ceed; 84 CeedQFunctionGetCeed(qf, &ceed); 85 CeedQFunction_Hip *data; 86 CeedCallBackend(CeedCalloc(1, &data)); 87 CeedCallBackend(CeedQFunctionSetData(qf, data)); 88 CeedInt num_input_fields, num_output_fields; 89 CeedCallBackend(CeedQFunctionGetNumArgs(qf, &num_input_fields, &num_output_fields)); 90 91 // Read QFunction source 92 CeedCallBackend(CeedQFunctionGetKernelName(qf, &data->qfunction_name)); 93 CeedDebug256(ceed, 2, "----- Loading QFunction User Source -----\n"); 94 CeedCallBackend(CeedQFunctionLoadSourceToBuffer(qf, &data->qfunction_source)); 95 CeedDebug256(ceed, 2, "----- Loading QFunction User Source Complete! -----\n"); 96 97 // Register backend functions 98 CeedCallBackend(CeedSetBackendFunction(ceed, "QFunction", qf, "Apply", CeedQFunctionApply_Hip)); 99 CeedCallBackend(CeedSetBackendFunction(ceed, "QFunction", qf, "Destroy", CeedQFunctionDestroy_Hip)); 100 return CEED_ERROR_SUCCESS; 101 } 102 103 //------------------------------------------------------------------------------ 104