You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
93 lines
5.1 KiB
93 lines
5.1 KiB
#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
|
|
|