3 * Implementation of GiNaC's clifford objects.
4 * No real implementation yet, to be done. */
7 * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany
9 * This program is free software; you can redistribute it and/or modify
10 * it under the terms of the GNU General Public License as published by
11 * the Free Software Foundation; either version 2 of the License, or
12 * (at your option) any later version.
14 * This program is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 * GNU General Public License for more details.
19 * You should have received a copy of the GNU General Public License
20 * along with this program; if not, write to the Free Software
21 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
33 #ifndef NO_NAMESPACE_GINAC
35 #endif // ndef NO_NAMESPACE_GINAC
37 GINAC_IMPLEMENT_REGISTERED_CLASS(clifford, lortensor)
40 // default constructor, destructor, copy constructor assignment operator and helpers
47 debugmsg("clifford default constructor",LOGLEVEL_CONSTRUCT);
48 tinfo_key = TINFO_clifford;
53 debugmsg("clifford destructor",LOGLEVEL_DESTRUCT);
57 clifford::clifford(const clifford & other)
59 debugmsg("clifford copy constructor",LOGLEVEL_CONSTRUCT);
63 const clifford & clifford::operator=(const clifford & other)
65 debugmsg("clifford operator=",LOGLEVEL_ASSIGNMENT);
75 void clifford::copy(const clifford & other)
77 inherited::copy(other);
80 void clifford::destroy(bool call_parent)
82 if (call_parent) inherited::destroy(call_parent);
91 clifford::clifford(const std::string & n, const ex & mu) : inherited(lortensor_symbolic, n, mu)
93 debugmsg("clifford constructor from string,ex",LOGLEVEL_CONSTRUCT);
94 tinfo_key=TINFO_clifford;
97 clifford::clifford(const std::string & n, const exvector & iv) : inherited(lortensor_symbolic, n, iv)
99 debugmsg("clifford constructor from string,exvector", LOGLEVEL_CONSTRUCT);
100 GINAC_ASSERT(all_of_type_lorentzidx());
101 tinfo_key=TINFO_clifford;
104 clifford::clifford(const std::string & n, exvector *ivp) : inherited(lortensor_symbolic, n, *ivp)
106 debugmsg("clifford constructor from string,exvector", LOGLEVEL_CONSTRUCT);
107 GINAC_ASSERT(all_of_type_lorentzidx());
108 tinfo_key=TINFO_clifford;
115 /** Construct object from archive_node. */
116 clifford::clifford(const archive_node &n, const lst &sym_lst) : inherited(n, sym_lst)
118 debugmsg("clifford constructor from archive_node", LOGLEVEL_CONSTRUCT);
121 /** Unarchive the object. */
122 ex clifford::unarchive(const archive_node &n, const lst &sym_lst)
124 return (new clifford(n, sym_lst))->setflag(status_flags::dynallocated);
127 /** Archive the object. */
128 void clifford::archive(archive_node &n) const
130 inherited::archive(n);
134 // functions overriding virtual functions from bases classes
139 basic * clifford::duplicate() const
141 debugmsg("clifford duplicate",LOGLEVEL_DUPLICATE);
142 return new clifford(*this);
145 void clifford::printraw(std::ostream & os) const
147 debugmsg("clifford printraw",LOGLEVEL_PRINT);
148 os << "clifford(" << "indices=";
150 os << ",hash=" << hashvalue << ",flags=" << flags << ")";
153 void clifford::printtree(std::ostream & os, unsigned indent) const
155 debugmsg("clifford printtree",LOGLEVEL_PRINT);
156 os << std::string(indent,' ') << " (clifford): "
157 << seq.size() << "indices=";
158 printtreeindices(os, indent);
159 os << ", hash=" << hashvalue
160 << " (0x" << std::hex << hashvalue << std::dec << ")"
161 << ", flags=" << flags << std::endl;
164 void clifford::print(std::ostream & os, unsigned upper_precedence) const
166 debugmsg("clifford print",LOGLEVEL_PRINT);
171 bool clifford::info(unsigned inf) const
173 return inherited::info(inf);
178 int clifford::compare_same_type(const basic & other) const
180 GINAC_ASSERT(is_of_type(other,clifford));
181 // only compare indices
182 return exprseq::compare_same_type(other);
185 bool clifford::is_equal_same_type(const basic & other) const
187 GINAC_ASSERT(is_of_type(other,clifford));
188 // only compare indices
189 return exprseq::is_equal_same_type(other);
192 ex clifford::thisexprseq(const exvector & v) const
194 return clifford(name, v);
197 ex clifford::thisexprseq(exvector *vp) const
199 return clifford(name, vp);
202 ex clifford::simplify_ncmul(const exvector & v) const
204 return simplified_ncmul(v);
211 /** Construct an object representing a Dirac gamma matrix. The index must
212 * be of class lorentzidx.
215 * @return newly constructed object */
216 clifford clifford_gamma(const ex & mu)
218 return clifford("gamma", mu);
221 #ifndef NO_NAMESPACE_GINAC
223 #endif // ndef NO_NAMESPACE_GINAC