/** @file lortensor.h
*
- * Interface to GiNaC´s lortensor objects. */
+ * Interface to GiNaC's Lorentz tensors. */
/*
- * GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2001 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
#include <string>
#include <vector>
-#include "config.h"
#include <iostream>
-#include <ginac/indexed.h>
-#include <ginac/lorentzidx.h>
+#include "indexed.h"
+#include "lorentzidx.h"
-#ifndef NO_GINAC_NAMESPACE
namespace GiNaC {
-#endif // ndef NO_GINAC_NAMESPACE
-/** Base class for lortensor object */
+
+/** This class holds an object carrying Lorentz indices (of class
+ * lorentzidx). It can represent a general (symbolic) tensor of type
+ * (p,q), or one of the constant tensors g (the metric), delta (unity
+ * matrix) or epsilon (4-dimensional totally antisymmetric tensor). */
class lortensor : public indexed
{
- // friends
- friend lortensor lortensor_g(ex const & mu, ex const & nu);
- // friend lortensor lortensor_delta(ex const & mu, ex const & nu);
- friend lortensor lortensor_epsilon(ex const & mu, ex const & nu,
- ex const & rho, ex const & sigma);
- friend lortensor lortensor_rankn(string const & n, exvector const & iv);
- friend lortensor lortensor_rank1(string const & n, ex const & mu);
- friend lortensor lortensor_rank2(string const & n, ex const & mu, ex const & nu);
- friend ex simplify_lortensor_mul(ex const & m);
- friend ex simplify_lortensor(ex const & e);
-
- // types
+ GINAC_DECLARE_REGISTERED_CLASS(lortensor, indexed)
+
+// friends
+
+ friend lortensor lortensor_g(const ex & mu, const ex & nu);
+ friend lortensor lortensor_delta(const ex & mu, const ex & nu);
+ friend lortensor lortensor_epsilon(const ex & mu, const ex & nu,
+ const ex & rho, const ex & sigma);
+ friend lortensor lortensor_vector(const std::string & n, const ex & mu);
+ friend lortensor lortensor_symbolic(const std::string & name, const exvector & iv);
+
+ friend ex simplify_lortensor_mul(const ex & m);
+ friend ex simplify_lortensor(const ex & e);
+
+// types
public:
- typedef enum { invalid,
- lortensor_g,
- lortensor_rankn,
- lortensor_rank1,
- lortensor_rank2,
- // lortensor_delta,
- lortensor_epsilon
- } lortensor_types;
-
- // member functions
-
- // default constructor, destructor, copy constructor assignment operator and helpers
-public:
- lortensor();
- ~lortensor();
- lortensor(lortensor const & other);
- lortensor const & operator=(lortensor const & other);
-protected:
- void copy(lortensor const & other);
- void destroy(bool call_parent);
+ typedef enum {
+ invalid, /**< not properly constructed */
+ lortensor_g, /**< metric tensor */
+ lortensor_delta, /**< unity matrix */
+ lortensor_epsilon, /**< four-dimensional totally antisymmetric tensor */
+ lortensor_symbolic /**< general symbolic Lorentz tensor */
+ } lortensor_types;
- // other constructors
+// member functions
+
+ // other constructors
protected:
- lortensor(lortensor_types const lt, string const & n);
- lortensor(lortensor_types const lt, string const & n, ex const & mu);
- lortensor(lortensor_types const lt, string const & n, ex const & mu, ex const & nu);
- lortensor(lortensor_types const lt, string const & n, ex const & mu, ex const & nu,
- ex const & rho);
- lortensor(lortensor_types const lt, string const & n, ex const & mu, ex const & nu, ex const & rho, ex const & sigma);
- lortensor(lortensor_types const lt, string const & n, exvector const & iv);
- lortensor(lortensor_types const lt, string const & n, unsigned s, exvector const & iv);
- lortensor(lortensor_types const lt, string const & n, unsigned s, exvector * ivp);
-
- //functions overriding virtual functions from base classes
+ lortensor(lortensor_types const lt, const std::string & n);
+ lortensor(lortensor_types const lt, const std::string & n, const ex & mu);
+ lortensor(lortensor_types const lt, const std::string & n, const ex & mu, const ex & nu);
+ lortensor(lortensor_types const lt, const std::string & n, const ex & mu, const ex & nu, const ex & rho);
+ lortensor(lortensor_types const lt, const std::string & n, const ex & mu, const ex & nu, const ex & rho, const ex & sigma);
+ lortensor(lortensor_types const lt, const std::string & n, const exvector & iv);
+ lortensor(lortensor_types const lt, const std::string & n, unsigned s, const exvector & iv);
+ lortensor(lortensor_types const lt, const std::string & n, unsigned s, exvector * ivp);
+
+ //functions overriding virtual functions from base classes
public:
- basic * duplicate() const;
- void printraw(ostream & os) const;
- void printtree(ostream & os, unsigned indent) const;
- void print(ostream & os, unsigned upper_precedence=0) const;
- void printcsrc(ostream & os, unsigned type, unsigned upper_precedence=0) const;
- bool info(unsigned inf) const;
- ex eval(int level=0) const;
+ void printraw(std::ostream & os) const;
+ void printtree(std::ostream & os, unsigned indent) const;
+ void print(std::ostream & os, unsigned upper_precedence=0) const;
+ bool info(unsigned inf) const;
+ ex eval(int level=0) const;
protected:
- int compare_same_type(basic const & other) const;
- bool is_equal_same_type(basic const & other) const;
- unsigned return_type(void) const;
- unsigned return_type_tinfo(void) const;
- ex thisexprseq(exvector const & v) const;
- ex thisexprseq(exvector * vp) const;
+ bool is_equal_same_type(const basic & other) const;
+ unsigned return_type(void) const;
+ unsigned return_type_tinfo(void) const;
+ ex thisexprseq(const exvector & v) const;
+ ex thisexprseq(exvector * vp) const;
- // new virtual functions which can be overridden by derived classes
- // none
+ // new virtual functions which can be overridden by derived classes
+ // none
- //non virtual functions in this class
+ //non virtual functions in this class
public:
- void setname(string const & n);
- string getname(void) const {return name;}
+ void setname(const std::string & n) {name = n;}
+ std::string getname(void) const {return name;}
protected:
- bool all_of_type_lorentzidx(void) const;
+ bool all_of_type_lorentzidx(void) const;
private:
- string & autoname_prefix(void);
+ std::string & autoname_prefix(void);
- //member variables
+ //member variables
protected:
- lortensor_types type;
- string name;
- unsigned serial;
+ lortensor_types type; /**< Type of object */
+ std::string name; /**< Name of symbolic tensor */
+ unsigned serial; /**< Unique serial number for comparing symbolic tensors */
private:
- static unsigned next_serial;
+ static unsigned next_serial;
};
-// global constants
-
- extern const lortensor some_lortensor;
- extern type_info const & typeid_lortensor;
-
-// utility functions !!!!!!!!!!!!!!!!!! (to do)
-
-
+// utility functions
+
inline const lortensor & ex_to_lortensor(const ex &e)
{
return static_cast<const lortensor &>(*e.bp);
return static_cast<lortensor &>(*e.bp);
}
-lortensor lortensor_g(ex const & mu, ex const & nu);
-ex simplify_lortensor_mul(ex const & m);
-ex simplify_lortensor(ex const & e);
-ex Dim(void);
+lortensor lortensor_g(const ex & mu, const ex & nu);
+lortensor lortensor_delta(const ex & mu, const ex & nu);
+lortensor lortensor_epsilon(const ex & mu, const ex & nu,
+ const ex & rho, const ex & sigma);
+lortensor lortensor_vector(const std::string & n, const ex & mu);
+lortensor lortensor_symbolic(const std::string & name, const exvector & iv);
+
+ex simplify_lortensor_mul(const ex & m);
+ex simplify_lortensor(const ex & e);
-#ifndef NO_GINAC_NAMESPACE
} // namespace GiNaC
-#endif // ndef NO_GINAC_NAMESPACE
#endif // ndef __GINAC_LORTENSOR_H__