#ifndef SOR_H
#define SOR_H
#include "main.h"

#define F(v, a, b, c) *((v) + ((a)*NPTS*NPTS) + ((b)*NPTS) + (c))
#define Dd(v, a, b) *(v + (a*NPTS) + b)
/* v[a][b][c] */
/* (double) *(v + (a*NPTS*NPTS) + (b*NPTS) + c) */
/* *(*((*(v+a))+b)+c) */


fpt relax(double a, fpt *f, int *b);
fpt charge(fpt *f, int *b, int o);

#endif