Nestedsums library: root_of_unity.h Source File

00001 00008 /* 00009 * Copyright (C) 2001 Stefan Weinzierl 00010 * 00011 * This program is free software; you can redistribute it and/or modify 00012 * it under the terms of the GNU General Public License as published by 00013 * the Free Software Foundation; either version 2 of the License, or 00014 * (at your option) any later version. 00015 * 00016 * This program is distributed in the hope that it will be useful, 00017 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00018 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00019 * GNU General Public License for more details. 00020 * 00021 * You should have received a copy of the GNU General Public License 00022 * along with this program; if not, write to the Free Software 00023 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00024 */ 00025 00026 #ifndef __NESTEDSUMS_ROOT_OF_UNITY_H__ 00027 #define __NESTEDSUMS_ROOT_OF_UNITY_H__ 00028 00029 #include "ginac/ginac.h" 00030 00031 namespace nestedsums { 00032 00041 class root_of_unity : public GiNaC::basic 00042 { 00043 GINAC_DECLARE_REGISTERED_CLASS(root_of_unity, GiNaC::basic) 00044 00045 // ctors 00046 public: 00047 root_of_unity(unsigned l, unsigned k); 00048 00049 // functions overriding virtual functions from base classes 00050 public: 00051 void archive(GiNaC::archive_node &node) const; 00052 void read_archive(const GiNaC::archive_node &node, GiNaC::lst &sym_lst); 00053 00054 unsigned return_type(void) const { return GiNaC::return_types::noncommutative; } 00055 00056 void print(const GiNaC::print_context & c, unsigned level = 0) const; 00057 unsigned precedence(void) const {return 50;} 00058 GiNaC::ex eval(int level = 0) const; 00059 GiNaC::ex subs(const GiNaC::exmap & m, unsigned options = 0) const; 00060 00061 protected: 00062 GiNaC::ex eval_ncmul(const GiNaC::exvector & v) const; 00063 GiNaC::ex derivative(const GiNaC::symbol & s) const; 00064 unsigned calchash(void) const; 00065 00066 // new virtual functions which can be overridden by derived classes 00067 00068 // non-virtual functions 00069 00070 // friends : 00071 00072 // member variables : 00073 00074 protected: 00075 00076 unsigned l; 00077 unsigned k; 00078 }; 00079 00080 GINAC_DECLARE_UNARCHIVER(root_of_unity); 00081 00082 // utility functions 00083 00089 inline GiNaC::ex create_root_of_unity(unsigned l ,unsigned k) 00090 { 00091 return (new root_of_unity(l,k))->setflag(GiNaC::status_flags::dynallocated); 00092 } 00093 00094 } // namespace nestedsums 00095 00096 #endif // ndef __NESTEDSUMS_ROOT_OF_UNITY_H__