common.h

#ifndef __COMMON_H
#define __COMMON_H
#include <limits>
#include "trottersuzuki.h"

void print_matrix(string filename, double * matrix, size_t stride, size_t width, size_t height);
void stamp(Lattice *grid, State *state, string fileprefix);
void stamp_matrix(Lattice *grid, double *matrix, string filename);

void calculate_borders(int coord, int dim, int * start, int *end, int *inner_start, int *inner_end, int length, int halo, int periodic_bound);
void my_abort(string err);
void memcpy2D(void * dst, size_t dstride, const void * src, size_t sstride, size_t width, size_t height);
double bessel_j_zeros(int l, int x);

#endif