X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Ftensor.cpp;h=42056ed150d61668c5a5923fd26fa159bebd340d;hp=bae0cc1af8308632ed0cb6ba332ac1b67f613051;hb=10365850aa3803337bfea1fc201b81b6752096d4;hpb=695f6ae955ec530cded8f21efd5569df39447f76 diff --git a/ginac/tensor.cpp b/ginac/tensor.cpp index bae0cc1a..42056ed1 100644 --- a/ginac/tensor.cpp +++ b/ginac/tensor.cpp @@ -3,7 +3,7 @@ * Implementation of GiNaC's special tensors. */ /* - * GiNaC Copyright (C) 1999-2005 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2011 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 @@ -17,13 +17,9 @@ * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ -#include -#include -#include - #include "tensor.h" #include "idx.h" #include "indexed.h" @@ -36,6 +32,10 @@ #include "archive.h" #include "utils.h" +#include +#include +#include + namespace GiNaC { GINAC_IMPLEMENT_REGISTERED_CLASS(tensor, basic) @@ -64,7 +64,7 @@ GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(tensepsilon, tensor, // constructors ////////// -tensor::tensor() : inherited(TINFO_tensor) +tensor::tensor() { setflag(status_flags::evaluated | status_flags::expanded); } @@ -74,44 +74,34 @@ DEFAULT_CTOR(tensmetric) minkmetric::minkmetric() : pos_sig(false) { - tinfo_key = TINFO_minkmetric; } spinmetric::spinmetric() { - tinfo_key = TINFO_spinmetric; } minkmetric::minkmetric(bool ps) : pos_sig(ps) { - tinfo_key = TINFO_minkmetric; } tensepsilon::tensepsilon() : minkowski(false), pos_sig(false) { - tinfo_key = TINFO_tensepsilon; } tensepsilon::tensepsilon(bool mink, bool ps) : minkowski(mink), pos_sig(ps) { - tinfo_key = TINFO_tensepsilon; } ////////// // archiving ////////// -DEFAULT_ARCHIVING(tensor) -DEFAULT_ARCHIVING(tensdelta) -DEFAULT_ARCHIVING(tensmetric) -DEFAULT_ARCHIVING(spinmetric) -DEFAULT_UNARCHIVE(minkmetric) -DEFAULT_UNARCHIVE(tensepsilon) - -minkmetric::minkmetric(const archive_node &n, lst &sym_lst) : inherited(n, sym_lst) +void minkmetric::read_archive(const archive_node& n, lst& sym_lst) { + inherited::read_archive(n, sym_lst); n.find_bool("pos_sig", pos_sig); } +GINAC_BIND_UNARCHIVER(minkmetric); void minkmetric::archive(archive_node &n) const { @@ -119,11 +109,13 @@ void minkmetric::archive(archive_node &n) const n.add_bool("pos_sig", pos_sig); } -tensepsilon::tensepsilon(const archive_node &n, lst &sym_lst) : inherited(n, sym_lst) +void tensepsilon::read_archive(const archive_node& n, lst& sym_lst) { + inherited::read_archive(n, sym_lst); n.find_bool("minkowski", minkowski); n.find_bool("pos_sig", pos_sig); } +GINAC_BIND_UNARCHIVER(tensepsilon); void tensepsilon::archive(archive_node &n) const { @@ -132,6 +124,10 @@ void tensepsilon::archive(archive_node &n) const n.add_bool("pos_sig", pos_sig); } +GINAC_BIND_UNARCHIVER(tensdelta); +GINAC_BIND_UNARCHIVER(tensmetric); +GINAC_BIND_UNARCHIVER(spinmetric); + ////////// // functions overriding virtual functions from base classes ////////// @@ -141,6 +137,22 @@ DEFAULT_COMPARE(tensdelta) DEFAULT_COMPARE(tensmetric) DEFAULT_COMPARE(spinmetric) +bool tensdelta::info(unsigned inf) const +{ + if(inf == info_flags::real) + return true; + + return false; +} + +bool tensmetric::info(unsigned inf) const +{ + if(inf == info_flags::real) + return true; + + return false; +} + int minkmetric::compare_same_type(const basic & other) const { GINAC_ASSERT(is_a(other)); @@ -152,6 +164,14 @@ int minkmetric::compare_same_type(const basic & other) const return inherited::compare_same_type(other); } +bool minkmetric::info(unsigned inf) const +{ + if(inf == info_flags::real) + return true; + + return false; +} + int tensepsilon::compare_same_type(const basic & other) const { GINAC_ASSERT(is_a(other)); @@ -165,6 +185,22 @@ int tensepsilon::compare_same_type(const basic & other) const return inherited::compare_same_type(other); } +bool tensepsilon::info(unsigned inf) const +{ + if(inf == info_flags::real) + return true; + + return false; +} + +bool spinmetric::info(unsigned inf) const +{ + if(inf == info_flags::real) + return true; + + return false; +} + DEFAULT_PRINT_LATEX(tensdelta, "delta", "\\delta") DEFAULT_PRINT(tensmetric, "g") DEFAULT_PRINT_LATEX(minkmetric, "eta", "\\eta") @@ -328,13 +364,17 @@ ex tensepsilon::eval_indexed(const basic & i) const if (minkowski) { for (size_t j=1; j(x)) + if (!is_a(x)) { throw(std::runtime_error("indices of epsilon tensor in Minkowski space must be of type varidx")); - if (ex_to(x).is_covariant()) - if (ex_to(x).get_value().is_zero()) + } + if (ex_to(x).is_covariant()) { + if (ex_to(x).get_value().is_zero()) { sign = (pos_sig ? -sign : sign); - else + } + else { sign = (pos_sig ? sign : -sign); + } + } } } @@ -355,6 +395,8 @@ bool tensor::replace_contr_index(exvector::iterator self, exvector::iterator oth again: if (self_idx->is_symbolic()) { for (size_t i=1; inops(); i++) { + if (! is_a(other->op(i))) + continue; const idx &other_idx = ex_to(other->op(i)); if (is_dummy_pair(*self_idx, other_idx)) { @@ -588,6 +630,9 @@ ex epsilon_tensor(const ex & i1, const ex & i2) if (!ex_to(i1).get_dim().is_equal(_ex2)) throw(std::runtime_error("index dimension of epsilon tensor must match number of indices")); + if(is_a(i1.op(0))||is_a(i2.op(0))) + return indexed(epsilon, antisymmetric2(), i1, i2).hold(); + return indexed(epsilon, antisymmetric2(), i1, i2); } @@ -604,6 +649,9 @@ ex epsilon_tensor(const ex & i1, const ex & i2, const ex & i3) if (!ex_to(i1).get_dim().is_equal(_ex3)) throw(std::runtime_error("index dimension of epsilon tensor must match number of indices")); + if(is_a(i1.op(0))||is_a(i2.op(0))||is_a(i3.op(0))) + return indexed(epsilon, antisymmetric3(), i1, i2, i3).hold(); + return indexed(epsilon, antisymmetric3(), i1, i2, i3); } @@ -621,6 +669,9 @@ ex lorentz_eps(const ex & i1, const ex & i2, const ex & i3, const ex & i4, bool if (!ex_to(i1).get_dim().is_equal(_ex4)) throw(std::runtime_error("index dimension of epsilon tensor must match number of indices")); + if(is_a(i1.op(0))||is_a(i2.op(0))||is_a(i3.op(0))||is_a(i4.op(0))) + return indexed(pos_sig ? epsilon_pos : epsilon_neg, antisymmetric4(), i1, i2, i3, i4).hold(); + return indexed(pos_sig ? epsilon_pos : epsilon_neg, antisymmetric4(), i1, i2, i3, i4); }