blob: 48e490e760904611807ddc6a73fbb51229f68888 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
// fd_grid.cc -- grid with finite differencing operations
// $Header$
//
// fd_grid::dx_coeff
// fd_grid::dxx_coeff
//
#include <stdio.h>
#include <assert.h>
#include <math.h>
#ifdef STANDALONE_TEST
#include "fake_cctk.h"
#else
#include "cctk.h"
#endif
#include "config.h"
#include "stdc.h"
#include "../jtutil/util.hh"
#include "../jtutil/array.hh"
#include "../jtutil/linear_map.hh"
using jtutil::error_exit;
#include "coords.hh"
#include "grid.hh"
#include "fd_grid.hh"
// all the code in this file is inside this namespace
namespace AHFinderDirect
{
//*****************************************************************************
//
// This function computes a single coefficient of a 1st derivative
// molecule, for unit grid spacing.
//
//static
fp fd_grid::dx_coeff(int m)
{
#ifndef FINITE_DIFF_ORDER
#error "must define FINITE_DIFF_ORDER!"
#endif
switch (m)
{
#if (FINITE_DIFF_ORDER == 2)
case -1: return FD_GRID__ORDER2__DX__COEFF_M1;
case 0: return FD_GRID__ORDER2__DX__COEFF_0;
case +1: return FD_GRID__ORDER2__DX__COEFF_P1;
#elif (FINITE_DIFF_ORDER == 4)
case -2: return FD_GRID__ORDER4__DX__COEFF_M2;
case -1: return FD_GRID__ORDER4__DX__COEFF_M1;
case 0: return FD_GRID__ORDER4__DX__COEFF_0;
case +1: return FD_GRID__ORDER4__DX__COEFF_P1;
case +2: return FD_GRID__ORDER4__DX__COEFF_P2;
#else
#error "unknown value " FINITE_DIFF_ORDER " for FINITE_DIFF_ORDER!"
#endif
default:
error_exit(ERROR_EXIT,
"***** fd_grid::dx_coeff(): m=%d is outside order=%d molecule radius=%d\n",
m, FINITE_DIFF_ORDER, FD_GRID__MOL_RADIUS); /*NOTREACHED*/
}
}
//*****************************************************************************
//
// This function computes a single coefficient of a 2nd derivative
// molecule, for unit grid spacing.
//
//static
fp fd_grid::dxx_coeff(int m)
{
#ifndef FINITE_DIFF_ORDER
#error "must define FINITE_DIFF_ORDER!"
#endif
switch (m)
{
#if (FINITE_DIFF_ORDER == 2)
case -1: return FD_GRID__ORDER2__DXX__COEFF_M1;
case 0: return FD_GRID__ORDER2__DXX__COEFF_0;
case +1: return FD_GRID__ORDER2__DXX__COEFF_P1;
#elif (FINITE_DIFF_ORDER == 4)
case -2: return FD_GRID__ORDER4__DXX__COEFF_M2;
case -1: return FD_GRID__ORDER4__DXX__COEFF_M1;
case 0: return FD_GRID__ORDER4__DXX__COEFF_0;
case +1: return FD_GRID__ORDER4__DXX__COEFF_P1;
case +2: return FD_GRID__ORDER4__DXX__COEFF_P2;
#else
#error "unknown value " FINITE_DIFF_ORDER " for FINITE_DIFF_ORDER!"
#endif
default:
error_exit(ERROR_EXIT,
"***** fd_grid::dxx_coeff(): m=%d is outside order=%d molecule radius=%d\n",
m, FINITE_DIFF_ORDER, FD_GRID__MOL_RADIUS); /*NOTREACHED*/
}
}
//******************************************************************************
} // namespace AHFinderDirect
|
