X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Finifcns_nstdsums.cpp;h=4312377f10cf2d868972db41a56dac127973855d;hp=9dd590d35c6eb88626c8b8cab0c981fdf295d445;hb=cc94094751459129e1392a93dfa0264adac789a5;hpb=5d367f84b43b0636f943fe239f7a58b4d147f002

diff --git a/ginac/inifcns_nstdsums.cpp b/ginac/inifcns_nstdsums.cpp
index 9dd590d3..4312377f 100644
--- a/ginac/inifcns_nstdsums.cpp
+++ b/ginac/inifcns_nstdsums.cpp
@@ -47,7 +47,7 @@
  */
 
 /*
- *  GiNaC Copyright (C) 1999-2008 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2010 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
@@ -64,11 +64,6 @@
  *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
  */
 
-#include <sstream>
-#include <stdexcept>
-#include <vector>
-#include <cln/cln.h>
-
 #include "inifcns.h"
 
 #include "add.h"
@@ -84,6 +79,10 @@
 #include "utils.h"
 #include "wildcard.h"
 
+#include <cln/cln.h>
+#include <sstream>
+#include <stdexcept>
+#include <vector>
 
 namespace GiNaC {
 
@@ -210,7 +209,7 @@ void double_Xn()
 			}
 		}
 		// X_n
