void sine_gordon_exact ( double a, double x, double y, double *u, double *uxy );
double sine_gordon_residual ( double u, double uxy );