Actual source code: petsc_p4est_package.c

petsc-3.9.4 2018-09-11
Report Typos and Errors
  1: #include <petscdmforest.h>
  2: #include <petsc/private/petscimpl.h>
  3: #include <p4est_base.h>
  4: #include "petsc_p4est_package.h"

  6: static const char*const SCLogTypes[] = {"DEFAULT","ALWAYS","TRACE","DEBUG","VERBOSE","INFO","STATISTICS","PRODUCTION","ESSENTIAL","ERROR","SILENT","SCLogTypes","SC_LP_",0};

  8: static PetscBool    PetscP4estInitialized = PETSC_FALSE;
  9: static PetscBool    PetscBeganSc          = PETSC_FALSE;
 10: static PetscClassId P4ESTLOGGING_CLASSID;

 12: PetscObject P4estLoggingObject; /* Just a vehicle for its classid */

 14: static void PetscScLogHandler(FILE *log_stream, const char *filename, int lineno,int package, int category,int priority, const char *msg)
 15: {
 16:   PetscInfo_Private(filename,P4estLoggingObject,":%d{%s} %s",lineno,package == sc_package_id ? "sc" : package == p4est_package_id ? "p4est" : "",msg);
 17: }

 19: /* p4est tries to abort: if possible, use setjmp to enable at least a little unwinding */
 20: #if defined(PETSC_HAVE_SETJMP_H) && defined(PETSC_USE_ERRORCHECKING)
 21: #include <setjmp.h>
 22: PETSC_VISIBILITY_INTERNAL jmp_buf PetscScJumpBuf;
 23: PETSC_INTERN void PetscScAbort_longjmp(void)
 24: {
 25:   PetscError(PETSC_COMM_SELF,-1,"p4est function","p4est file",PETSC_ERR_LIB,PETSC_ERROR_INITIAL,"Error in p4est stack call\n");
 26:   longjmp(PetscScJumpBuf,1);
 27:   return;
 28: }

 30: #define PetscScAbort PetscScAbort_longjmp
 31: #else
 32: #define PetscScAbort NULL
 33: #endif

 35: static PetscErrorCode PetscP4estFinalize(void)
 36: {

 40:   if (PetscBeganSc) {
 41:     /* We do not want libsc to abort on a mismatched allocation and prevent further Petsc unwinding */
 42:     PetscStackCallP4est(sc_package_set_abort_alloc_mismatch,(sc_package_id,0));
 43:     PetscStackCallP4est(sc_package_set_abort_alloc_mismatch,(p4est_package_id,0));
 44:     PetscStackCallP4est(sc_package_set_abort_alloc_mismatch,(-1,0));
 45:     PetscStackCallP4est(sc_finalize,());
 46:   }
 47:   PetscHeaderDestroy(&P4estLoggingObject);
 48:   return(0);
 49: }

 51: PetscErrorCode PetscP4estInitialize(void)
 52: {
 53:   PetscBool      psc_catch_signals    = PETSC_FALSE;
 54:   PetscBool      psc_print_backtrace  = PETSC_TRUE;
 55:   int            psc_log_threshold    = SC_LP_DEFAULT;
 56:   int            pp4est_log_threshold = SC_LP_DEFAULT;
 57:   char           logList[256];
 58:   PetscBool      opt,pkg;

 62:   if (PetscP4estInitialized) return(0);
 63:   PetscP4estInitialized = PETSC_TRUE;

 65:   /* Register Classes */
 66:   PetscClassIdRegister("p4est logging",&P4ESTLOGGING_CLASSID);
 67:   /* Process info exclusions */
 68:   PetscOptionsGetString(NULL,NULL,"-info_exclude",logList,sizeof(logList),&opt);
 69:   if (opt) {
 70:     PetscStrInList("p4est",logList,',',&pkg);
 71:     if (pkg) {PetscInfoDeactivateClass(P4ESTLOGGING_CLASSID);}
 72:   }
 73:   /* Process summary exclusions */
 74:   PetscOptionsGetString(NULL,NULL,"-log_exclude",logList,sizeof(logList),&opt);
 75:   if (opt) {
 76:     PetscStrInList("p4est",logList,',',&pkg);
 77:     if (pkg) {PetscLogEventDeactivateClass(P4ESTLOGGING_CLASSID);}
 78:   }
 79:   PetscHeaderCreate(P4estLoggingObject,P4ESTLOGGING_CLASSID,"p4est","p4est logging","DM",PETSC_COMM_WORLD,NULL,PetscObjectView);
 80:   if (sc_package_id == -1) {
 81:     int       log_threshold_shifted = psc_log_threshold + 1;
 82:     PetscBool set;

 84:     PetscBeganSc = PETSC_TRUE;
 85:     PetscOptionsGetBool(NULL,NULL,"-petsc_sc_catch_signals",&psc_catch_signals,NULL);
 86:     PetscOptionsGetBool(NULL,NULL,"-petsc_sc_print_backtrace",&psc_print_backtrace,NULL);
 87:     PetscOptionsGetEnum(NULL,NULL,"-petsc_sc_log_threshold",SCLogTypes,(PetscEnum*)&log_threshold_shifted,&set);
 88:     if (set) psc_log_threshold = log_threshold_shifted - 1;
 89:     sc_init(PETSC_COMM_WORLD,(int)psc_catch_signals,(int)psc_print_backtrace,PetscScLogHandler,psc_log_threshold);
 90:     if (sc_package_id == -1) SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_LIB,"Could not initialize libsc package used by p4est");
 91:     sc_set_abort_handler(PetscScAbort);
 92:   }
 93:   if (p4est_package_id == -1) {
 94:     int       log_threshold_shifted = pp4est_log_threshold + 1;
 95:     PetscBool set;

 97:     PetscOptionsGetEnum(NULL,NULL,"-petsc_p4est_log_threshold",SCLogTypes,(PetscEnum*)&log_threshold_shifted,&set);
 98:     if (set) pp4est_log_threshold = log_threshold_shifted - 1;
 99:     PetscStackCallP4est(p4est_init,(PetscScLogHandler,pp4est_log_threshold));
100:     if (p4est_package_id == -1) SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_LIB,"Could not initialize p4est");
101:   }
102:   DMForestRegisterType(DMP4EST);
103:   DMForestRegisterType(DMP8EST);
104:   PetscRegisterFinalize(PetscP4estFinalize);
105:   return(0);
106: }