#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_odeiv.h>
#include "pism/verification/tests/exactTestM.h"

Macros
#define	MAX(a, b) (a > b ? a : b)

#define	MIN(a, b) (a > b ? b : a)

#define	SperA 31556926.0 /* seconds per year; 365.2422 days */

#define	g 9.81

#define	rho 910.0 /* ice density; kg/m^3 */

#define	rhow 1028.0 /* sea water density; kg/m^3 */

#define	n 3.0 /* Glen power */

#define	barB

#define	H0 500.0 /* m */

#define	Rg 300.0e3 /* m; 300 km */

#define	Rc 600.0e3 /* m; 600 km */

#define	NOT_DONE 8966

#define	INVALID_METHOD 8968

#define	NEGATIVE_R 8969

Functions
double	F_M (double x, double alpha, double r, double Q)

double	dF_M (double x, double alpha, double r, double Q)

int	funcM_ode_G (double r, const double alpha[], double f[], void *params)

int	exactM_old (double r, double alpha, double Drr, const double EPS_ABS, const double EPS_REL, const int ode_method)

struct TestMParameters	exactM (double r, double EPS_ABS, double EPS_REL, int ode_method)

Macros

Functions