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

#ifndef TRI_D_H
#define TRI_D_H

static fp_t DDNUM = 720.000000000000000; 
 
const fp_t T2D00_DP1[6][6] = {
{30.0000000000000000 ,30.0000000000000000 ,30.0000000000000000 ,120.0000000000000000 ,0.0000000000000000 ,-120.0000000000000000},
{30.0000000000000000 ,90.0000000000000000 ,30.0000000000000000 ,120.0000000000000000 ,0.0000000000000000 ,-120.0000000000000000},
{30.0000000000000000 ,30.0000000000000000 ,30.0000000000000000 ,120.0000000000000000 ,0.0000000000000000 ,-120.0000000000000000},
{120.0000000000000000 ,120.0000000000000000 ,120.0000000000000000 ,480.0000000000000000 ,0.0000000000000000 ,-480.0000000000000000},
{0.0000000000000000 ,0.0000000000000000 ,0.0000000000000000 ,0.0000000000000000 ,480.0000000000000000 ,0.0000000000000000},
{-120.0000000000000000 ,-120.0000000000000000 ,-120.0000000000000000 ,-480.0000000000000000 ,0.0000000000000000 ,480.0000000000000000}
};

// const fp_t T2D00_DP1[6][6] = {
// {1.0, 1.0, 5.0, 1.0, 1.0, 1.0},
// {1.0, 1.0, 1.0, 1.0, 1.0, 1.0},
// {1.0, 1.0, 1.0, 1.0, 1.0, 1.0},
// {1.0, 1.0, 1.0, 1.0, 1.0, 1.0},
// {1.0, 1.0, 1.0, 1.0, 1.0, 1.0},
// {1.0, 1.0, 1.0, 1.0, 1.0, 1.0}
// };

const fp_t T2D01_DP1[6][6] = { 
{45.0000000000000000 ,15.0000000000000000 ,-15.0000000000000000 ,60.0000000000000000 ,60.0000000000000000 ,-60.0000000000000000},
{75.0000000000000000 ,45.0000000000000000 ,-45.0000000000000000 ,-60.0000000000000000 ,180.0000000000000000 ,-180.0000000000000000},
{45.0000000000000000 ,15.0000000000000000 ,-15.0000000000000000 ,60.0000000000000000 ,60.0000000000000000 ,-60.0000000000000000},
{180.0000000000000000 ,60.0000000000000000 ,-60.0000000000000000 ,240.0000000000000000 ,240.0000000000000000 ,-240.0000000000000000},
{-60.0000000000000000 ,60.0000000000000000 ,-60.0000000000000000 ,240.0000000000000000 ,240.0000000000000000 ,-240.0000000000000000},
{-180.0000000000000000 ,-60.0000000000000000 ,60.0000000000000000 ,-240.0000000000000000 ,-240.0000000000000000 ,240.0000000000000000}
};

const fp_t T2D10_DP1[6][6] = { 
{45.0000000000000000 ,75.0000000000000000 ,45.0000000000000000 ,180.0000000000000000 ,-60.0000000000000000 ,-180.0000000000000000},
{15.0000000000000000 ,45.0000000000000000 ,15.0000000000000000 ,60.0000000000000000 ,60.0000000000000000 ,-60.0000000000000000},
{-15.0000000000000000 ,-45.0000000000000000 ,-15.0000000000000000 ,-60.0000000000000000 ,-60.0000000000000000 ,60.0000000000000000},
{60.0000000000000000 ,-60.0000000000000000 ,60.0000000000000000 ,240.0000000000000000 ,240.0000000000000000 ,-240.0000000000000000},
{60.0000000000000000 ,180.0000000000000000 ,60.0000000000000000 ,240.0000000000000000 ,240.0000000000000000 ,-240.0000000000000000},
{-60.0000000000000000 ,-180.0000000000000000 ,-60.0000000000000000 ,-240.0000000000000000 ,-240.0000000000000000 ,240.0000000000000000}
};

const fp_t T2D11_DP1[6][6] = {
{90.0000000000000000 ,30.0000000000000000 ,-30.0000000000000000 ,0.0000000000000000 ,120.0000000000000000 ,-120.0000000000000000},
{30.0000000000000000 ,30.0000000000000000 ,-30.0000000000000000 ,0.0000000000000000 ,120.0000000000000000 ,-120.0000000000000000},
{-30.0000000000000000 ,-30.0000000000000000 ,30.0000000000000000 ,0.0000000000000000 ,-120.0000000000000000 ,120.0000000000000000},
{0.0000000000000000 ,0.0000000000000000 ,0.0000000000000000 ,480.0000000000000000 ,0.0000000000000000 ,0.0000000000000000},
{120.0000000000000000 ,120.0000000000000000 ,-120.0000000000000000 ,0.0000000000000000 ,480.0000000000000000 ,-480.0000000000000000},
{-120.0000000000000000 ,-120.0000000000000000 ,120.0000000000000000 ,0.0000000000000000 ,-480.0000000000000000 ,480.0000000000000000}
};


const fp_t w_dot_nxw[6][6] 
= {-8.03966578630366e-18 ,    -0.0833333333333333 ,    -0.0833333333333333  ,    -0.166666666666667  ,   0.166666666666667  ,   0.166666666666667,
0.0833333333333333  , 2.2351612531664e-18  ,   -0.0833333333333333  ,    -0.166666666666667  ,   0.166666666666667  ,    -0.166666666666667,
0.0833333333333333 ,   0.0833333333333333 ,  -9.90790307867125e-18  ,   0.166666666666667 ,    0.166666666666667  ,    -0.166666666666667,
0.166666666666667  ,   0.166666666666667  ,    -0.166666666666667 ,  -5.87078337278095e-18 , 7.63278329429795e-17 ,  -7.63278329429795e-17,
-0.166666666666667  ,    -0.166666666666667   ,   -0.166666666666667 , 7.63278329429795e-17  , -5.87078337278095e-18 , 5.87078337278095e-18,
-0.166666666666667  ,   0.166666666666667   ,  0.166666666666667  , -7.63278329429795e-17 , 5.87078337278095e-18 ,  -5.87078337278095e-18};



// /* 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)}
// };

#endif