#include "../dgtd-performance.hpp"

static int DNUM = 720;
static int T2D00[8][8] = 
{
	{30,30,30,120,0,-120,48,24},
	{30,90,30,120,0,-120,72,36},
	{30,30,30,120,0,-120,48,24},
	{120,120,120,480,0,-480,192,96},
	{0,0,0,0,480,0,0,-48},
	{-120,-120,-120,-480,0,480,-192,-96},
	{48,72,48,192,0,-192,96,48},
	{24,36,24,96,-48,-96,48,32}
};

static int T2D01[8][8] = 
{
	{45,15,-15,60,60,-60,24,-12},
	{75,45,-45,-60,180,-180,36,-48},
	{45,15,-15,60,60,-60,24,-12},
	{180,60,-60,240,240,-240,96,-48},
	{-60,60,-60,240,240,-240,48,0},
	{-180,-60,60,-240,-240,240,-96,48},
	{84,36,-36,48,144,-144,48,-32},
	{48,12,-12,0,48,-48,16,-16}
};

static int T2D10[8][8] = 
{
	{45,75,45,180,-60,-180,84,48},
	{15,45,15,60,60,-60,36,12},
	{-15,-45,-15,-60,-60,60,-36,-12},
	{60,-60,60,240,240,-240,48,0},
	{60,180,60,240,240,-240,144,48},
	{-60,-180,-60,-240,-240,240,-144,-48},
	{24,36,24,96,48,-96,48,16},
	{-12,-48,-12,-48,0,48,-32,-16}
};

static int T2D11[8][8] = 
{
	{90,30,-30,0,120,-120,36,-36},
	{30,30,-30,0,120,-120,24,-24},
	{-30,-30,30,0,-120,120,-24,24},
	{0,0,0,480,0,0,48,48},
	{120,120,-120,0,480,-480,96,-96},
	{-120,-120,120,0,-480,480,-96,96},
	{36,24,-24,48,96,-96,32,-16},
	{-36,-24,24,48,-96,96,-16,32}
};

/* the divergence constant matrices */
static int DIV2D0[3][8] = 
{
	{0, 0, 0, -2, 0, 0, -1, -1},
	{0, 0, 0, -2, 0, 0,  0,  0},
	{0, 0, 0, -2, 0, 0,  0,  0}
};

static int DIV2D1[3][8] = 
{
	{0, 0, 0, 0, -2, 0, 0, 0},
	{0, 0, 0, 0, -2, 0, 0, 1},
	{0, 0, 0, 0, -2, 0, 0, 0}
};

static int DIV2D2[3][8] = 
{
	{0, 0, 0, 0, 0, -2, 0, 0},
	{0, 0, 0, 0, 0, -2, 0, 0},
	{0, 0, 0, 0, 0, -2, 1, 0}
};

static fp_t T2D[3][3] = 
{
	{ (1.0 / 6.0), (1.0 / 12.0), (1.0 / 12.0)},
	{ (1.0 / 12.0), (1.0 / 6.0), (1.0 / 12.0)},
	{ (1.0 / 12.0), (1.0 / 12.0), (1.0 / 6.0)}
};