27 unsigned nz =
N,
nr = 0, nv = 0;
30 data = std::vector<poly2d_term>((N + 2) / 2, v3);
36 data[nv].coeff = -
data[nv - 1].coeff * (nz + 1) * (nz + 2) / (nr *
nr);
120 for (
unsigned jp = 2; jp <=
MaxM; ++jp)
126 out <<
"TrigArr: TASize = " <<
TASize <<
"\tMaxM = " <<
MaxM << std::endl;
130 std::streamsize old_prec = out.precision(), wdt = prec + 7;
133 for (jm = 0; jm <=
MaxM; ++jm) {
134 out << std::setw(wdt) << std::left << jm;
137 for (jm = 0; jm <=
MaxM; ++jm) {
141 for (jm = 0; jm <=
MaxM; ++jm) {
144 out <<
"|" << std::endl;
145 out.precision(old_prec);
147 out <<
"\tTrigArr size is not adjusted." << std::endl;
150 out <<
"\tTrigArr is not allocated." << std::endl;
159 p_out.
data.reserve(2 *
L);
164 for (it = 0; it <
data.size(); ++it) {
170 p_out.
data.push_back(term);
173 for (it = 0; it <
data.size(); ++it) {
179 p_out.
data.push_back(term);
183 if (!p_out.
data.empty()) {
198 p_out.
data.reserve(2 *
L);
203 for (it = 0; it <
data.size(); ++it) {
209 p_out.
data.push_back(term);
212 for (it = 0; it <
data.size(); ++it) {
218 p_out.
data.push_back(term);
222 if (!p_out.
data.empty()) {
static void FillTrigArr(const double phi)
std::vector< poly2d_term > data
static std::set< poly2d_base * > poly2d_base_set
static void SetPhi(const double phi)
void Scale(const double C)
static trig_pair * TrigArr
Abs< T >::type abs(const T &t)
static void PrintTrigArr(std::ostream &out=std::cout, const std::streamsize prec=5)
static void SetTrigArrSize(const unsigned N)