1*d636dbe3SBarry Smith /* $Id: petsc.h,v 1.56 1995/08/22 02:23:57 bsmith Exp bsmith $ */ 22eac72dbSBarry Smith 32eac72dbSBarry Smith #if !defined(__PETSC_PACKAGE) 42eac72dbSBarry Smith #define __PETSC_PACKAGE 52eac72dbSBarry Smith 6e0dd00f9SBarry Smith #define PETSC_VERSION_NUMBER "PETSc Version 2.0.Beta.7 Released ?, 1995." 77e59f0e8SBarry Smith 8f0479e8cSBarry Smith #include <stdio.h> 928988994SBarry Smith #if defined(PARCH_sun4) 106abc6512SBarry Smith int fprintf(FILE*,char*,...); 116abc6512SBarry Smith int printf(char*,...); 126abc6512SBarry Smith int fflush(FILE*); 136abc6512SBarry Smith int fclose(FILE*); 1428988994SBarry Smith #endif 15f0479e8cSBarry Smith 16f0479e8cSBarry Smith /* MPI interface */ 17f0479e8cSBarry Smith #include "mpi.h" 187f813858SBarry Smith #include "mpiu.h" 19f0479e8cSBarry Smith 20c6a45a97SBarry Smith #if defined(PETSC_COMPLEX) 21c6a45a97SBarry Smith /* work around for bug in alpha g++ compiler */ 22c6a45a97SBarry Smith #if defined(PARCH_alpha) 23c6a45a97SBarry Smith #define hypot(a,b) (double) sqrt((a)*(a)+(b)*(b)) 24c6a45a97SBarry Smith /* extern double hypot(double,double); */ 25c6a45a97SBarry Smith #endif 268ed539a5SBarry Smith #include <complex.h> 2720563c6bSBarry Smith #define PETSCREAL(a) real(a) 288ed539a5SBarry Smith #define Scalar complex 298ed539a5SBarry Smith #else 3020563c6bSBarry Smith #define PETSCREAL(a) a 318ed539a5SBarry Smith #define Scalar double 328ed539a5SBarry Smith #endif 338ed539a5SBarry Smith 34d90ea2a1SBarry Smith extern void *(*PetscMalloc)(unsigned int,int,char*); 35d90ea2a1SBarry Smith extern int (*PetscFree)(void *,int,char*); 3678b31e54SBarry Smith #define PETSCMALLOC(a) (*PetscMalloc)(a,__LINE__,__FILE__) 3778b31e54SBarry Smith #define PETSCFREE(a) (*PetscFree)(a,__LINE__,__FILE__) 387f813858SBarry Smith extern int PetscSetMalloc(void *(*)(unsigned int,int,char*), 397f813858SBarry Smith int (*)(void *,int,char*)); 40c7485abaSBarry Smith extern int TrDump(FILE *); 41e9f47e53SBarry Smith extern int TrGetMaximumAllocated(double*); 428ed539a5SBarry Smith 4378b31e54SBarry Smith #define PETSCNEW(A) (A*) PETSCMALLOC(sizeof(A)) 4478b31e54SBarry Smith #define PETSCMEMCPY(a,b,n) memcpy((char*)(a),(char*)(b),n) 4578b31e54SBarry Smith #define PETSCMEMSET(a,b,n) memset((char*)(a),(int)(b),n) 4620563c6bSBarry Smith #include <memory.h> 472eac72dbSBarry Smith 48145fe58dSBarry Smith #define PETSCMIN(a,b) ( ((a)<(b)) ? (a) : (b) ) 49145fe58dSBarry Smith #define PETSCMAX(a,b) ( ((a)<(b)) ? (b) : (a) ) 5081745bc7SBarry Smith #define PETSCABS(a) ( ((a)<0) ? -(a) : (a) ) 51145fe58dSBarry Smith 522eac72dbSBarry Smith /* Macros for error checking */ 53f0479e8cSBarry Smith #if !defined(__DIR__) 54f0479e8cSBarry Smith #define __DIR__ 0 55f0479e8cSBarry Smith #endif 56*d636dbe3SBarry Smith 57*d636dbe3SBarry Smith /* 58*d636dbe3SBarry Smith Unable to malloc error and no supported function 59*d636dbe3SBarry Smith */ 60*d636dbe3SBarry Smith #define PETSC_ERR_MEM 55 61*d636dbe3SBarry Smith #define PETSC_ERR_SUP 56 62*d636dbe3SBarry Smith 6328988994SBarry Smith #if defined(PETSC_DEBUG) 6478b31e54SBarry Smith #define SETERRQ(n,s) {return PetscError(__LINE__,__DIR__,__FILE__,s,n);} 65d6dfbf8fSBarry Smith #define SETERRA(n,s) \ 66d6dfbf8fSBarry Smith {int _ierr = PetscError(__LINE__,__DIR__,__FILE__,s,n);\ 67d6dfbf8fSBarry Smith MPI_Abort(MPI_COMM_WORLD,_ierr);} 6878b31e54SBarry Smith #define CHKERRQ(n) {if (n) SETERRQ(n,(char *)0);} 69d6dfbf8fSBarry Smith #define CHKERRA(n) {if (n) SETERRA(n,(char *)0);} 70*d636dbe3SBarry Smith #define CHKPTRQ(p) if (!p) SETERRQ(PETSC_ERR_MEM,(char*)0); 71*d636dbe3SBarry Smith #define CHKPTRA(p) if (!p) SETERRA(PETSC_ERR_MEM,(char*)0); 7228988994SBarry Smith #else 7378b31e54SBarry Smith #define SETERRQ(n,s) {return PetscError(__LINE__,__DIR__,__FILE__,s,n);} 7428988994SBarry Smith #define SETERRA(n,s) \ 7528988994SBarry Smith {int _ierr = PetscError(__LINE__,__DIR__,__FILE__,s,n);\ 7628988994SBarry Smith MPI_Abort(MPI_COMM_WORLD,_ierr);} 7778b31e54SBarry Smith #define CHKERRQ(n) {if (n) SETERRQ(n,(char *)0);} 7828988994SBarry Smith #define CHKERRA(n) {if (n) SETERRA(n,(char *)0);} 79*d636dbe3SBarry Smith #define CHKPTRQ(p) if (!p) SETERRQ(PETSC_ERR_MEM,(char*)0); 80*d636dbe3SBarry Smith #define CHKPTRA(p) if (!p) SETERRA(PETSC_ERR_MEM,(char*)0); 8128988994SBarry Smith #endif 822eac72dbSBarry Smith 838ed539a5SBarry Smith typedef struct _PetscObject* PetscObject; 841a941147SBarry Smith #define PETSC_COOKIE 1211211 856b5873e3SBarry Smith #define PETSC_DECIDE -1 8655b5a45fSLois Curfman McInnes #define PETSC_DEFAULT -2 87da69df5fSBarry Smith 88ca9b4cbeSLois Curfman McInnes typedef enum { PETSC_FALSE, PETSC_TRUE } PetscTruth; 89ca9b4cbeSLois Curfman McInnes 90da69df5fSBarry Smith #include "viewer.h" 918c8d4905SLois Curfman McInnes #include "options.h" 922eac72dbSBarry Smith 93a0a59b22SBarry Smith extern int PetscInitialize(int*,char***,char*,char*); 94a0a59b22SBarry Smith extern int PetscFinalize(); 95a0a59b22SBarry Smith 96c60448a5SBarry Smith extern int PetscObjectDestroy(PetscObject); 978f4c47bcSLois Curfman McInnes extern int PetscObjectExists(PetscObject,int*); 987f223b93SBarry Smith extern int PetscObjectGetComm(PetscObject,MPI_Comm *comm); 99a5a9c739SBarry Smith extern int PetscObjectSetName(PetscObject,char*); 100a5a9c739SBarry Smith extern int PetscObjectGetName(PetscObject,char**); 1012eac72dbSBarry Smith 102f0479e8cSBarry Smith extern int PetscDefaultErrorHandler(int,char*,char*,char*,int,void*); 103f0479e8cSBarry Smith extern int PetscAbortErrorHandler(int,char*,char*,char*,int,void* ); 104f0479e8cSBarry Smith extern int PetscAttachDebuggerErrorHandler(int,char*,char*,char*,int,void*); 105f0479e8cSBarry Smith extern int PetscError(int,char*,char*,char*,int); 106f0479e8cSBarry Smith extern int PetscPushErrorHandler(int 107f0479e8cSBarry Smith (*handler)(int,char*,char*,char*,int,void*),void* ); 1088ed539a5SBarry Smith extern int PetscPopErrorHandler(); 1098ed539a5SBarry Smith 1108ed539a5SBarry Smith extern int PetscSetDebugger(char *,int,char *); 1118ed539a5SBarry Smith extern int PetscAttachDebugger(); 1128ed539a5SBarry Smith 113f0479e8cSBarry Smith extern int PetscDefaultSignalHandler(int,void*); 114f0479e8cSBarry Smith extern int PetscPushSignalHandler(int (*)(int,void *),void*); 115f0479e8cSBarry Smith extern int PetscPopSignalHandler(); 116f0479e8cSBarry Smith extern int PetscSetFPTrap(int); 1176b5873e3SBarry Smith #define FP_TRAP_OFF 0 1186b5873e3SBarry Smith #define FP_TRAP_ON 1 119fba13676SBarry Smith #define FP_TRAP_ALWAYS 2 1208ed539a5SBarry Smith 12131b41497SLois Curfman McInnes /* 12231b41497SLois Curfman McInnes Definitions used for the Fortran interface: 12331b41497SLois Curfman McInnes FORTRANCAPS: Names are uppercase, no trailing underscore 12431b41497SLois Curfman McInnes FORTRANUNDERSCORE: Names are lowercase, trailing underscore 12531b41497SLois Curfman McInnes */ 126602c2b74SBarry Smith #if defined(PARCH_cray) || defined(PARCH_NCUBE) || defined(PARCH_t3d) 1278ed539a5SBarry Smith #define FORTRANCAPS 128d90ea2a1SBarry Smith #elif !defined(PARCH_rs6000) && !defined(PARCH_NeXT) && !defined(PARCH_hpux) 1298ed539a5SBarry Smith #define FORTRANUNDERSCORE 1302eac72dbSBarry Smith #endif 1312eac72dbSBarry Smith 13219b02663SBarry Smith #include "phead.h" 13319b02663SBarry Smith #include "plog.h" 1347857610eSBarry Smith 1357ec12d72SBarry Smith extern void PetscSleep(int); 1367ec12d72SBarry Smith 1372eac72dbSBarry Smith #endif 138