#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)

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

#endif