/*
 * self-defined petsc.h header for provesa examples: ex2a.c, ex2b.c,ex2c.c,ex2d.c.
 * https://repo.anl-external.org/repos/provesa/codes/mxm/
 */

#ifndef _PETSC_
#define _PETSC_

#include <math.h>
#define PetscExpScalar(x)     exp(x)

/* -------Types------- */

/* PetscErrorCode - datatype used for return error code from almost all PETSc functions */
typedef int PetscErrorCode;

/* PETSc type that represents a PetscReal.
 * This is the same as a PetscReal except in code that is automatically differentiated it 
 * is treated as a constant (not an indendent or dependent variable) 
 */
typedef double PassiveReal;

/* PETSc type that represents integer - used primarily to represent size of arrays and indexing
 * into arrays. Its size can be * configured with the option-with-64-bit-indices - to be either 
 * 32bit or 64bit [default 32 bit ints]
 */
typedef int PetscInt;

/* PetscReal - PETSc type that represents a real number version of PetscScalar */
typedef double PetscReal;

/* PetscScalar - PETSc type that represents either a double precision real number, a double precision
 * complex number, a single precision real number, a long double or an int - if the code is configured
 * with --with-scalar-type=real, complex --with-precision=single,double,__float128
 */
typedef double PetscScalar;

/* C struct that contains information about a structured grid and a processors logical location in it. */
typedef struct DMDALocalInfo {
  // PetscInt         dim,dof,sw;
  PetscInt         mx,my,mz;    /* global number of grid points in each direction */
  PetscInt         xs,ys,zs;    /* starting point of this processor, excluding ghosts */
  PetscInt         xm,ym,zm;    /* number of grid points on this processor, excluding ghosts */
  //PetscInt         gxs,gys,gzs;    /* starting point of this processor including ghosts */
  //PetscInt         gxm,gym,gzm;    /* number of grid points on this processor including ghosts */
} DMDALocalInfo;

/* ----------Functions---------- */

/* First executable line of each PETSc function, used for error handling. Final line of PETSc 
 * functions should be PetscFunctionReturn(0);
 */
#define PetscFunctionBegin

/* Last executable line of each PETSc function used for error handling. Replaces return() */
//void PetscFunctionReturn(0);
#define PetscFunctionReturn(x)    return(x)

/* Adds floating point operations to the global counter. */
PetscErrorCode PetscLogFlops(double );

/* Checks error code, if non-zero it calls the error handler and then returns */
PetscErrorCode CHKERRQ(PetscErrorCode );

/*
#undef __FUNCT__
#define __FUNCT__ "PetscLogFlops"
PETSC_STATIC_INLINE PetscErrorCode PetscLogFlops(PetscLogDouble n)
{
  PetscFunctionBegin;
#if defined(PETSC_USE_DEBUG)
  if (n < 0) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Cannot log negative flops");
#endif
  petsc_TotalFlops += PETSC_FLOPS_PER_OP*n;
  PetscFunctionReturn(0);
}
*/

#endif