* functions. */
/*
- * GiNaC Copyright (C) 1999-2015 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2018 Johannes Gutenberg University Mainz, Germany
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
// in this case n more (or less) terms are needed
// (sadly, to generate them, we have to start from the beginning)
if (n == 0 && coeff == 1) {
- epvector epv;
- ex acc = dynallocate<pseries>(rel, epv);
- epv.reserve(2);
- epv.push_back(expair(-1, _ex0));
- epv.push_back(expair(Order(_ex1), order));
- ex rest = pseries(rel, std::move(epv)).add_series(argser);
+ ex rest = pseries(rel, epvector{expair(-1, _ex0), expair(Order(_ex1), order)}).add_series(argser);
+ ex acc = dynallocate<pseries>(rel, epvector());
for (int i = order-1; i>0; --i) {
epvector cterm { expair(i%2 ? _ex1/i : _ex_1/i, _ex0) };
acc = pseries(rel, std::move(cterm)).add_series(ex_to<pseries>(acc));
const ex &point = rel.rhs();
const symbol foo;
const ex replarg = series(log(arg), s==foo, order).subs(foo==point, subs_options::no_pattern);
- epvector seq { expair(-I*csgn(arg*I)*Pi, _ex0),
- expair(Order(_ex1), order) };
+ epvector seq;
+ if (order > 0) {
+ seq.reserve(2);
+ seq.push_back(expair(-I*csgn(arg*I)*Pi, _ex0));
+ }
+ seq.push_back(expair(Order(_ex1), order));
return series(replarg - I*Pi + pseries(rel, std::move(seq)), rel, order);
}
throw do_taylor(); // caught by function::series()
Order0correction += log((I*arg_pt+_ex_1)/(I*arg_pt+_ex1))*I*_ex_1_2;
else
Order0correction += log((I*arg_pt+_ex1)/(I*arg_pt+_ex_1))*I*_ex1_2;
- epvector seq { expair(Order0correction, _ex0),
- expair(Order(_ex1), order) };
+ epvector seq;
+ if (order > 0) {
+ seq.reserve(2);
+ seq.push_back(expair(Order0correction, _ex0));
+ }
+ seq.push_back(expair(Order(_ex1), order));
return series(replarg - pseries(rel, std::move(seq)), rel, order);
}
throw do_taylor();
Order0correction += log((arg_pt+_ex_1)/(arg_pt+_ex1))*_ex1_2;
else
Order0correction += log((arg_pt+_ex1)/(arg_pt+_ex_1))*_ex_1_2;
- epvector seq { expair(Order0correction, _ex0),
- expair(Order(_ex1), order) };
+ epvector seq;
+ if (order > 0) {
+ seq.reserve(2);
+ seq.push_back(expair(Order0correction, _ex0));
+ }
+ seq.push_back(expair(Order(_ex1), order));
return series(replarg - pseries(rel, std::move(seq)), rel, order);
}
throw do_taylor();