11447629fSBarry Smith 21447629fSBarry Smith #include <../src/vec/is/ao/aoimpl.h> /*I "petscao.h" I*/ 31447629fSBarry Smith 4*da8c939bSJacob Faibussowitsch static PetscBool AOPackageInitialized = PETSC_FALSE; 5*da8c939bSJacob Faibussowitsch static PetscBool AORegisterAllCalled = PETSC_FALSE; 6*da8c939bSJacob Faibussowitsch 7*da8c939bSJacob Faibussowitsch /*@C 8*da8c939bSJacob Faibussowitsch AOFinalizePackage - This function finalizes everything in the `AO` package. It is called 9*da8c939bSJacob Faibussowitsch from `PetscFinalize()`. 10*da8c939bSJacob Faibussowitsch 11*da8c939bSJacob Faibussowitsch Level: developer 12*da8c939bSJacob Faibussowitsch 13*da8c939bSJacob Faibussowitsch .seealso: `AOInitializePackage()`, `PetscInitialize()` 14*da8c939bSJacob Faibussowitsch @*/ 15*da8c939bSJacob Faibussowitsch PetscErrorCode AOFinalizePackage(void) 16*da8c939bSJacob Faibussowitsch { 17*da8c939bSJacob Faibussowitsch PetscFunctionBegin; 18*da8c939bSJacob Faibussowitsch PetscCall(PetscFunctionListDestroy(&AOList)); 19*da8c939bSJacob Faibussowitsch AOPackageInitialized = PETSC_FALSE; 20*da8c939bSJacob Faibussowitsch AORegisterAllCalled = PETSC_FALSE; 21*da8c939bSJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 22*da8c939bSJacob Faibussowitsch } 23*da8c939bSJacob Faibussowitsch 24*da8c939bSJacob Faibussowitsch /*@C 25*da8c939bSJacob Faibussowitsch AOInitializePackage - This function initializes everything in the `AO` package. It is called 26*da8c939bSJacob Faibussowitsch from `PetscDLLibraryRegister_petscvec()` when using dynamic libraries, and on the first call to `AOCreate()` 27*da8c939bSJacob Faibussowitsch when using static or shared libraries. 28*da8c939bSJacob Faibussowitsch 29*da8c939bSJacob Faibussowitsch Level: developer 30*da8c939bSJacob Faibussowitsch 31*da8c939bSJacob Faibussowitsch .seealso: `AOFinalizePackage()`, `PetscInitialize()` 32*da8c939bSJacob Faibussowitsch @*/ 33*da8c939bSJacob Faibussowitsch PetscErrorCode AOInitializePackage(void) 34*da8c939bSJacob Faibussowitsch { 35*da8c939bSJacob Faibussowitsch char logList[256]; 36*da8c939bSJacob Faibussowitsch PetscBool opt, pkg; 37*da8c939bSJacob Faibussowitsch 38*da8c939bSJacob Faibussowitsch PetscFunctionBegin; 39*da8c939bSJacob Faibussowitsch if (AOPackageInitialized) PetscFunctionReturn(PETSC_SUCCESS); 40*da8c939bSJacob Faibussowitsch AOPackageInitialized = PETSC_TRUE; 41*da8c939bSJacob Faibussowitsch /* Register Classes */ 42*da8c939bSJacob Faibussowitsch PetscCall(PetscClassIdRegister("Application Order", &AO_CLASSID)); 43*da8c939bSJacob Faibussowitsch /* Register Constructors */ 44*da8c939bSJacob Faibussowitsch PetscCall(AORegisterAll()); 45*da8c939bSJacob Faibussowitsch /* Register Events */ 46*da8c939bSJacob Faibussowitsch PetscCall(PetscLogEventRegister("AOPetscToApplication", AO_CLASSID, &AO_PetscToApplication)); 47*da8c939bSJacob Faibussowitsch PetscCall(PetscLogEventRegister("AOApplicationToPetsc", AO_CLASSID, &AO_ApplicationToPetsc)); 48*da8c939bSJacob Faibussowitsch /* Process Info */ 49*da8c939bSJacob Faibussowitsch { 50*da8c939bSJacob Faibussowitsch PetscClassId classids[1]; 51*da8c939bSJacob Faibussowitsch 52*da8c939bSJacob Faibussowitsch classids[0] = AO_CLASSID; 53*da8c939bSJacob Faibussowitsch PetscCall(PetscInfoProcessClass("ao", 1, classids)); 54*da8c939bSJacob Faibussowitsch } 55*da8c939bSJacob Faibussowitsch /* Process summary exclusions */ 56*da8c939bSJacob Faibussowitsch PetscCall(PetscOptionsGetString(NULL, NULL, "-log_exclude", logList, sizeof(logList), &opt)); 57*da8c939bSJacob Faibussowitsch if (opt) { 58*da8c939bSJacob Faibussowitsch PetscCall(PetscStrInList("ao", logList, ',', &pkg)); 59*da8c939bSJacob Faibussowitsch if (pkg) PetscCall(PetscLogEventExcludeClass(AO_CLASSID)); 60*da8c939bSJacob Faibussowitsch } 61*da8c939bSJacob Faibussowitsch /* Register package finalizer */ 62*da8c939bSJacob Faibussowitsch PetscCall(PetscRegisterFinalize(AOFinalizePackage)); 63*da8c939bSJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 64*da8c939bSJacob Faibussowitsch } 651447629fSBarry Smith 661447629fSBarry Smith /*@C 67cab54364SBarry Smith AOSetType - Builds an application ordering for a particular `AOType` 681447629fSBarry Smith 69c3339decSBarry Smith Collective 701447629fSBarry Smith 711447629fSBarry Smith Input Parameters: 72cab54364SBarry Smith + ao - The `AO` object 731447629fSBarry Smith - method - The name of the AO type 741447629fSBarry Smith 751447629fSBarry Smith Options Database Key: 76cab54364SBarry Smith . -ao_type <type> - Sets the `AO` type; use -help for a list of available types 771447629fSBarry Smith 781447629fSBarry Smith Level: intermediate 791447629fSBarry Smith 80cab54364SBarry Smith Notes: 81cab54364SBarry Smith See "petsc/include/petscao.h" for available AO types (for instance, `AOBASIC` and `AOMEMORYSCALABLE`). 82cab54364SBarry Smith 83cab54364SBarry Smith `AO` are usually created via the convenience routines such as `AOCreateBasic()` or `AOCreateMemoryScalable()` 84cab54364SBarry Smith 85cab54364SBarry Smith .seealso: `AO`, `AOType`, `AOCreateBasic()`, `AOCreateMemoryScalable()`, `AOGetType()`, `AOCreate()` 861447629fSBarry Smith @*/ 87d71ae5a4SJacob Faibussowitsch PetscErrorCode AOSetType(AO ao, AOType method) 88d71ae5a4SJacob Faibussowitsch { 891447629fSBarry Smith PetscErrorCode (*r)(AO); 901447629fSBarry Smith PetscBool match; 911447629fSBarry Smith 921447629fSBarry Smith PetscFunctionBegin; 931447629fSBarry Smith PetscValidHeaderSpecific(ao, AO_CLASSID, 1); 949566063dSJacob Faibussowitsch PetscCall(PetscObjectTypeCompare((PetscObject)ao, method, &match)); 953ba16761SJacob Faibussowitsch if (match) PetscFunctionReturn(PETSC_SUCCESS); 961447629fSBarry Smith 979566063dSJacob Faibussowitsch PetscCall(AORegisterAll()); 989566063dSJacob Faibussowitsch PetscCall(PetscFunctionListFind(AOList, method, &r)); 996adde796SStefano Zampini PetscCheck(r, PetscObjectComm((PetscObject)ao), PETSC_ERR_ARG_UNKNOWN_TYPE, "Unknown AO type: %s", method); 100dbbe0bcdSBarry Smith PetscTryTypeMethod(ao, destroy); 1011447629fSBarry Smith ao->ops->destroy = NULL; 1021447629fSBarry Smith 1039566063dSJacob Faibussowitsch PetscCall((*r)(ao)); 1043ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 1051447629fSBarry Smith } 1061447629fSBarry Smith 1071447629fSBarry Smith /*@C 108cab54364SBarry Smith AOGetType - Gets the `AO` type name (as a string) from the AO. 1091447629fSBarry Smith 1101447629fSBarry Smith Not Collective 1111447629fSBarry Smith 1121447629fSBarry Smith Input Parameter: 1131447629fSBarry Smith . ao - The vector 1141447629fSBarry Smith 1151447629fSBarry Smith Output Parameter: 116cab54364SBarry Smith . type - The `AO` type name 1171447629fSBarry Smith 1181447629fSBarry Smith Level: intermediate 1191447629fSBarry Smith 120cab54364SBarry Smith .seealso: `AO`, `AOType`, `AOSetType()`, `AOCreate()` 1211447629fSBarry Smith @*/ 122d71ae5a4SJacob Faibussowitsch PetscErrorCode AOGetType(AO ao, AOType *type) 123d71ae5a4SJacob Faibussowitsch { 1241447629fSBarry Smith PetscFunctionBegin; 1251447629fSBarry Smith PetscValidHeaderSpecific(ao, AO_CLASSID, 1); 1264f572ea9SToby Isaac PetscAssertPointer(type, 2); 1279566063dSJacob Faibussowitsch PetscCall(AORegisterAll()); 1281447629fSBarry Smith *type = ((PetscObject)ao)->type_name; 1293ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 1301447629fSBarry Smith } 1311447629fSBarry Smith 1321447629fSBarry Smith /*--------------------------------------------------------------------------------------------------------------------*/ 1331447629fSBarry Smith 134*da8c939bSJacob Faibussowitsch PetscFunctionList AOList = NULL; 135*da8c939bSJacob Faibussowitsch 1361447629fSBarry Smith /*@C 1371d36bdfdSBarry Smith AORegister - Register an application ordering method 1381d36bdfdSBarry Smith 1391d36bdfdSBarry Smith Not Collective 1401d36bdfdSBarry Smith 1411d36bdfdSBarry Smith Input Parameters: 142cab54364SBarry Smith + sname - the name (`AOType`) of the `AO` scheme 1431d36bdfdSBarry Smith - function - the create routine for the application ordering method 1441447629fSBarry Smith 1451447629fSBarry Smith Level: advanced 1461d36bdfdSBarry Smith 147cab54364SBarry Smith .seealso: `AO`, `AOType`, `AOCreate()`, `AORegisterAll()`, `AOBASIC`, `AOADVANCED`, `AOMAPPING`, `AOMEMORYSCALABLE` 1481447629fSBarry Smith @*/ 149d71ae5a4SJacob Faibussowitsch PetscErrorCode AORegister(const char sname[], PetscErrorCode (*function)(AO)) 150d71ae5a4SJacob Faibussowitsch { 1511447629fSBarry Smith PetscFunctionBegin; 1529566063dSJacob Faibussowitsch PetscCall(AOInitializePackage()); 1539566063dSJacob Faibussowitsch PetscCall(PetscFunctionListAdd(&AOList, sname, function)); 1543ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 1551447629fSBarry Smith } 156*da8c939bSJacob Faibussowitsch 157*da8c939bSJacob Faibussowitsch PETSC_INTERN PetscErrorCode AOCreate_Basic(AO ao); 158*da8c939bSJacob Faibussowitsch PETSC_INTERN PetscErrorCode AOCreate_MemoryScalable(AO ao); 159*da8c939bSJacob Faibussowitsch 160*da8c939bSJacob Faibussowitsch /*@C 161*da8c939bSJacob Faibussowitsch AORegisterAll - Registers all of the application ordering components in the `AO` package. 162*da8c939bSJacob Faibussowitsch 163*da8c939bSJacob Faibussowitsch Not Collective 164*da8c939bSJacob Faibussowitsch 165*da8c939bSJacob Faibussowitsch Level: advanced 166*da8c939bSJacob Faibussowitsch 167*da8c939bSJacob Faibussowitsch .seealso: `AO`, `AOType`, `AORegister()`, `AORegisterDestroy()` 168*da8c939bSJacob Faibussowitsch @*/ 169*da8c939bSJacob Faibussowitsch PetscErrorCode AORegisterAll(void) 170*da8c939bSJacob Faibussowitsch { 171*da8c939bSJacob Faibussowitsch PetscFunctionBegin; 172*da8c939bSJacob Faibussowitsch if (AORegisterAllCalled) PetscFunctionReturn(PETSC_SUCCESS); 173*da8c939bSJacob Faibussowitsch AORegisterAllCalled = PETSC_TRUE; 174*da8c939bSJacob Faibussowitsch 175*da8c939bSJacob Faibussowitsch PetscCall(AORegister(AOBASIC, AOCreate_Basic)); 176*da8c939bSJacob Faibussowitsch PetscCall(AORegister(AOMEMORYSCALABLE, AOCreate_MemoryScalable)); 177*da8c939bSJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 178*da8c939bSJacob Faibussowitsch } 179