-		for (int n=2; n<Xn.size(); ++n) {
+		for (size_t n=2; n<Xn.size(); ++n) {
 			for (int i=xninitsize+1; i<=xend; ++i) {
 				if (i & 1) {
 					result = 0; // k == 0
@@ -570,11 +569,12 @@ ex G_eval(const Gparameter& a, int scale, const exvector& gsyms)
 
 			Gparameter newa;
 			Gparameter::const_iterator it2 = short_a.begin();
-			for (--it2; it2 != it;) {
-				++it2;
+			for (; it2 != it; ++it2) {
 				newa.push_back(*it2);
 			}
+			newa.push_back(*it);
 			newa.push_back(a[0]);
+			it2 = it;
 			++it2;
 			for (; it2 != short_a.end(); ++it2) {
 				newa.push_back(*it2);	
@@ -659,6 +659,8 @@ Gparameter::const_iterator check_parameter_G(const Gparameter& a, int scale,
 			++trailing_zeros;
 		}
 	}
+	if (lastnonzero == a.end())
+		return a.end();
 	return ++lastnonzero;
 }
 
@@ -1091,7 +1093,7 @@ G_do_trafo(const std::vector<cln::cl_N>& x, const std::vector<int>& s,
 	Gparameter a(x.size());
 	exmap subslst;
 	std::size_t pos = 1;
-	int scale;
+	int scale = pos;
 	for (sortmap_t::const_iterator it = sortmap.begin(); it != sortmap.end(); ++it) {
 		if (it->second < x.size()) {
 			if (s[it->second] > 0) {
@@ -1684,7 +1686,7 @@ static void Li_print_latex(const ex& m_, const ex& x_, const print_context& c)
 	} else {
 		x = lst(x_);
 	}
-	c.s << "\\mbox{Li}_{";
+	c.s << "\\mathrm{Li}_{";
 	lst::const_iterator itm = m.begin();
 	(*itm).print(c);
 	itm++;
@@ -2200,7 +2202,7 @@ static ex S_deriv(const ex& n, const ex& p, const ex& x, unsigned deriv_param)
 
 static void S_print_latex(const ex& n, const ex& p, const ex& x, const print_context& c)
 {
-	c.s << "\\mbox{S}_{";
+	c.s << "\\mathrm{S}_{";
 	n.print(c);
 	c.s << ",";
 	p.print(c);
@@ -2281,7 +2283,7 @@ bool convert_parameter_H_to_Li(const lst& l, lst& m, lst& s, ex& pf)
 		}
 	}
 	if (has_negative_parameters) {
-		for (int i=0; i<m.nops(); i++) {
+		for (std::size_t i=0; i<m.nops(); i++) {
 			if (m.op(i) < 0) {
 				m.let_op(i) = -m.op(i);
 				s.append(-1);
@@ -2321,7 +2323,7 @@ struct map_trafo_H_convert_to_Li : public map_function
 					s.let_op(0) = s.op(0) * arg;
 					return pf * Li(m, s).hold();
 				} else {
-					for (int i=0; i<m.nops(); i++) {
+					for (std::size_t i=0; i<m.nops(); i++) {
 						s.append(1);
 					}
 					s.let_op(0) = s.op(0) * arg;
@@ -2403,7 +2405,7 @@ struct map_trafo_H_reduce_trailing_zeros : public map_function
 					
 					//
 					parameter.remove_last();
-					int lastentry = parameter.nops();
+					std::size_t lastentry = parameter.nops();
 					while ((lastentry > 0) && (parameter[lastentry-1] == 0)) {
 						lastentry--;
 					}
@@ -2492,9 +2494,9 @@ ex trafo_H_mult(const ex& h1, const ex& h2)
 			hlong = h2.op(0).op(0);
 		}
 	}
-	for (int i=0; i<=hlong.nops(); i++) {
+	for (std::size_t i=0; i<=hlong.nops(); i++) {
 		lst newparameter;
-		int j=0;
+		std::size_t j=0;
 		for (; j<i; j++) {
 			newparameter.append(hlong[j]);
 		}
@@ -2522,7 +2524,7 @@ struct map_trafo_H_mult : public map_function
 			ex result = 1;
 			ex firstH;
 			lst Hlst;
-			for (int pos=0; pos<e.nops(); pos++) {
+			for (std::size_t pos=0; pos<e.nops(); pos++) {
 				if (is_a<power>(e.op(pos)) && is_a<function>(e.op(pos).op(0))) {
 					std::string name = ex_to<function>(e.op(pos).op(0)).get_name();
 					if (name == "H") {
@@ -2556,7 +2558,7 @@ struct map_trafo_H_mult : public map_function
 			if (Hlst.nops() > 0) {
 				ex buffer = trafo_H_mult(firstH, Hlst.op(0));
 				result *= buffer;
-				for (int i=1; i<Hlst.nops(); i++) {
+				for (std::size_t i=1; i<Hlst.nops(); i++) {
 					result *= Hlst.op(i);
 				}
 				result = result.expand();
@@ -2584,7 +2586,7 @@ ex trafo_H_1tx_prepend_zero(const ex& e, const ex& arg)
 	if (name == "H") {
 		h = e;
 	} else {
-		for (int i=0; i<e.nops(); i++) {
+		for (std::size_t i=0; i<e.nops(); i++) {
 			if (is_a<function>(e.op(i))) {
 				std::string name = ex_to<function>(e.op(i)).get_name();
 				if (name == "H") {
@@ -2616,7 +2618,7 @@ ex trafo_H_prepend_one(const ex& e, const ex& arg)
 	if (name == "H") {
 		h = e;
 	} else {
-		for (int i=0; i<e.nops(); i++) {
+		for (std::size_t i=0; i<e.nops(); i++) {
 			if (is_a<function>(e.op(i))) {
 				std::string name = ex_to<function>(e.op(i)).get_name();
 				if (name == "H") {
@@ -2647,7 +2649,7 @@ ex trafo_H_1tx_prepend_minusone(const ex& e, const ex& arg)
 	if (name == "H") {
 		h = e;
 	} else {
-		for (int i=0; i<e.nops(); i++) {
+		for (std::size_t i=0; i<e.nops(); i++) {
 			if (is_a<function>(e.op(i))) {
 				std::string name = ex_to<function>(e.op(i)).get_name();
 				if (name == "H") {
@@ -2680,7 +2682,7 @@ ex trafo_H_1mxt1px_prepend_minusone(const ex& e, const ex& arg)
 	if (name == "H") {
 		h = e;
 	} else {
-		for (int i=0; i<e.nops(); i++) {
+		for (std::size_t i = 0; i < e.nops(); i++) {
 			if (is_a<function>(e.op(i))) {
 				std::string name = ex_to<function>(e.op(i)).get_name();
 				if (name == "H") {
@@ -2711,7 +2713,7 @@ ex trafo_H_1mxt1px_prepend_one(const ex& e, const ex& arg)
 	if (name == "H") {
 		h = e;
 	} else {
-		for (int i=0; i<e.nops(); i++) {
+		for (std::size_t i = 0; i < e.nops(); i++) {
 			if (is_a<function>(e.op(i))) {
 				std::string name = ex_to<function>(e.op(i)).get_name();
 				if (name == "H") {
@@ -2749,7 +2751,7 @@ struct map_trafo_H_1mx : public map_function
 				// special cases if all parameters are either 0, 1 or -1
 				bool allthesame = true;
 				if (parameter.op(0) == 0) {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != 0) {
 							allthesame = false;
 							break;
@@ -2765,7 +2767,7 @@ struct map_trafo_H_1mx : public map_function
 				} else if (parameter.op(0) == -1) {
 					throw std::runtime_error("map_trafo_H_1mx: cannot handle weights equal -1!");
 				} else {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != 1) {
 							allthesame = false;
 							break;
@@ -2791,7 +2793,7 @@ struct map_trafo_H_1mx : public map_function
 					map_trafo_H_1mx recursion;
 					ex buffer = recursion(H(newparameter, arg).hold());
 					if (is_a<add>(buffer)) {
-						for (int i=0; i<buffer.nops(); i++) {
+						for (std::size_t i = 0; i < buffer.nops(); i++) {
 							res -= trafo_H_prepend_one(buffer.op(i), arg);
 						}
 					} else {
@@ -2805,13 +2807,13 @@ struct map_trafo_H_1mx : public map_function
 					map_trafo_H_1mx recursion;
 					map_trafo_H_mult unify;
 					ex res = H(lst(1), arg).hold() * H(newparameter, arg).hold();
-					int firstzero = 0;
+					std::size_t firstzero = 0;
 					while (parameter.op(firstzero) == 1) {
 						firstzero++;
 					}
-					for (int i=firstzero-1; i<parameter.nops()-1; i++) {
+					for (std::size_t i = firstzero-1; i < parameter.nops()-1; i++) {
 						lst newparameter;
-						int j=0;
+						std::size_t j=0;
 						for (; j<=i; j++) {
 							newparameter.append(parameter[j+1]);
 						}
@@ -2850,7 +2852,7 @@ struct map_trafo_H_1overx : public map_function
 				// special cases if all parameters are either 0, 1 or -1
 				bool allthesame = true;
 				if (parameter.op(0) == 0) {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != 0) {
 							allthesame = false;
 							break;
@@ -2860,7 +2862,7 @@ struct map_trafo_H_1overx : public map_function
 						return pow(-1, parameter.nops()) * H(parameter, 1/arg).hold();
 					}
 				} else if (parameter.op(0) == -1) {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != -1) {
 							allthesame = false;
 							break;
@@ -2872,7 +2874,7 @@ struct map_trafo_H_1overx : public map_function
 						       / factorial(parameter.nops())).expand());
 					}
 				} else {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != 1) {
 							allthesame = false;
 							break;
@@ -2895,7 +2897,7 @@ struct map_trafo_H_1overx : public map_function
 					map_trafo_H_1overx recursion;
 					ex buffer = recursion(H(newparameter, arg).hold());
 					if (is_a<add>(buffer)) {
-						for (int i=0; i<buffer.nops(); i++) {
+						for (std::size_t i = 0; i < buffer.nops(); i++) {
 							res += trafo_H_1tx_prepend_zero(buffer.op(i), arg);
 						}
 					} else {
@@ -2910,7 +2912,7 @@ struct map_trafo_H_1overx : public map_function
 					map_trafo_H_1overx recursion;
 					ex buffer = recursion(H(newparameter, arg).hold());
 					if (is_a<add>(buffer)) {
-						for (int i=0; i<buffer.nops(); i++) {
+						for (std::size_t i = 0; i < buffer.nops(); i++) {
 							res += trafo_H_1tx_prepend_zero(buffer.op(i), arg) - trafo_H_1tx_prepend_minusone(buffer.op(i), arg);
 						}
 					} else {
@@ -2924,13 +2926,13 @@ struct map_trafo_H_1overx : public map_function
 					map_trafo_H_1overx recursion;
 					map_trafo_H_mult unify;
 					ex res = H(lst(1), arg).hold() * H(newparameter, arg).hold();
-					int firstzero = 0;
+					std::size_t firstzero = 0;
 					while (parameter.op(firstzero) == 1) {
 						firstzero++;
 					}
-					for (int i=firstzero-1; i<parameter.nops()-1; i++) {
+					for (std::size_t i = firstzero-1; i < parameter.nops() - 1; i++) {
 						lst newparameter;
-						int j=0;
+						std::size_t j = 0;
 						for (; j<=i; j++) {
 							newparameter.append(parameter[j+1]);
 						}
@@ -2971,7 +2973,7 @@ struct map_trafo_H_1mxt1px : public map_function
 				// special cases if all parameters are either 0, 1 or -1
 				bool allthesame = true;
 				if (parameter.op(0) == 0) {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != 0) {
 							allthesame = false;
 							break;
@@ -2983,7 +2985,7 @@ struct map_trafo_H_1mxt1px : public map_function
 						       / factorial(parameter.nops())).expand());
 					}
 				} else if (parameter.op(0) == -1) {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != -1) {
 							allthesame = false;
 							break;
@@ -2995,7 +2997,7 @@ struct map_trafo_H_1mxt1px : public map_function
 						       / factorial(parameter.nops())).expand());
 					}
 				} else {
-					for (int i=1; i<parameter.nops(); i++) {
+					for (std::size_t i = 1; i < parameter.nops(); i++) {
 						if (parameter.op(i) != 1) {
 							allthesame = false;
 							break;
@@ -3018,7 +3020,7 @@ struct map_trafo_H_1mxt1px : public map_function
 					map_trafo_H_1mxt1px recursion;
 					ex buffer = recursion(H(newparameter, arg).hold());
 					if (is_a<add>(buffer)) {
-						for (int i=0; i<buffer.nops(); i++) {
+						for (std::size_t i = 0; i < buffer.nops(); i++) {
 							res -= trafo_H_1mxt1px_prepend_one(buffer.op(i), arg) + trafo_H_1mxt1px_prepend_minusone(buffer.op(i), arg);
 						}
 					} else {
@@ -3033,7 +3035,7 @@ struct map_trafo_H_1mxt1px : public map_function
 					map_trafo_H_1mxt1px recursion;
 					ex buffer = recursion(H(newparameter, arg).hold());
 					if (is_a<add>(buffer)) {
-						for (int i=0; i<buffer.nops(); i++) {
+						for (std::size_t i = 0; i < buffer.nops(); i++) {
 							res -= trafo_H_1mxt1px_prepend_minusone(buffer.op(i), arg);
 						}
 					} else {
@@ -3047,13 +3049,13 @@ struct map_trafo_H_1mxt1px : public map_function
 					map_trafo_H_1mxt1px recursion;
 					map_trafo_H_mult unify;
 					ex res = H(lst(1), arg).hold() * H(newparameter, arg).hold();
-					int firstzero = 0;
+					std::size_t firstzero = 0;
 					while (parameter.op(firstzero) == 1) {
 						firstzero++;
 					}
-					for (int i=firstzero-1; i<parameter.nops()-1; i++) {
+					for (std::size_t i = firstzero - 1; i < parameter.nops() - 1; i++) {
 						lst newparameter;
-						int j=0;
+						std::size_t j=0;
 						for (; j<=i; j++) {
 							newparameter.append(parameter[j+1]);
 						}
@@ -3127,7 +3129,7 @@ static ex H_evalf(const ex& x1, const ex& x2)
 			}
 		}
 
-		for (int i=0; i<x1.nops(); i++) {
+		for (std::size_t i = 0; i < x1.nops(); i++) {
 			if (!x1.op(i).info(info_flags::integer)) {
 				return H(x1, x2).hold();
 			}
@@ -3199,7 +3201,7 @@ static ex H_evalf(const ex& x1, const ex& x2)
 		// ensure that the realpart of the argument is positive
 		if (cln::realpart(x) < 0) {
 			x = -x;
-			for (int i=0; i<m.nops(); i++) {
+			for (std::size_t i = 0; i < m.nops(); i++) {
 				if (m.op(i) != 0) {
 					m.let_op(i) = -m.op(i);
 					res *= -1;
@@ -3406,7 +3408,7 @@ static void H_print_latex(const ex& m_, const ex& x, const print_context& c)
 	} else {
 		m = lst(m_);
 	}
-	c.s << "\\mbox{H}_{";
+	c.s << "\\mathrm{H}_{";
 	lst::const_iterator itm = m.begin();
 	(*itm).print(c);
 	itm++;
@@ -3482,7 +3484,7 @@ static void initcX(std::vector<cln::cl_N>& crX,
 	int Sm = 0;
 	int Smp1 = 0;
 	std::vector<std::vector<cln::cl_N> > crG(s.size() - 1, std::vector<cln::cl_N>(L2 + 1));
-	for (int m=0; m < s.size() - 1; m++) {
+	for (int m=0; m < (int)s.size() - 1; m++) {
 		Sm += s[m];
 		Smp1 = Sm + s[m+1];
 		for (int i = 0; i <= L2; i++)
@@ -3708,7 +3710,7 @@ cln::cl_N zeta_do_Hoelder_convolution(const std::vector<int>& m_, const std::vec
 	s_p[0] = s_p[0] * cln::cl_N("1/2");
 	// convert notations
 	int sig = 1;
-	for (int i=0; i<s_.size(); i++) {
+	for (std::size_t i = 0; i < s_.size(); i++) {
 		if (s_[i] < 0) {
 			sig = -sig;
 			s_p[i] = -s_p[i];