Serial md code

/* a few global setting for all of pardy */
#ifndef GLOBALH
#define GLOBALH
/* global_rank and global_size should be set, once, at the start of main. */
#define DIM 3
typedef signed char tiny_int;
static inline double max3(double x1, double x2, double x3)
{
    /* Helper function to find maximum of three double precision numbers */
    double m12 = (x1>x2) ? x1 : x2;
    return (m12>x3) ? m12 : x3;
}
#endif